老逗的gcd 莫比乌斯反演

老逗的gcd 莫比乌斯反演_第1张图片

这个题一看就是莫比乌斯反演,怎么处理无平方因子数呢?由莫比乌斯函数的性质,不难写出:


简单反演得:

那么令

重点在于处理的这个函数的值,用传统筛的话可以加一个小优化就是底数和倍数都只枚举Mob不为0的数

即设Cur数组保存所有Mob不为0的数

#include
#include
using namespace std;
typedef long long LL;
const LL Mn=10000005;
LL vst[Mn+5],p[Mn+5],Mob[Mn+5],Cur[Mn+5],f[Mn+5];
void LineShaker(){
	    LL i,j;
	    Mob[1]=1;Cur[++Cur[0]]=1;
	    for(i=2;i<=Mn;i++){
			if(!vst[i]){p[++p[0]]=i;Mob[i]=-1;}
			if(Mob[i])Cur[++Cur[0]]=i;
			for(j=1;j<=p[0]&&p[j]*i<=Mn;j++){
				vst[p[j]*i]=1;
				if(i%p[j]==0){Mob[i*p[j]]=0;break;}
				Mob[i*p[j]]=-Mob[i];
		    }
		}
		for(i=1;i<=Cur[0];i++)
			for(j=1;j<=Cur[0]&&Cur[j]*Cur[i]<=Mn;j++)
				f[Cur[i]*Cur[j]]+=Cur[i]*Mob[Cur[i]]*Mob[Cur[i]]*Mob[Cur[j]];
		for(i=1;i<=Mn;i++)f[i]+=f[i-1];
}
LL Cal(LL x,LL y){
	 LL i,pos,sum=0;
	 if(x>y)swap(x,y);
	 for(i=1;i<=x;i=pos+1){
		 pos=min(x/(x/i),y/(y/i));
		 sum+=(f[pos]-f[i-1])*(x/i)*(y/i);
	 }
	 return sum;
}
void solve(){
	    LL T,x,y;
	    scanf("%lld",&T);
	    for(;T--;){
			scanf("%lld%lld",&x,&y);
			printf("%lld\n",Cal(x,y));
		}
}
int main(){
	  LineShaker();
	  solve();
      return 0;
}

可以缓慢地AC,但姿势不优,考虑分离质因子效应,由Mob的性质不难发现:


那么用线性筛筛出每个数最小因子的幂,最小因子则可以利用积性求解

#include
#include
using namespace std;
typedef long long LL;
const LL Mn=10000005;
LL f[Mn+5];
int vst[Mn+5],p[Mn+5],Mob[Mn+5],l[Mn+5],Mi[Mn+5];
void LineShaker(){
	    LL i,j;
	    Mob[1]=1;
	    for(i=2;i<=Mn;i++){
			if(!vst[i]){p[++p[0]]=i;Mob[i]=-1;l[i]=i;Mi[i]=i;}
			for(j=1;j<=p[0]&&p[j]*i<=Mn;j++){
				vst[p[j]*i]=1;Mi[p[j]*i]=p[j];
				if(i%p[j]==0){
					Mob[i*p[j]]=0;l[i*p[j]]=l[i]*p[j];
					break;
				}
				Mob[i*p[j]]=-Mob[i];l[i*p[j]]=p[j];
		    }
		}
		f[1]=1;
		for(i=2;i<=Mn;i++){
			if(i==l[i]){
				if(!vst[i])f[i]=i-1;
				else if(i==Mi[i]*Mi[i])f[i]=-Mi[i];
				else f[i]=0;
			}
			f[i]=f[i/l[i]]*f[l[i]];
		}
		for(i=1;i<=Mn;i++)f[i]+=f[i-1];
}
LL Cal(LL x,LL y){
	 LL i,pos,sum=0;
	 if(x>y)swap(x,y);
	 for(i=1;i<=x;i=pos+1){
		 pos=min(x/(x/i),y/(y/i));
		 sum+=(f[pos]-f[i-1])*(x/i)*(y/i);
	 }
	 return sum;
}
void solve(){
	    LL T,x,y;
	    scanf("%lld",&T);
	    for(;T--;){
			scanf("%lld%lld",&x,&y);
			printf("%lld\n",Cal(x,y));
		}
}
int main(){
	  LineShaker();
	  solve();
      return 0;
}
/*
5 
4 5
6 7
3 6
8 9
986 349
*/

这样理论达到线性预处理,显然已经最优。

你可能感兴趣的:(莫比乌斯反演)