筛法求素数 函数模板

const int MAX=20000000;
//an存储素数数组，下标从0开始，个数为num
int an[2000000],num;
bool bn[MAX+10];

void prime()
{
num = 0;
memset(bn,1,sizeof(bn));
bn[0] = bn[1] = 0;
for(int i = 2;i<MAX;i++)
{
if(bn[i]) an[num++] = i;
for(int j = 0;j<num&&an[j]*i<MAX;j++)
{
bn[an[j]*i] = 0;
if(i%an[j] == 0) break;
}
}
}