【OI学习笔记】数论-欧拉函数

“读读欧拉，他是所有人的老师”——法国数学家拉普拉斯

upd:
2022/8/29 修正文章中的 Latex 用法

板块：基础数论
前置知识：唯一分解定理
难度：较易
前置知识一览:

• 唯一分解定理：对于一个数 $n\in N$ 且 $n$ 不为质数，则 $n$ 可以被唯一分解为 $n=\prod p_i^{k_i}$（$p_i$ 为质数）。最早的证明由欧几里得给出。

什么是欧拉函数-定义

计算 $1\sim n$ 中与 $n$ 成互质数的个数的称为欧拉函数，记作 $\phi (n)$.形式化地，设 $n=\prod _{i=1}^n{p}_i^{k_i}$，则 $\phi (n)=n\prod _{i=1}^n(1-\frac{1}{p_i})=\prod _{i=1}^n{p}_i^{k_i-1}\times (p_i-1)$。

欧拉函数为什么对-证明

已知欧拉函数是一个积性函数，当 $m,n$ 互质时，$\phi (nm)=\phi (m)\times \phi (n)$，
由唯一分解定理可知$n=\prod _{i=1}^x{p}_i^{k_i}$，
所以 $\phi (n)=\prod _{i=1}^x\phi (p_i^{k_i})$，
对于 $\forall s,\phi (p_s^{k_s})=p_s^{k_s}-p_s^{k_s-1}$，
从定义出发 $\phi (p_s^{k_s})$ 等于小于或等于 $p_s^{k_s}$ 的正整数中与 $p_s^{k_s}$ 互质的数的总数，
从 1 到 $p_s^{k_s}$ 中共有 $p_s^{k_s}$ 个数字，其中与 $p_s^{k_s}$ 不互质的数有 $p_s,2p_s,...,p_s^{k_s-1}\times p_s$，共 $p_s^{k_s-1}$ 项，
所以 $\phi (p_s)=p_s^{k_s}-p_s^{k_s-1}=p_s^{a_s}\times (1-\frac{1}{p_s})$.
因此有：
$$\begin{gathered} \phi (n)=\prod _{i=1}^n\phi (p_i^{k_i}) \\ =\prod _{i=1}^n(p_i^{k_i}-p_i^{k_i-1}) \\ =\prod _{i=1}^n{p}_i^{k_i}\times (1-\frac{1}{p_i}) \\ =\prod _{i=1}^n{p}_i^{k_i}\prod _{i=1}^n(1-\frac{1}{p_i}) \\ =n\times \prod _{i=1}^n(1-\frac{1}{p_i}). \end{gathered}$$
至于我给出的另外一个形式，是这样得来的，由上式的倒数第二行得出：
$$\begin{gathered} \prod _{i=1}^n{p}_i^{k_i}\times (1-\frac{1}{p_i}) \\ =\prod _{i=1}^n(p_i^{k_i}-p_i^{k_i-1}) \\ =\prod _{i=1}^n(p_i^{k_i-1}\times p_i-p_i^{k_i-1}) \\ =\prod p_i^{k_i-1}\times (p_i-1) \end{gathered}$$

欧拉函数怎么实现-代码

公式法求 $\phi (n)$

这里用的公式是我给出的两个公式的前者。
因为每次乘的是 $1-\frac{1}{p_i}$，那么把 1 看做 $\frac{p_i}{p_i}$，也就是$\frac{p_i}{p_i}-\frac{1}{p_i}=\frac{p_i-1}{p_i}$.

	    scanf("%d", &n);
        int res = n;
        for (int i = 2; i <= n / i; i++)
        {
            if (n % i == 0)
            {
                res = res / i * (i - 1); //做了一点点小小的优化
                while (n % i == 0) n /= i;
            }
        }
        if (n > 1) res = res / n * (n - 1);
        printf("%d\n", res);

用欧拉筛求所有数的欧拉函数和

#include <iostream>
#include <cstdio>

using namespace std;

const int N = 1e6 + 10;

typedef long long LL;

int n;
int primes[N], phi[N], cnt;
bool st[N];

LL get_eulers(int x)
{
    phi[1] = 1;
    for (int i = 2; i <= x; i++)
    {
        if (!st[i])
        {
            primes[cnt++] = i;
            phi[i] = i - 1;
        }
        for (int j = 0; primes[j] <= x / i; j++)
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
            {
                phi[primes[j] * i] = phi[i] * primes[j];
                break;
            }
            phi[primes[j] * i] = phi[i] * (primes[j] - 1);
        }
    }
    
    LL res = 0;
    for (int i = 1; i <= x; i++) res += phi[i];
    
    return res;
}

int main()
{
    scanf("%d", &n);
    printf("%lld", get_eulers(n));
    
    return;
}

*扩展：

• 欧拉定理：$\forall a\perp n,a^{\phi (n)}\equiv 1(\mathrm{mod}n)$。对于任意一个与 $n$ 互质的数 $a$，$a^{\phi (n)}$ 与 1 对模 $n$ 同余。