Description

設d(x)為x的約數個數，給定N、M，求

Input

輸入文件包含多組測試數據。

第一行，一個整數T，表示測試數據的組數。
接下來的T行，每行兩個整數N、M。

Output

T行，每行一個整數，表示你所求的答案。

Sample Input

2
7 4
5 6

Sample Output

110
121

HINT

1<=N, M<=50000
1<=T<=50000

Solution

莫比烏斯反演
但這題更多的是套路
首先，一個神奇的東東：\(d(nm)= \sum_{i|n}\sum_{j|m}[gcd(i,j)=1]\)
這個東西是個套路，證明的話可以百度，用的確實多

#include<bits/stdc++.h>
#define ll long long
const int MAXN=50000+10;
int T,cnt,prime[MAXN],mu[MAXN],s[MAXN],f[MAXN];
bool vis[MAXN];
template<typename T> inline void read(T &x)
{
    T data=0,w=1;
    char ch=0;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')w=-1,ch=getchar();
    while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
    x=data*w;
}
template<typename T> inline void write(T x,char c='\0')
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
    if(c!='\0')putchar(c);
}
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline void init()
{
    memset(vis,1,sizeof(vis));
    vis[0]=vis[1]=0;
    mu[1]=1;
    for(register int i=2;i<MAXN;++i)
    {
        if(vis[i])
        {
            prime[++cnt]=i;
            mu[i]=-1;
        }
        for(register int j=1;j<=cnt&&i*prime[j]<MAXN;++j)
        {
            vis[i*prime[j]]=0;
            if(i%prime[j])mu[i*prime[j]]=-mu[i];
            else break;
        }
    }
    for(register int i=1;i<MAXN;++i)s[i]=s[i-1]+mu[i];
    for(register int k=1;k<MAXN;++k)
        for(register int i=1;;)
        {
            if(i>k)break;
            int j=k/(k/i);
            f[k]+=(k/i)*(j-i+1);
            i=j+1;
        }
}
inline ll solve(int N,int M)
{
    ll res=0;
    for(register int i=1;;)
    {
        if(i>min(N,M))break;
        int j=min(N/(N/i),M/(M/i));
        res+=(ll)f[N/i]*(ll)f[M/i]*(ll)(s[j]-s[i-1]);
        i=j+1;
    }
    return res;
}
int main()
{
    read(T);
    init();
    while(T--)
    {
        int N,M;
        read(N);read(M);
        write(solve(N,M),'\n');
    }
    return 0;
}

【刷題】BZOJ 3994 [SDOI2015]約數個數和