3505: [Cqoi2014]数三角形

Description

给定一个nxm的网格，请计算三点都在格点上的三角形共有多少个。下图为4x4的网格上的一个三角形。

注意三角形的三点不能共线。

Input

输入一行，包含两个空格分隔的正整数m和n。

Output

输出一个正整数，为所求三角形数量。

Sample Input

2 2

Sample Output

76

数据范围

1<=m,n<=1000

思路：

数三角形，我们首先看看， n 行 m 列， 所以一共有 (n+1)(m+1) 个点。设为 t
根据排列组合，所以我们可以找到
t(t-1)(t-2)/6 种方案。
然后我们在找三点共线的方案数。
假设两个点， (a,b), (c,d) 这两个点分别是三个点的两个端点，那么中间一共有 gcd(c-a,d-b)-1 个点满足条件。
很多线段都是重复的，所以我们可以从(0,0) 进行平移就可以了。
假设两个点是 （0,0）（a，b） 那么他可以平移的得到的线段就是，，
（n-a+1）（m-j+1）
如果（a，b）不在坐标轴上，那么所得到的方案数要 乘上 2；、

#include<bits/stdc++.h>
using namespace std; 
int g[2000][2000];
int n,m;
int getgcd(int a, int b){
	if (g[a][b]) return g[a][b];
	if (!a) return g[a][b] = b;
	if (!b) return g[a][b] = a;
	return g[a][b] = getgcd(b, a%b);
}
void cal(){
	for (int i = 0; i <= n; i++)
		g[i][0] = i;
	for (int i = 0; i <= m; i++)
		g[0][i] = i;
	for (int i = 1; i <= n; ++i)
		for (int j = 1; j <= m; ++j)
			g[i][j] = getgcd(i,j);
	return;
}
int main(){
	scanf("%d%d",&n,&m);
	cal();
	long long ans,t;
	t = 1ll*(n+1)*(m+1);
	ans = t*(t-1)*(t-2)/6;
	for (int i = 0; i <= n; ++i)
	for (int j = 0; j <= m; ++j)
	if (i || j){
		if (!i || !j) ans -= (1ll*(g[i][j]-1)*(n-i+1)*(m-j+1)); else
		ans -= (1ll*2*(g[i][j]-1)*(n-i+1)*(m-j+1));
	}	
	printf("%lld\n",ans);
	return 0;
}