UVALive 3720 Highways(组合计数、递推)

题意：n行m列的点阵，问有多少非垂直并且非水平的直线至少穿过点阵中的两个点。

思路：第一眼看的时候感觉和一道平面坐标中能看到点的题挺像的，于是按着那题的思路推了个用欧拉公式算的。写完以后发现不对，欧拉公式推的并不能包含所有情况，而且会有重复统计的问题。换个思路，想到了用几个矩阵组合出答案的方法，首先只考虑一个方向的直线(如斜率小于0的)，则a[i][j]=a[i-1][j]+a[i][j-1]-a[i-1][j-1]+c[i][j]，a[i][j]表示(i,j)表示的矩阵的方案数，c[i][j]表示经过(i,j)点的直线和前面的点的不重复的组合数。可以发现gcd(i,j)==1则有这样一条直线，但这样统计同样也会有重复，计算的时候要把这部分减去。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=300+10;
int a[maxn][maxn],sum[maxn][maxn];
int gcd(int a,int b) { return b==0?a:gcd(b,a%b); }
void Init()
{
    memset(a,0,sizeof(a));
    memset(sum,0,sizeof(sum));
    for(int i=1;i<maxn;++i)
        for(int j=1;j<maxn;++j)
        {
            int tmp=gcd(i,j);
            if(tmp==1)
                a[i][j]=1;
            else if(tmp==2)
                a[i][j]=-1;
        }
    for(int i=1;i<maxn;++i)
        for(int j=1;j<maxn;++j)
            a[i][j]+=a[i][j-1];
    for(int i=1;i<maxn;++i)
        for(int j=1;j<maxn;++j)
            a[i][j]+=a[i-1][j];
    for(int i=1;i<maxn;++i)
        for(int j=1;j<maxn;++j)
            sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+a[i][j];
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    Init();
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0) break;
        n--;m--;
        printf("%d\n",sum[n][m]*2);
    }
    return 0;
}