容斥原理，莫比乌斯反演

莫比乌斯反演

本文写的主要是有关于组合计数的莫比乌斯反演，关于纯数学处理的莫比乌斯反演请看莫比乌斯反演。

莫比乌斯反演原理：
$[n=m]=\sum _{d|\frac{n}{m}}\mu (d)$

莫比乌斯反演：
$f(n)=\sum _{d|n}g(d)\iff g(n)=\sum _{d|n}\mu (d)f(\frac{n}{d})$

$f(n)=\stackrel{m}{\sum _{n|d}}g(d)\iff g(n)=\stackrel{m}{\sum _{n|d}}\mu \left(\frac{d}{n}\right)f(d)$

URAL-1091-Tmutarakan Exams

给定$n$个数，其数值取值范围为$[1,m]$，请问有多少种情况使得这$n$个数的最大公约数为$1$.
$n,m=10^5$

这题修正法比较好想：

如果不限制公约数关系，则答案为$m^n$，这里面算重了不互质的情况，每种情况只多算了一次。

如果设$f(k)$表示这$n$个数的最大公约数为$k$的倍数的情况数，设$g(k)$表示这$n$个数的最大公约数为$k$的情况数，则：
${(m/k)}^n=f(k)=\stackrel{m}{\sum _{k|d}}g(d)$

则：$g(k)=\stackrel{m}{\sum _{k|d}}\mu \left(\frac{d}{n}\right)f(d)$

答案为$g(1)$

P9118

这个题可以考虑修正：

$\lfloor \sqrt[i]{n}\rfloor -1$表示$[2,n]$当中完全$i$次方数的个数（因为$1^k=1$，特判），会发现算重了，具体是这样算重的：
设$g(x)$表示$[2,n]$当中能够表示为$a^x$并且不能表示为$b^y(y>x)$的数的个数，则$\lfloor \sqrt[x]{n}\rfloor -1=f(x)=\stackrel{\log {}_{2}n}{\sum _{x|d}}g(d)$

反演可以求出$g$，求出$g$的后缀和$+1$即为答案。

时间复杂度$O(k\log k)$

还有另一种可以求出系数的修正：

考虑如何求出$g(x)$，设容斥系数为$h_x(d)$，答案应为$\sum _{x|d}{h}_x(d)\cdot f(d)$

考虑一个实际上能够表示为$a^K$并且不能表示为$b^y(y>K)$的数的贡献为：
$\sum _{x|d|K}{h}_x(d)$

我们希望它等于$1$，可以指定它在$x=K$时做一次贡献，即：
$\sum _{x|d|K}{h}_x(d)=[x=K]$
$\sum _{d|\frac{K}{x}}{h}_x(dx)=\left[K=x\right]$

我们知道有：$[n=m]=\sum _{d|\frac{n}{m}}\mu (d)$，也就是希望$h_x(dx)=\mu (d)$
因此：$h_x(i)=\mu \left(\frac{i}{x}\right)$

答案应为$\sum _{x|d}\mu \left(\frac{d}{x}\right)\cdot f(d)$

同类题目

CF1884D

首先补集转化为求有多少二元组是不好的。

考虑批量处理本质相同的情况。一个二元组的本质就是他们的$gcd$。

这个题我赛时的思路是考虑修正。

$f(x)$表示$x|gcd(i,j)$的二元组$(i,j)$数量，设序列中有$y$个数是$x$的倍数，则$f(x)=\frac{y(y-1)}{2}$。

那么$f(x)=\sum _{x|d}g(d)$
可以反演求出$g(x)$表示$x=gcd$的二元数量。

答案为$\frac{n(n-1)}{2}-\sum _{i}g(i)[i被至少一个数整除]$

一点小细节是对于枚举倍数的时候记得对序列去重，不然对于一个全$1$序列来说，枚举每个数的倍数，复杂度就不是调和数了。

P3172

容易想到修正，把倍数和恰好的关系式列出来：
${\left(\lfloor \frac{R}{x}\rfloor -\lfloor \frac{L-1}{x}\rfloor \right)}^n=f(x)=\sum _{x|d}g(d)$，其中$f(x)$表示$x|gcd$的情况数，$g(x)$表示$gcd=x$的情况。

则：$ans=g(K)=\sum _{x|d}\mu \left(\frac{d}{x}\right)f(d)$

杜教筛预处理莫比乌斯函数前缀和即可，可以做到$O\left(H^{\frac{2}{3}}\right)$。

不过当时好像是2015年，可能杜教筛还不是很火，注意到$x>R-L+1$时${\left(\lfloor \frac{R}{x}\rfloor -\lfloor \frac{L-1}{x}\rfloor \right)}^n=0或1$，即${\left(\lfloor \frac{R}{x}\rfloor -\lfloor \frac{L-1}{x}\rfloor \right)}^1$，于是还可以纯数论分块。

后记

于是皆大欢喜。