-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathsieve.cpp
44 lines (38 loc) · 1.07 KB
/
sieve.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
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
// TC: O(n) // not the traditional sieve of eratosthenes (O(n*log(log(n)))), requires memory
///////////////////// start yanking ////////////////////
const int N = 10000000;
vector<int> primes;
vector<bool> sieve(N+1, true);
int mpf[N+1]; // minimum prime factor
void build_sieve() { // call me (precomputation)
sieve[0] = 0, sieve[1] = 0;
for (int i = 2; i <= N; i++) {
if (sieve[i]) {
primes.push_back(i);
mpf[i] = i;
}
for (int j = 0; j < primes.size() && i*primes[j] <= N && primes[j] <= mpf[i]; j++) {
mpf[i*primes[j]] = primes[j];
sieve[i*primes[j]] = false;
}
}
}
bool is_prime(int n) { // call me (result)
assert(n < N && n > 0);
return sieve[n];
}
///////////////////// stop yanking /////////////////////
int main() {
build_sieve();
for (int i = 0; i < primes.size(); i++) {
cout << primes[i] << " ";
}
cout << "\n";
int n;
cin >> n;
cout << is_prime(n);
return 0;
}