莫比乌斯反演--懵逼反演系列

反演的定义

设有数论函数$f(n),g(n)$，其中$g(n)$已知
且$f(n),g(n)$满足关系$g(n)={\sum }_{i=0}^n{a}_{n,i}f(i)$

反演就是利用$g(n)$和已知关系推导出$f(n)$关于$g(i)$的表达
$f(n)={\sum }_{i=0}^n{b}_{n,i}g(i)$

---

莫比乌斯函数

了解莫比乌斯反演之前肯定要先了解莫比乌斯函数
设正整数$N$按算术基本定理分解质因数为$N={\prod }_{i=1}^m{p}_i^{c_i}$

$$\mu (N)=\{\begin{array}{r}0,\exists i\in [1,m],c_i>1\\ 1,m\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}2),\forall i\in [1,m]c_i=1\\ -1,m\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}2),\forall i\in [1,m]c_i=1\end{array}$$

即若$N$含相同质因子，则$\mu (N)=0$
若$N$所有质因子各不相等且质因子个数为偶数，则$\mu (N)=1$
若$N$所有质因子各不相等且质因子个数为奇数，则$\mu (N)=-1$

特别的$\mu (1)=1$

莫比乌斯函数的性质

• ${\sum }_{d|n}\mu (d)=[n=1]$

这个性质是接下来证明的重要论据

• 莫比乌斯函数是积性函数

这个意味着莫比乌斯函数可以用埃式筛或线性筛求解

• 对任意正整数$N$,${\sum }_{d|n}\frac{\mu (d)}{d}=\frac{\phi (n)}{n}$

这个挺好玩的，但是没啥用

---

莫比乌斯反演

常见的莫比乌斯反演有两种形式

• $$g(n)=\sum _{d|n}f(d)\Leftarrow \Rightarrow f(n)=\sum _{d|n}\mu (d)g(\frac{n}{d})$$
• $$g(n)=\sum _{n|d}f(d)\Leftarrow \Rightarrow f(n)=\sum _{n|d}\mu (\frac{d}{n})g(d)$$

形式一证明

$$f(n)=\sum _{d|n}\mu (d)g(\frac{n}{d})$$
$$=\sum _{d|n}\mu (d)\sum _{k|\frac{n}{d}}f(k)=\sum _{d|n}\sum _{k|\frac{n}{d}}\mu (d)f(k)$$
注意到$d,k$满足$d|n,k|(n/d)$，所以容易得出$k|n,d|(n/k)$，继续变换

$$\sum _{d|n}\sum _{k|\frac{n}{d}}\mu (d)f(k)=\sum _{k|n}\sum _{d|\frac{n}{k}}f(k)\mu (d)=\sum _{k|n}f(k)\sum _{d|\frac{n}{k}}\mu (d)$$

根据前文所述莫比乌斯函数性质一
有${\sum }_{d|\frac{n}{k}}\mu (d)=[n=k]$，带入后得$f(n)=f(n)\ast 1$成立

证毕

形式二证明

证明过程与前面证明类似
$$f(n)=\sum _{k=1}^{+\infty }\mu (k)F(nk)$$
$$=\sum _{k=1}^{+\infty }\mu (k)\sum _{nk|t}f(t)$$
$$=\sum _{n|t}f(t)\sum _{k|\frac{t}{n}}\mu (k)$$

同样带入性质一得证

---

莫比乌斯反演的应用

洛谷P3455 [POI2007]ZAP-Queries
对于给定的整数a,b和d，有多少正整数对x,y，满足x<=a，y<=b，并且gcd(x,y)=d
$1\le n\le 50000,1\le d\le a,b\le 50000$

莫比乌斯反演的经典例题
设$f(k)$表示满足条件的$(x,y)$数量，即$f(k)={\sum }_{i=1}^a{\sum }_{j=1}^b[gcd(i,j)=k]$
$g(k)$表示满足$gcd(x,y)\mathrm{m}\mathrm{o}\mathrm{d}k==0$得$(x,y)$数量$(x\le a,y\le b)$

那么有$$g(n)=\lfloor \frac{a}{n}\rfloor \lfloor \frac{b}{n}\rfloor =\sum _{n|k}f(k)$$

于是根据上述反演公式得$$ans=f(d)=\sum _{d|k}\mu (\lfloor \frac{k}{d}\rfloor )g(k)$$

令$t=k/d$
$$\sum _{t=1}^{min(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor )}\mu (t)\lfloor \frac{a}{td}\rfloor \lfloor \frac{b}{td}\rfloor$$

到这里已经可以$O(n)$处理每次询问，但由于多组询问复杂度还是不太乐观
对于$\lfloor \frac{a}{td}\rfloor \lfloor \frac{b}{td}\rfloor$这一部分用整除分块处理，处理一次询问的时间复杂度为$O(\sqrt{n})$

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long lt;
  
lt read()
{
    lt f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}
  
const int maxn=50010;
int T;
int miu[maxn],vis[maxn];
int prim[maxn],cnt;
lt sum[maxn];

void Miu(int n)
{
    miu[1]=1;
    for(int i=2;i<=n;++i)
    {
        if(!vis[i]){ miu[i]=-1; prim[++cnt]=i;}
        for(int j=1;j<=cnt;++j)
        {
            if(prim[j]*i>n) break;
            vis[prim[j]*i]=1;
            if(i%prim[j]==0) break;
            else miu[i*prim[j]]=-miu[i];
        }
    } 
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+miu[i];
}

int main()
{
    T=read(); Miu(maxn-10);
    while(T--)
    {
        lt a=read(),b=read(),d=read();
        
        lt ans=0,lim=min(a/d,b/d);
        int ad=a/d,bd=b/d;
        for(int ll=1,rr;ll<=lim;ll=rr+1)
        {
            rr=min(ad/(ad/ll),bd/(bd/ll));
            ans+=(a/(ll*d))*(b/(ll*d))*(sum[rr]-sum[ll-1]);
        }
        
        printf("%lld\n",ans);
    }
    return 0;
}

---

换个角度看反演

其实大多时候解反演题并不一定要设出$f(n),g(n)$
把莫比乌斯反演公式的推导过程铺开
我们会发现其实就是在想办法构造${\sum }_{d|n}\mu (d)=[n==1]$

那上面这题为例换个方法反演，首先$ans={\sum }_{i=1}^a{\sum }_{j=1}^b[gcd(i,j)=d]$
若$gcd(x,y)=1$，则有$gcd(x\ast d,y\ast d)=d$，于是上述式子变形得

$$\sum _{x=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{y=1}^{\lfloor \frac{b}{d}\rfloor }[gcd(x,y)=1]$$

现在我们惊喜的发现$[gcd(x,y)=1]$可以用莫比乌斯函数性质一做个替换

$$\sum _{x=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{y=1}^{\lfloor \frac{b}{d}\rfloor }\sum _{t|gcd(x,y)}\mu (t)$$

把枚举约数提前

$$\sum _{t=1}^{min(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor )}\sum _{x=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{y=1}^{\lfloor \frac{b}{d}\rfloor }\mu (t)[t|gcd(x,y)]$$

若$[t|gcd(x,y)]=1$成立，则$t|x,t|y$，所以把枚举$x,y$变成枚举$t$的倍数
最后把$\mu (t)$提前

$$\sum _{t=1}^{min(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor )}\mu (t)\sum _{t_1=1}^{\lfloor \frac{a}{dt}\rfloor }\sum _{t_2=1}^{\lfloor \frac{b}{dt}\rfloor }1=\sum _{t=1}^{min(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor )}\mu (t)\lfloor \frac{a}{td}\rfloor \lfloor \frac{b}{td}\rfloor$$

竟然得到了完全一样的答案！！

所以这说明不是遇上反演的题就一定要凑出$f,g$和公式所述的关系
而是应该思考如何将表达式构造出能利用${\sum }_{d|n}\mu (d)=[n==1]$的形式

---

继续懵逼反演的应用

BZOJ2301 [HAOI2011]Problem b
对于给定的整数a,b,c,d,k，有多少正整数对x,y，满足a<=x<=b，c<=y<=d，并且gcd(x,y)=k
1≤n≤50000，1≤a≤b≤50000，1≤c≤d≤50000，1≤k≤50000

和上一题大同小异，加一个容斥即可
如果被卡常的话可以考虑把一些没必要开long long的变量换成int
以及$a/k,b/k$这种重复计算的可以先处理，可以快很多

lt query(int a,int b)
{
    lt res=0;
    int ak=a/k,bk=b/k,lim=min(a/k,b/k);
    for(int ll=1,rr;ll<=lim;ll=rr+1)
    {
        rr=min(ak/(ak/ll),bk/(bk/ll));
        res+=(lt)(ak/ll)*(bk/ll)*(sum[rr]-sum[ll-1]);
    }
    return res;
}

lt ans=query(b,d)-query(b,c-1)-query(a-1,d)+query(a-1,c-1);

HDU - 1695 GCD
对于给定的整数a,b,c,d,k，有多少正整数对x,y，满足1<=x<=b，1<=y<=d，并且gcd(x,y)=k
(x=i,y=j)与(x=j,y=i)等价
T<=3000,0 < a <= b <= 100,000, 0 < c <= d <= 100,000, 0 <= k <= 100,000

和第一题区别在于去重，由于只有$x<=b,y<=b(b<d)$这部分会有重复，所以减去这部分总数的一半
以及注意特判k==0的情况

if(b>d) swap(b,d);
lt ans1=query(b,d,k),ans2=query(b,b,k);
printf("Case %d: %lld\n",cs,ans1-ans2/2);

---

洛谷P2257 YY的GCD
给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对
T = 10000，N, M <= 10000000

首先这题$f(k),g(k)$定义依然同上 (一般看到gcd的题都会想到这么设函数)
回顾上面反演后的式子
$$f(d)=\sum _{d|k}\mu (\lfloor \frac{k}{d}\rfloor )g(k)=\sum _{t=1}^{min(\lfloor \frac{N}{d}\rfloor ,\lfloor \frac{M}{d}\rfloor )}\mu (t)\lfloor \frac{N}{td}\rfloor \lfloor \frac{M}{td}\rfloor$$

现在考虑将这个式子应用到此题
$$ans=\sum _{p\in prime}f(p)=\sum _{p\in prime}\sum _{p|k}\mu (\lfloor \frac{k}{p}\rfloor )g(k)=\sum _{p\in prime}\sum _{t=1}^{min(\lfloor \frac{N}{p}\rfloor ,\lfloor \frac{M}{p}\rfloor )}\mu (t)\lfloor \frac{N}{tp}\rfloor \lfloor \frac{M}{tp}\rfloor$$

令$T=tp$，继续变形得

$$\sum _{T=1}^{min(N,M)}\sum _{p|T,p\in prime}\mu (\frac{T}{p})\lfloor \frac{N}{T}\rfloor \lfloor \frac{M}{T}\rfloor =\sum _{T=1}^{min(N,M)}\lfloor \frac{N}{T}\rfloor \lfloor \frac{M}{T}\rfloor \sum _{p|T,p\in prime}\mu (\frac{T}{p})$$
于是后面那一坨可以预处理前缀和，每次询问用整除分块处理即可

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long lt;
  
int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}
  
const int maxn=10000010;
int T;
int miu[maxn],vis[maxn];
int prim[maxn],cnt;
lt sum[maxn];

void Miu(int n)
{
    miu[1]=1;
    for(int i=2;i<=n;++i)
    {
        if(!vis[i]){ miu[i]=-1; prim[++cnt]=i;}
        for(int j=1;j<=cnt;++j)
        {
            if(prim[j]*i>n) break;
            vis[prim[j]*i]=1;
            if(i%prim[j]==0) break;
            else miu[i*prim[j]]=-miu[i];
        }
    } 
    for(int i=1;i<=cnt;++i)
    for(int j=1;prim[i]*j<=n;++j) sum[prim[i]*j]+=miu[j];
    for(int i=1;i<=n;++i) sum[i]=sum[i-1]+sum[i];
}

lt query(int n,int m)
{
    lt res=0; int lim=min(n,m);
    for(int ll=1,rr;ll<=lim;ll=rr+1)
    {
        rr=min(n/(n/ll),m/(m/ll));
        res+=(lt)(n/ll)*(m/ll)*(sum[rr]-sum[ll-1]);
    }
    return res;
}

int main()
{
    T=read(); Miu(maxn-5);
    while(T--)
    {
        int n=read(),m=read();
        printf("%lld\n",query(n,m));
    }
    return 0;
}

---

其他懵逼反演应用推导比较长，戳下面链接看噢
BZOJ2154 || 洛谷P1829 [国家集训队]Crash的数字表格 题解
BZOJ3529 || 洛谷P3312 [SDOI2014]数表 题解
BZOJ3994 || P3327 [SDOI2015]约数个数和 题解