【LG-P2257】 YY的GCD

最新推荐文章于 2024-11-12 12:29:51 发布

伍叁壹_

最新推荐文章于 2024-11-12 12:29:51 发布

阅读量63

点赞数

分类专栏：题解数学文章标签：概率论算法图论

本文链接：https://blog.csdn.net/ez_gsn/article/details/123547681

版权

题解同时被 2 个专栏收录

63 篇文章 1 订阅

订阅专栏

数学

16 篇文章 0 订阅

订阅专栏

P2257 YY的GCD

给定 $N$ , $M$ ，求 $\leq x \leq N$ ， $\leq y \leq M$ 且 $\gcd(x, y)$ 为质数的 $(x, y)$ 有多少对。

推柿子：（ $n\leq m$ ）

$\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=k]\ k\in prime$

（常规操作，同除 $k$ ：）

$=\sum_{k=1,\ k\in prime}^n\sum_{i=1}^{\left\lfloor\frac{n}{k}\right\rfloor}\sum_{j=1}^{\left\lfloor\frac{m}{k}\right\rfloor}[\gcd(i,j)=1]$

（根据莫比乌斯函数的性质，有：）

$=\sum_{k=2,\ k\in prime}^n\sum_{i=1}^{\left\lfloor\frac{n}{k}\right\rfloor}\sum_{j=1}^{\left\lfloor\frac{m}{k}\right\rfloor}\sum_{d|\gcd(i,j)} \mu(d)$

（因为 $d$ 同为 $i$ 和 $j$ 的因数，所以反过来枚举 $d$ 的倍数，即直接计算其倍数的数量即可：）

$=\sum_{k=1,\ k\in prime}^n\sum_{d=1}^{\left\lfloor\frac{n}{k}\right\rfloor} \mu(d) \times\left\lfloor\frac{n}{kd}\right\rfloor\times\left\lfloor\frac{m}{kd}\right\rfloor$

（此时设 $T = k d$ ，有：）

$=\sum_{T=1}^n\left\lfloor\frac{n}{T}\right\rfloor\times\left\lfloor\frac{m}{T}\right\rfloor\sum_{k=1,\ k\in prime}^n \mu(\frac{T}{k})$

到这里之后，对于 $\mu$ 的值可以线性筛出来，然后处理出一个数组 $f_i=\sum\limits_{k=1,\ k\in prime}^n \mu(\frac{i}{k})$ 。

再求出 $f$ 的前缀和，最后柿子的前半部分可以用数论分块做，二者结合即可。

CODE：

#include<bits/stdc++.h>
using namespace std;

#define ll long long
#define maxn 10000005

int t, n, m;
ll mu[maxn + 5], tot, prm[maxn + 5]; 
ll f[maxn + 5], sum[maxn + 5];
bool vis[maxn + 5];

inline void pre()
{
	mu[1] = 1;
	for(int i = 2; i <= maxn; ++i)
	{
		if(!vis[i])
			prm[++tot] = i, mu[i] = -1, vis[i] = 1;
		for(int j = 1; j <= tot and i * prm[j] <= maxn; ++j)
		{
			vis[i * prm[j]] = 1;
			if(i % prm[j] == 0) break;
			mu[i * prm[j]] = mu[i] * mu[prm[j]];
		}
	}
	for(int i = 1; i <= maxn; ++i)
		for(int j = 1; j <= tot and i * prm[j] <= maxn; ++j)
			f[i * prm[j]] += mu[i];
	for(int i = 1; i <= maxn; ++i)
		sum[i] = sum[i - 1] + f[i];
}

inline ll slv(int a, int b)
{
	long long ans = 0;
	for(int l = 1, r; l <= a; l = r + 1)
	{
		r = min(a / (a / l), b / (b / l));
		ans += (long long)(sum[r] - sum[l - 1]) * (long long)(a / l) * (long long)(b / l);
	}
	return ans;
}

signed main()
{
	scanf("%d", &t);
	pre();
	while(t--)
	{
		scanf("%d %d", &n, &m);
		if(n > m) swap(n, m);
		printf("%lld\n", slv(n, m));
	}
	return 0;
}