显示代码纯文本
#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;
}