題目連結
題意：
對於一個n行m列的方陣 第i行第j列的值是gcd(i, j) * 2 - 1 求方陣和 n, m <= 1e5

這道題很妙了
f[x] 為滿足 gcd(i, j) = x 的數對(i, j)的個數
g[x] 為i，j因數包括x的數對(i, j)的個數
顯然g[x] = ( n / x ) * ( m / x )
而f[x] + f[2x] + f[3x] + … = g[x]
所以f[x] = g[x] - f[2x] - f[3x] - …
然後從n到1計算f即可
複雜度O(n ln n)

#include<iostream>
#include
 
<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int N = 1e5 + 5;
long long f[N], ans;
int n, m;
int main()
{
 scanf("%d%d", &n, &m);
 if(n > m) swap(n, m);
 for(int i = n; i >= 1; --i){
  f[i] = 1ll * (n / i) * (m / i);
  for
 
(int j = i << 1; j <= n; j += i) f[i] -= f[j];
  ans +=  (f[i] * ((i << 1) - 1));
 }
 printf("%lld", ans);
    return 0;
}