#include #include #include using namespace std; const int MAX_PRIMES=100000000; bitset nums; // true for each prime vector primes; // list of primes void sieve(int n) { nums.set(); // set all values to true nums[0] = nums[1] = false; for(int i=2; i<=sqrt(n)+1; i++) // filter out non-primes if(nums[i]) for(int j=i*i; j