学习笔记第五节：莫比乌斯反演

最新推荐文章于 2024-09-02 13:47:18 发布

Deep_Kevin

最新推荐文章于 2024-09-02 13:47:18 发布

阅读量329

点赞数 1

分类专栏：学习笔记莫比乌斯反演

本文链接：https://blog.csdn.net/Deep_Kevin/article/details/80224492

版权

学习笔记同时被 2 个专栏收录

61 篇文章 3 订阅

订阅专栏

莫比乌斯反演

12 篇文章 0 订阅

订阅专栏

正题

现在在这里再来讲述一遍莫比乌斯反演。

莫比乌斯反演基于一条公式 $\sum_{d\mid n} \mu(d) =[n=1]$ 。

不懂的小伙伴先去找数论函数的介绍看看。

假如现在要求 $\sum_{i=1}^n\sum_{j=1}^m [gcd(i,j)=1]$ , $n<=m<=10^{7}$ .~~开到更大可以用杜教筛来做~~。

那么我们可以把后面的 [gcd(i,j)=1] ，换成 $\sum_{d\mid gcd(i,j)} \mu(d)$

就变成了

$\sum_{i=1}^{n}\sum_{j=1}^m\sum_{d\mid gcd(i,j)}\mu(d) \\=\sum_{i=1}^{n}\sum_{j=1}^m\sum_{d\mid i\ \Lambda \ d\mid j}\mu(d) \\=\sum_{d=1}^{n}\mu(d)\ floor(n/i)\ floor(m/i)$

那么我们又可以预处理 $\mu(d)$ 的前缀和，又可以对后面的东西进行整除分块。

时间就是 $O(\sqrt n)$ 了。（认为n,m同阶）。

程序中就是直接按照定义欧拉筛出 mu(d) ，然后做一遍前缀和，后面整除分块本人喜欢这种写法。

mu[1]=1;
for(int i=2;i<=maxn;i++){
	if(!vis[i]) {p[++p[0]]=i;mu[i]=-1;}
	for(int j=1;j<=p[0] && (temp=i*p[j])<maxn;j++){
		vis[temp]=true;
		if(i%p[j]==0) break;
		mu[temp]=-mu[i];
	}
	mu[i]+=mu[i-1];
}
scanf("%d",&T);
while(T--){
	scanf("%d %d",&n,&m);
	l=1;
	while(l<=n){
		r=min(n/(n/l),m/(m/l));
		ans+=(long long)(n/l)*(m/l)*(mu[r]-mu[l-1]);
		l=r+1;
	}
	printf("%lld\n",ans);
}

其实莫比乌斯反演是这样的：

有 h=f*1 ，*为狄利克雷卷积。