Технологии программирования, Решето Эратосфена

Пайтон? Да ну, он же медленный для работы с таким количеством информации. А вот С/С++ для этого - самое то! Информацию о всех простых числах на отрезке [0; 4294967295] получить легко - у меня даже на телефоне Самсунг Гэлэкси это происходит за чуть больше чем шесть с половиной минут, а на нормальном "железе" ещё быстрее. А что дальше делать с этой информацией? Код внизу просто заполняет битовый массив информацией о простоте или непростоте всех неотрицательных чисел до 4294967295 включительно, а потом по введённому с консоли числy из рассматриваемого диапазона говорит да, если введено простое число, и нет, если наоборот. Для хранения данных, как уже́ было сказано, надо чуть более полгигабайта. Число 409 в начале скрина - это время работы сита Эратосфена в секундах для данного диапазона чисел.
#include <iostream>
#include <cstdlib>
#include <ctime>
#include <cmath>
using namespace std;
unsigned long *a, b[32], c[32];
void fal(unsigned long x)
{ a[x >> 5] &= c[x & 31]; }
unsigned bit(unsigned long x)
{ return a[x >> 5] & b[x & 31]; }
int main()
{ unsigned long i, j, k, l, m, n = 4294967295, t;
unsigned long long N; i = 1; b[0] = i; c[0] = ~i;
t = time(NULL); for (j = 1; j < 32; j++) { i <<= 1;
b[j] = i; c[j] = ~i; } m = ceil((n + 1.) / 32);
a = new unsigned long [m]; for (i = 0; i <= m; i++)
a[i] = 4294967295; fal(0); fal(1); i = ceil(sqrt(n));
for (j = 2; j <= i; j++) if (bit(j)) { k = n / j; for (l = 2;
l <= k; l++) fal(l * j); } cout << time(NULL) - t << endl;
for(;;) { cout << "N » "; cin >> N; t = time(NULL);
if (N <= n ) { if (bit(N)) cout << "yes" << endl;
else cout << "no" << endl; continue; } } }