莫比乌斯教我变魔术

---

莫比乌斯反演

详细的不赘述了，都在这里
若$f(n)$是积性函数，则如下$F(x)$也是积性函数

$$F(x)=\sum_{d\mid x}f(d)$$ 则莫比乌斯反演公式为： $$f(x)=\sum_{d\mid x}\mu(d)F({x\over d})=\sum_{d\mid x}\mu({x\over d})F(d)$$

---

容斥原理

来考虑最简单的问题：1~n中的x的倍数$(1\le x\le n)$有多少个？

用下图来表示容斥的情况。每个圆内均有一个数字x，圆表示x的倍数的集合$F(x)$。
若$x\mid y$，则$F(y)\in F(x)$

定义$f(x)$为1~n中最大约数为x的数集，则$f(x)$就不是完整的圆了，可得

$$F(x)=\sum_{x\mid y,y\mid n}f(y)$$

其实$F(x)$很好求。如果反过来要求$f(x)$应该怎么办？

举例：n=30， 求$f(2)$
$f(2)=F(2)-F(4)-F(6)-F(10)+F(12)-F(14)+F(20)-F(22)-F(26)+F(28)+F(30)$

只有在$y|x时,F(y)$才会被考虑。设$F(y)$前的系数为$a_y$设p为素数：
当${y\over x}=p$时，$F(y)$只属于$F(x)$，故$a_y=-1$；
当${y\over x}=p_1p_2$时，$F(y)$属于$F(x),F(xp_1),F(xp_2),又a_{xp_1}=a_{xp_2}=-1$，故$a_y=1$；
…
当${y\over x}=p_1p_2...p_s$时，改变状态的有$-C_s^1+C_s^2-C_s^3...+(-1)^{s-1}C_s^{s-1}=-1-(-1)^s，故a_y=(-1)^s$；
若${x\over y_i}=p_1^{v_1}p_2^{v_2}...p_s^{v_s}且至少有一个v大于1$，则已有标记为

$$-C_s^1+C_s^2-C_s^3...+(-1)^sC_s^s=\sum_{i=0}^s{C_s^i(-1)^i}-C_s^0=(1-1)^s-1=-1$$ 则确定 $a_i=0，故a_y=0$.
这证明过程和莫比乌斯函数的 证明是一样的，得到的系数即 $\mu({y\over x})$。得到 $$f(x)=\sum_{x\mid y}\mu({y\over x})F(y)$$

感受一下差别

莫比乌斯	容斥
d是x的约数	y是x的倍数
$\mu({x\over d})$自变量是跳跃的	$\mu({y\over x})$自变量是连续的
$f(d)$为有限项	$f(y)$为无限项（除非规定了一个范围）

因此要根据不同的已知条件来确定是哪一种。

校赛C题

求区间[a,b]中与n互质的数的个数。

分析：求[1,a]中与n互质的个数。
设F(y)是[1,a]中y的倍数的个数，则F(y)=a/y；设f(x)是[1,a]中与n的最大公约数是x的数的个数。注意，这里和上面容斥原理定义的函数略有不同。$f(x)$的条件=旧有条件&新条件，造成$F(y)$的条件范围也缩小了。

$$F(y)=\sum_{y\mid x,y\le a}^{x\mid n}f(x)~~~\Rightarrow~~~ f(x)=\sum_{x\mid y}^{x\mid n}\mu({y\over x})F(y)$$ 即要求 $$f(1)=\sum_{1\le y\le a}^{y\mid n}\mu(y){a\over y}$$

POJ 3904

给n个数，求其中四个数的gcd是1的情况有多少种。

校赛D题

给n个数，求其中三个数的gcd是m的情况有多少种。
令$f(x)=三个数的gcd是x的情况数$
$F(x)=三个数的gcd是x的倍数的情况数$
n个数中有约数为x的数有y个，则$F(x)=C_y^3$
求$f(m)$

---

超级变变变

gcd

1. 求$\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}gcd(a,b)$
   预备：$\sum\limits_{d\mid n}\phi(d)=n$（证明）
   若$d\mid gcd(a,b)$，则$\sum\limits_{d\mid gcd(a,b)}\phi(d)=gcd(a,b)$，且$d\mid a,~d\mid b$,
   

   $$\begin{align*}\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}gcd(a,b)&=\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}\sum\limits_{d\mid a,d\mid b}\phi(d)\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d)\sum\limits_{a=1,d\mid a}^{N}\sum\limits_{b=1,d\mid b}^{M}1\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d)(\sum\limits_{a=1,d\mid a}^{N}1)(\sum\limits_{b=1,d\mid b}^{M}1)\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d){N\over d}{M\over d} \end{align*}$$

2. 求$\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}[gcd(a,b)=1]$
   预备：$\sum\limits_{d\mid n}\mu(d)=[n=1]$（证明）
   若$d\mid gcd(a,b)$，则$\sum\limits_{d\mid gcd(a,b)}\mu(d)=[gcd(a,b)=1]$，且$d\mid a,~d\mid b$
   

   $$\begin{align*} \sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}[gcd(a,b)=1]&=\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}\sum\limits_{d\mid a,d\mid b}\mu(d)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\mu(d)\sum\limits_{a=1,d\mid a}^{N}\sum\limits_{b=1,d\mid b}^{M}1\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\mu(d)(\sum\limits_{a=1,d\mid a}^{N}1)(\sum\limits_{b=1,d\mid b}^{M}1)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\mu(d){N\over d}{M\over d} \end{align*}$$

3. 求$\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}gcd(a,b)$（另一种证法）
   

   $$\begin{align*} \sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}gcd(a,b)&=\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}g·[gcd(a,b)=g]\\\ &=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{a=1}^{N/g}\sum\limits_{b=1}^{M/g}[gcd(a,b)=1]\\\ ^{(2式结果)}&=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{d=1}^{min\{N/g,M/g\}}\mu(d){N/g\over d}{M/g\over d}\\\ ^{(x=d·g)}&=\sum\limits_{x=1}^{min\{N,M\}}\sum\limits_{g\mid x}g·\mu({x\over g}){N\over x}{M\over x}\\\ &=\sum\limits_{x=1}^{min\{N,M\}}{N\over x}{M\over x}\sum\limits_{g\mid x}g·\mu({x\over g}) \end{align*}$$
   与1的结果对比 $$\sum\limits_{g\mid x}g·\mu({x\over g})=\phi(x)$$ 两者相等。

4. 求$\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}lcm(a,b)$
   

   $$\begin{align*}\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}lcm(a,b)&=\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}{ab\over gcd(a,b)}\\ &=\sum\limits_{g=1}^{min\{N,M\}}{1\over g}\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}ab[gcd(a,b)=g]\\ &=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{a=1}^{N/g}\sum\limits_{b=1}^{M/g}ab[gcd(a,b)=1]\\ ^{由2得到}&=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{d=1}^{min\{N/g,M/g\}}\mu(d)(\sum\limits_{a=1,d\mid a}^{N/g}a)(\sum\limits_{b=1,d\mid b}^{M/g}b)\\ &=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{d=1}^{min\{N/g,M/g\}}\mu(d)(d\sum\limits_{a=1}^{N/g/d}a)(d\sum\limits_{b=1}^{M/g/d}b)\\ &=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{d=1}^{min\{N/g,M/g\}}\mu(d)(d\sum\limits_{a=1}^{N/g/d}a)(d\sum\limits_{b=1}^{M/g/d}b)\\ &=\sum\limits_{g=1}^{min\{N,M\}}g\sum\limits_{d=1}^{min\{N/g,M/g\}}\mu(d)·{1\over 4}d^2·{N/g\over d}({N/g\over d}+1){M/g\over d}({M/g\over d}+1)\\ ^{x=dg}&={1\over 4}\sum\limits_{x=1}^{min\{N,M\}}{N\over x}{M\over x}({N\over x}+1)({M\over x}+1)x\sum\limits_{d\mid x}d·\mu(d) \end{align*}$$
   这里$x\sum\limits_{d\mid x}d·\mu(d)$可以先预处理出来。

   ---

$O(2\sqrt N+2\sqrt M)$优化

当遇到类似于$\sum\limits_{x=1}^{min\{N,M\}}\phi(x){N\over x}{M\over x}$的通项公式时，一项一项求得速度很慢。观察${N\over x}和{M\over x}$的变化曲线，可以发现，当x比较大的时候，${N\over x}和{M\over x}$将保持在一些区间内不变。
如N=100时，x在区间[51,100]内${N\over x}$保持为1；M=56，x在区间[29,56]内${M\over x}$保持为1。因此x
在区间[51,56]内${N\over x}{M\over x}$保持不变为1。

对于N，x每逢一个约数都使得$N\over x$发生改变，故这样的结点大约有$2\sqrt N$个。
同时，要对欧拉函数求前缀和，使查询更快。

LL gao(int n, int m){
    LL re = 0;
    if(n > m) swap(n, m);
    int last;
    for(int i = 1;i <= n;i = last+1){
        last = min(n/(n/i), m/(m/i));
        re += 1LL*(n/i)*(m/i)*(sum[last]-sum[i-1]);
    }
    return re;
}

---

HDU 5528

欧拉函数变换

---

南大省赛

求$G(N,M)=\sum\limits_{a=1}^N\sum\limits_{b=1}^Ma^2b^2gcd(a,b)$
若$d\mid gcd(a,b)$，则$\sum\limits_{d\mid gcd(a,b)}\phi(d)=gcd(a,b)$，且$d\mid a,~d\mid b$

$$\begin{align*} \sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}a^2b^2gcd(a,b)&=\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}a^2b^2\sum\limits_{d\mid a,d\mid b}\phi(d)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d)\sum\limits_{a=1,d\mid a}^{N}\sum\limits_{b=1,d\mid b}^{M}a^2b^2\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d)(\sum\limits_{a=1,d\mid a}^{N}a^2)(\sum\limits_{b=1,d\mid b}^{M}b^2)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d)·d^4·(\sum\limits_{a=1}^{N/d}a^2)(\sum\limits_{b=1}^{M/d}b^2)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\phi(d)·d^4·{{N\over d}({N\over d}+1)(2{N\over d}+1)\over 6}·{{M\over d}({M\over d}+1)(2{M\over d}+1)\over 6} \end{align*}$$
另一种解法：
$$\begin{align*} \sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}a^2b^2[gcd(a,b)=1]&=\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}a^2b^2\sum\limits_{d\mid a,d\mid b}\mu(d)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\mu(d)\sum\limits_{a=1,d\mid a}^{N}\sum\limits_{b=1,d\mid b}^{M}a^2b^2\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\mu(d)(\sum\limits_{a=1,d\mid a}^{N}a^2)(\sum\limits_{b=1,d\mid b}^{M}b^2)\\\ &=\sum\limits_{d=1}^{min\{N,M\}}\mu(d)·d^4·(\sum\limits_{a=1}^{N/d}a^2)(\sum\limits_{b=1}^{M/d}b^2) \end{align*}$$

$$\begin{align*} \sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}a^2b^2gcd(a,b) =&\sum\limits_{a=1}^{N}\sum\limits_{b=1}^{M}a^2b^2g·[gcd(a,b)=g]\\\ =&\sum\limits_{g=1}^{min\{N,M\}}g^5\sum\limits_{a=1}^{N/g}\sum\limits_{b=1}^{M/g}a^2b^2[gcd(a,b)=1]\\\ ^{(上式结果)}=&\sum\limits_{g=1}^{min\{N,M\}}g^5\sum\limits_{d=1}^{min\{N/g,M/g\}}\mu(d)·d^4·(\sum\limits_{a=1}^{N/g/d}a^2)(\sum\limits_{b=1}^{M/g/d}b^2)\\\ ^{(x=d·g)}=&\sum\limits_{x=1}^{min\{N,M\}}\sum\limits_{g\mid x}g^5·\mu({x\over g})·({x\over g})^4(\sum\limits_{a=1}^{N/x}a^2)(\sum\limits_{b=1}^{M/x}b^2)\\\ =&\sum\limits_{x=1}^{min\{N,M\}}(\sum\limits_{a=1}^{N/x}a^2)(\sum\limits_{b=1}^{M/x}b^2)\sum\limits_{g\mid x}g^5·\mu({x\over g})·({x\over g})^4\\\ =&\sum\limits_{x=1}^{min\{N,M\}}(\sum\limits_{a=1}^{N/x}a^2)(\sum\limits_{b=1}^{M/x}b^2)\sum\limits_{g\mid x}g.x^4·\mu({x\over g})\\\ =&\sum\limits_{x=1}^{min\{N,M\}}x^4(\sum\limits_{a=1}^{N/x}a^2)(\sum\limits_{b=1}^{M/x}b^2)\sum\limits_{g\mid x}g·\mu({x\over g})\\\ =&\sum\limits_{x=1}^{min\{N,M\}}x^4(\sum\limits_{a=1}^{N/x}a^2)(\sum\limits_{b=1}^{M/x}b^2)\phi(x)\\\ \end{align*}$$

---

计蒜客初赛第二场D题

在n*m个格点中选择两个格点作出直线。求所有不重合直线的条数。

同一斜率的直线有如图(n-1-a+m-1-b-1)条，其中$1<a<n,1<b<m,gcd(a,b)=1$，故所有的直线条数为

$$\sum_{a=1}^{n-1}\sum_{b=1}^{m-1}(n+m-1-a-b)[gcd(a,b)=1]$$
用莫比乌斯函数展开慢慢推。

---

公式打到手抽筋…