[线性筛][筛素数/筛约数个数]

线性筛中每个数只会被最小的素因子筛一次

1)线性筛素数

void getprime(int siz){
    memset(isprime,1,sizeof(isprime));
    isprime[1]=0;
    for(int i=2;i<=siz;i++){
        if(isprime[i])prime[++tot]=i;
        for(int j=1;j<=tot&&i*prime[j]<=siz;j++){
            isprime[i*prime[j]]=0;
            if(i%prime[j]==0)break;
        }
    }
}

2)线性筛约数个数

每个数的约数个数为(a1+1)*(a2+1)*...*(ak+1)，首先如果数i是素数，可以直接得到约数个数w[i]=2，又由于线性筛的时候，每个数是被最小素因子筛掉，所以如果i%prime[j]!=0，那么prime[j]将是i*prime[j]的最小素数，并且个数为1，所以w[i*prime[j]]=w[i]*(1+1)，如果i%prime[j]==0，那么prime[j]既是i*prime[j]的最小素数，又是i的最小素数(因为prime[j]是i*prime[j]的最小素数嘛，所以不可能存在比prime[j]更小的素数)，所以此时w[i*prime[j]]=w[i]/(n[i]+1)*(n[i*prime[j]]+1)，其中n[i]表示i的最小素数的个数，那么显然n[i*prime[j]]=n[i]+1

void getprime(int siz){
    memset(isprime,1,sizeof(isprime));
    isprime[1]=0;
    for(int i=2;i<=siz;i++){
        if(isprime[i]){prime[++tot]=i;w[i]=2;ww[i]=1;}
        for(int j=1;j<=tot&&i*prime[j]<=siz;j++){
            isprime[i*prime[j]]=0;
            if(i%prime[j]!=0){
                w[i*prime[j]]=w[i]*2;
                ww[i*prime[j]]=1;
            }
            else{
                ww[i*prime[j]]=ww[i]+1;
                w[i*prime[j]]=w[i]/(ww[i]+1)*(ww[i*prime[j]]+1);
                break;
            }
        }
    }
}