欧拉函数有一个简单的定义:对于任意一个正整数 n,欧拉函数 φ(n) 表示小于等于 n 的正整数中与 n 互质的数的个数。
互质又叫做互素,指的是两个数的最大公约数为 1。例如,2 和 3 就是一对互质的数,因为它们没有其他公共因子。
因此,当 n 是质数的时候,φ(n) 就是 n-1,因为小于等于 n 的正整数中只有 1 到 n-1 都和 n 互质。
当 n 不是质数时,我们可以将 n 分解成多个质数的乘积,即 n = p1^k1 * p2^k2 * ... * pn^kn。其中,p1, p2, ..., pn 是 n 的所有质因数。
因为任意一个小于等于 n 的正整数 x,要么和 n 相互质,要么和 n 的某个质因数 pi 有公共因子。因此,我们可以考虑计算出所有和 n 有公共因子的正整数的个数,然后用 n 减去这个数就是小于等于 n 的正整数中与 n 互质的数的个数了。
一个小于等于 n 的正整数 x 与 n 有公共因子的充分必要条件是 x 是 n 的某个质因数 pi 的倍数。具体而言,x 要么是 p1, 2p1, 3p1, ..., k1p1 中的一个;要么是 p2, 2p2, 3p2, ..., k2p2 中的一个;以此类推。其中,k1, k2, ..., kn 是满足 pk^k <= n 的最大整数。
因此,我们可以计算出每个质因数 pi 对应的有公共因子的正整数的个数,然后使用容斥原理将它们相加减去重复计算的部分。
对于一个质数 p,小于等于 n 的正整数中与之互质的数的个数就是 p-1,因为除了 p 本身,其它小于等于 p 的正整数都与之互质。
因此,我们可以计算出小于等于 n 的所有正整数的欧拉函数值 φ(n),然后用 n 减去 φ(n) 就是小于等于 n 且与 n 互质的数的个数了。
以下是欧拉函数实现:
int phi(int n) { // 定义一个函数 phi,它接受一个正整数 n 作为参数,返回 n 的欧拉函数值。
int ans = n; // 初始化 ans 为 n,这是因为如果 n 是质数,则 φ(n) 的值就为 n - 1。
for (int i = 2; i * i <= n; i++) { // 遍历 2 至 sqrt(n) 的所有数 i。
if (n % i == 0) { // 如果 n 能够被 i 整除,即 i 是 n 的一个质因数。
ans = ans / i * (i - 1); // 则计算出 i 对 φ(n) 的贡献并将其从 ans 中除去。
while (n % i == 0) n /= i; // 不断地将 n 除以 i,直到 n 不能再被 i 整除。
}
}
if (n > 1) ans = ans / n * (n - 1); // 处理 n 的最后一个质因数,并将其贡献加入 ans 中。
return ans; // 返回 n 的欧拉函数值。
}
综上所述,这个代码的主要思路是遍历 n 的所有小于等于 sqrt(n) 的质因数 i,计算出 i 对 φ(n) 的贡献并将其从 ans 中除去,然后将剩余的质因数的贡献加入 ans 中。其中,当处理每个质因数 i 的时候,我们使用了一个 while 循环将 n 不断地除以 i,以便处理 i 的多次出现。最后,如果 n 是大于 1 的数,即 n 还有至少一个大于 sqrt(n) 的质因数,则需要将其贡献加入 ans 中。