莫比乌斯反演+线性筛--bzoj3309: DZY Loves Math

传送门

某个很厉害的老师讲的题quq

很值得一提的是：

线性筛不止可以筛积性函数，它的本质是用最小质因子来完成不重复的线性筛出函数

因此满足和最小质因子有关的即使不是积性函数的函数都可以用它来筛

---

然后说这道题

首先可以把式子写出来，然后运用各种推莫比乌斯的套路推推推

推到Σ(k)[a/k][b/k]Σ(d|k)f(d)μ(k/d) 的时候发现后面那个就是反演式子

可以设g(n)=f*μ(n) 乘号表示卷积，然后我们就只要考虑g的取值就好了

但此时陷入了僵局，虽然发现了只有当n中所有质因子的次幂都相等的时候g才不为0，并且g的取值就是(-1)^(x+1)

（其中x是n的质因子个数）

但是光想着如果在线性筛中筛出每个数的质因子次幂，这显然是不可做的

其实就只要筛出miu函数就好了，然后我们可以枚举每个miu不为0的数

他们的质因子次幂一定都是相等的，然后不断累乘就可以算出所有g不为0的数了

 

然后就是老套路，前缀和+数论分块

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 10000001
#define LL long long
using namespace std;
int t,n,m,cnt,pri[700000],g[maxn],miu[maxn];
bool vis[maxn];

inline int rd(){
	int x=0,f=1;char c=' ';
	while(c<'0' || c>'9') f=c=='-'?-1:1,c=getchar();
	while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
	return x*f;
}

inline void get_miu(){
	miu[1]=1;
	for(int i=2;i<maxn;i++){
		if(!vis[i]) pri[++cnt]=i,miu[i]=-1;
		for(int j=1;j<=cnt && i*pri[j]<maxn;j++){
			vis[i*pri[j]]=1;
			if(i%pri[j]==0) {
				miu[i*pri[j]]=0;
				break;
			}
			miu[i*pri[j]]=-miu[i];
		}
	}
	for(int i=2;i<maxn;i++)
		if(miu[i]!=0)
			for(LL j=i;j<maxn;j*=i)
				g[j]=-miu[i];
	for(int i=2;i<maxn;i++) g[i]+=g[i-1];
}

int main(){
	get_miu();
	t=rd();
	while(t--){
		n=rd(); m=rd(); LL ans=0;
		if(n>m) swap(n,m);
		int l,r;
		for(l=1;l<=n;l=r+1){
			int x=n/l,y=m/l; r=min(n/x,m/y);
			ans+=1LL*(g[r]-g[l-1])*x*y;
		}
		printf("%lld\n",ans);
	}
	return 0;
}