利用容斥解决欧拉函数相关问题&&[NOI2010] 能量采集

例题 P1447 [NOI2010] 能量采集

Label

利用容斥原理求最大公约数为$k$的数对个数（或：求${\sum }_{i=1}^{n}gcd(i,n)$）

Description

P1447 [NOI2010] 能量采集

Solution

首先，我们应该推得：对于某株坐标$(x,y)$的植物，若$gcd(x,y)=d$，则挡在此植物前面的植物应有$(n-1)$株（$(\frac{x}{d},\frac{y}{d}),(\frac{2x}{d},\frac{2y}{d})...(\frac{(n-1)x}{d},\frac{(n-1)y}{d})$）。

所以，设$f_d={\sum }_{x=1}^n{\sum }_{y=1}^m[gcd(x,y)==d]$，则此题答案即为${\sum }_{d=1}^n(2d+1)f_d$。

由于$x,y$范围不一致，故此处$f$序列无法用欧拉函数求解。

根据容斥原理，$f_d$代表的所有元素首先取自于公约数集合内含有$d$的有序数对的集合，在此题中，这个集合的元素个数是$\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$。之后从此集合内去掉公约数大于$d$的元素，剩下的即为最大公约数为$d$的有序数对。又由于若两个数的公约数为$d$，那么它们的公约数显然是$d$的倍数，所以我们可以得到如下求$f_d$的公式：
​ $$f_d=\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor -\sum _{i=2}^{id\le min(m,n)}f[id]$$

程序复杂度$O(nlogn)$

显然，这个过程经稍微变动后可以用来解决求${\sum }_{i=1}^{n}gcd(i,n)$的问题。

利用容斥原理推导欧拉函数表达式$\phi (n)$

我们考虑将$n$改写为唯一分解定理的形式：
​ $$n=\prod p_i^{c_i}$$

考虑公式$$|\bigcap _{i=1}^n{S}_i|=|U|-|\bigcup _{i=1}^n{\complement }_U{S}_i|$$，只不过此处我们定义 S d i = { x ∣ d i ∤ x , x ∈ N + , x < n } ( d i 为 n 的 素 因 子 ) S_{d_i}=\{x|d_i\nmid x,x\in\N^+,x<n \}(d_i 为 n 的素因子) Sdi​​={x∣di​∤x,x∈N+,x<n}(di​为n的素因子)。

这样一来，我们便可直接导出欧拉函数任意项$\phi (n)$的公式：
$$\phi (n)=|\bigcap _{i=1}^n{S}_{d_i}|=|U|-|\bigcup _{i=1}^n{\complement }_U{S}_{d_i}|=n-\sum _{i=1}^n(-1{)}^{i-1}\frac{n}{{\prod }_{j=1}^i{p}_j}=n\prod _{i=1}^k(1-\frac{1}{p_i})$$

Code(P1447)