线性筛d&&sd（约数个数&&约数和）

$d$函数

$d(n)$表示n的约数个数和

百度一下，我们就可以得到一个计算d函数的暴力做法

设$n=p_1^{a_1}p_2^{a_2}...p_k^{a_k}$
$d(n)=(a_1+1)(a_2+1)...(a_k+1)$

显然d是一个积性函数，我们可以用线性筛求出

设$prime[i]$表示第$i$个质数
$num[i]$表示$i$的最小质因子出现次数

首先我们可以观察一下线性筛素数：

void prepare() {
    for (int i=2;i<N;i++) {
        if (!no[i]) sshu[++tot]=i;
        for (int j=1;j<=tot&&sshu[j]*i<N;j++) {
            no[sshu[j]*i]=1;
            if (i%sshu[j]==0) {
                no[sshu[j]*i]=1;
                break;
            }
        }
    }
}

我们之所以可以保证上述算法是线性的
就是因为：if (i%sshu[j]==0) break;
也就是说，每个数只会被ta的最小质因子筛掉

① 当前数是质数

很显然，$d[n]=2,num[n]=1$

②i%prime[j]!=0

这种情况下，$i$没有$prime[j]$这个质因子
然而$i*prime[j]$中包含了一个$prime[j]$
前面我们已经得到了$d[i]$
$d[i]=(a_1+1)(a_2+1)...(a_k+1)$

然后我们在这之后补上一个$prime[j]$
$d[i*prime[j]]=(a_1+1)(a_2+1)...(a_k+1)(1+1)$
$d[i*prime[j]]=d[i]*d[prime[j]]=d[i]*2$

而且由于当前的$prime[j]$必然是$i*prime[j]$的最小素因子(因为从小到大枚举)， 我们要记录下这个最小素因子的个数

$num[i*prime[j]] = 1$

③i%prime[j]=0

这种情况下，i中必然包含了至少一个$prime[j]$
而且$prime[j]$也必定是$i$的最小质因数

而$i*prime[j]$比起$i$则是多了一个最小质因子

$d[i]=(a_1+1)(a_2+1)...(a_k+1)$
$d[i*prime[j]]=(a_1+1+1)(a_2+1)...(a_k+1)$

那么$d[i*prime[j]]=d[i]/(a_1+1)*(a_1+2)$
$num[i*prime[j]]=num[i]+1$

void prepare() {
    d[1]=1; num[1]=1;
    for (int i=2;i<N;i++) {
        if (!no[i]) {
            sshu[++tot]=i;
            d[i]=2; num[i]=1;
        }
        for (int j=1;j<=tot&&sshu[j]*i<N;j++) {
            int v=sshu[j]*i;
            no[v]=1;
            if (i%sshu[j]==0) {
                num[v]=num[i]+1;
                d[v]=d[i]/num[v]*(num[v]+1);
                break;
            }
            d[v]=d[i]<<1; num[v]=1;
        }
    }
    //for (int i=1;i<=10;i++) printf("%d\n",d[i]);
}

$sd$函数

$sd(n)$表示$n$的所有约数之和

再次百度一下，我们就可以得到一种比较暴力的$sd$计算方法

设$n=p_1^{a_1}p_2^{a_2}...p_k^{a_k}$
$sd(n)=(1+p_1+p_1^2+p_1^3+...+p_1^{a_1})(1+p_2+p_2^2+...+p_2^{a_2})...(1+p_3+p_3^2+...+p_3^{a_3})$

又是一个积性函数，我们可以考虑用线性筛筛出sd函数

设$prime[i]$表示第$i$个质数
$sp[i]$表示$i$的$(1+p_{1}+p_{1}^2+...+p_{1}^{a_1})$

① 当前数是质数

显然$sd[n]=n+1,sq[n]=n+1$

②i%prime[j]!=0

这种情况，$i$中不包含$prime[j]$这个质因子
而且$prime[j]$一定是$i*prime[j]$的最小质因子

$sd[i]=(1+p_1+p_1^2+...+p_1^{a_1})(1+p_2+p_2^2+...+p_2^{a_2})...(1+p_3+p_3^2+...+p_3^{a_3})$
$sd[i*prime[j]]=(1+p_1+p_1^2+...+p_1^{a_1})...(1+p_3+p_3^2+...+p_3^{a_3})(1+prime[j])$

$sd[i*prime[j]]=sd[i]*sd[prime[j]]$
$sp[i*prime[j]]=sp[prime[j]]=prime[j]+1$

③i%prime[j]=0

这种情况下，$i$中至少包含了一个$prime[j]$

那么：
$sd[i]=(1+p_1+p_1^2+...+p_1^{a_1})(1+p_2+p_2^2+...+p_2^{a_2})...(1+p_3+p_3^2+...+p_3^{a_3})$
$sd[i*prime[j]]=(1+p_1+p_1^2+...+p_1^{a_1}+p_1^{a_1+1})...(1+p_3+p_3^2+...+p_3^{a_3})$

我们观察一下两式唯一的不同：
$(1+p_1+p_1^2+...+p_1^{a_1})$
$(1+p_1+p_1^2+...+p_1^{a_1}+p_1^{a_1+1})$

$sd[i*prime[j]]=sd[i]/(1+p_1+p_1^2+...+p_1^{a_1})*(1+p_1+p_1^2+...+p_1^{a_1}+p_1^{a_1+1})$

那么
$sd[i*prime[j]]=sd[i]/sp[i]*(sp[i]*prime[j]+1)$
$sp[i*prime[j]]=sp[i]*prime[j]+1$

void prepare() {
    sd[1]=1; sp[1]=1;
    for (int i=2;i<N;i+=) {
        if (!no[i]) {
            sshu[++tot]=i;
            sd[i]=i+1;
            sp[i]=i+1;
        }
        for (int j=1;j<=tot&&sshu[j]*i<N;j++) {
            int v=sshu[j]*i;
            no[v]=1;
            if (i%sshu[j]==0) {
                sp[v]=sp[i]*sshu[j]+1;
                sd[v]=sd[i]/sp[i]*sp[v];
                break;
            }
            sd[v]=sd[i]*sd[sshu[j]];
            sp[v]=sshu[j]+1;
        }
    }
}