记录编号 286029 评测结果 AAAAAAAAAA
题目名称 [BZOJ 2820] YY的GCD 最终得分 100
用户昵称 GravatarTenderRun 是否通过 通过
代码语言 C++ 运行时间 18.072 s
提交时间 2016-07-26 09:05:11 内存使用 124.29 MiB
显示代码纯文本
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxn=10000010;
int prime[maxn],cnt;
int mu[maxn],sum[maxn];
bool check[maxn];

void Prepare(){
	mu[1]=1;
	for(int i=2;i<=10000000;i++){
		if(!check[i]){
			prime[++cnt]=i;
			mu[i]=-1;
		}
		for(int j=1;j<=cnt;j++){
			if(prime[j]*i>10000000)break;
			check[prime[j]*i]=true;
			if(i%prime[j]==0){
				mu[prime[j]*i]=0;
				break;
			}
			mu[prime[j]*i]=mu[i]*-1;
		}
	}
	for(int i=1;i<=10000000;i++)
		for(int j=1;j<=cnt&&prime[j]*i<=10000000;j++)
			sum[i*prime[j]]+=mu[i];
	for(int i=1;i<=10000000;i++)sum[i]+=sum[i-1];
}

int T;
int a,b;
long long C(int n,int m){
	if(n>m)swap(n,m);
	long long ret=0;int p;
	for(int i=1;i<=n;i=p+1){
		p=min(n/(n/i),m/(m/i));
		ret+=1ll*(sum[p]-sum[i-1])*(n/i)*(m/i);
	}
	return ret;	
}

int main(){
	freopen("YYnoGCD.in","r",stdin);
	freopen("YYnoGCD.out","w",stdout);
	Prepare();
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&a,&b);
		printf("%lld\n",C(a,b));
	}
	return 0;
}