POJ 3034 DP

#include "stdio.h"
#include "string.h"
#include "stdlib.h"

struct node
{
    int x,y;
}map[11][1010];

int temp[31][31],dp[2][31][31],time[11];

int Max(int a,int b)
{
    if(a<b) return b;else return a;
}
int gcd(int a,int b)
{
    if (b==0) return a; else return gcd(b,a%b);
}
int cal(int x,int y,int tx,int ty) // 记录在某时间从(x,y)走到(tx,ty） 所产生的价值
{
    int dx,dy,tp,ans;
    if (x==tx && y==ty) return temp[x][y];
    dx=tx-x;
    dy=ty-y;
    tp=gcd(abs(dx),abs(dy));
    dx/=tp;
    dy/=tp;

    ans=0;
    for (x,y;x!=tx || y!=ty;x+=dx,y+=dy)
        ans+=temp[x][y];
    ans+=temp[tx][ty];
    return ans;
}
int main()
{
    int n,d,m,x,y,t,i,j,l,ans,a,b;
    while (scanf("%d%d%d",&n,&d,&m)!=EOF)
    {
        if (n+d+m==0) break;
        memset(map,0,sizeof(map));
        memset(time,0,sizeof(time));
        while (m--)
        {
            scanf("%d%d%d",&x,&y,&t);
            time[t]++; //记录时间ｔ一共出现了多少地鼠
            map[t][time[t]].x=x+d; // 记录时间t出现的每一个地鼠
            map[t][time[t]].y=y+d;
        }
        n+=d+d;

        memset(dp,0,sizeof(dp));

        for (i=1;i<=10;i++)
        {
            memset(temp,0,sizeof(temp)); 
            for (j=1;j<=time[i];j++)
                temp[map[i][j].x][map[i][j].y]=1;// 记录时间i出现的地鼠图
            for (j=0;j<n;j++) 
                for (l=0;l<n;l++) // 枚举位置
                {
                    dp[i%2][j][l]=dp[1-i%2][j][l];
                    for (a=-d;a<=d;a++)
                        for (b=-d;b<=d;b++) //枚举移动方式
                        {
                            if (a*a+b*b>d*d) continue;
                            x=a+j;
                            y=b+l;
                            if (x<0 || x>=n || y<0 || y>=n) continue;
                            dp[i%2][j][l]=Max(dp[i%2][j][l],dp[1-i%2][x][y]+cal(x,y,j,l));
                        }
                }
        }
        ans=0;
        for (i=0;i<n;i++)
            for (j=0;j<n;j++)
            ans=Max(ans,dp[0][i][j]);

        printf("%d\n",ans);
    }
    return 0;
}