题意:给定多组n,m,求一个边长为n*m的网格表中的网格顶点能组成多少个三角形。
分析:网格点组成三角形,直接点求就是C((n+1)*(m+1),3)的总的取三个顶点的情况数减去三点共线的情况,我们能直接减去在平行于x轴和y轴的直线上的情况(n+1)*C(m+1,3)+(m+1)*C(n+1,3),还有就是斜线上的三点共线,我们只算斜率为正的直线上的情况,最后*2即可。我们设dp[i][j]表示在边长为i*j的矩形中斜线上三点共线的方案数,那么有dp[i][j]=do[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+gcd(i,j)-1。预处理好这些量就能直接处理啦。
代码:
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=1010;
const int MAX=151;
const int MOD1=1000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const ll INF=10000000010;
typedef unsigned long long ull;
int gcd(int a,int b) {
return b ? gcd(b,a%b):a;
}
ll dp[N][N];
void deal() {
int i,j;
for (i=2;i<=1000;i++)
for (j=2;j<=1000;j++)
dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+gcd(i,j)-1;
for (i=2;i<=1000;i++)
for (j=2;j<=1000;j++)
dp[i][j]+=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];
}
ll C(int n,int m) {
ll ret=1;
for (int i=0;i<m;i++) ret*=n-i;
for (int i=1;i<=m;i++) ret/=i;
return ret;
}
int main()
{
int n,m,ca=0;
deal();
while (scanf("%d%d", &n, &m)&&(n+m)) {
n++;m++;
printf("Case %d: %lld\n", ++ca, C(n*m,3)-n*C(m,3)-m*C(n,3)-2*dp[n-1][m-1]);
}
return 0;
}