线性筛板子

求素数

#include<cstdio>
#include<cstring>
using namespace std;
const int N=10000000;
bool noprime[N+5];
int prime[N+5],size;
void euler(int n)
{
    for(int i=2;i<=n;i++){
        if(!noprime[i]) prime[++size]=i;
        for(int j=1;j<=size&&i*prime[j]<=n;j++){
            noprime[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
    noprime[1]=1;//这里要记得特判 
}
int main()
{
    int n,m;
    scanf("%d %d",&n,&m);
    euler(n);
    while(m--){
        int x;
        scanf("%d",&x);
        if(noprime[x]) printf("No\n");
        else printf("Yes\n");
    }
    return 0;
}

欧拉函数——phi[ ]
莫比乌斯函数——mob[ ]
正约数个数——facnum[ ]
质因数个数——d[ ]

int prime[N],phi[N],mob[N],/**/facnum[N],d[N];bool noprime[N];
/*
莫比乌斯函数mob[i]
若i为奇数个不同素数之积mob[i] = -1
若i为偶数个不同素数之积mob[i] = 1
若i有平方因子则mob[i] = 0。
虽然在这里不用莫比乌斯函数 
*/
void euler()
{
    int size=0;
    phi[1]=1;mob[1]=1;facnum[1]=1;
    for(int i=2;i<=N-5;i++){
        if(!noprime[i]){
            prime[++size]=i;
            phi[i]=i-1;
            mob[i]=-1;
            facnum[i]=2;
            d[i]=1;
        } 
        for(int j=1;j<=size&&i*prime[j]<=N-5;j++){
            noprime[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=prime[j]*phi[i];//!!
                facnum[i*prime[j]]=facnum[i]/(1+d[i])*(2+d[i]);
                d[i*prime[j]]=d[i]+1;
                mob[i*prime[j]]=0;
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
            facnum[i*prime[j]]=facnum[i]*2;
            d[i*prime[j]]=1;
            mob[i*prime[j]]=mob[i]*(-1);
        }
    }
}