gcd求和的几种方法

最大公约数 : Greatest Common Divisor, 简称 : gcd
定义符号 : $d|n$: $n\%d=0$
定义函数 : $\varphi(n)$: 欧拉函数
定义函数 : $\mu(n)$: 莫比乌斯函数

---

一维形式: $\sum\limits_{i=1}^{n}gcd(n,i)$

1. $1\le n\le10^6$

相信大家都会

for(int i=1;i<=n;++i)
    ans+=gcd(n,i);

---

2. $1\le n\le10^{9}$

设: $gcd(n,i)=d$
显然: $n\%d=0$
原问题转换为:
$\sum\limits_{d|n}d\cdotp\sum\limits_{i=1}^n\big[ gcd(n,i)=d\big]$
$\sum\limits_{d|n}d\cdotp\sum\limits_{i=1}^{\frac{n}{d}}\big[ gcd(\frac{n}{d},i)=1\big]$
$\sum\limits_{d|n}d\cdotp \varphi(\frac{n}{d})$
$O(\sqrt n)$枚举$n$的因子, 再$O(\sqrt\frac{n}{d})$计算$\varphi(\frac{n}{d})$就好了

---

3. $1\le n\le10^{12}$

不难发现$\sum\limits_{d|n}d\cdotp \varphi(\frac{n}{d})$是狄利克雷卷积形式
$f(n)=n$ 与 $\varphi(n)$均为积性函数

所以$g(n)=\sum\limits_{d|n}d\cdotp \varphi(\frac{n}{d})$为积性函数
若: $n=p_1^{e_1}\cdotp p_2^{e_2}\cdots p_k^{e_k}\qquad p_1, p_2,\cdots ,p_k$均为质数

那么: $g(n)=g(p_1^{e_1})\cdot g(p_2^{e_2})\times\cdots \times g(p_k^{e_k})$

考虑: $g(p^k)=\sum\limits_{d|p^k}d\cdotp\varphi(\frac{p^k}{d})\qquad p$为质数
即: $g(p^k)=\sum\limits_{i=0}^k p^{k-i}\cdotp \varphi(p^i)=\sum\limits_{i=1}^k p^{k-i}\cdotp \varphi(p^i)+p^k\cdotp\varphi(1)$
$g(p^k)=p^k+\sum\limits_{i=1}^k p^{k-i}\cdotp(p^i-p^{i-1})=p^k+\sum\limits_{i=1}^k (p^{k}-p^{k-1})$

$g(p^k)=(k+1)\cdotp p^k-k\cdotp p^{k-1}$

复杂度: $O(\sqrt n)$

---

二维形式: $\sum\limits_{i=1}^n \sum\limits_{j=1}^mgcd(i,j)$

1. $1\le n,m\le10^3$

相信大家都会

for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
        ans+=gcd(i,j);

---

2. $1\le n=m\le10^5\qquad q$组询问

设: $g(n)=\sum\limits_{i=1}^{n}gcd(n,i)$
我们已经可以在$O(\sqrt n)$内得到$g(n)$

$\begin{align} \sum\limits_{i=1}^n\sum\limits_{j=1}^ngcd(i,j)&=2\cdotp\sum\limits_{i=1}^n\sum\limits_{j=1}^igcd(i,j)-\sum\limits_{i=1}^ngcd(i,i)\\ &=2\cdotp\sum\limits_{i=1}^n\sum\limits_{j=1}^igcd(i,j)-\frac{(n+1)\cdotp n}{2}\\ &=2\cdotp \sum\limits_{i=1}^ng(i)-\frac{(n+1)\cdotp n}{2}\\ \end{align}$
复杂度: $O(n\sqrt n+q)$

---

3. $1\le n=m\le 10^7\qquad q$组询问

解法1:

枚举: $gcd(i,j)=d\in[1,n]$
原问题转换为:

$\sum\limits_{d=1}^nd\sum\limits_{i=1}^n\sum\limits_{j=1}^n[gcd(i,j)=d]$

$\sum\limits_{d=1}^nd\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor\frac{n}{d}\rfloor}[gcd(i,j)=1]$

设: $f(n)=\sum\limits_{i=1}^n\sum\limits_{j=1}^n[gcd(i,j)=1]$
设: $g(n)=\sum\limits_{i=1}^n\sum\limits_{j=1}^i[gcd(i,j)=1]=\sum\limits_{i=1}^n\varphi(i)$
不难发现: $f(n)=2\cdotp g(n)-1$

$\sum\limits_{d=1}^nd\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor\frac{n}{d}\rfloor}[gcd(i,j)=1]=\sum\limits_{d=1}^nd\cdotp f(\lfloor\frac{n}{d}\rfloor)=\sum\limits_{d=1}^nd\cdotp(2\cdotp g(\lfloor\frac{n}{d}\rfloor)-1)$

显然: $\sum\limits_{d=1}^nd\cdotp(2\cdotp g(\lfloor\frac{n}{d}\rfloor)-1)$可以分块

复杂度: $O(n+q\cdotp\sqrt n)$

解法2:

$\sum\limits_{i=1}^n\sum\limits_{j=1}^ngcd(i,j)=2\cdotp\sum\limits_{i=1}^n\sum\limits_{j=1}^igcd(i,j)-\frac{(n+1)\cdotp n}{2}$

设: $g(n)=\sum\limits_{i=1}^n\sum\limits_{j=1}^igcd(i,j)=\sum\limits_{i=1}^n\sum\limits_{d|i}d\cdotp\varphi(\frac{i}{d})$

$\sum\limits_{d|n}d\cdotp\varphi(\frac{n}{d})$是积性函数, 可以线性筛, 同时维护它的前缀和就可以$O(1)$回答询问

复杂度: $O(n+q)$

---

4. $1\le n\le m\le10^5$

枚举: $gcd(i,j)=d\in[1,n]$

原问题转换为:
$\sum\limits_{d=1}^nd\sum\limits_{i=1}^n\sum\limits_{j=1}^m[gcd(i,j)=d]$

设: $g(x)=\sum\limits_{i=1}^n\sum\limits_{j=1}^m[gcd(i,j)=x]$

设: $f(x)=\sum\limits_{i=1}^n\sum\limits_{j=1}^m[x\big|gcd(i,j)]=\lfloor\frac{n}{x}\rfloor\cdotp\lfloor\frac{m}{x}\rfloor$

莫比乌斯反演: $f(x)=\sum\limits_{x|d}g(d)\quad\Rightarrow\quad g(x)=\sum\limits_{x|d}\mu(\frac{d}{x})f(d)$

$g(x)=\sum\limits_{x|d}\mu(\frac{d}{x})\cdotp\lfloor\frac{n}{d}\rfloor\cdotp\lfloor\frac{m}{d}\rfloor$

$\sum\limits_{d=1}^nd\sum\limits_{d|i}\mu(\frac{i}{d})\cdotp\lfloor\frac{n}{i}\rfloor\cdotp\lfloor\frac{m}{i}\rfloor$

复杂度: $O(nlogn)$

---

5. $1\le n\le m\le10^6$

$\sum\limits_{d=1}^nd\sum\limits_{i=1}^n\sum\limits_{j=1}^m[gcd(i,j)=d]\quad\Rightarrow\quad\sum\limits_{d=1}^nd\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor \frac{m}{d}\rfloor}[gcd(i,j)=1]$

莫比乌斯反演得: $\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor \frac{m}{d}\rfloor}[gcd(i,j)=1] =\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)\cdotp\lfloor\frac{\lfloor\frac{n}{d}\rfloor}{i}\rfloor\cdotp\lfloor\frac{\lfloor\frac{m}{d}\rfloor}{i}\rfloor$

得: $\sum\limits_{d=1}^nd\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)\cdotp\lfloor\frac{\lfloor\frac{n}{d}\rfloor}{i}\rfloor\cdotp\lfloor\frac{\lfloor\frac{m}{d}\rfloor}{i}\rfloor$

令: $f(x,y)=\sum\limits_{i=1}^{x}\mu(i)\cdotp\lfloor\frac{x}{i}\rfloor\cdotp\lfloor\frac{y}{i}\rfloor$

利用数论分块可以在$O(\sqrt x+\sqrt y)$得到函数值

同样的, 已知$f(x,y)$的值, $\sum\limits_{d=1}^nd\cdotp f(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor)$可以在$O(\sqrt n+\sqrt m)$内得到答案

所以可以分块套分块, 复杂度: $O\big((\sqrt n+ \sqrt m)^2\big)$

---

6. $1\le n\le m\le10^7\qquad q$组询问

$\sum\limits_{d=1}^nd\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)\cdotp\lfloor\frac{\lfloor\frac{n}{d}\rfloor}{i}\rfloor\cdotp\lfloor\frac{\lfloor\frac{m}{d}\rfloor}{i}\rfloor=\sum\limits_{d=1}^nd\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)\cdotp\lfloor\frac{n}{i\cdotp d}\rfloor\cdotp\lfloor\frac{m}{i\cdotp d}\rfloor$

令: $k=i\cdotp d\qquad k\in[1,n]$

枚举 $k$: $\sum\limits_{k=1}^n\sum\limits_{d|k}d\cdotp\mu(\frac{k}{d})\cdotp\lfloor\frac{n}{k}\rfloor\cdotp\lfloor\frac{m}{k}\rfloor=\sum\limits_{k=1}^n\lfloor\frac{n}{k}\rfloor\cdotp\lfloor\frac{m}{k}\rfloor\sum\limits_{d|k}d\cdotp\mu(\frac{k}{d})$

不难发现: $\sum\limits_{d|k}d\cdotp\mu(\frac{k}{d})$是狄利克雷卷积形式
$f(n)=n$ 与 $\mu(n)$均为积性函数

所以$g(n)=\sum\limits_{d|n}d\cdotp \mu(\frac{n}{d})$为积性函数
若: $n=p_1^{e_1}\cdotp p_2^{e_2}\cdots p_k^{e_k}\qquad p_1, p_2,\cdots ,p_k$均为质数

那么: $g(n)=g(p_1^{e_1})\cdot g(p_2^{e_2})\times\cdots \times g(p_k^{e_k})$

考虑: $g(p^k)=\sum\limits_{d|p^k}d\cdotp\mu(\frac{p^k}{d})\qquad p$为质数
$g(p^k)=\sum\limits_{i=0}^k p^{k-i}\cdotp \mu(p^i)=p^k-p^{k-1}=\varphi(p^k)$

所以: $g(n)=\varphi(n)$

$\sum\limits_{k=1}^n\lfloor\frac{n}{k}\rfloor\cdotp\lfloor\frac{m}{k}\rfloor\sum\limits_{d|k}d\cdotp\mu(\frac{k}{d})=\sum\limits_{k=1}^n\lfloor\frac{n}{k}\rfloor\cdotp\lfloor\frac{m}{k}\rfloor\cdotp \varphi(k)$

复杂度: $O\big(n+q\cdotp (\sqrt n+\sqrt m)\big)$

---

7. $1\le n=m\le 10^{11}$

用杜教筛来求欧拉函数$\varphi(x)$前缀和$s(n)$

$s(n)=\frac{(n+1)\cdotp n}{2}-\sum\limits_{i=2}^ns(\lfloor\frac{n}{i}\rfloor)$

$\sum\limits_{d=1}^nd\cdotp(2\cdotp s(\lfloor\frac{n}{d}\rfloor)-1)\qquad$(3)
或
$\sum\limits_{k=1}^n\lfloor\frac{n}{k}\rfloor^2\cdotp \varphi(k)\qquad$(6)

复杂度: $O(n^{\frac{2}{3}})$

---

总结:

对于 $m=n$的情况
我们得到了$2$个不同的等式:
$\sum\limits_{i=1}^n\sum\limits_{j=1}^ngcd(i,j)= \begin{cases} \sum\limits_{d=1}^nd\cdotp(2\cdotp s(\lfloor\frac{n}{d}\rfloor)-1)\qquad\Big(s(n)=\sum\limits_{i=1}^n\varphi(i)\Big)\\\ \sum\limits_{k=1}^n\lfloor\frac{n}{k}\rfloor^2\cdotp \varphi(k)\\ \end{cases}$

它们之间是不是有一些联系呢?
$\begin{align} &\quad\sum\limits_{d=1}^nd\cdotp(2\cdotp s(\lfloor\frac{n} {d}\rfloor)-1)\\ &=\Big(\sum\limits_{d=1}^n2\cdotp d\cdotp s(\lfloor\frac{n} {d}\rfloor)\Big)-\frac{(n+1)\cdotp n}{2}\\ &=\Big(\sum\limits_{d=1}^n2\cdotp d\cdotp s(\lfloor\frac{n} {d}\rfloor)\Big)-\sum\limits_{d=1}^n s(\lfloor\frac{n} {d}\rfloor)\\ &=\sum\limits_{d=1}^n(2\cdotp d-1)\cdotp s(\lfloor\frac{n} {d}\rfloor)\\ &=\sum\limits_{d=1}^n(2\cdotp d-1)\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\varphi(i)\\ &=\sum\limits_{d=1}^n\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}(2\cdotp d-1)\cdotp \varphi(i)\\ &=\sum\limits_{d=1}^n\sum\limits_{d|i}(2\cdotp d-1)\cdotp \varphi(\frac{i}{d})\\ &=\sum\limits_{i=1}^n\sum\limits_{d|i}(2\cdotp d-1)\cdotp \varphi(\frac{i}{d})\\ &=\sum\limits_{i=1}^n(f\times g)(i)\\ \end{align}$

其中$\begin{cases}f(x)=\varphi(x)\\ g(x)=2\cdotp x-1 \end{cases}$

$\begin{align} &\quad\sum\limits_{k=1}^n\lfloor\frac{n}{k}\rfloor^2\cdotp \varphi(k)\\ &=\sum\limits_{k=1}^n\sum_{i=1}^{\lfloor\frac{n}{k}\rfloor}(2\cdotp i-1)\cdotp \varphi(k)\\ &=\sum\limits_{k=1}^n\sum\limits_{k|i}(2\cdotp \frac{i}{k}-1)\cdotp \varphi(k)\\ &=\sum\limits_{i=1}^n\sum\limits_{k|i}(2\cdotp \frac{i}{k}-1)\cdotp \varphi(k)\\ &=\sum\limits_{i=1}^n(g\times f)(i)\\ \end{align}$
所以这两个等式都表达的是同一个狄利克雷卷积的前缀和, 自然也就相等了