BZOJ 3505 [Cqoi2014]数三角形（组合数学）

【题目链接】 http://www.lydsy.com/JudgeOnline/problem.php?id=3505

【题目大意】

　　给定一个nxm的网格，请计算三点都在格点上的三角形共有多少个。
　　注意三角形的三点不能共线。

【题解】

　　我们计算三个点组合的情况，去除横竖三共线，以及斜着三点共线的情况即可。
　　一个矩形斜对角上的整点数为其长宽的最大公约数+1.

【代码】

#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long LL;
LL c[1001000][4],ans;
int n,m;
void init(){
    c[0][0]=1;
    for(int i=1;i<=n*m;i++){
        c[i][0]=1;
        for(int j=1;j<=3;j++)c[i][j]=c[i-1][j-1]+c[i-1][j];
    }
}
void solve(){
    ans=c[n*m][3]-n*c[m][3]-m*c[n][3];
    for(int i=1;i<n;i++)for(int j=1;j<m;j++){
        LL t=__gcd(i,j)+1;
        if(t>2)ans-=(t-2)*2*(n-i)*(m-j);
    }
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        ans=0;n++;m++;
        init(); solve();
        printf("%lld\n",ans);
    }return 0;
}