筛法求素数 函数模板

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;
        }
    }
}