poj 3034 DP

这道题一开始思路是对的，但是就是感觉会超时。后来看了一下别人的思路，处理点(i,j)到点（x，y）上有几个老鼠时用到了gcd();感觉很巧妙，学到了一个技巧。

设 dp[t][i][j]表示时间为t是，锤子在(i, yj坐标的最大值。枚举x,y，求t + 1时间到(i, j)的最优值， dis(i, j, x, y) <= d。在计算(x, y) -> (i, j)的时候不能一个单位一个单位的

加，否则会超时。dx表示x -> i的偏移量，dy表示y -> j的偏移量。然后求出dx, dy的最大公约数dd。 dx /= dd, dy /= dd。这样从(x, y) 到(i, j)时没次x += dx, y += dy;

另外最重要的一点是，这道题有神坑。如果想不到的话，就等着WA到死吧。

中转点可能在区域外，如下数据：　　　　//相当重要！！！

20 5 4
1 0 1
0 1 1
0 5 2
1 6 2
答案应该是4而不是3，这就需要扩展区域为n+2*d就够啦。

就是说，你可能把锤子放到map外面，然后沿着某个斜线打到老鼠更多（因为在map外，某个点可能对应着两条有很多老鼠的斜线，这样就可以打到更多）。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX 40
int n,d,m;
int map[MAX][MAX][25];
int dp[25][MAX][MAX];
int x,y,t,xx,yy,dx,dy,MAX_t;
int ans;
bool dis(int x1,int y1,int x2,int y2)
{
    int tempx=(x1-x2);
    int tempy=(y1-y2);
    return (tempx*tempx+tempy*tempy)<=d*d;
}
int gcd(int a,int b)
{
    if(a==0) return b;
    return gcd(b%a,a);
}
void work(int N)
{
    int i,j,k;
    ans=0;
    for(k=1;k<=MAX_t;k++){
        for(i=0;i<N;i++)
        {
            for(j=0;j<N;j++)
            {
                for(x=0;x<N;x++)
                {
                    for(y=0;y<N;y++)
                    {
                        if(dis(i,j,x,y))
                        {
                            if(i==x&&j==y)
                            {
                                dp[k][i][j]=max(dp[k][i][j],dp[k-1][i][j]+map[i][j][k]);
                            }
                            else
                            {
                                dx=x-i;dy=y-j;
                                int dd=gcd(abs(dx),abs(dy));
                                dx/=dd;dy/=dd;
                                xx=i;yy=j;
                                int num=0;
                                while(1)
                                {
                                    num+=map[xx][yy][k];
                                    if(xx==x&&yy==y) break;
                                    xx+=dx;
                                    yy+=dy;
                                }
                                dp[k][i][j]=max(dp[k][i][j],dp[k-1][x][y]+num);
                            }
                        }
                    }
                }
            }
        }
    }
    for(i=d;i<N;i++)
    {
        for(j=d;j<N;j++)
        {
            ans=max(ans,dp[MAX_t][i][j]);
        }
    }
}
int main()
{
    int i,j;
    while(scanf("%d%d%d",&n,&d,&m),n)
    {
        memset(map,0,sizeof(map));
        memset(dp,0,sizeof(dp));
        MAX_t=0;
        for(i=0;i<m;i++)
        {
            scanf("%d%d%d",&x,&y,&t);
            MAX_t=max(MAX_t,t);
            map[x+d][y+d][t]=1;
        }
        work(n+2*d);
        printf("%d\n",ans);
    }
    return 0;
}