【莫比乌斯反演】BZOJ2920-YY的GCD

【题目大意】

给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对。

【思路】

太神了这道题……蒟蒻只能放放题解：戳，明早再过来看看还会不会推导过程……

实用的结论：

嗯……

/**************************************************************
    Problem: 2820
    Language: C++
    Result: Accepted
    Time:4164 ms
    Memory:196600 kb
****************************************************************/

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int INF=0x7fffffff;
const int MAXN=+;
typedef long long ll;
int miu[MAXN],g[MAXN],prime[MAXN],pnum=;
ll sum[MAXN];
int N,M;

void get_miu(int maxn)
{
    miu[]=;
    g[]=;
    sum[]=sum[]=;
    for (int i=;i<maxn;i++) miu[i]=-INF;
    for (int i=;i<maxn;i++)
    {
        if (miu[i]==-INF)
        {
            miu[i]=-;
            prime[++pnum]=i;
            g[i]=;
        }
        for (int j=;j<=pnum;j++)
        {
            if (i*prime[j]>=maxn) break;
            if (i%prime[j]==)
            {
                miu[i*prime[j]]=;
                g[i*prime[j]]=miu[i];
            }
                else
                {
                    miu[i*prime[j]]=-miu[i];
                    g[i*prime[j]]=miu[i]-g[i];
                }
        }
        sum[i]=sum[i-]+g[i];
    }
} 

void get_ans()
{
    ll ans=;
    scanf("%d%d",&N,&M);
    if (N>M) swap(N,M);
    int pos;
    for (int t=;t<=N;t=pos+)
    {
        pos=min(N/(N/t),M/(M/t));
        ans+=(ll)(sum[pos]-sum[t-])*(N/t)*(M/t);
    }
    printf("%lld\n",ans);
}

int main()
{
    get_miu(MAXN);
    int T;
    scanf("%d",&T);
    while (T--) get_ans();
    return ;
}