【BZOJ 2820】 YY的GCD

最新推荐文章于 2022-01-28 15:38:33 发布

Regina8023

最新推荐文章于 2022-01-28 15:38:33 发布

阅读量1.3k

点赞数 1

分类专栏： ~Mathematic~ 文章标签： OI BZOJ 莫比乌斯函数

本文链接：https://blog.csdn.net/Regina8023/article/details/43877249

版权

~Mathematic~ 专栏收录该内容

20 篇文章 0 订阅

订阅专栏

2820: YY的GCD

Time Limit: 10 Sec Memory Limit: 512 MB
Submit: 807 Solved: 404
[ Submit][ Status]

Description

神犇YY虐完数论后给傻×kAc出了一题

给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对

kAc 这种傻×必然不会了，于是向你来请教 ……

多组输入

Input

第一行一个整数T 表述数据组数

接下来T行，每行两个正整数，表示N, M

Output

T 行，每行一个整数表示第i组数据的结果

Sample Input

2
10 10
100 100

Sample Output

30
2791

HINT

T = 10000

N, M <= 10000000

莫比乌斯函数。

iwtwiioi的题解非常详细orz

大概说一下做法，求gcd(x,y)=k的(x,y)的对数就是【BZOJ 1101】，那么这道题可以枚举k来求，但是会TLE，那么我们考虑公式的变形。（详见iwtwiioi的题解。。）

接下来的关键就是求g[x]函数：根据线性筛的原理，为了方便，把g[x]变成求g[kp]。

k是线性筛的时候当前要处理的数，p是枚举已经筛出来的质数(一定<k)。

就像求莫比乌斯函数一样，分三种情况：

1.k为质数，g[k]=1

2.k%p=0，说明p的指数>=2，然后根据g[x]函数的计算方法来计算此时的g[kp]，分p'=p和p'!=p两种情况。

得出g[kp]=mu[k]

3.k%p!=0，说明p的指数=1，然后根据g[x]函数的计算方法来计算此时的g[kp]，分p'=p和p'!=p两种情况。

得出g[kp]=-g[k]

于是在求莫比乌斯函数的同时就把g数组求好了，按照原来的分块方法计算即可~

#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#define LL long long
#define M 10000005
using namespace std;
int check[M],mu[M],p[M],sum[M],g[M];
void Getmobius()
{
	mu[1]=1;
	int tot=0;
	for (int i=2;i<=M;i++)
	{
		if (!check[i])
		{
			p[++tot]=i;
			mu[i]=-1;
			g[i]=1;
		}
		for (int j=1;j<=tot;j++)
		{
			if (p[j]*i>M) break;
			check[p[j]*i]=1;
			if (i%p[j]==0)
			{
				mu[p[j]*i]=0;
				g[p[j]*i]=mu[i];
				break;
			}
			mu[p[j]*i]=-mu[i];
			g[p[j]*i]=mu[i]-g[i];
		}
	}
	for (int i=1;i<=M;i++)
		sum[i]=sum[i-1]+g[i];
}
LL Calc(int a,int b)
{	
	LL ans=0LL;
	int pos;
	for (int d=1;d<=min(a,b);d=pos+1)
	{
		pos=min(a/(a/d),b/(b/d));
		ans+=(LL)(sum[pos]-sum[d-1])*(a/d)*(b/d);
	}
	return ans;
}
int main()
{
        Getmobius();
	int T;
	scanf("%d",&T);
	while (T--)
	{
		int n,m;
		scanf("%d%d",&n,&m);
		printf("%lld\n",Calc(n,m));
	}
	return 0;
}