【数论】欧拉函数

欧拉函数

一、定义

对于正整数 $n$，在小于 $n$ 的正整数中与 $n$ 互质的数的数目是欧拉函数，记作 $\phi (n)$。

二、性质

性质一

根据定义，我们可以得出公式：

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

性质二

$$\phi (n)=n\prod _{\begin{array}{c}p|n\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)$$

其中，$\mathbb{P}$ 为质数集。

下面证明该性质：

如果 $p$ 是 $n$ 质因数，那么有 $\frac{1}{p}$ 的数是 $p$ 的倍数， 即有 $\left(1-\frac{1}{p}\right)$ 的数与 $n$ 互质，那么与 $n$ 的质因数均不相等的数的个数就是 $n\left(1-\frac{1}{p_1}\right)\left(1-\frac{1}{p_2}\right)...\left(1-\frac{1}{p_k}\right)$，得证。

性质三

如果 $n$ 是质数，小于等于 $n$ 的数均与 $n$ 互质，那么有 $\phi (n)=n-1$。

性质四

如果 $n$ 与 $m$ 互质，那么有

$$\phi (nm)=\phi (n)\phi (m)$$

下面证明该性质：

$$\phi (n)\phi (m)=n\prod _{\begin{array}{c}p|n\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)m\prod _{\begin{array}{c}p|m\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)=nm\prod _{\begin{array}{c}p|n\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)\prod _{\begin{array}{c}p|m\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)$$

因为 $n$ 和 $m$ 互质，所以 $nm$ 的质因数就是 $n$ 的质因数与 $m$ 的质因数的交集，所以

$$\phi (n)\phi (m)=nm\prod _{\begin{array}{c}p|nm\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)=\phi (nm)$$

性质五

当 $n>2$ 时，$\phi (n)$ 是偶数。

下面证明该性质：

当 $\gcd (n,i)=1$ 时，$\gcd (n,n-i)=1$。

若 $i=n-i$，则 $n=2i$，$n$ 与 $i$ 不互质。

所以 $\phi (n)$ 总为偶数。（$n>2$）

性质六

在区间 $[1,n]$ 中所有与 $n$ 互质的数的和为 $\frac{\phi (n)n}{2}$

下面证明该性质：

根据性质五，如果 $\gcd (n,i)=1$，那么有 $\gcd (n,n-i)=1$，$i$ 和 $n-i$ 的和是 $n$，这样的数对的个数是 $\frac{\phi (n)}{2}$ 个，所以和为 $\frac{\phi (n)n}{2}$。

性质七

$$\phi (n^2)=n\phi (n)$$

下面证明该性质：

因为 $n^2$ 的质因数和 $n$ 的质因数相同，因此

$$\phi (n^2)=n^2\prod _{\begin{array}{c}p|n^2\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)=n^2\prod _{\begin{array}{c}p|n\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)=n\times n\prod _{\begin{array}{c}p|n\\ p\in \mathbb{P}\end{array}}\left(1-\frac{1}{p}\right)=n\phi (n)$$

三、求欧拉函数

快速求欧拉函数

根据性质二，我们可以在 $O(\sqrt{n})$ 的时间复杂度内完成。

代码如下：

void phi(int 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);
	return ans;
}

线性求欧拉函数

根据性质三、性质四和性质七，我们就可以线性求欧拉函数。

• 如果 $n$ 为质数，那么 $\phi (n)=n-1$。
• 如果 $m$ 为 $n$ 的最小平方因子，那么 $\phi (n)=m\phi (\frac{n}{m})$
• 如果 $m$ 为 $n$ 的最小质因数，那么 $\phi (n)=\phi (m)\phi (\frac{n}{m})=(m-1)\phi (\frac{n}{m})$。

以下是一个可以参考的代码：

phi[1] = 1;
for (int i = 2; i <= n; i ++ )
{
	if (!vis[i])
	{
		prime[ ++ cnt] = i;
		phi[i] = i - 1;
	}
	for (int j = 1; prime[j] <= n / i; j ++ )
	{
		vis[i * prime[j]] = 1;
		if (i % prime[j] == 0)
		{
			phi[i * prime[j]] = prime[j] * phi[i];
			break;
		}
		phi[i * prime[j]] = (prime[j] - 1) * phi[i];
	}
}

四、例题

P2158 [SDOI2008] 仪仗队

我们假设坐标的下标从 $0$ 开始。

设一个人的位置是 $(x,y)$ ，若 $gcd(x,y)=k\ne 1$，那么在坐标 $(\frac{x}{k},\frac{y}{k})$ 的人会将这个人遮挡，所以这道题其实就是在求 $\gcd (x,y)=1$ 的个数，公式如下：

$$\sum _{i=0}^{n-1}\sum _{j=0}^{n-1}[\gcd (i,j)=1]$$

通过观察图片可以发现，图片延 $y=x$ 对称，所以原式

$$=\sum _{i=0}^{n-1}\sum _{j=0}^{i-1}[\gcd (i,j)=1]+\sum _{i=0}^{n-1}\sum _{j=i}^i[\gcd (i,j)=1]+\sum _{i=0}^{n-1}\sum _{j=i+1}^{n-1}[\gcd (i,j)=1]$$

$$=2\sum _{i=0}^{n-1}\sum _{j=0}^{i-1}[\gcd (i,j)=1]+1$$

$$=2\sum _{i=1}^{n-1}\phi (i)+1$$

需要注意，当 $n=1$ 时，答案为 $0$。

代码如下：

#include <bits/stdc++.h>

using namespace std;

const int N = 1000005;

int n;
long long prime[N], vis[N], phi[N], cnt;

int main()
{
	scanf("%d", &n);
	if (n == 1)
	{
	    puts("0");
	    return 0;
	}
	phi[1] = 1;
	for (int i = 2; i <= n; i ++ )
	{
		if (!vis[i])
		{
			prime[ ++ cnt] = i;
			phi[i] = i - 1;
		}
		for (int j = 1; prime[j] <= n / i; j ++ )
		{
			vis[i * prime[j]] = 1;
			if (i % prime[j] == 0)
			{
				phi[i * prime[j]] = prime[j] * phi[i];
				break;
			}
			phi[i * prime[j]] = (prime[j] - 1) * phi[i];
		}
	}
	long long ans = 0;
	for (int i = 1; i < n; i ++ )
		ans += phi[i];
	printf("%d\n", 2 * ans + 1);
	
    return 0;
}

P2568 GCD

枚举每一个质数 $p$，求 $\gcd (x,y)=p$ 的情况，此时一定有 $p|x$，$p|y$，那么 $\gcd (x,y)=p$ 等价于 $\gcd (\frac{x}{p},\frac{y}{p})=1$，于是就有公式：

$$\sum _{p\in \mathbb{P}}\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)=p]$$

$$=\sum _{p\in \mathbb{P}}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{p}\rfloor }[\gcd (i,j)=1]$$

$$=\sum _{p\in \mathbb{P}}\left(\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{i-1}[\gcd (i,j)=1]+\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=i}^i[\gcd (i,j)=1]+\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=i+1}^{\lfloor \frac{n}{p}\rfloor }[\gcd (i,j)=1]\right)$$

$$=\sum _{p\in \mathbb{P}}\left(2\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{i-1}[\gcd (i,j)=1]+1\right)$$

$$=\sum _{p\in \mathbb{P}}\left[2\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }(\phi (i)-1)+1\right]$$

$$=\sum _{p\in \mathbb{P}}\left(2\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\phi (i)-1\right)$$

对于倒数第 $2$ 个式子，由于与 $1$ 互质的数的个数 $0$ 但 $\phi (1)=1$，需要减 $1$。

#include <bits/stdc++.h>

using namespace std;

const int N = 10000005;

int n;
long long prime[N], vis[N], phi[N], Phi[N], cnt;

int main()
{
	scanf("%d",n);
	phi[1] = 1;
	for(int i = 2; i <= n; i ++ )
	{
		if(!vis[i])
		{
			prime[ ++ cnt] = i;
			phi[i] = i - 1;
		}
		for(int j = 1; prime[j] <= n / i; j ++ )
		{
			vis[i * prime[j]] = 1;
			if(i % prime[j] == 0)
			{
				phi[i * prime[j]] = prime[j] * phi[i];
				break;
			}
			phi[i * prime[j]] = (prime[j] - 1) * phi[i];
		}
	}
	long long ans = 0;
	for(int i = 1; i < n; i ++ )
	{
		Phi[i] = Phi[i - 1] + phi[i];
	}
	for(int i = 1; i <= cnt; i ++ )
	{
		ans += 2 * Phi[n / prime[i]] - 1;
	}
	printf("%d\n", ans);
    return 0;
}