BZOJ 2820 YY的GCD 莫比乌斯反演

由于这道题比较鬼畜，而且公式巨难打出，所以我粘了两页PoPoQQQ的PPT orz，并对其中一些部分解释一下

最下面的公式就是换了一种枚举的方式，其中u括号里的东西就是把d变成T/p，其中p|T的条件只要带回到原表达式中即可理解。

这里注意这里的每个前缀和代表的意义是对于每一个T满足条件的Σu之和。
下面附上我的代码

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cstring>
#include<string>
#include<iostream>
#include<iomanip>
#include<algorithm>
using namespace std;
#define M 10001000
bool pd[M];
int u[M];
long long zs[1000000];
int top=0;
int sum[M];
void shai()
{
    u[1]=1;
    for(int i=2;i<M;i++)
    {
        if(!pd[i])
        {
            zs[++top]=i,
            u[i]=-1;
        }
        for(int j=1;(long long)zs[j]*i<M;j++)
        {
            pd[zs[j]*i]=true;
            if(i%zs[j]==0) break;
            u[i*zs[j]]=-u[i];
        }
    }
    for(int i=1;i<=top;i++)
        for(int j=zs[i];j<M;j+=zs[i])
            sum[j]+=u[j/zs[i]];
    for(int i=1;i<M;i++) sum[i]+=sum[i-1];
}
long long get_ans(int n,int m)
{
    long long re=0;
    int last;
    for(int i=1;i<=n && i<=m;i=last+1)
    {
        last=min(n/(n/i),m/(m/i));
        re+=(long long)(sum[last]-sum[i-1])*(n/i)*(m/i);
    }
    return re;
}
int main()
{
    int T;
    scanf("%d",&T);
    shai();
    while(T--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        printf("%lld\n",get_ans(n,m));
    }
    return 0;
}