3309: DZY Loves Math 莫比乌斯反演

令$F(i)$表示$i$所含质因子最大幂指数，$f(i)$表示$gcd(x,y)=i$的数对$(x,y)$的数量，然后莫比乌斯反演得到下式：

$$Ans=\sum_{T=1}^n\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\sum_{d\mid T}F(d)\mu(\frac{T}{d})$$
令 $G(T)=\sum_{d\mid T}F(d)\mu(\frac{T}{d})$，然后我们大力讨论一下。
首先设 $T=P_1^{x_1}\times P_2^{x_2}\times..\times P_k^{x_k},d==P_1^{y_1}\times P_2^{y_2}\times..\times P_k^{y_k}$。
要使 $\mu(\frac{T}{d})\neq0$，必须满足 $y_i=x_i$或者 $y_i=x_i-1$。
令 $a=Max_{1\leq i\leq k}(x_i)$，则 $F(d)=a$或 $a-1$。
假设 $T$中一共有$q$个质因子的指数为 $a$。

一、当$q=k$时，
① $f(d)=a$时

$$\sum_{f(d)=a}f(d)\mu(\frac{T}{d})=a\times\sum_{f(d)=a}\mu(\frac{T}{d})=a\times\sum_{i=1}^k(-1)^{k-i}C_k^i=a\times (-1)^{k+1}$$
② $f(d)=a-1$时
$$\sum_{f(d)=a-1}f(d)\mu(\frac{T}{d})=(a-1)\times\sum_{f(d)=a-1}\mu(\frac{T}{d})=(a-1)\times(-1)^k$$
所以当 $q=k$时， $G(T)=a\times(-1)^{k+1}+(a-1)\times(-1)^k=(-1)^{k+1}$
二、当 $q<k$时
① $f(d)=a$时
$$\sum_{f(d)=a}f(d)\mu(\frac{T}{d})=a\times\sum_{f(d)=a}\mu(\frac{T}{d})=a\times\sum_{i=1}^q(-1)^{q-i}C_q^i\times\sum_{j=0}^{k-q}(-1)^jC_{k-q}^j=0$$
② $f(d)=a-1$时
$$\sum_{f(d)=a-1}f(d)\mu(\frac{T}{d})=(a-1)\times\sum_{f(d)=a-1}\mu(\frac{T}{d})=(a-1)\times(-1)^q\times\sum_{j=0}^{k-q}(-1)^jC_{k-q}^j=0$$
所以当 $q<k$时， $G(T)=0$
知道这些，我们可以线筛的时候记录一下最小质因子出现的次数和除去最小质因子 $p_1^{x_1}$后的数，然后判断幂指数是否相等。

#include<iostream>
#include<cstdio>
#define ll long long 
#define MAXN 10000005
using namespace std;
int prime[MAXN],t[MAXN],last[MAXN];
ll g[MAXN];
bool flag[MAXN];
inline ll read()
{
    ll a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void prepare()
{
    for (int i=2;i<MAXN;i++)
    {
        if (!flag[i]) prime[++prime[0]]=i,last[i]=t[i]=g[i]=1;
        for (int j=1;j<=prime[0]&&i*prime[j]<MAXN;j++)
        {
            flag[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                last[i*prime[j]]=last[i];
                t[i*prime[j]]=t[i]+1;
                if (last[i*prime[j]]==1) g[i*prime[j]]=1;
                else g[i*prime[j]]=(t[last[i*prime[j]]]==t[i*prime[j]]?-g[last[i*prime[j]]]:0);
                break;
            }
            last[i*prime[j]]=i;
            t[i*prime[j]]=1;
            g[i*prime[j]]=(t[i]==1?-g[i]:0);
        }
    }
    for (int i=1;i<MAXN;i++) g[i]+=g[i-1];
}
inline ll solve(ll n,ll m)
{
    ll ans=0;
    for (int i=1,pos;i<=n;i=pos+1)
    {
        pos=min(n/(n/i),m/(m/i));
        ans+=(g[pos]-g[i-1])*(n/i)*(m/i);
    }
    return ans;
}
int main()
{
    prepare();
    int testcase=read();
    while (testcase--)
    {
        ll n=read(),m=read();
        if (n>m) swap(n,m);
        printf("%lld\n",solve(n,m));
    }
    return 0;
}