C++ 代码

#include<iostream>
using namespace std;

typedef long long LL;

LL C(int a){
    return (LL)a * (a - 1) * (a - 2) / 6;
}

int gcd(int a, int b){
    return b ? gcd(b, a % b) : a;
}

int main(){
    int n, m;
    scanf("%d %d", &n, &m);
    n++,m++;
    LL res = C(n*m) - n * C(m) - m * C(n);
    for(int i = 1; i < n; i++)
        for(int j = 1; j < m; j++){
            res -= (LL)2 * (gcd(i, j) - 1) * (n - i) * (m - j);
        }
    printf("%lld", res);
    return 0;
}