【详解】欧拉函数

互质与欧拉函数

互质

任意自然数 $a,b$，若 $gcd(a,b)=1$，则称 $a,b$ 互质。

对于三个数或更多个数，将 $gcd(a,b,c)=1$ 的情况称为 $a,b,c$ 互质。把 $gcd(a,b)=gcd(b,c)=gcd(a,c)=1$ 称为两两互质，两者有明显的区别，显然，后者要求更高。

欧拉函数

$1\sim N$ 中与 $N$ 互质的个数被称为欧拉函数，记为 $\phi (N)$ 。

在算术基本定理中，$N=p_1^{c_1}\times p_2^{c_2}\times ...\times p_m^{c_m}$，$p_1,p_2...p_m$ 为质数，则：

$\phi (N)=N\times \frac{p_1-1}{p_1}\times \frac{p_2-1}{p_2}\times ...\times \frac{p_m-1}{p_m}$

证明：

若 $p$ 是 $N$ 的质因子，则 $1\sim N$ 中 $p$ 的倍数有：$p,2p,3p...(\frac{N}{p})\times p$，共 $\frac{N}{p}$ (整除)个。同理，若 $q$ 也是 $N$ 的质因子，则 $1\sim N$ 中 $q$ 的倍数有：$q,2q,3q...(\frac{N}{q})\times q$，共 $\frac{N}{q}$ (整除)个。同时 $p\times q$ 的倍数即被 $p$ 排除了，又被 $q$ 排除了，所有需要加回来。因此，$1\sim N$ 中不与 $N$ 含有共同质因子 $p$ 和 $q$ 的个数为：

$N-\frac{N}{q}-\frac{N}{p}+\frac{N}{pq}=N\times (1-\frac{1}{p}-\frac{1}{q}+\frac{1}{pq})=N\times (1-\frac{1}{p})\times (1-\frac{1}{q})$

这种思想被称为容斥思想（后续会学习）。同理，可以得到 $1\sim N$ 中不与 $N$ 含有共同因子的数的个数，也就是与 $N$ 互质的个数。

根据欧拉函数的计算公式，就可以求解出欧拉函数。

【参考程序】

//求解n的欧拉函数
int phi(int n)
{
	int ans=n;
	for(int i=2;i<=sqrt(n);i++)
	{
		if(n%i==0)
		{
			ans=ans/i*(i-1);
			while(n%i==0) n/=i;
		}
	}
	if(n>1) ans=ans/n*(n-1);
	return ans;
}

时间复杂度 $O(\sqrt{N})$ 。

性质

性质 $1$：$n>1$，$1\sim n$ 中与 $n$ 互质的数的和为 $\frac{n\times \phi (n)}{2}$ 。

性质 $2$： 若 $a,b$ 互质，则 $\phi (ab)=\phi (a)\times \phi (b)$。

证明：

因为 $gcd(n,x)=gcd(n,n-x)$ (九章算术.更相减损术)，$n$ 和 $x、n-x$ 不互质，因此与 $n$ 不互质的数是成对出现的，和为 $n$ ，平均值为 $\frac{n}{2}$，则与 $n$ 互质的数的平均值为 $\frac{n}{2}$ ，$1\sim n$ 中与 $n$ 互质的数的和为 $\frac{n\times \phi (n)}{2}$ ，性质 $1$ 成立。

根据欧拉函数的计算公式，对 $a,b$ 分解质因数，直接可以得到性质 $2$ 。

积性函数

如果 $a,b$ 互质，有 $f(a,b)=f(a)\times f(b)$，那么称函数 $f$ 为积性函数。

性质

性质 $3$： 若 $f$ 为积性函数，且在算术定理中，$n={\prod }_{i=1}^m{p}_i^{c_i}$ ，则 $f(n)={\prod }_{i=1}^{m}f(p_i^{c_i})$

性质 $4$： 设 $p$ 是质数，若 $p|n$ 且 $p^2|n$，则 $\phi (n)=\phi (\frac{n}{p})\times p$

性质 $5$： 设 $p$ 是质数，若 $p|n$ 且 $p^2$ 不能整除 $n$ ，则 $\phi (n)=\phi (\frac{n}{p})\times (p-1)$

证明： 略，可参考《算法竞赛进阶指南》$0x32$ 约数。

例 1：「SDOI2008」仪仗队

【题目描述】

作为体育委员， $C$ 君负责这次运动会仪仗队的训练。仪仗队是由学生组成的 $N\ast N$的方阵，为了保证队伍在行进中整齐划一， $C$ 君会跟在仪仗队的左后方，根据其视线所及的学生人数来判断队伍是否整齐(如下图)。

现在， $C$ 君希望你告诉他队伍整齐时能看到的学生人数。

【输入】

共一个数 $N$ 。

【输出】

共一个数，即 $C$ 君应看到的学生人数。

【样例输入】

4

【样例输出】

9

【数据范围】

对于 $100\%$ 的数据，$1\le N\le 40000$ 。

【算法分析】

设左下角为 $(0,0)$ ，右上角为$(N-1,N-1)$ 。

对于$(x,y)$ ，和 $(0,0)$ 组成的线的斜率为 $\frac{y}{x}$ ，若 $x,y$ 存在最大公约数 $p(p>1)$，那么$(\frac{x}{p},\frac{y}{p})$ 和 $(0,0)$ 斜率也为 $\frac{y}{x}$，$(\frac{x}{p},\frac{y}{p})$和 $(x,y)$、$(0,0)$ 在同一条线上，且在 $(x,y)$ 之前，$(0,0)$ 就只能看到 $(\frac{x}{p},\frac{y}{p})$。

因此，除了 $(1,0),(0,1),(1,1)$，$(x,y) $ 能被看到必须满足 $gcd(x,y)=1，1\le x，y\le N-1，x\ne y$ 。

同时，能看到的人数关于 $(0,0)$ 和 $(N-1,N-1)$ 的直线对称，因此可以只需要考虑一半即可，即 $1\le x<y\le N-1$ 。对于任意 $y(2\le y\le N-1)$ 有多少 $x$ 满足 $gcd(x,y)=1$，且 $x<y$ 。这样的 $x$ 个数就是 $\phi (y)$ 。

因此，最后的答案：$ans=3+2\times (\phi (2)+\phi (3)+...+\phi (N-1))$ 。

求 $2\sim (N-1)$ 的欧拉函数可以使用朴素做法，逐个求解，时间复杂度为 $O(N\sqrt{N})$ 。

利用埃氏筛法，按照欧拉函数的公式，在 $O(N\log N)$ 的时间求出，

int phi[100010];
void euler(int n)
{
	for(int i=2;i<=n;i++)	phi[i]=i;
	for(int i=2;i<=n;i++)
	{
		if(phi[i]==i)				//质数
		{
			for(int j=i;j<=n;j+=i)	//i，i+i,i+2i...都包含i的质因子
				phi[j]=phi[j]/i*(i-1);
		}
	}
}

线性筛法中，每个合数 $n$ 只会被它最小的质因子 $p$ 筛一次，可以利用性质 $4$ 和性质 $5$ ，用 $\phi (\frac{n}{p})$ 递推 $\phi (n)$ ，事件复杂度为 $O(N)$ 。

方法一：朴素算法

#include <bits/stdc++.h>
using namespace std;
int phi(int n)  //求解n的欧拉函数
{
    int ans=n;
    for(int i=2;i<=sqrt(n);i++)
    {
        if(n%i==0)
        {
            ans=ans/i*(i-1);
            while(n%i==0) n/=i;
        }
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}
 
int main()
{
    int n,ans=0;
    cin>>n;
    if(n==1) {cout<<0<<endl;return 0;}
    for(int i=2;i<=n-1;i++)
        ans+=phi(i);
    cout<<3+2*ans<<endl;
    return 0;
} 

方法二：线性筛法

#include <bits/stdc++.h>
using namespace std;
const int N=40010;
int  v[N],prime[N],phi[N];
void euler(int n)
{
	memset(v,0,sizeof(v));	//保存最小质因子 
	int k=0;				//质数个数 
	for(int i=2;i<=n;i++)
	{
		if(v[i]==0)
		{
			v[i]=i;
			prime[++k]=i;
			phi[i]=i-1;		//质数i的欧拉函数为i-1 
		}
		//给当前的数i乘以更小的质因子
		for(int j=1;j<=k;j++)
		{
			if(prime[j]>v[i]||prime[j]>n/i) break;
			v[i*prime[j]]=prime[j];					//更新i*prime[j]的最小质因子
			phi[i*prime[j]]=phi[i]*(i%prime[j]?prime[j]-1:prime[j]);	//使用性质4和性质5
		} 
	 } 
}
int main()
{
	int n,ans=0;
	cin>>n;
	if(n==1) {cout<<0<<endl;return 0;}
	euler(n);
	for(int i=2;i<=n-1;i++)
		ans+=phi[i];
	cout<<3+2*ans<<endl; 
    return 0;
}