ZOJ 3647 Gao the Grid

直接是算出所有的三角形，再减去共线的三角形。GCD求线上有多少个第三点，就是当前这种线有多少个不符合的三角形

#include <iostream>
#include <cstdio>
#include <cmath>
#include <vector>
#include <cstring>
#include <algorithm>
#include <string>
#include <set>
#include <functional>
#include <numeric>
#include <sstream>
#include <stack>
#include <map>
#include <queue>
using namespace std;
long long c(long long n,long long m)
{
	long long sum = 1;
	for(long long i=1;i<=m;i++)
	{
		sum = sum*(n-m+i)/i; 
	} 
	return sum; 
} 
long long gcd(long long a,long long b)
{
	return b==0? a:gcd(b,a%b); 
} 
int main()
{
	long long n,m;
	while(cin >> n >> m)
	{
		long long ans = c((n+1)*(m+1),3);
	//	cout<<ans<<endl; 
		if(n>=2)ans-=(m+1)*c(n+1,3);
		if(m>=2) ans-=(n+1)*c(m+1,3);
	//	cout<<ans<<endl; 
		for(int i=2;i<=n;i++)
		{
			for(int j=2;j<=m;j++)
			{
				long long cnt = gcd(i,j)-1;
				//if(cnt>=3) 
				ans -= (n-i+1)*(m-j+1)*2*cnt; 
			} 
		} 
		cout<<ans<<endl; 
	} 
	return 0; 
}