-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathPRIME1.cpp
70 lines (63 loc) · 1.47 KB
/
PRIME1.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
/**
Rasul Kerimov (CoderINusE)
PRIME1 - Prime Generator
Segmented sieve of Eratosthenes
*/
#include <bits/stdc++.h>
#define MAXN 100005
using namespace std;
typedef pair<int,int> p;
typedef long long ll;
pair <int, int> a[MAXN + 2];
vector <int> primes;
int pr[MAXN + 2];
void simple_eratosphen(ll limit) {
int n = (int) sqrt(limit) + 1;
int d[n + 3];
memset(d, 0, sizeof(d));
d[0] = d[1] = 1;
for(int i = 2; i <= n; ++i) {
if(!d[i]) {
if(1ll * i * i <= n)
for(int j = i * i; j <= n; j += i)
d[j] = 1;
}
}
for(int i = 2; i <= n; ++i) {
if(!d[i]) primes.push_back(i);
}
}
void erat(ll n) {
simple_eratosphen(n);
}
void segment(ll l, ll r) {
int len = r - l + 1;
memset(pr, 0, sizeof(pr));
for(int i = 0; i < primes.size(); ++i) {
ll left = (l / primes[i]) * primes[i];
if(left < l) left += primes[i];
if(left == primes[i]) left += primes[i];
if(left > r) break;
while(left <= r) {
pr[(int)(left - l)] = 1;
left += primes[i];
}
}
if(l == 1) pr[0] = 1;
for(int i = 0; i < r - l + 1; ++i) {
if(pr[i] == 0) {
printf("%lld\n", l + i);
}
}
printf("\n");
}
int main() {
int n; scanf("%d", &n);
ll l, r;
erat(1e9 + 1);
for(int i = 0; i < n; ++i) {
scanf("%lld%lld", &l, &r);
assert(l <= r);
segment(l, r);
}
}