欧拉函数常见题型

<<<书接上回  带来一些简单欧拉函数运用 (本文所用性质均在上文中提及)

欧拉函数题主要在于怎样将问题转化为与欧拉函数有关的问题

基础题

求约数个数

问题 ：f(i) 表示约数个数，求， 

由性质，f(i) 为积性函数，所以考虑线性筛

令 i = ab ，a 为 i 最小质因数的最高次幂，b = i / a , 这样保证了 gcd(a,b)=1 , 同时 i 的最小质因数的最高次幂 hp[i] 也可以用线性筛得到

具体细节见代码注释

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+100;
int prime[N],tot,n;
bool ip[N];
int f[N];//约数个数
int hp[N];//最小质因数最高次幂
int mi[N];//最小质因数最高次方，24=2^3*3,mi[24]=3
int main(){
	cin>>n;
	f[1]=1;
	for(int i=2;i<=n;i++) {
		if(!ip[i]){
			prime[++tot]=i;
			hp[i]=i;
			f[i]=2;
			mi[i]=1;//质数只有自己一个质因数，最高次幂=1
		}
		for(int j=1;j<=tot&&i*prime[j]<=n;j++){
			int x=i*prime[j];
			ip[x]=1;
			if(i%prime[j]!=0) {
				hp[x]=prime[j];
				mi[x]=1;
				int a=hp[x],b=x/a;
				f[x]=f[a]*f[b];
			}
			else{
				hp[x]=hp[i]*hp[prime[j]];
				mi[x]=mi[i]+mi[prime[j]];
				int a=hp[x],b=x/a;
				if(b==1) f[x]=mi[x]+1;//如果x=p^k,f[x]=k+1;
				else f[x]=f[a]*f[b];
				break;
			}
		}
	}
	long long ans=0;
	for(int i=1;i<=n;i++) ans+=f[i];
	cout<<ans;
	
	return 0;
}

最大公约数的和 洛谷P1390

问题 ： 

正整数N，求1至N这N个数中每任意两个不同的数的最大公约数的和。1 <= N <= 2000000

解决：

考虑枚举最大公约数 d  ∈ [1,n] ，则需要知道有多少个数对 (x,y) , 满足 gcd(x,y) = d 。

得到 x=a*d ，y=b*d ， gcd(a,b) = 1

不妨设 a>b , 考虑枚举 a ，易得 a 范围 ∈ [2,n/d] ，b 的个数即为

所以总方案数 =  ， 维护一个 phi[x] 的前缀和即可

代码不多赘述

GCD(x,y)是质数 洛谷P2568

问题：

给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对。

先说一个发现：关于欧拉函数与最大公约数的题，很多都是通过枚举最大公约数解决的

解决：枚举 质数 p ，令 x = ap , y = bp  ， 得到 gcd(a,b)=1 。

同样的，不妨设 a<b ，则对于 b ∈ [1 , N/p] ，满足条件的 a 个数为 

由于 (a,b) , (b,a) 不同，所以要将上面柿子 *2 ，但 i=1 的时候不 *2

s[i]表示phi[i] 前缀和

GCD(X,N)>=M

问题：

T组测试数据 , 1<=m<=n <= 1e9

给定整数N和M，问题是求有多少个整数X，满足1 <= X <= N且gcd（X，N）> = M。

解决：

同样的，枚举 n 的因数 d 作为 gcd  (因数个数小于 sqrt(n) ) ， 就是求 

注意由于 n 太大不能线性筛预处理，要每次根据计算式得出欧拉函数值

好题：

 一、非互质数求和

问题：

如果一个正整数x小于N，而且x与N不互质，那么整数x称为“好数”。

给定一个正整数N，求所有的“好数”的和。

答案 %1e9+7

解决：

正难则反。考虑非好数总和

易发现非好数个数 为 ，但要求的是数的和。

怎么办？使用 欧几里得定理  if gcd(n,i)==1 ，则有 gcd(n,n-i)==1  ，这意味着将非好数两两一组，每组的和为 n。

所以非好数的和：      ，  再用总和减去即可。 tips：注意模数处理

代码：

#include<bits/stdc++.h>
using namespace std;
#define int long long
int t,n,m;
int mod=1000000007;
int ksm(int a,int b){
	int ans=1,base=a;
	while(b){
		if(b&1) ans*=base;
		base*=base;
		ans%=mod;
		base%=mod;
		b>>=1;
	}
	return ans%mod;
}
int ny(int x){
	return ksm(x,mod-2);
}
int eul(int x){
	if(x==1) return 1;
	int ans=x;
	for(int i=2;i*i<=x;i++) {
		if(x%i==0) {
			ans*=(i-1);
			ans%=mod;
			ans*=ny(i);
			ans%=mod;
			while(x%i==0) x/=i;
		}
	}
	if(x>1) ans=ans*(x-1)%mod*ny(x)%mod;
	return ans;
}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>t;
	while(t!=0){
		cout<<((t*(t-1)%mod*ny(2)%mod)-(t*eul(t)%mod*ny(2)%mod)+mod)%mod<<endl;
		cin>>t;
	}
	
	
	return 0;
}

二、USB的数学题

题意：

t 组测试数据 ， 每次输入 n ( n 总共不超过 1e7 ) 。每次求 

解决：

枚举 d = gcd（i,n）, 算出有多少个 i 满足 gcd(i,n)=d ，将 i 的个数 × n/d 贡献进ans

所以本质上是求 ， 等价于 

易发现上式为积性函数，所以用线性筛求出即可

代码：

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e7+10;
int n,t,x,prime[N],m,ans,g[N],highPow[N],mi[N],w[2][N];
bool b[N];
signed main()
{
	scanf("%lld",&t);
	n=N-10;
	for(int i=2;i<=n;i++)
	{
		if(!b[i])
		{
			g[i]=i*(i-1)+1;
			highPow[i]=i;
			mi[i]=1;
			prime[++m]=i;
		}
		for(int j=1;j<=m&&i*prime[j]<=n;j++)
		{
			int x=i*prime[j]; 
			b[x]=true; 
			if(i%prime[j]) 
			{
				int A=i,B=prime[j];
				g[x]=g[A]*g[B];
				highPow[x]=B;
				mi[x]=1;
			}
			else
			{
				highPow[x]=highPow[i]*prime[j];
				mi[x]=mi[i]+1;
				int A=highPow[x];
				int B=i*prime[j]/A;
				if(B==1)
				{
					w[0][0]=w[1][0]=1;
					w[0][1]=prime[j];
					w[1][1]=prime[j]-1;
					for(int k=2;k<=mi[x];k++)
					{
						w[0][k]=w[0][k-1]*prime[j];
						w[1][k]=w[1][k-1]*prime[j];
					}
					for(int k=0;k<=mi[x];k++)
					{
						g[x]+=w[0][k]*w[1][k];
					}
				}
				else
				{
					g[x]=g[A]*g[B];
				}
				break; 
			}
		}
	}
	while(t--)
	{
		scanf("%lld",&x);
		if(x!=1)printf("%lld\n",g[x]); 
		else printf("1\n");
	}
	return 0;
}

————撒花————