数论知识总结-莫比乌斯反演

数论知识总结-莫比乌斯反演

---

NOIP爆零の蒟蒻又来学数论辣

用yyb的博客学了一下莫比乌斯反演

---

莫比乌斯函数

玄学

[ mu(x)=left{ egin{aligned} 1 & (x=1) \ 0 & (x=p^2 imes k,kin {mathbb{Z}^+}) \ (-1)^k & (x=Pi_{i=1}^{n}{p_i^{k_i}},k_i=1) \ end{aligned} ight. ]

x=1时，(mu(x)=1)
x包含质数的平方时，(mu(x)=0)
否则令k=x的质因子个数，(mu(x)=(-1)^k)

---

莫比乌斯反演

有两种形式？？

如果有两个函数(f(x),g(x))：

若满足(g(x)=sum_{d|x}f(d))，就能推出(f(x)=sum_{d|x}mu(frac{x}{d})g(d))

若满足(g(x)=sum_{x|d,dle n}f(d))，就能推出(f(x)=sum_{x|d,dle n}mu(frac{d}{x})g(d))

如果求不出(f(x))，但能快速求出(g(x))，任意一个(f(x))就可以快速算出来。

证明不会嗷嗷嗷嗷嗷嗷，能背就好辣

---

例题

[POI2007]ZAP-Queries

将a/=d,b/=d，问题转化为求gcd(i,j)=1的对数

设(f(x))为满足(gcd(i,j)=x)的有序数对((i,j))数量，

(F(x))为(gcd(i,j))为(x)倍数的有序数对((i,j))数量。

关系式：(F(x)=sum_{x|d,dle n}f(d))。

现在要求(f(1))

(F(x))很好求，(F(x)=lfloorfrac{a}{x} floor imeslfloorfrac{b}{x} floor)

so (f(x)=sum_{x|d,dle n}mu(frac{d}{x})F(d))

这样枚举可以获得70分的好成绩

好像还有个东西叫数论分块

(F(x))可以分成很多值相等的段，而且最多(O(sqrt{n}))段

于是就可以做了

如果求无序数对的话可以写个solve函数，ans=solve(a,b)-solve(min(a,b),min(a,b))/2

// It is made by XZZ
#include<cstdio>
#include<algorithm>
#define il inline
#define rg register
#define vd void
#define sta static
#define int long long
using std::min;
il int gi(){
    rg int x=0,f=1;rg char ch=getchar();
    while(ch<'0'||ch>'9')f=ch=='-'?-1:f,ch=getchar();
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int maxn=50001;
int mu[maxn],Mu[maxn],pr[maxn],g[maxn];
bool flg[maxn];
main(){
    mu[1]=1;
    for(rg int i=2;i<maxn;++i){
        if(!flg[i])pr[++pr[0]]=i,mu[i]=-1;
        for(rg int j=1;i*pr[j]<maxn&&j<=pr[0];++j){
            flg[i*pr[j]]=1;
            if(i%pr[j]==0){mu[i*pr[j]]=0;break;}
            mu[i*pr[j]]=-mu[i];
        }
    }
    for(rg int i=1;i<maxn;++i)Mu[i]=mu[i]+Mu[i-1];
    int T=gi(),a,b,k,ans,i,j;
    for(rg int zsy=1;zsy<=T;++zsy){
        a=gi(),b=gi(),k=gi();
        if(k==0){puts("0");continue;}
        a/=k,b/=k;ans=0;i=1;
        while(i<=min(a,b)){
            j=min(a/(a/i),b/(b/i));
            ans+=(Mu[j]-Mu[i-1])*(a/j)*(b/j);
            i=j+1;
        }
        printf("%lld
",ans);
    }
    return 0;
}

---

[蛤OI2011]Problem B

和上面一样，简单利用容斥。

// It is made by XZZ
#include<cstdio>
#include<algorithm>
#define il inline
#define rg register
#define vd void
#define sta static
using namespace std;
typedef long long ll;
il int gi(){
    rg int x=0,f=1;rg char ch=getchar();
    while(ch<'0'||ch>'9')f=ch=='-'?-1:f,ch=getchar();
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int maxn=50001;
int mu[maxn],Mu[maxn],pr[maxn];
bool yes[maxn];
il ll solve(int x,int y,int k){
    if(x==0||y==0)return 0;
    if(x>y)swap(x,y);
    x/=k,y/=k;
    int i=1,j;
    ll ret=0;
    while(i<=x){
        j=min(x/(x/i),y/(y/i));
        ret+=1ll*(Mu[j]-Mu[i-1])*(x/j)*(y/j);
        i=j+1;
    }
    return ret;
}
main(){
    mu[1]=1;
    for(rg int i=2;i<maxn;++i){
        if(!yes[i])pr[++pr[0]]=i,mu[i]=-1;
        for(rg int j=1;i*pr[j]<maxn&&j<=pr[0];++j){
            yes[i*pr[j]]=1;
            if(i%pr[j]==0){mu[i*pr[j]]=0;break;}
            mu[i*pr[j]]=-mu[i];
        }
    }
    for(rg int i=1;i<maxn;++i)Mu[i]=Mu[i-1]+mu[i];
    int T=gi(),a,b,c,d,k;
    while(T--){
        a=gi(),b=gi(),c=gi(),d=gi(),k=gi();
        printf("%lld
",solve(b,d,k)-solve(b,c-1,k)-solve(a-1,d,k)+solve(a-1,c-1,k));
    }
    return 0;
}

---

长期更新