2019-08-13  2024-09-15    1282 字  3 分钟

给定任意正整数 $math_inline$n$math_inline$ ,那么在小于等于 $math_inline$n$math_inline$ 的所有正整数之中,有多少个与 $math_inline$n$math_inline$ 构成互质关系?

计算这个值的方法就叫做欧拉函数 $math_inline$\phi(n)$math_inline$ 表示:在 $math_inline$1$math_inline$ 到 $math_inline$n$math_inline$ 之中,与n构成互质关系的数的数量。

分析

情况一

如果 $math_inline$n=1$math_inline$ ,则 $math_inline$\phi(1)=1$math_inline$ 。因为1与任何数(包括自身)都构成互质关系

情况二

如果 $math_inline$n$math_inline$ 是质数,则 $math_inline$\phi(n)=n-1$math_inline$ 。因为质数与小于它的每一个数,都构成互质关系。比如 $math_inline$5$math_inline$ 与 $math_inline$1、2、3、4$math_inline$ 都构成互质关系。

情况三

如果 $math_inline$n$math_inline$ 是质数的某一个次方,即 $math_inline$n=p^k$math_inline$ ( $math_inline$p$math_inline$ 为质数, $math_inline$k$math_inline$ 为大于等于 $math_inline$1$math_inline$ 的整数),则 $math_inline$\phi(p^k)=p^k-p^{k-1}$math_inline$

比如 $math_inline$\phi(8)=\phi(2^3)=2^3-2^2=8-4=4$math_inline$

这是因为只有当一个数不包含质数 $math_inline$p$math_inline$ ,才有可能与 $math_inline$n$math_inline$ 互质。而包含质数 $math_inline$p$math_inline$ 的数一共有 $math_inline$p^{k-1}$math_inline$ 个。

即 $math_inline$1\times p、2\times p、3\times p、... 、p^{k-1}\times p$math_inline$ ,把他们去除,剩下的就是与 $math_inline$n$math_inline$ 互质的数。

上面的式子还可以写成下面的形式:

$math_inline$\phi(p^k)=p^k-p^{k-1}=p^k(1-\frac{1}{p})$math_inline$

可以看出,上面的第二种情况是 $math_inline$k=1$math_inline$ 时的特例

情况四

如果n可以分解成两个互质的整数之积: $math_inline$n=p_1\times p_2$math_inline$

则: $math_inline$\phi(n)=\phi(p_1p_2)=\phi(p_1)\phi(p_2)$math_inline$

即积的欧拉函数等于各个因子的欧拉函数之积。比如 $math_inline$\phi(56)=\phi(8\times 7)=\phi(8)\phi(7)=4\times 6=24$math_inline$

对于素数 $math_inline$p$math_inline$ , $math_inline$\phi(p)=p-1$math_inline$ ,对于两个素数 $math_inline$p,q$math_inline$ , $math_inline$\phi(pq)=pq-1$math_inline$

欧拉函数是积性函数,但不是完全积性函数

欧拉函数的积性

若 $math_inline$m,n$math_inline$ 互质,则 $math_inline$\phi(mn)=\phi(m)\phi(n)$math_inline$ 。

由“ $math_inline$m,n$math_inline$ 互质”可知 $math_inline$m,n$math_inline$ 无公因数,所以:

$math_inline$\phi(m)\phi(n)=m(1-\frac{1}{p_1})(1-\frac{1}{p_2})(1-\frac{1}{p_3})...(1-\frac{1}{p_n})\cdot n(1-\frac{1}{p_1'})(1-\frac{1}{p_2'})(1-\frac{1}{p_3'})...(1-\frac{1}{p_n'})$math_inline$

其中 $math_inline$p_1、p_2、p_3...p_n$math_inline$ 为 $math_inline$m$math_inline$ 的质因数, $math_inline$p_1'、p_2'、p_3'...p_n'$math_inline$ 为 $math_inline$n$math_inline$ 的质因数,而 $math_inline$m,n$math_inline$ 无公因数。

所以 $math_inline$p_1、p_2、p_3...p_n、p_1'、p_2'、p_3'...p_n$math_inline$ 互不相同,所以 $math_inline$p_1、p_2、p_3...p_n、p_1'、p_2'、p_3'...p_n$math_inline$ 均为 $math_inline$mn$math_inline$ 的质因数且为 $math_inline$mn$math_inline$ 的质因数的全集。

所以 $math_inline$\phi(mn)=mn(1-\frac{1}{p_1})(1-\frac{1}{p_2})(1-\frac{1}{p_3})...(1-\frac{1}{p_n})\cdot (1-\frac{1}{p_1'})(1-\frac{1}{p_2'})(1-\frac{1}{p_3'})...(1-\frac{1}{p_n'})$math_inline$

所以 $math_inline$\phi(mn)=\phi(m)\phi(n)$math_inline$ 。

以上部分涉及到“中国剩余定理”

情况五

因为任意一个大于1的正整数,都可以写成一系列质数的积。(分解质因数)

$math_inline$n=p_1^{k_1}p_2^{k_2}...p_r^{k_r}$math_inline$

根据结论四,得到

$math_inline$\phi(n)=\phi(p_1^{k_1})\phi(p_2^{k_2})...\phi(p_r^{k_r})$math_inline$

再根据结论三,得到

$math_inline$\phi(n)=p_1^{k_1}p_2^{k_2}...p_r^{k_r}(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_r})$math_inline$

也就等于

$math_inline$\phi(n)=n(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_r})$math_inline$

结论

比如 $math_inline$\phi(1323)=(3^3\times 7^2)=1323(1-\frac{1}{3})(1-\frac{1}{7})=756$math_inline$

即 $math_inline$\phi(mn)=\phi(n)\phi(m)$math_inline$ 只在 $math_inline$(m,n)=1$math_inline$ 时成立。

对于一个正整数 $math_inline$N$math_inline$ 的素数幂分解 $math_inline$N=p_1^{q_1}p_2^{q_2}...p_n^{q_n}$math_inline$

$math_inline$\phi(N)=N(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_n})$math_inline$

除了 $math_inline$N=2, \phi(N)$math_inline$ 都是偶数

设 $math_inline$N$math_inline$ 为正整数, $math_inline$\sum{\phi(d)}=N(d|N)$math_inline$

欧拉函数

根据性质二,我们可以在 $math_inline$O(\sqrt{n})$math_inline$ 的时间复杂度内求出一个数的欧拉函数值。

Code
//直接求解欧拉函数  
int euler(int n) { //返回euler(n)   
    int res = n, a = n;
    for (int i = 2; i * i <= a; i++) {
        if (a % i == 0) {
            res = res / i * (i - 1);//先进行除法是为了防止中间数据的溢出   
            while (a % i == 0) a /= i;
        }
    }
    if (a > 1) res = res / a * (a - 1);
    return res;
}

欧拉函数线性筛

如果要求 $math_inline$N$math_inline$ 以内的所有数的欧拉函数 $math_inline$O(n)$math_inline$

Code
const int MAXN = 1e6;
//欧拉线性筛:在线性时间内筛素数的同时求出所有数的欧拉函数
int tot;
int phi[MAXN]; //保存各个数字的欧拉函数
int prime[MAXN]; //按顺序保存素数
bool mark[MAXN]; //判断是否是素数
void get_phi(int N){
    phi[1] = 1;
    for(int i = 2; i <= MAXN; i++){ //相当于分解质因数的逆过程
        if(!mark[i]){
            prime[++tot] = i;
            phi[i] = i-1;
        }
        for(int j = 1; j <= tot; j++){
            if(i * prime[j] > N) break;
            mark[i * prime[j]] = 1; //确定i*prime[j]不是素数
            if(i % prime[j] == 0){ //判断prime[j] 是否为 i的约数
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
            else{ //prime[j] - 1 就是 phi[prime[j]],利用了欧拉函数的积性
                phi[i * prime[j]] = phi[i] * (prime[j] - 1);
            }
        }
    }
}