传送门
反正就是推一大堆式子,然后就推出来了。。。
http://hzwer.com/6142.html
膜hzwer
代码:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define ll long long
#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
using namespace std;
inline int read(){
int x=0,f=1;char ch;
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
return x*f;
}
const int N=1e7+5;
int prime[N>>1],vis[N],cnt,mu[N];
ll f[N];
inline void init(int n){
mu[1]=1;
for(int i=2;i<=n;++i){
if(!vis[i])prime[++cnt]=i,mu[i]=-1;
for(int j=1;j<=cnt&&i*prime[j]<=n;++j){
vis[i*prime[j]]=1;
if(i%prime[j])mu[i*prime[j]]=-mu[i];
else{
mu[i*prime[j]]=0;
break;
}
}
}
for(int i=1;i<=cnt;++i){
int p=prime[i];
for(int j=1;j*p<=n;j++)f[j*p]+=mu[j];
}
for(int i=1;i<=n;++i)f[i]+=f[i-1];
}
int T,n,m;
int main(){
init(1e7);
T=read();
while(T--){
n=read();m=read();
if(n>m)swap(n,m);
ll ans=0;
for(int i=1,j;i<=n;i=j+1){
j=min(n/(n/i),m/(m/i));
ans+=(f[j]-f[i-1])*(n/i)*(m/i);
}
printf("%lld\n",ans);
}
return 0;
}