欧拉函数

欧拉函数 F(n)    [1,n]中和n互质的数的个数 

 求一个数的欧拉函数

//分解质因数后 N=p1^a1*p2^a2*.....*pk^ak
//则F(N) = N*(1-1/p1)*(1-1/p2)*...*(1-1/pk) 
int phi(int n)
{
	int res=n;
	for(int i=2;i*i<=n;i++){
		if(n%i==0){
			res=res/i*(i-1);
			while(n%i==0) n/=i;
		}
	}
	if(n>1) res=res/n*(n-1);
	return res;
}

求[1,n]的欧拉函数（线性筛）

int phi[N];
bool vis[N];
int prime[N];
int cnt;
//筛法求1~n的欧拉 
void get_euler(int n){
	//和1互质的只有1 
	phi[1]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			//与质数N互质的数个数为N-1 
			prime[cnt++]=i;
			phi[i]=i-1;
		}	
		//线性筛 
		for(int j=0;prime[j]*i<=n;j++){
			int t=prime[j]*i;
			vis[t]=true;
			//1.prime[j]是i的质因子 
			//  prime[j]*i的质因子都出现在i的质因子中
			//  prime[j]*i的欧拉就等于i的欧拉*prime[j]
			//2.prime[j]不是i的质因子 
			//prime[j]*i的欧拉等于i的欧拉乘prime[j]*(1-1/prime[j]) 
			if(i%prime[j]==0){
				phi[t]=phi[i]*prime[j];
				break;
			}else{
				phi[t]=phi[i]*(prime[j]-1);
			}
		}
	}
}

---

可见的点

201. 可见的点 - AcWing题库

        考虑左上半边（y>=x），满足gcd(x,y)==1 即x y 互质的点为所求，所以y的欧拉即为在左上半边可以匹配的x的数量，两边对称*2，再加上特殊的边界点即可

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=1e3+5;
int phi[N];
bool vis[N];
int prime[N];
int cnt;
void get_euler(int n){
	phi[1]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			prime[cnt++]=i;
			phi[i]=i-1;
		}	
		for(int j=0;prime[j]*i<=n;j++){
			int t=prime[j]*i;
			vis[t]=true;
			if(i%prime[j]==0){
				phi[t]=phi[i]*prime[j];
				break;
			}else{
				phi[t]=phi[i]*(prime[j]-1);
			}
		}
	}
}
int main(){
	int t;cin>>t;
	get_euler(N-1);
	for(int j=1;j<=t;j++){
		ll res=0;
		int n;cin>>n;
		if(n!=0){
			for(int i=2;i<=n;i++){
				res+=phi[i];
			}
			res*=2;
			res+=3;
		}
		cout<<j<<" "<<n<<" "<<res<<endl;
	}
} 

最大公约数

220. 最大公约数 - AcWing题库

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=1e7+5;
int phi[N];
bool vis[N];
int prime[N];
int cnt;
ll phi_sum[N];
void get_euler(int n){
	phi[1]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			prime[cnt++]=i;
			phi[i]=i-1;
		}	
		for(int j=0;prime[j]*i<=n;j++){
			int t=prime[j]*i;
			vis[t]=true;
			if(i%prime[j]==0){
				phi[t]=phi[i]*prime[j];
				break;
			}else{
				phi[t]=phi[i]*(prime[j]-1);
			}
		}
	}
	for(int i=1;i<=n;i++){
		phi_sum[i]=phi_sum[i-1]+phi[i];
	}
}
int main(){
	int n;cin>>n;
	get_euler(n);
	ll res=0;
	for(int i=0;i<cnt;i++){
		res+=2*(phi_sum[n/prime[i]]-phi_sum[1])+1;
	}
	cout<<res;
}