forked from neiljaviya/algos-1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
segmentedSieve.cpp
46 lines (42 loc) · 1.21 KB
/
segmentedSieve.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
//Counting or primes upto a number using segmented sieve
#include<bits/stdc++.h>
using namespace std;
int count_primes(int n) {
const int S = 10000;
vector<int> primes;
int nsqrt = sqrt(n);
vector<char> is_prime(nsqrt + 2, true);
for (int i = 2; i <= nsqrt; i++) {
if (is_prime[i]) {
primes.push_back(i);
for (int j = i * i; j <= nsqrt; j += i)
is_prime[j] = false;
}
}
int result = 0;
vector<char> block(S);
for (int k = 0; k * S <= n; k++) {
fill(block.begin(), block.end(), true);
int start = k * S;
for (int p : primes) {
int start_idx = (start + p - 1) / p;
int j = max(start_idx, p) * p - start;
for (; j < S; j += p)
block[j] = false;
}
if (k == 0)
block[0] = block[1] = false;
for (int i = 0; i < S && start + i <= n; i++) {
if (block[i])
result++;
}
}
return result;
}
int main()
{
cout<<"Enter n upto which you want count of prime numbers"<<endl;
int n;
cin>>n;
cout<<"Number of primes upto "<<n<<" are "<<count_primes(n)<<endl;
}