Advertisement
Guest User

Sieve

a guest
May 25th, 2016
55
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 0.52 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2.  
  3. using namespace std;
  4. const int sz = 10000000;
  5. bool isnt_prime[sz + 1];
  6. void fix() {
  7.     isnt_prime[0] = isnt_prime[1] = true;
  8.     for (int i = 4; i <= sz; i += 2) isnt_prime[i] = true;
  9.     for (int i = 3; i <= sqrt(sz); i++) {
  10.         if (!isnt_prime[i]) {
  11.             for (int j = 3; j * i <= sz; j += 2) isnt_prime[i * j] = true;
  12.         }
  13.     }
  14.  
  15. }
  16. int main() {
  17.     fix();
  18.     for (int i = 0; i <= 100; i++) {
  19.         if (!isnt_prime[i]) cout << i << " ";
  20.     }
  21.     return 0;
  22. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement