枚举倍数的一种灵活的变形：$g(d)={\sum }_{d|i}^{n}f(i)={\sum }_{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i\cdot d)$
很显然，但有时能发挥大作用。
其实本质还是要理解西格玛究竟是在算什么

解析

一个重要的前置结论：
$$d(i\cdot j)=\sum _{x|i,y|j}[\gcd (x,y)=1]$$
证明：
尝试把所有互质的数对 $(x,y)$ 和 $i\cdot j$ 的所有因子建立双射关系。首先，对于只有 $i,j$ 一方有的质因子，直接归给有该项质因子的那一方对应的 $x$ 或 $y$。对于双方都有的质因子 $p$，若 $i$ 的对应次数为 $a$，$j$ 的对应次数为 $b$，如果 $i\cdot j$ 的某个因子 $p$ 的次数 $c$ 满足 $c\le a$，就当成 $x$ 含有 $p^c$ 的因子；若 $c>a$，就当成 $y$ 含有 $p^{c-a}$ 次方的因子。这样我们就建立起了双射关系，且所有的 $x,y$ 都是互质的，证毕。

我们再化一下这个 $d$ 函数：
$$\begin{gathered} d(i,j)=\sum _{x|i}\sum _{y|j}[\gcd (i,j)=1] \\ =\sum _{x|i}\sum _{y|j}\sum _{k|\gcd (x,y)}\mu (k) \\ =\sum _{k|i,k|j}\mu (k)\sum _{k|x,x|i}\sum _{k|y,y|j}1 \\ =\sum _{k|i,k|j}\mu (k)\sum _{x|\frac{i}{k}}^{\frac{i}{k}}\sum _{y|\frac{j}{k}}^{\frac{j}{k}}1 \\ =\sum _{k|i,k|j}\mu (k)d(\frac{i}{k})d(\frac{j}{k}) \end{gathered}$$
带回原式：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}d(ij)=\sum _{i=1}^n\sum _{j=1}^m\sum _{k|i,k|j}\mu (k)d(\frac{i}{k})d(\frac{j}{k}) \\ =\sum _k\mu (k)\sum _{k|i}^n\sum _{k|j}^{m}d(\frac{i}{k})d(\frac{j}{k}) \\ =\sum _k\mu (k)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }d(i)d(j) \\ =\sum _k\mu (k)(\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }d(i))\cdot (\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }d(j)) \end{gathered}$$
线性筛出 $\mu ,d$，预处理出前缀和，整除分块即可。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define debug(...) fprintf(stderr,__VA_ARGS__)
inline ll read(){ll x(0),f(1);char c=getchar();while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}while(isdigit(c)) {x=(x<<1)+(x<<3)+c-'0';c=getchar();}return x*f;
}
const int N=1e5+100;
int n,m;int p[N],tot,vis[N],mu[N],msum[N],d[N];
ll dsum[N];
void init(int n){mu[1]=1;d[1]=1;for(int i=2;i<=n;i++){if(!vis[i]) p[++tot]=i,mu[i]=-1,d[i]=2;for(int j=1;j<=tot&&p[j]<=n/i;j++){vis[i*p[j]]=1;if(i%p[j]==0){mu[i*p[j]]=0;d[i*p[j]]=d[i]+d[i]-d[i/p[j]];break;}mu[i*p[j]]=-mu[i];d[i*p[j]]=d[i]*2;}}for(int i=1;i<=n;i++){msum[i]=msum[i-1]+mu[i];dsum[i]=dsum[i-1]+d[i];}return;
}signed main() {
#ifndef ONLINE_JUDGE//freopen("a.in","r",stdin);//freopen("a.out","w",stdout);
#endifinit(5e4);int T=read();while(T--){n=read();m=read();ll ans(0);for(int l=1,r;l<=min(n,m);l=r+1){r=min(n/(n/l),m/(m/l));ans+=(msum[r]-msum[l-1])*dsum[n/l]*dsum[m/l];}printf("%lld\n",ans);}return 0;
}
/*
5
5
8
1
46
3564
*/