洛谷P2257 YY的GCD

题目链接：https://www.luogu.org/problemnew/show/P2257

之前就听队友说过莫比乌斯极其阔怕，然后就找了些莫比乌斯的博客来看，看了半天发现这个东西很强大，可以把很多难以处理的问题简化，就是数学公式有点难推，这题算是自己做的第一道莫比乌斯吧，虽然还不是很明白(dalao都是些什么显然已经可得，看的有点云里雾里)，留篇博客放这，以后再来看看吧

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define N 10000000
#define LL long long
bool vis[N+10];
int a[N],s[N+10],mu[N+10];
void init()
{
    int i,j,tot=0;
    memset(vis,false,sizeof(vis));
    mu[1]=1;
    for(i=2;i<=N;i++)
    {
        if(!vis[i])
        {
            a[tot++]=i;
            mu[i]=-1;
        }
        for(j=0;j<tot&&i*a[j]<=N;j++)
        {
            vis[i*a[j]]=true;
            if(i%a[j]==0)
                break;
            mu[i*a[j]]=-mu[i];
        }
    }
    for(i=0;i<tot;i++)
    for(j=1;j*a[i]<=N;j++)
        s[j*a[i]]+=mu[j];
    for(i=1;i<=N;i++)
        s[i]+=s[i-1];
}
int main()
{
    int i,j,k,t,n,m,l,r;
    LL ans;
    init();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        if(n>m)
            swap(n,m);
        ans=0;
        for(l=1;l<=n;l=r+1)
        {
            r=min(n/(n/l),m/(m/l));
            ans=ans+(LL)(n/l)*(m/l)*(s[r]-s[l-1]);
            l=r+1;
        }
        printf("%lld\n",ans);
    }
    return 0;
}