除数函数d(n)与约数和

除数函数 $d(n)$

定义：除数函数 $d(n)$ 是一个算术函数，用于计算正整数 $n$ 的正约数个数。

例如 $d(6)=4$, 因为 $1,2,3,6$ 都是它的约数。

---

公式

① 若 $n=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_k^{a_k}$， 其中 $p_i$ 为质数，则 $d(n)=(a_1+1)\ast (a_2+1)\ast ...\ast (a_k+1)$

这一点很好理解，将 $n$ 分解质因数后，每个质因子都有 $a_i+1$ 种取法（可以不取，即 $0$ 个），做一个乘法原理，即为 $n$ 的约数个数

---

求解 d(n)

一、求解单个值

设 $n=a\ast b,a<=b$, 则有 $a\le \sqrt{n}$ , 所以我们只需要枚举 $\sqrt{n}$ 内的约数 $a$， 大于 $\sqrt{n}$ 的约数可以通过 $\frac{n}{a}$ 得到。注意一下，当 $a=b=\sqrt{n}$ 时，别重复计数。时间复杂度 $O(\sqrt{n})$

int dn = 0;
for (int i = 1; i*i <= n; i++) {
	if (n % i == 0) {
		dn += 2;
		if (i * i == n) dn--;
	}
}

二、求解$d(1)\sim d(n)$

可以使用线性筛求解， 时间复杂度 $O(n)$

// 用 d[i] 表示 i 的约数个数，num[i] 表示 i 的最小质因子出现次数
void calc_dn()
{
	d[1] = 1;
	for (int i = 2; i < N; i++) {
		if (!flag[i]) prime[++lp] = i, d[i] = 2, num[i] = 1;
		for (int j = 1; j <= lp && i * prime[j] < N; j++) {
			flag[prime[j] * i] = 1;
			if (i % prime[j] == 0) {
				num[i * prime[j]] = num[i] + 1;
				d[i * prime[j]] = d[i] / num[i * prime[j]] * (num[i * prime[j]] + 1);
				break;
			}
			else {
				num[i * prime[j]] = 1;
				d[i * prime[j]] = d[i] * 2;
			}
		}
	}
}

除数函数的上界

参考论文

感兴趣可以自己看上述论文，这里给出结论：

$\log d(n)<1.066\frac{\log n}{\log \log n},n\ge 3$

大多数人关心的应该是具体的数值，下边给出：

① $10^9$ 以内，$d(n)\le 1344$

② $2^{31}-1$ 以内，$d(n)\le 1600$

③ $10^{18}$ 以内，$d(n)\le 103680$

④ $2^{63}-1$ 以内，$d(n)\le 161280$

参考：除数函数的渐近上界？

约数和

咕咕咕