算法基础课——第四章 数学知识（二）

---

第四章 数学知识（二）

---

如无特殊说明，所有数均为正整数.

---

欧拉函数

---

• 互质：若 $gcd(a,b)=1$​ 则称 $a,b$​​​ 互质.

• 欧拉函数的定义：

$1\sim N$​ 中与 $N$​ 互质的数的个数被称为欧拉函数，记为 $\phi (N)$​.
若在算数基本定理中，$N=P_1^{a_1}{P}_2^{a_2}\cdots P_m^{a_m}$​，则：$\phi (N)=N\times \frac{P_1-1}{P_1}\times \frac{P_2-1}{P_2}\times \cdots \times \frac{P_m-1}{P_m}$.​​

• $\phi (N)=N\times \frac{P_1-1}{P_1}\times \frac{P_2-1}{P_2}\times \cdots \times \frac{P_m-1}{P_m}$.

例：$1\sim 6$ 中与 $6$ 互质的数有 $1,5$，则 $\phi (6)=6(1-\frac{1}{2})(1-\frac{1}{3})=2$.​

证明：利用 容斥原理，$N$ 的质因子只有 $P_1\sim P_k$ 这 $k$ 个.

$①$ 从 $1\sim N$ 中去掉 $P_1\sim P_k$ 的所有倍数，即 $N-\frac{N}{P_1}-\frac{N}{P_2}-\cdots -\frac{N}{P_k}$，这里面就会多去掉了一部分，比如一些数即是 $P_1$ 的倍数，也是 $P_2$ 的倍数；
$②$ 加上所有 $P_i{P}_j$ 的倍数（$i,j$ 都是 $1\sim k$ 中任选的两个数，$i\ne j$），即 $+\frac{N}{P_1{P}_2}+\frac{N}{P_1{P}_3}+\cdots$，但如果一个数是 $P_1,P_2,P_3$ 的倍数，它会被去掉三次，而它会被 $P_1{P}_2,P_1{P}_3,P_2{P}_3$ 加上三次，效果是没有减没有加，但其实这个数应该被去掉；
$③$ 减去所有$P_i{P}_j{P}_k$ 的倍数，即 $-\frac{N}{P_1{P}_2{P}_3}-\frac{N}{P_1{P}_2{P}_4}-\cdots$，随后依次类推 $+\frac{N}{P_1{P}_2{P}_3{P}_4}+\cdots -\frac{N}{P_1{P}_2{P}_3{P}_4{P}_5}-\cdots$.​

这样，实际上公式的展开就是证明的过程中得到的各项，每项都可以提出来一个 $N$；​

通过验证各项系数来推导正确性：
$\frac{N}{P_1}$​ 的系数为第一个括号内选 $-\frac{1}{P_1}$​，其余括号内选 $1$​，得 $-\frac{N}{P_1}$；​
$\frac{N}{P_1{P}_2}$​ 的系数为第一个括号内选 $-\frac{1}{P_1}$​，第二个括号内选 $-\frac{N}{P_2}$​，其余括号内选 $1$​，可以得到 $+\frac{N}{P_1{P}_2}$；​
以此类推，即可得证.

AcWing 873. 欧拉函数

原题链接

给定 $n$ 个正整数 $a_i$，请你求出每个数的欧拉函数。

输入格式

第一行包含整数 $n$。

接下来 $n$ 行，每行包含一个正整数 $a_i$​。

输出格式

输出共 $n$ 行，每行输出一个正整数 $a_i$ 的欧拉函数。

数据范围

$1\le n\le 100$,
$1\le a_i\le 2\times 10^9$

输入样例：

3
3
6
8

输出样例：

2
2
4

时/空限制： 1s / 64MB
来源： 模板题
算法标签：数学知识 欧拉函数

yxc’s Solution

• 直接使用公式来求.
• 时间复杂度为 $O(\sqrt{n})$.

因为想要得到的是整除的结果，所以将 $(1-\frac{1}{P})$​ 通分为 $(P-1)/P$​，先进行乘法运算会导致溢出，所以代码中先进行除法运算.

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int main(){
   
	int n,a;
	scanf("%d",&n);
	while(n--){
   
		scanf("%d",&a);
		int res=a;
		for(int i=2;i<=a/i;++i)
			if(a%i==0){
   
				res=res/i*(i-1);
				while(a%i==0) a/=i;
			}
		if(a>1) res=res/a*(a-1);
		printf("%d\n",res);
	}
	return 0;
}

AcWing 874. 筛法求欧拉函数

题目链接

给定一个正整数 $n$，求 $1\sim n$ 中每个数的欧拉函数之和。

输入格式

共一行，包含一个整数 $n$。

输出格式

共一行，包含一个整数，表示 $1\sim n$ 中每个数的欧拉函数之和。

数据范围

$1\le n\le 10^6$

输入样例：

6

输出样例：

12

时/空限制： 1s / 64MB
来源： 模板题
算法标签：数学知识 欧拉函数 线性筛法

yxc’s Solution

• 如果每个数都用公式来做，则时间复杂度会达到 $O(n\sqrt{n})$，是不理想的.
• 借用线性筛法求素数的代码，可以达到 $O(n)$ 的时间复杂度求每个数的欧拉函数的效果.

for(int i=2;i<=n;++i){
   
	if(!st[i]) primes[cnt++]=i;
	for(int j=0;primes[j]<=n/i;++j){
   
		st[primes[j]*i]=true;
		if(i%primes[j]==0) break;
}

线性筛法这个代码是很好用的，它能顺便求出来很多东西，不光可以求欧拉函数.

• 根据数据范围 $1\le n\le 10^6$，则欧拉函数的和可能会很大（达到越 $n^2$ 的数量级），所以答案应该用 long long 存储；

• 如果一个数 $P$​ 是质数，则它的欧拉函数应为 $P-1$，即 $1\sim P-1$ 都与其互质；

• 当 i%primes[j]==0 时，$\varphi ($​​i*primes[j]$)$​ 应该为多少？

  • 从公式的角度有 primes[j] 是 i 的一个质因子，所以 $\varphi ($i$)$ 已经乘过 $(1-\frac{1}{primes[j]})$ 了，所以 $\phi ($i$)$ 中形如的 $(1-\frac{1}{P})$ 的式子和 $\phi ($i*primes[j]$)$ 中是一模一样的.

  或者从分解质因数的角度，i*primes[j] 分解质因数后的结果只是比 i 多乘了一个 primes[j] ，故不考虑指数的情况下，i*primes[j] 和 i 的质因子是相同的.

  • 即，假设 $\phi ($i ) = i × ( 1 − 1 P 1 ) × ⋯ × ( 1 − 1 P k ) )=i\times (1-\frac{1}{P_1})\times \cdots \tim