欧拉函数试除法与筛法求解

欧拉函数的定义
1∼N 中与 N 互质的数的个数被称为欧拉函数，记为 ϕ(N)。
若在算数基本定理中，N=p^a1₁p^a2₂…p^am_m，则：
ϕ(N) = N×（p1−1）/p1×（p2−1）/p2×…×（p_m−1）/pm
与n互质的数的和为：n*ϕ(n)/2
证明如下：

则根据该公式可得求一个数的欧拉函数的代码为：

试除法求解
即在质因数分解的过程中将欧拉函数求出，复杂度为$\sqrt{n}$

#include <bits/stdc++.h>
using namespace std;
int main(){
		int k,ans;
		cin>>k;
		ans=k;
		for(int i=2;i<=k/i;i++){
			if(k%i==0){
				while(k%i==0)k=k/i;
				ans=ans*(i-1)/i;
			}
		}
		if(k>1)ans=ans*(k-1)/k;
		cout<<ans<<endl;
	return 0;
} 

线性筛法
代码解释

代码如下：

#include<bits/stdc++.h>
using namespace std;
int n,b[3000],phi[30000],pre[30000],t;
int main(){
	cin>>n;
	for(int i=2;i<n;i++){
		if(!b[i]){
			pre[++t]=i;
			phi[i]=i-1;
		}
		for(int j=1;j<=t;j++){
			if(pre[j]*i>n)break;
			b[i*pre[j]]=1;
			if(i%pre[j]==0){
				phi[i*pre[j]]=phi[i]*pre[j];
				//pre[j]为i的因子，故phi[i*pre[j]]=phi[i]*pre[j],因子不变，变得仅仅是N.
				break;
			}
			phi[i*pre[j]]=phi[i]/pre[j]*(1-pre[j]);
			//因为i和pre[j]互质，所以phi[i*pre[j]]=phi[i]*phi[pre[j]]=phi[i]*(pre[j]-1)
		}
	}
	cout<<phi[n]<<endl;
}