数三角形 (Counting Triangles,Dhaka 2005,LA 3295)

#include <iostream>
#include <string.h>
#include <string>
#include <fstream>
#include <algorithm>
#include <stdio.h>
#include <vector>
#include <queue>
#include <set>
#include <cmath>
using namespace std;
const double eps = 1e-8;
const int INF=0x7fffffff;
unsigned long long uINF = ~0LL;
#define MAXN 10000007
typedef long long LL;
LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}
LL g[1001][1001];
void init()
{
    for(int i=1;i<=1000;i++)
        for(int j=1;j<=1000;j++)
        g[i][j]=gcd(i,j);
}
int main()
{
    LL n,m;int t=1;
    init();
    while(scanf("%lld%lld",&n,&m),n+m)
    {
        LL sum=n*m+n+m+1;
        sum=sum*(sum-1)*(sum-2)/6;
        LL x=(n+1)*n*(n-1)/6*(m+1);
        LL y=(m+1)*m*(m-1)/6*(n+1);
        LL ans=0;
        for(int i=2;i<=n;i++)
            for(int j=2;j<=m;j++)
            ans+=(g[i][j]-1)*(n-i+1)*(m-j+1);

            ans*=2;
        printf("Case %d: %lld
",t++,sum-x-y-ans);
    }

    return 0;
}

网格点数 所有三元组 然后排除共线的三元组

和 LA3720 一样 枚举斜率 ~