欧拉函数（C++）

欧拉函数

定义

\qquad 一组数 $x_1,x_2,\dots ,x_s$ 称为是模 $m$ 的既约剩余系（简称缩系），如果对 $\forall j,s\in Z,i\le j\le s,(x_j,m)=1$ ，并定义 $\varphi (m)=s$ ，即 { 1 , 2 , 3 , ⋯   , m } \{1,2,3,\cdots,m\} {1,2,3,⋯,m} 中与 $m$ 互素的个数，$\varphi (n)$称为欧拉函数。

\qquad 通俗地讲，欧拉函数 $\varphi (n)$ 表示1~n中与n互质的数的个数。

\qquad 举个例子：

\qquad 在1~6中，只有 $1$ 和 $5$ 两个数字与6互质，故 $\varphi (6)=2$ 。

\qquad 在1~18中，有 $1，2，5，7，11，13，17$ 一共6个数字与18互质，故 $\varphi (18)=6$ 。

性质

1. 如果一个数 $p$ 是质数，那么它的欧拉函数为 $p-1$，也就是说1~p-1的所有数字均与 $p$ 互质。

2. 设 $m=m_1{m}_2(m\in {\mathbb{N}}^{\ast })$，则有：若 $m_1$ 与 $m_2$ 互质，则 $\varphi (m)=\varphi (m_1)\varphi (m_2)$；若 $m_1$ 与 $m_2$ 不互质，则 $\varphi (m)=m_1\varphi (m_2)$，其中 $m_1\le m_2$.

3. 对于 $\forall m\in {\mathbb{N}}^{\ast }$，有
   $$\sum _{d|m}\varphi (d)=\sum _{d|m}\varphi \left(\frac{m}{d}\right)$$

4. 欧拉定理：
   $$若a与m互质，则a^{\varphi (m)}\equiv 1(modm).$$

公式

\qquad 根据算术基本定理，一个数N可以被唯一分解为有限个质数之积的形式：
$$N=P_1^{{\alpha }_1}{P}_2^{{\alpha }_2}{P}_3^{{\alpha }_3}\cdots P_k^{{\alpha }_k}$$
\qquad 在此基础上， $\varphi (N)$的计算公式为：
$$\varphi (N)=N\left(1-\frac{1}{P_1}\right)\left(1-\frac{1}{P_2}\right)\cdots \left(1-\frac{1}{P_k}\right)$$
\qquad 举个例子：
$$6=2^1\ast 3^1\Rightarrow \varphi (6)=6\ast \left(1-\frac{1}{2}\right)\left(1-\frac{1}{3}\right)=6\ast \frac{1}{2}\ast \frac{2}{3}=2$$

$$18=2^1\ast 3^2\Rightarrow \varphi (18)=18\ast \left(1-\frac{1}{2}\right)\left(1-\frac{1}{3}\right)=18\ast \frac{1}{2}\ast \frac{2}{3}=6$$

证明

\qquad 根据容斥原理：

\qquad 1.从 1~N 中去掉 $P_1,P_2,...P_k$ 的所有倍数。（共减 $C_k^1$ 项）
$$N-\frac{N}{P_1}-\frac{N}{P_2}-\cdots \frac{N}{P_k}$$
\qquad 2.再加上所有 $P_i\ast P_j$ 的倍数。（这些数字被重复去掉了，共加 $C_k^2$ 项）
$$+\frac{N}{P_1{P}_2}+\frac{N}{P_1{P}_3}+\frac{N}{P_1{P}_4}+\cdots +\frac{N}{P_{k-1}{P}_k}$$
\qquad 3.再减去所有 $P_i\ast P_j\ast P_h$ 的倍数。（共减 $C_k^3$ 项）
$$-\frac{N}{P_1{P}_2{P}_3}-\frac{N}{P_1{P}_2{P}_4}-\cdots -\frac{N}{P_{k-2}{P}_{k-1}{P}_k}$$
\qquad 以此类推…

\qquad 最终所得式子即为上述公式的展开式。

求一个数的欧拉函数

时间复杂度

\qquad 直接用公式求解，决速步为分解质因数，因此时间复杂度为 $O(\sqrt{n})$.

代码

由于C++中 / 为整除运算，我们也不希望出现小数，因此将公式稍微变形，将 $res\times \left(1-\frac{1}{P}\right)$ 的操作改为等价的 $\frac{res}{P}\times (P-1)$ 操作。

int phi(int n) {
	int res = n, p = 2;
	while (n != 1) {
		if (n % p == 0) {
			res = res / p * (p - 1);
		while (n % p == 0)	n /= p;
		}
		p++;
	}
	return res;
}

求n以内的每个数的欧拉函数

思路：

\qquad 利用素数线性筛（欧氏筛法），结合欧拉函数的一些性质，改写为一个线性的筛法。

时间复杂度：

$O(n)$

代码：

#include <iostream>
using namespace std;
const int N = 1e6 + 5;
int p[N], t;		//质数及下标
int phi[N];			//欧拉函数
bool s[N];			//是否为合数

void get_euler(int n) {		//计算n以内的欧拉函数
	phi[1] = 1;
	for (int i = 2; i <= n; i++) {
		if (!s[i]) {
			p[++t] = i;
			phi[i] = i - 1;
		}
		for (int j = 1; p[j] <= n / i; j++) {
			s[i * p[j]] = true;
			if (i % p[j] == 0) {
				phi[i * p[j]] = phi[i] * p[j];
				break;
			}
			phi[i * p[j]] = phi[i] * (p[j] - 1);
		}
	}
}

代码解释：

\qquad 第13行： 若 $i$ 为质数，则 $\varphi (i)=i-1$.

\qquad 第18行： 此时 $p[j]$ 为 $i$ 的最小质因数，若：
$$\begin{gathered} i=P_1^{{\alpha }_1}{P}_2^{{\alpha }_2}{P}_3^{{\alpha }_3}\cdots P_k^{{\alpha }_k} \\ \varphi (i)=i\left(1-\frac{1}{P_1}\right)\left(1-\frac{1}{P_2}\right)\cdots \left(1-\frac{1}{P_k}\right) \end{gathered}$$
\qquad 则：
$$\begin{gathered} i\ast p[j]=P_1^{{\alpha }_1+1}{P}_2^{{\alpha }_2}{P}_3^{{\alpha }_3}\cdots P_k^{{\alpha }_k} \\ \varphi (i\ast p[j])=i\times p[j]\times \left(1-\frac{1}{P_1}\right)\left(1-\frac{1}{P_2}\right)\cdots \left(1-\frac{1}{P_k}\right) \end{gathered}$$
\qquad 因此：
$$\varphi (i\ast p[j])=p[j]\times \varphi (i)$$
\qquad 第21行： 此时 $p[j]$ 小于 $i$ 的最小质因数，若：
$$\begin{gathered} i=P_1^{{\alpha }_1}{P}_2^{{\alpha }_2}{P}_3^{{\alpha }_3}\cdots P_k^{{\alpha }_k} \\ \varphi (i)=i\left(1-\frac{1}{P_1}\right)\left(1-\frac{1}{P_2}\right)\cdots \left(1-\frac{1}{P_k}\right) \end{gathered}$$
\qquad 则：
$$\begin{gathered} i\ast p[j]=P_1^{{\alpha }_1}{P}_2^{{\alpha }_2}{P}_3^{{\alpha }_3}\cdots P_k^{{\alpha }_k}\cdot p[j] \\ \varphi (i\ast p[j])=i\times p[j]\times \left(1-\frac{1}{P_1}\right)\left(1-\frac{1}{P_2}\right)\cdots \left(1-\frac{1}{P_k}\right)\cdot \left(1-\frac{1}{p[j]}\right) \end{gathered}$$
\qquad 因此：
$$\varphi (i\ast p[j])=\varphi (i)\times p[j]\times \left(1-\frac{1}{p[j]}\right)=\varphi (i)\times (p[j]-1)$$