【学习笔记】欧拉函数

0.概述

一个函数。只有毒瘤出题人会用到它。

1.定义

$$\phi (n)=\sum _{x=1}^n[\gcd (x,n)=1]$$

即，所有小于等于 $n$ 的正整数，与 $n$ 互质的有多少个。顺带一说，互质也可写为 $a\perp b$ 。

2.计算

2.1.积性函数

如果 $a$ 和 $b$ 互质，那么
$$\phi (ab)=\phi (a)\cdot \phi (b)$$

证明：发现 $\gcd (x,ab)=1$ 的充要条件是 $\gcd (x,a)=1\wedge \gcd (x,b)=1$ 。

那么设 $x=kb+t(0⩽k<a,0⩽t<b)$，由 $x\perp b\iff t\perp b$ 可知，不同的 $t$ 的数量是 $\phi (b)$ 个。

另一方面，对于某个 $t$，肯定有 { k b + t   m o d   a    ∣    k ∈ [ 0 , a ) } = { 0 , 1 , 2 , … , a − 1 } \{kb+t\bmod a\;|\;k\in[0,a)\}=\{0,1,2,\dots,a{\rm-}1\} {kb+tmoda∣k∈[0,a)}={0,1,2,…,a−1}，所以可行的 $k$ 只有 $\phi (a)$ 个。二者相乘就是 $\phi (ab)=\phi (a)\cdot \phi (b)$ 。

2.2.一般的计算

如果 $p$ 是一个质数，$k$ 是正整数，那么
$$\phi (p^k)=p^{k-1}(p-1)$$

证明：$x\perp p^k\iff x\perp p\iff x=kp(1⩽k⩽p^{k-1})$，故 $\phi (p^k)=p^k-p^{k-1}$，证毕。

再结合 积性函数，如果 $n$ 的质因数分解是 $n=p_1^{t_1}{p}_2^{t_2}{p}_3^{t_3}\dots p_m^{t_m}$ ，那么
$$\phi (n)=\prod _{x=1}^m{p}_x^{t_x-1}(p_x-1)$$

2.3.特殊的计算

其实就是 一般的计算 变了形。
$$\phi (n)=n\prod _{p\mid n}\left(1-\frac{1}{p}\right)$$

也可以理解为容斥，即考虑是否去掉某个质因数的倍数。注意这里隐式规定 $p\in \text{Prime}$ 。

3.性质

3.0.积性函数

前面已经讲过 积性函数，这里就略过。

3.1.与 $n$ 的关系

还有一个很美妙的结论（使得杜教筛在 $\phi$ 上可以大显身手）：
$$\sum _{d|n}\phi (d)=n$$

快速证明：由容斥可知 $id\ast \mu =\phi$，两边同乘 $I$ 得 $id=\phi \ast I$，证毕。其中 $\ast$ 为狄利克雷卷积。

初等证明：将 { 1 , 2 , 3 , … , n } \{1,2,3,\dots,n\} {1,2,3,…,n} 按如下规则划分： S d = { x ∣ gcd ⁡ ( x , n ) = d } S_d=\{x\mid\gcd(x,n)=d\} Sd​={x∣gcd(x,n)=d} 。显然 S 1 ∪ S 2 ∪ ⋯ ∪ S n = { 1 , 2 , 3 , … , n } S_1\cup S_2\cup\cdots\cup S_n=\{1,2,3,\dots,n\} S1​∪S2​∪⋯∪Sn​={1,2,3,…,n} 。

对于任意 $x\in S_d$，可以写成 $x=kd(1⩽k⩽\frac{n}{d}$ 但不一定取遍 $)$ 。

由定义，$\gcd (x,n)=\gcd (kd,n)=d\Rightarrow \gcd (k,\frac{n}{d})=1$ 。

这样一来，$|S_d|=\phi (\frac{n}{d})$ 。于是 $n={\sum }_{d|n}|S_d|={\sum }_{d|n}\phi (\frac{n}{d})$，得证。

3.2.与莫比乌斯挂钩

$$\sum _{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$$

比较显然且用处不大，此处从略。

3.3.“偶性”

$$2\mid \phi (n)(n⩾3)$$

证明：如果 $n$ 中含有任意一个奇素因数 $p$ ，那么根据 一般的计算，可知 $2\mid (p-1)\mid \phi (n)$ 。

而不包含奇素因数的数只有 $1$ 和 $2$ 。

3.4.“收敛” 速度快

接下来的这一点，是很多题复杂度的保障：如果 $2\mid n$ ，那么 $\phi (n)⩽\frac{n}{2}$ 。

证明：$n$ 中含有 $p=2$ 这一质因子。根据 特殊的计算，$\phi (n)⩽n(1-\frac{1}{p})=\frac{n}{2}$ 。

再结合 偶性 可知，每次令 $n_i=\phi (n_{i-1})$，那么 $i⩾\lceil {\log }_{2}n\rceil$ 即可让 $n_i=1$ 。

3.5.欧拉降幂

我另开了一篇博客写这个……

4.代码实现

4.1.求单点值

直接用 一般的计算，暴搜质因数分解，复杂度 $\mathcal{O}(\sqrt{n})$ 。

int getPhi(int n){
	int res = 1;
	for(int x=2; x<=n/x; ++x){
		if(n%x) continue;
		n /= x, res *= (x-1);
		while(n%x == 0)
			n /= x, res *= x;
	}
	if(n != 1) res *= (n-1);
	return res;
}

我个人更喜欢写法二，感觉常数小一点。

long long __phi(long long n){
	long long res = n;
	for(int i=2; i<=n/i; ++i)
		if(n%i == 0){
			res -= res/i, n /= i;
			while(n%i == 0) n /= i;
		}
	if(n != 1) res -= res/n;
	return res;
}

4.2.奇怪法求单点值

int getPhi(int n){
	vector<int> ys; // 质因数 
	while(n != 1){
		d = findDivisor(n); // 找到一个质因数 
		ys.push_back(d);
		while(n%d == 0) n /= d;
	}
	sort(ys.begin(),ys.end());
	n = ys.size(); int res = n;
	for(int S=1; S<(1<<n); ++S){
		int sgn = 1, cnt = n;
// 是质因子集合S中每一个质因子的倍数的个数 
		for(int i=0; i<n; ++i)
			if(S>>i&1){
				sgn = -sgn;
				cnt /= ys[i];
			}
		res += sgn*cnt;
	}
	return res;
}

好处是什么呢？——可以求 ${\sum }_{x=1}^n[\gcd (x,m)=1]$ 。只需要稍稍作一点改动。

4.3.奇怪法

int getCoprime(int n,int m){
	vector<int> ys; // 质因数 
	while(m != 1){
		d = findDivisor(m); // 找到一个质因数 
		ys.push_back(d);
		while(m%d == 0) m /= d;
	}
	sort(ys.begin(),ys.end());
	m = ys.size(); int res = n;
	for(int S=1; S<(1<<m); ++S){
		int sgn = 1, cnt = n;
// 是质因子集合S中每一个质因子的倍数的个数 
		for(int i=0; i<n; ++i)
			if(S>>i&1)
				sgn = -sgn, cnt /= ys[i];
		res += sgn*cnt;
	}
	return res;
}

4.4.所有因数的单点值

能够做到 $\sqrt{n}$ 的复杂度——也就是说，均摊 $\mathcal{O}(1)$ 。

首先线性筛求出 $x\in [1,\sqrt{n}]$ 的 $\phi (x)$，然后 $\mathcal{O}(\sqrt{n})$ 的求出 $\phi (n)$ 。接下来要做的，就是研究 $\phi (x),\phi (\frac{n}{x}),\phi (n)$ 的关系。

观察 特殊的计算 中的右半部分，也就是 ${\prod }_{p|n}(1-\frac{1}{p})$ ，将其记为 $\lambda (n)$ 。于是就有 $n\lambda (n)=\phi (n)$ 。

考虑这样一个恒等式：记 $d=\gcd (x,\frac{n}{x})$ ，那么
$$\lambda (x)\cdot \lambda \left(\frac{n}{x}\right)=\lambda (n)\lambda (d)$$

证明：如果 $p\mid d$，在两边都会被计算两次；如果 $p\nmid n$，那么两边都一次也不会计算；其他的情况，两边都只计算一次。证毕。

把两边同时乘一个 $n$ 可以得到：
$$\begin{gathered} x\lambda (x)\cdot \frac{n}{x}\cdot \lambda \left(\frac{n}{x}\right)=n\lambda (n)\lambda (d) \\ \Rightarrow \phi (x)\phi \left(\frac{n}{x}\right)=\phi (n)\frac{\phi (d)}{d} \\ \Rightarrow \phi \left(\frac{n}{x}\right)=\frac{\phi (n)\phi (d)}{d\phi (x)} \end{gathered}$$

更形式化地，
$$\phi (a)\phi (b)=\phi (ab)\frac{\phi [\gcd (a,b)]}{\gcd (a,b)}$$

可以发现，这是自洽的——在 $\gcd (a,b)=1$ 的情况下，自然满足 积性函数 咯。

求 $\gcd (x,\frac{n}{x})$ 怎样避免 $\mathcal{O}(\log n)$ 辗转相除呢？一种方法是，由于 $x⩽\sqrt{n}$，我们可以用快速 $\gcd$ 算法。这样不费脑子。

另一种方法是，求出 $x\in [1,\sqrt{n}]$ 的质因数表示。具体来说，存储最小质因数与对应的次数，以及除以它得到的值（即记录 $x=p^{k}v$ 其中 $p$ 为最小质因子且 $p\perp v$ ）。我们已经用 $\mathcal{O}(\sqrt{n})$ 搜出了 $n$ 的质因数分解（用 $\text{unordered map}$ 桶计数就行，因为不同 $p$ 的数量甚至低于 $\log n$，可视为常数），那么 $\gcd (x,\frac{n}{x})$ 中 $p$ 的指数易得，最后乘上 $\gcd (v,\frac{n}{v})(v=\frac{x}{p^k})$ 即可。

5.例题

$\text{Naive Great Common Divisor}$

电磁看题。应该算比较常规的方法。

$\text{The Luckiest Number}$

点此跳转。用欧拉公式解决关于 $x$ 的方程 $a^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$ 。

$\text{POAIBE}$

电磁跳转。$\text{p.s.}$ 这个名字是缩写。

杜教筛求 $\phi$ 的前缀和

我们已经知道了 $\phi \ast I=id$ ，直接用即可。

放公式：令 $S_n={\sum }_{i=1}^n\phi (i)$ ，一定有
$$S_n=\frac{n(n+1)}{2}-\sum _{d=1}^{n}S(\lfloor \frac{n}{d}\rfloor )$$

map<int,long long> wk; /* WK(freopen) or2 */
long long djs(int n){
	if(n < SqrtN)
		return phi[n]; /* 线性筛中得到的结果 */
	if(wk.count(n)) return wk[n]; /* 记忆化 */
	long long res = 1ll*n*(n+1)>>1;
	for(int l=1,r; l<=n; l=r){
		r = n/(n/l)+1; /* 整除分块 */
		res -= (r-l)*djs(n/l);
	}
	return wk[n] = res;
}

复杂度是 $\mathcal{O}(n^{\frac{2}{3}})$ 的。计算过程我放在别的博客的不显眼之处了，不知道也罢。