Advertisement
Rapptz

bitset prime sieve

Dec 31st, 2012
396
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 0.46 KB | None | 0 0
  1. #include <bitset>
  2. #include <iostream>
  3.  
  4. int main () {
  5.     using ulong = unsigned long long;
  6.     ulong sum = 0;
  7.     ulong j = 0;
  8.     const ulong bound = 2000000;
  9.     std::bitset<bound> sieve;
  10.     sieve.set(0);
  11.     sieve.set(1);
  12.     for (ulong i = 2; i < bound; ++i) {
  13.         if (sieve[i] == false) {
  14.             sum += i;
  15.             for (j = 2*i; j < bound; j += i) {
  16.                 sieve.set(j);
  17.             }
  18.         }
  19.     }
  20.     std::cout << sum;
  21. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement