欧拉函数模板

欧拉函数：$\phi (n)$表示$1\sim n$中与$n$互素的个数。
互素：两个整数的公因子只有$1$。

假设 $n$ 可以表示为 $n=P_1^{{\alpha }_1}\ast P_2^{{\alpha }_2}\ast ...\ast P_k^{{\alpha }_k}$，其中 $P_i$ 均为 $n$ 的质因子，那么存在结论 $\phi (n)=n\ast (1-\frac{1}{P_1})\ast (1-\frac{1}{P_2})\ast ...\ast (1-\frac{1}{P_k})$。

---

计算 $\phi (n)$ ：

#include<bits/stdc++.h>
using namespace std;

int main()
{
	int n;
	cin >> n;
	int ans = n;
	for (int i = 2;i <= n / i;i ++) {
		if (n % i == 0) {
			ans = ans / i * (i - 1); // 避免出现小数 
			while (n % i == 0) n /= i;
		}
	}
	if (n > 1) ans = ans / n * (n - 1);
	cout << ans << endl;

	return 0;
}

---

计算 $1\sim n$ 的欧拉函数之和：

可由线性筛法（欧拉筛 / 线性筛）同步求得：

1. 初始化 $\phi (1)=1$；

2. 当枚举到的 $i$ 是素数时。一个素数的与之互素的个数为 $1\sim i-1$，即 $i-1$ 个，即

   $\phi (i)=i-1$；

3. 当枚举到的 $i$ 满足 $i\%P_j=0$，其中 $P_j$ 为线性筛法中第二层循环中枚举到的已保存的素数。$i$ 可以表示为 $P_1^{{\alpha }_1}\ast ...\ast P_t^{{\alpha }_t}$，其中 $P$ 为 $i$ 的质因子，由于 $i\%P_j==0$，所以 $P_j$ 也是 $i$ 的质因子，因此 $P_j\le P_t$。

   $P_j\ast i=P_1^{{\alpha }_1}\ast ...\ast P_j^{{\alpha }_j+1}\ast ...\ast P_t^{{\alpha }_t}$，

   $\phi (i)=i\ast (1-\frac{1}{P_1})\ast (1-\frac{1}{P_2})\ast ...\ast (1-\frac{1}{P_t})$，

   由于欧拉函数只与质因子的种类有关而与每个质因子的个数无关，故

   $\phi (P_j\ast i)=P_j\ast i\ast (1-\frac{1}{P_1})\ast (1-\frac{1}{P_2})\ast ...\ast (1-\frac{1}{P_t})$，

   进一步表示为

   $\phi (P_j\ast i)=\phi (i)\ast P_j$；

4. 当枚举到的 $i$ 满足 $i\%P_j\ne 0$。由于 $P_j$ 是素数，所以 $P_j\ast i$ 比 $i$ 多了一个素数 $P_j$，

   $P_j\ast i=P_1^{{\alpha }_1}\ast ...\ast P_t^{{\alpha }_t}\ast P_j$，故

   $\phi (P_j\ast i)=P_j\ast i\ast (1-\frac{1}{P_1})\ast (1-\frac{1}{P_2})\ast ...\ast (1-\frac{1}{P_t})\ast (1-\frac{1}{P_j})$，

   进一步表示为

   $\phi (P_j\ast i)=P_j\ast \phi (i)\ast (1-\frac{1}{P_j})=\phi (i)\ast (P_j-1)$.

以上就是我们要加的四行代码。

---

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1e6+10;

int phi[N]; // phi[i] 表示 i 的欧拉函数值
int p[N], cnt; // 统计素数
int vis[N];
LL ans, n; 

int main()
{
	cin >> n;
	phi[1] = 1; // 情况1，注意初始化！ 
	for (int i = 2;i <= n;i ++) { // 线性筛 2~n 
		if (!vis[i]) { // 说明i是素数 
			phi[i] = i - 1; // 情况2
			p[cnt++] = i;
		}
		for (int j = 0;j < cnt;j ++) {
			if (i * p[j] > n) break;
			vis[i * p[j]] = 1;
			if (i % p[j] == 0) {
				phi[i * p[j]] = phi[i] * p[j]; // 情况3
				break;
			}
			phi[i * p[j]] = phi[i] * (p[j] - 1); // 情况4 
		}
	} 
	
	for (int i = 1;i <= n;i ++)
		ans += phi[i]; 

	cout << ans << endl;
	return 0;
}