素数筛选模板

const int N=1e6+11;
int p[N],pr[N],cnt;
void init()
{
    for(int i=2;i<N;i++)
    {
        if(!p[i])    pr[++cnt]=i;
        for(int j=1;j<=cnt&&i*pr[j]<N;j++)
        {
            p[i*pr[j]]=1;
            if(i%pr[j]==0)    break;
        }
    }
}

 此时大素数的判定：

int pd(int x){
    for(int i=1;pr[i]*pr[i]<=x;i++){
        if(x%pr[i]==0)
            return 0;
    }
    return 1;
}