我在 spoj.com 上遇到了这个问题
http://www.spoj.com/problems/PRIME1/
在阅读了一些关于素数生成算法的信息后，我发现“阿特金斯筛子”是最快的素数生成算法。在本教程http://www.geeksforgeeks.org/sieve-of-atkin/的帮助下，我已经部分实现了算法
当我运行代码时，它给了我分段错误。调试后我才知道这段代码不是最优的，它使用了更多的内存和存储空间。有人可以告诉我如何优化程序。这是我的代码片段。是的，它是不完整的。
#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int sieveOfAtkin(long int limit)
{
// if limit is either 2 or 3, print them
if (limit >= 2)
{
cout << 2 << endl;
}
if (limit >= 3)
{
cout << 3 << endl;
}
bool sieve[limit];
for (long int i=0; i<=limit; i++)
sieve[i] = false;
for (long int x=1; x<limit; x++)
{
for (long int y=1; y<limit; y++)
{
long int n = (4*x*x) + (y*y);
if (n<=limit && (n%12==1 || n%12==5))
sieve[n] ^= true;
n = (3*x*x) + (y*y);
if (n<=limit && (n%12==7))
sieve[n] ^= true;
n = (3*x*x) - (y*y);
if (x>y && n<=limit && (n%12==11))
sieve[n] ^= true;
}
}
for (long int r=5; r*r<limit; r++)
{
if(sieve[r])
{
for(long int i=r*r; i<limit; i+=r*r)
sieve[r] = false;
}
}
for (long int i=5; i<limit; i++)
if(sieve[i])
cout << i << endl;
return 0;
}
int main(void)
{
long int x;
cin >> x;
sieveOfAtkin(x);
return 0;
}