【数论】欧拉函数

1 数论函数

函数： 函数具备两个要素——值域和定义域。
$$\underset{值域}{\underbrace{f(x)}}=\underset{定义域}{\underbrace{x-1}}$$

• 定义域：变量 $x$ 的取值范围
• 值域：对应的函数 $f(x)$ 的取值范围

数论函数： 非特殊要求情况下，定义域为正整数的函数。

对于数论函数，有一下三种类型：

1. 普通的函数
2. 完全积性函数： 对于两个整数 $n$，$m$，函数满足条件

$$f(mn)=f(m)f(n)$$

3. 积性函数： 对于两个素数 $n$，$m$，函数满足条件

$$f(mn)=f(m)f(n)$$

常数函数：
$$1(n)=1$$
单位函数：

$[表达式]$：如果表达式为真，则等于 $1$ ；反之等于 $0$。
$$\epsilon (n)=[n=1]$$
单位函数也是一种完全积性函数。
$$\epsilon (mn)=\{\begin{array}{ll}1& m=1,n=1\\ 0& m\ne 1\vee n\ne 1\end{array}$$

2 欧拉函数

$\phi (n)$：$[1,n]$ 中与 $n$ 互质的数的个数$(n\in {\mathbb{N}}^{\ast })$。

定义式：
$$\phi (n)=\sum _{i=1}^n[(i,n)=1]$$
欧拉函数式积性函数（非完全）。

欧拉函数公式：
$$\phi (n)=n\prod _{\begin{array}{rl}& p\mid n\\ & p为素数\end{array}}\left(1-\frac{1}{p}\right)$$
补充：容斥原理

规律为“奇加偶减”，例如三个集合的容斥：$A+B+C-A\cap B-A\cap C-B\cap C+A\cap B\cap C$

证明：

$讨论n由两种质因子p,q的情况:$

$[1,n]中p的倍数：$
$$p,2p,3p,\dots ,\frac{n}{p}p$$
$[1,n]中q的倍数：$
$$q,2q,3q,\dots ,\frac{n}{q}q$$
$则一共有：$
$$\frac{n}{p}+\frac{n}{q}-\frac{n}{pq}$$
$则[1,n]中剩余的数即为互质的数：$
$$\begin{array}{rl}\phi (n)& =n-\frac{n}{p}-\frac{n}{q}+\frac{n}{pq}\\ & =n\left(n-\frac{1}{p}\right)\left(n-\frac{1}{q}\right)\end{array}$$
$同理，得到公式：$
$$\phi (n)=n\prod _{\begin{array}{rl}& p\mid n\\ & p为素数\end{array}}\left(1-\frac{1}{p}\right)$$
性质：
$$\phi (n)=n-1(n为素数)\text{\hspace{1em}(1)}$$

$$若n=p^k,\phi (n)=p^k-p^{k-1}(p为素数)\text{\hspace{1em}(2)}$$

$$\phi (2n)=\phi (n)(n为素数)\text{\hspace{1em}(3)}$$

$$\sum _{d\mid n}\phi (d)=n\text{\hspace{1em}(4)}$$

$$若p\mid n,且p^2\mid n,则\phi (n)=\phi \left(\frac{n}{p}\right)p\text{\hspace{1em}(5)}$$

C++求欧拉函数Euler_func()：

void Euler_func(int n) {
    is_prime[0] = is_prime[1] = true;
    size = 0;
    for (int i = 2; i <= n; i++) {
        if (!is_prime[i]) {
            prime[++size] = i;  
            phi[i] = i - 1;   //性质1
        }
        for (int j = 1; j <= size && i * prime[j] <= n; j++) {
            is_prime[i * prime[j]] = true;
            if (i % prime[j] == 0) {
                phi[i * prime[j]] = phi[i] * prime[j];  //性质5
                break;
            }
            phi[i * prime[j]] = phi[i] * phi[prime[j]];  //积性函数
        }
    }
}

3 莫比乌斯函数

莫比乌斯函数：
$$\mu (n)=\{\begin{array}{ll}1& n=1\\ (-1{)}^r& n={\prod }_{i=1}^r{p}_i(p_1\ne p_2\ne \dots p_r,p_{1\sim r}为质数)\\ 0& else\end{array}$$
$\mu (n)=(-1{)}^r$ 的条件也可以说成 $n$ 无平方因子。

莫比乌斯函数是积性函数（非完全）。

C++求莫比乌斯函数Mobius_func()：

void Mobius_func(int n) {
    is_prime[0] = is_prime[1] = true;
    size = 0;
    mu[1] = 1;
    for (int i = 2; i <= n; i++) {
        if (!is_prime[i]) {
            prime[++size] = i;  
            mu[i] = -1;  //素数无平方因子，且只有一个质因子，r=1，mu(n) = -1
        }
        for (int j = 1; j <= size && i * prime[j] <= n; j++) {
            is_prime[i * prime[j]] = true;
            if (i % prime[j] == 0) {
                mu[i * prime[j]] = 0;  //有平方因子，mu(n) = 0
                break;
            }
            mu[i * prime[j]] = -mu[i];  //积性函数，mu(n) * mu(p) = -mu(n)
        }
    }
}

4 欧拉定理和欧拉降幂

欧拉定理：
$$若(a,m)=1,则a^{\phi (b)-1}\equiv 1(\mathrm{mod}m)$$
观察欧拉定理与费马小定理，可以发现，费马小定理是欧拉定理的特例。

根据欧拉定理，对于求 $x$ 的问题 $a^b\equiv x(\mathrm{mod}m)$，可以化简原式，将 $a$ 的指数 $b$ 减小，也就是欧拉降幂。

欧拉降幂：
$$a^b=\{\begin{array}{ll}{a}^b& b<\phi (m)\\ a^{(b\mathrm{mod}\phi (m))+\phi (m)}& b\ge \phi (m)\end{array}(\mathrm{mod}m)$$

5 例题：二元一次不定方程的解

题目描述：

需要求解的方程为：
$$ax+by=1$$
其中 $a$，$b$ 为已知的正整数。

找出 $x$ 的最小非负整数解和对应的 $y$，没有输出sorry。

题解：
$$\begin{array}{rl}ax+by& =1\\ ax& =1-by\\ ax& \equiv 1(\mathrm{mod}b)\end{array}$$
也就是说，$x$ 是 $a$ 模 $b$ 的逆元。

欧拉定理求逆元：
$$x=a^{\phi (b)-1}\mathrm{mod}b((a,b)=1)$$
在本题中，如果出现 $(a,b)\ne 1$ 的情况，则无符合题意的解，即输出sorry。