zoj 3647 Gao the Grid（矩阵对角线三点共线）

题意：给你一个矩阵，任选三点，问能形成多少个三角形？

如 n=3,m=4

解析：关键是求斜线上的三点个数： 枚举矩阵，有个公式，即一个（a*b）矩阵，在主对角线上，两端点+中间一点形成三点共线的个数为 (gcd（a，b）-1）*2

对于整x,y(可以看做以左下角为原点的坐标)
x/y=w/h
x*h=y*w
y=x*h/w
这儿可以化简一下
h=h1*d
w=w1*d
d=gcd(h,w)
y=x*h1/w1
那么x必须是w1 的倍数了。
而且y<=h
那么x*h1/w1<=h
x<=h/h1*w1=d*w1=w
也就是要求一下x<=w里面有几个 x%w1==0
那么x可以是w1,2*w1,3*w1,4*w1...,(d-1)*w1
d=gcd(w,h)
命题得证。

#include<iostream>
#include<cstdio>
#include<string.h>
using namespace std;

typedef long long LL;
LL gcd(LL a,LL b)
{
  return b==0?a:gcd(b,a%b);
}

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