欧拉函数c++

欧拉函数有一个简单的定义:对于任意一个正整数 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 中。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值