筛法求约数个数

Description

用$D_i$表示$i$的约数个数，$nu{m}_i$表示$i$的最小质因子出现次数。

Proof

引理1：

$$若n=\prod _{i=1}^m{p}_i^{c_i},则D_i=\prod _{i=1}^m{c}_i+1\text{\hspace{1em}(1.1)}$$

引理2：

$$若n=i\ast j(gcd(i,j)=1),则D_n=D_i\ast D_j\text{\hspace{1em}(1.2)}$$

$D$为积性函数

solution

线性筛

实现：

void pre() {
  	d[1] = 1;
  	for (int i = 2; i <= n; ++i) {
    	if (!v[i])
            v[i] = 1, p[++tot] = i, d[i] = 2, num[i] = 1;
    	for (int j = 1; j <= tot && i <= n / p[j]; ++j) {
      		v[p[j] * i] = 1;
      		if (i % p[j] == 0) {
        		num[i * p[j]] = num[i] + 1;
        		d[i * p[j]] = d[i] / num[i * p[j]] * (num[i * p[j]] + 1);
        		break;
      		}
        	else {
        		num[i * p[j]] = 1;
        		d[i * p[j]] = d[i] * 2;
      		}
    	}
  	}
}