记录编号 |
286029 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[BZOJ 2820] YY的GCD |
最终得分 |
100 |
用户昵称 |
TenderRun |
是否通过 |
通过 |
代码语言 |
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;
}