【笔记】莫比乌斯反演小归纳

莫比乌斯反演，又称懵逼钨丝繁衍，顾名思义就是一种让蒟蒻懵逼得像一根钨丝繁衍一般智障的算法
蒟蒻做了一天的懵逼钨丝繁衍，感觉十分智障

莫比乌斯反演

---

概念

若有两个函数，$f(x)$和$F(x)$，满足$$F(n)=\sum _{d|n}f(d)$$
则有$$f(x)=\sum _{d|n}\mu (d)F(\frac{n}{d})$$
其中$\mu$为莫比乌斯函数

证明

目前有两种证法 (蒟蒻只知道两种)：直接证和卷积

①直接证法：

$$\sum _{d|n}\mu (d)F(\frac{n}{d})=\sum _{d|n}\mu (d)\sum _{d^{\prime }|\frac{n}{d}}f(d^{\prime })=\sum _{d^{\prime }|n}f(d^{\prime })\sum _{d|\frac{n}{d^{\prime }}}\mu (d)=f(n)$$

②卷积证法：

引自知乎
定义俩函数，$F(x)$和$f(x)$，设卷积运算$$F\ast f=(F\ast f)(n)=\sum _{ij=n}F(i)f(j)$$
设卷积意义下的单位元$\iota$使得$\iota \ast f=f$
设乘法意义下的单位元$u$使得$(Ff)(n)=F(n)f(n)$，可得$u$是$\mu$在卷积意义下的逆元，即$u\ast \mu =\iota$
莫比乌斯反演公式：$$F(n)=\sum _{d|n}f(n)\iff f(n)=\sum _{d|n}\mu (d)F(\frac{n}{d})$$
变形：$$F=f\ast u\iff f=F\ast \mu$$
证明：$$F=f\ast u$$$$\iff F\ast \mu =f\ast u\ast \mu$$$$\iff F\ast \mu =f\ast \iota$$$$\iff F\ast \mu =f$$

线性筛莫比乌斯函数

具体蒟蒻不会证，贴上代码吧

	miu[1]=1;
    for(rg int i=2;i<=5e4;++i){
        if(!is[i])pri[++tot]=i,miu[i]=-1;
        for(rg int j=1;j<=tot;++j){
            rg int k=i*pri[j];if(k>5e4)break;
            is[k]=1;
            if(i%pri[j]==0){miu[k]=0;break;}
            else miu[k]-=miu[i];
        }
    }

应用

GCD HDU-1695

题意
给定$a$,$b$,$c$，求有序数对$(x,y)$满足$x\in [1,a],y\in [1,b],gcd(x,y)=c$的个数

题解
首先简化问题，设答案为$ans(a,b,c)$，可证答案为$ans(\frac{a}{c},\frac{b}{c},1)$

设$F(n)$表示$(x,y)$满足$x\in [1,a],y\in [1,b],c|gcd(x,y)$的个数
设$f(n)$表示$(x,y)$满足$x\in [1,a],y\in [1,b],gcd(x,y)=c$的个数
则可证$F(x)$与$f(x)$符合$F(n)={\sum }_{d|n}f(n)$，最后答案为$f(1)$
可证$F(d)=\frac{n}{d}\ast \frac{m}{d}$
答案为$f(1)={\sum }_d^{min(a,b)}\mu (d)F(d)$

代码

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
#define rg register

template <typename _Tp> inline void read(_Tp&x){
	char c11=getchar();x=0;bool booo=0;
	while(c11!='-'&&!isdigit(c11))c11=getchar();if(c11=='-'){c11=getchar();booo=1;}
	while(isdigit(c11)){x=x*10+c11-'0';c11=getchar();}if(booo)x=-x;return ;
}

const int N=100500;
ll miu[N],is[N],pri[N],tot=0,ans1,ans2;

void pre_miu(){
	miu[1]=1;
	for(rg int i=2;i<=1e5;++i){
		if(!is[i])pri[++tot]=i,miu[i]=-1;
		for(rg int j=1;j<=tot;++j){
			rg int k=pri[j]*i;if(k>1e5)break;
			is[k]=1;
			if(i%pri[j]==0){miu[k]=0;break;}
			else miu[k]-=miu[i];
		}
	}
	return ;
}

int main(){
	rg int T,a,b,c;read(T);
	pre_miu();
	for(rg int cases=1;cases<=T;++cases){
		ans1=ans2=0;
		read(a);read(a);read(b);read(b);read(c);
		if(!c){printf("Case %d: 0\n",cases);continue;}
		a/=c,b/=c;
		if(a>b)swap(a,b);
		for(rg int i=1;i<=a;++i)
			ans1+=(ll)miu[i]*(a/i)*(b/i),
			ans2+=(ll)miu[i]*(a/i)*(a/i);
		printf("Case %d: %lld\n",cases,ans1-(ans2>>1));
	}
	return 0;
}

[POI2007]ZAP-Queries

和楼上一样，不过数据大了一些，需要分块优化：由于有一段$\frac{n}{i}$是相同的，所以可以进行跳跃式加和，复杂度$O(\sqrt{min(n,m)})$

代码

		a/=c,b/=c,ans=0;
        if(a>b)swap(a,b);
        for(rg int i=1,j;i<=a;i=j+1){
            j=min(a/(a/i),b/(b/i));
            ans+=1ll*(sum[j]-sum[i-1])*(a/i)*(b/i);
        }

[HAOI2011]Problem b

又和楼上差不多，只不过加了下限，用上容斥

后面两题和第一题代码类似，加以修改即可