『题解』Luogu-P5176 公约数

P5176 公约数

Description

• 多测，$T$ 组数据。

• 给定 $n,m,p$，请求出
  $$\left[\sum _{i=1}^n\sum _{j=1}^m\sum _{k=1}^p\gcd (ij,ik,jk)\cdot \gcd (i,j,k)\cdot \left(\frac{\gcd (i,j)}{\gcd (i,k)\cdot \gcd (j,k)}+\frac{\gcd (i,k)}{\gcd (i,j)\cdot \gcd (j,k)}+\frac{\gcd (j,k)}{\gcd (i,j)\cdot \gcd (i,k)}\right)\right]\mathrm{m}\mathrm{o}\mathrm{d}(10^9+7)$$

• $T\le 10^3,n,m,p,\le 2\times 10^7$。

Solution

引理：
$$\gcd (ij,ik,jk)=\frac{\gcd (i,j)\gcd (i,k)\gcd (j,k)}{\gcd (i,j,k)}$$

证明：

只需证任意一种质因数 $p$ 的次数相等即可，因为若两个数的任意质因数的次数都相等，则这两个数相等。

设 $i,j,k$ 分别含 $a,b,c$ 个质因数 $p$，不妨设 $a\le b\le c$。

则 $\gcd (ik,ik,jk)$ 中的质因数 $p$ 的次数为 $\min (a+b,a+c,b+c)=a+b$；

$\frac{\gcd (i,j)\gcd (i,k)\gcd (j,k)}{\gcd (i,j,k)}$ 中的质因数 $p$ 的次数为 $\min (a,b)+\min (a,c)+\min (b,c)-\min (a,b,c)=a+a+b-a=a+b$。

所以这两个数的质因数 $p$ 的次数相等。

---

所以
$$\begin{array}{rl}ans& =\sum _{i=1}^n\sum _{j=1}^m\sum _{k=1}^p\frac{\gcd (i,j)\gcd (i,k)\gcd (j,k)}{\gcd (i,j,k)}\cdot \gcd (i,j,k)\cdot \left(\frac{\gcd (i,j)}{\gcd (i,k)\cdot \gcd (j,k)}+\frac{\gcd (i,k)}{\gcd (i,j)\cdot \gcd (j,k)}+\frac{\gcd (j,k)}{\gcd (i,j)\cdot \gcd (i,k)}\right)\\ & =\sum _{i=1}^n\sum _{j=1}^m\sum _{k=1}^p\gcd (i,j{)}^2+\gcd (i,k{)}^2+\gcd (j,k{)}^2\\ & =p\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j{)}^2+m\sum _{i=1}^n\sum _{k=1}^p\gcd (i,k{)}^2+n\sum _{j=1}^m\sum _{k=1}^p\gcd (j,k{)}^2\end{array}$$
设
$$f(n,m)=\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j{)}^2$$
则
$$ans=pf(n,m)+mf(n,p)+nf(m,p)$$
算 $f(n,m)$（不妨设 $n\le m$）
$$\begin{array}{rl}f(n,m)& =\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j{)}^2\\ & =\sum _{d=1}^n\sum _{i=1}^n\sum _{j=1}^m{d}^2[\gcd (i,j)=d]\\ & =\sum _{d=1}^n{d}^2\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }[\gcd (i,j)=1]\\ & =\sum _{d=1}^n{d}^2\sum _{k=1}^{\lfloor \frac{n}{d}\rfloor }\mu (k)\lfloor \frac{n}{dk}\rfloor \lfloor \frac{m}{dk}\rfloor \\ & =\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor \sum _{k\mid T}\mu (k){\left(\frac{T}{k}\right)}^2\\ & =\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor (\mu \ast {\mathrm{Id}}_2)(T)\end{array}$$
整除分块即可。

考虑对积性函数 $g=\mu \ast {\mathrm{Id}}_2$ 进行线性筛：

• $i\in \mathbb{P}$：$g(i)=p^2-1$；
• $p_j\nmid i$：$g(i\cdot p_j)=g(i)\cdot g(p_j)$；
• $p_j\mid i$：显然只有原来对 $i$ 有贡献的 $\mu (k)$ 现在还会有贡献，唯一不同的是后面的 ${\left(\frac{T}{k}\right)}^2$ 多乘了个 $p_j^2$，即 $g(i\cdot p_j)=g(i)\cdot p_j^2$。

时间复杂度为 $\mathrm{O}(n+T\sqrt{n})$。

Code

// 18 = 9 + 9 = 18.
#include <iostream>
#include <cstdio>
#define Debug(x) cout << #x << "=" << x << endl
typedef long long ll;
using namespace std;

const int MAXN = 2e7 + 5;
const int N = 2e7;
const int MOD = 1e9 + 7;
typedef int arr[MAXN];

int add(int a, int b) {return (a + b) % MOD;}
int sub(int a, int b) {return (a - b + MOD) % MOD;}
int mul(int a, int b) {return (ll)a * b % MOD;}

arr p, g, sum;
bool vis[MAXN];

void pre()
{
	g[1] = sum[1] = 1;
	for (int i = 2; i <= N; i++)
	{
		if (!vis[i])
		{
			p[++p[0]] = i;
			g[i] = mul(i, i) - 1;
		}
		sum[i] = add(sum[i - 1], g[i]);
		for (int j = 1; j <= p[0] && i * p[j] <= N; j++)
		{
			vis[i * p[j]] = true;
			if (i % p[j] == 0)
			{
				g[i * p[j]] = mul(g[i], mul(p[j], p[j]));
				break;
			}
			g[i * p[j]] = mul(g[i], g[p[j]]);
		}
	}
}

int GetSum(int l, int r)
{
	return sub(sum[r], sum[l - 1]);
}

int f(int n, int m)
{
	if (n > m)
	{
		swap(n, m);
	}
	int res = 0;
	for (int l = 1, r; l <= n; l = r + 1)
	{
		int k1 = n / l, k2 = m / l;
		r = min(n / k1, m / k2);
		res = add(res, mul(mul(k1, k2), GetSum(l, r)));
	}
	return res;
}

int main()
{
	pre();
	int t;
	scanf("%d", &t);
	while (t--)
	{
		int n, m, p;
		scanf("%d%d%d", &n, &m, &p);
		printf("%d\n", add(add(mul(p, f(n, m)), mul(m, f(n, p))), mul(n, f(m, p)))); 
	}
	return 0;
}