由於這道題比較鬼畜，而且公式巨難打出，所以我粘了兩頁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;
}