线性筛欧拉函数【模板】

最新推荐文章于 2023-08-18 20:25:24 发布

SSL_kyx

最新推荐文章于 2023-08-18 20:25:24 发布

阅读量130

点赞数

分类专栏： # 质数与同余数论

本文链接：https://blog.csdn.net/keyixi/article/details/117160038

版权

质数与同余同时被 2 个专栏收录

8 篇文章 0 订阅

订阅专栏

数论

8 篇文章 1 订阅

订阅专栏

该博客介绍了如何使用线性筛法来计算欧拉函数。欧拉函数φ(n)表示小于等于n且与n互质的正整数个数。文章详细阐述了欧拉函数的基本性质，如φ(p)=p-1（p为质数）和φ(pk)=pk-pk-1，并给出基于积性函数性质的计算公式。通过线性筛法，可以高效地计算出1到n的欧拉函数值，代码示例展示了具体实现过程。

摘要由CSDN通过智能技术生成

线性筛欧拉函数

欧拉函数

欧拉函数的定义：
若 $\gcd(a,b)=1$ ，则 $a$ 和 $b$ 互质
在自然数1到n中与n互质的个数称为欧拉函数，记作 $φ (n)$
~~下文皆以 $p$ 表示质数 (prime)~~

性质1： $φ (p) = p - 1$

易证明：显然，质数的因数只有1和它本身，且1与任何数互质，所以除了它本身，其它都互质。

性质2： $φ(p^k)=p^k-p^{k-1}$

证明：与性质1类似， $p^k$ 即总数减去它的因数数量 $p^{k-1}$ ，由于 $p^k$
根据容斥定理可以得到 $φ (n) = n * \prod (1 - p) / p$
由于φ是积性函数，所以 $φ(ab)=φ(a)*φ(b)(\gcd(a,b)=1)$
$\sum [d ∣ n] φ (d) = n ， \sum [g c d (n, i) = = 1] i = n * φ (n) / 2$
根据以上性质，可以用各种筛法解决

#include<cstdio>
#include<iostream>
#define Max 100000001
using namespace std;
int phi[Max],n,p[Max],pn; 
bool f[Max];
void Phi()//φ(欧拉函数)
{
	phi[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!f[i])
		{
			p[++pn]=i;
			phi[i]=i-1;//根据性质，i为质数， φ(i)=i-1
		}
		for(int j=1;j<=pn&&i*p[j]<=n;j++)
		{
			f[i*p[j]]=1;
			if(i%p[j]==0)
			{
				//根据通式，φ(i*p[j]) = (i*p[j])*∏(1-p)/p = (i*∏(1-p)/p)*p[j] = φ(i)*p[j]
				phi[i*p[j]]=phi[i]*p[j];
				break;
			}
			else //否则，i和 p[j]互质，因为 p[j]是质数 
			phi[i*p[j]]=phi[i]*(p[j]-1);//由积性函数性质得 φ(i*p[j]) = φ(i)*φ(p[j]) = φ(i)*(p[j]-1)
		}
	}
}
int main()
{
	printf("求1~n的欧拉函数\n输入n\n");
	scanf("%d",&n);
	Phi();
	for(int i=1;i<=n;i++)
	{
		printf("1~%d中与%d互质的个数:%d\n",i,i,phi[i]);
	}
}