很明显题目要求的东西可以写成$sum_{i=1}^{n}sum_{j=1}^m gcd(i,j)*2-1$（一点都不明显）
如果直接枚举肯定爆炸
那么我们设$f[i]$表示存在公因数$i$的数的对数
然而$i$并不一定是这几对数的最大公因数
那么怎么办呢？考虑容斥
以$i$为最大公因数的数的对数，就是有$i$为公因数的数，减去最大公因数为$2*i$的数，减去为$3*i$的数……
那么这个就可以一波容斥求出来了
时间复杂度为$O(nlogn)$
1 //minamoto 2 #include<cstdio> 3 #define ll long long 4 const int N=1e5+5; 5 ll f[N],ans;int n,m; 6 int main(){ 7 scanf("%d%d",&n,&m); 8 if(n>m) n^=m^=n^=m; 9 for(int i=n;i;--i){ 10 f[i]=1ll*(n/i)*(m/i); 11 for(int j=i<<1;j<=n;j+=i) f[i]-=f[j]; 12 ans+=((i<<1)-1)*f[i]; 13 } 14 printf("%lld ",ans); 15 return 0; 16 }