直接是算出所有的三角形,再减去共线的三角形。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;
}