[bzoj2820]：YY的GCD（莫比乌斯反演）

传送门
反正就是推一大堆式子，然后就推出来了。。。
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;
}