[SDOI2015]约数个数和 [莫比乌斯反演]

最新推荐文章于 2020-06-05 12:43:00 发布

FSYo

最新推荐文章于 2020-06-05 12:43:00 发布

阅读量132

点赞数

分类专栏：莫比乌斯反演

本文链接：https://blog.csdn.net/sslz_fsy/article/details/85031778

版权

莫比乌斯反演专栏收录该内容

28 篇文章 0 订阅

订阅专栏

传送门

$ans=f(1)=\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|i}[gcd(x,y)==1]=\sum _{x=1}^n\sum _{y=1}^m[n/x][m/y]=\sum _{d=1}^n \mu(d)F(d) =\sum _{d=1}^n \mu(d)\sum _{i=1}^n\sum _{j=1}^m[d|gcd(i,j)][n/i][m/j]=\sum _{d=1}^n \mu(d)\sum _{i=1}^{n/d}\sum _{j=1}^{m/d} [n/id][m/jd]$

用整除分块预处理 $\sum _{i=1}^n[n/i]$

发现对于连续的d , n/d的取值一样 , 也就是说 $\sum _{i=1}^{n/d}[n/id]$ 的取值一样 , 然后就可以整除分块做了

#include<bits/stdc++.h>
#define N 50050
#define LL long long
using namespace std;
int T,n,m,tot; LL s_mu[N],s_f[N];
int mu[N],p[N],isp[N]; 
void Solve(){
	scanf("%d%d",&n,&m); if(n>m) swap(n,m);
	LL ans = 0;
	for(int l=1,r;l<=n;l=r+1){
		int x1=n/l,x2=m/l; r = min(n/x1,m/x2);
		ans += (s_mu[r]-s_mu[l-1]) * s_f[x1] * s_f[x2];
	} printf("%lld\n",ans);
}
int main(){
	mu[1]=s_mu[1]=1;
	for(int i=2;i<=N-50;i++){
		if(!isp[i]) p[++tot]=i,mu[i]=-1;
		for(int j=1;j<=tot;j++){
			if(i*p[j]>N-50) break; isp[i*p[j]]=1; 
			if(i%p[j]==0) break; mu[i*p[j]] = -mu[i];
		}
		s_mu[i] = s_mu[i-1] + mu[i];
	}
	for(int i=1;i<=N-50;i++){ int ans=0;
		for(int l=1,r;l<=i;l=r+1){
			r = i/(i/l);
			ans += (r-l+1) * (i/l);
		} s_f[i] = ans;
	}
	scanf("%d",&T); while(T--) Solve(); 
	return 0;
}