3505: [Cqoi2014]数三角形

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
ll ans,t,c[4];
ll C(int n,int k){
    c[0]=1;
    for(int i=1;i<=k;i++)  c[i]=c[i-1]*(n-i+1)/i;
    return c[k];
}
int n,m;
void solve(){
    ans+=C(n*m,3);//总数 
    ans-=m*C(n,3);//减去同一行 
    ans-=n*C(m,3);//减去同一列 
    for(int i=1;i<n;i++){//再减去在斜直线上的 
        for(int j=1;j<m;j++){
            t=__gcd(i,j)+1;
            if(t>2) ans-=(t-2)*(n-i)*(m-j)*2;
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);n++;m++;
    solve();
    printf("%lld\n",ans);
    return 0;
}