BZOJ 2820: YY的GCD

2820: YY的GCD

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 1705  Solved: 908
[Submit][Status][Discuss]

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

Source

分析：

我们很容易就可以得出一下这个复杂度爆炸的式子：

$sum_ {prime(p)}^{n} sum_ {d=1}^{n/p} μ(d) leftlfloor frac{n}{pd} ight floor leftlfloor frac{m}{pd} ight floor$

但是有10000个case，这个复杂度明显爆炸...我们设$k=pd$

可以得到一下式子：

$sum_ {k=1}^{n} leftlfloor frac{n}{k} ight floor leftlfloor frac{m}{k} ight floor sum_ {prime(p)  pmid k} μ(frac {k}{p})$

然后对于后面的μ我们可以预处理前缀和，这样单个case复杂度大概可以降到$sqrt(n)+sqrt(m)$...

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
//by NeighThorn
#define int long long
using namespace std;

const int maxn=10000000+5;

int n,m,cas,cnt,f[maxn],mu[maxn],pri[maxn],vis[maxn];

inline void prework(void){
	mu[1]=1;
	for(int i=2;i<=10000000;i++){
		if(!vis[i])
			pri[++cnt]=i,mu[i]=-1;
		for(int j=1;j<=cnt&&i*pri[j]<=10000000;j++){
			vis[i*pri[j]]=1;
			if(i%pri[j]==0){
				mu[i*pri[j]]=0;
				break;
			}
			mu[i*pri[j]]=-mu[i];
		}
	}
	for(int i=1;i<=cnt;i++)
		for(int j=1;pri[i]*j<=10000000;j++)
			f[pri[i]*j]+=mu[j];
	for(int i=1;i<=10000000;i++)
		f[i]+=f[i-1];
}

signed main(void){
	scanf("%lld",&cas);prework();
	while(cas--){
		scanf("%lld%lld",&n,&m);
		if(n>m)
			swap(n,m);
		int ans=0,r;
		for(int i=1;i<=n;i=r+1){
			r=min(n/(n/i),m/(m/i));
			ans+=(f[r]-f[i-1])*(n/i)*(m/i);
		}
		printf("%lld
",ans);
	}
	return 0;
}

　　

---

By NeighThorn