线性筛获取约数个数

前提知识

首先
对于一个大于1正整数n可以分解质因数：
$$n=p_1^{d_1}+p_2^{d_2}+...+p_k^{d_k}，其中p_i为素数$$

那么根据排列组合，可以知道
$$约数的个数factor(n)=(1+d_1)(1+d_2)...(1+d_k)$$

开始筛选

首先我们引入一个数组$d(i)$表示$i$的最小质因子的次幂，(最小的原因是素数筛里每次都是用最小的质因子来筛合数的)

$$当i是素数时,显然有factor(i)=2,d[i]=1$$
记$p_j$表示一个素数

当i不是素数时，分两种情况;

$$当i\%p_j\ne 0且gcd(i,p_j)=1时,由积性函数的性质可得factor[i\ast p_j]=factor[i]\ast factor[p_j]=factor[i]\ast 2;d[i\ast p_j]=1$$
$$当i\%p_j=0时,说明p_j是i的最小质因子;因此factor[i\ast p_j]=factor[i]/(d[i]+1)\ast ((d[i]+1)+1);d[i\ast p_j]=d[i]+1;$$

Code

 	factor[1]=1;
    for(int i=2;i<=n;++i){
        if(!noprime[i]){
            d[i]=1;
            factor[i]=2;
            prime[++tot] = i;
        }
        for(int j=1;i*prime[j]<=n&&j<=tot;++j){
            noprime[i*prime[j]] = 1;
            if(i%prime[j]==0){
                factor[i*prime[j]]=factor[i]/(d[i]+1)*(d[i]+2);
                d[i*prime[j]] = d[i]+1;
                break;
            }
            factor[i*prime[j]]=factor[i]*2;
            d[i*prime[j]]=1;
        }
    }