-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy pathimprovedSieve.cpp
35 lines (32 loc) · 884 Bytes
/
improvedSieve.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
//Given N, find all primes less than N in O(N) time
#include <bits/stdc++.h>
#define ll long long
#define MAX_SIZE 100001
using namespace std;
vector<ll> prime; //to store the primes
vector<bool> isPrime(MAX_SIZE, true);
vector<ll> SPF(MAX_SIZE); //SPF[i] stores smallest prime factor of i
void improvedSieve(ll n) {
ll i,j;
isPrime[0] = isPrime[1] = false;
for(i=2;i<n;i++) {
if(isPrime[i]) {
prime.push_back(i);
SPF[i] = i; //since smallest prime factor of a prime is the prime itself
}
for(j=0;j<prime.size() && i*prime[j]<n && prime[j]<=SPF[i]; j++) {
isPrime[i*prime[j]] = false;
SPF[i*prime[j]] = prime[j];
}
}
}
int main() {
ll n;
cin>>n;
improvedSieve(n);
for(ll i=0;i<prime.size();i++ ) {
cout<<prime[i]<<" ";
}
cout<<endl;
return 0;
}