hdu2571_命运（DP，记忆化）

题目大意：

一个人在一个矩阵里头，要从起点[1，1]走到终点[n,m].而且规定只能向右或者向下走，向下一次只能走一格。但是如果向右走，则每次可以走一格或者走到该行的列数是当前所在列数倍数的格子，即：如果当前格子是（x,y），下一步可以是（x+1,y），(x,y+1)或者(x,y*k) 其中k>1。每一格都有一个整数，代表他能够获得的幸运值（有正有负）。问一种走法，使得他获得最大的幸运值。

题目分折：

如果用dp[i][j]表示他能够走到[i,j]，这个位置所能获得的最大幸运值，那么dp[i][j] = max(dp[i-1][j],dp[i][j-1],dp[i][j/k])+a[i][j].其中（2<k<=j && j%k==0).为了代码的方便，我们初始化所有的 dp[i][j] = -INF,只对dp[1][0]=dp[0][1] = 0;则可不考虑i==1||j==1的一系列麻烦的情况，这相当于我们在最外头加了一层，我们的值是从1开始。

如果我们想要用记忆化的思想，那么递归方程可以写为 ans = max(ans,DP(x+1,y),DP(x,y+1),DP(x,k*y))+a[x][y].结束条件是（x==n && y==m)return a[x][y];

题目链接： http://acm.hdu.edu.cn/showproblem.php?pid=2571

小乐一下：

我们每个人都有自己的命运，每个人都希望得到上天幸运之手的帮助，但生命之中，掌握自己命运的一定是自己。

我们对于可以用记忆化来写的程序，尽量用记忆化来定，因为它的代码比较简单，理解也比较容易。当然用DP来写是很能锻炼自己的思维能力的。

我的代码，更好更简洁的由你们来实现

DP：

#include<cstdio>
#include<cstring>
const int INF = 0x3fffffff;
int dp[25][1010];
int a[25][1010];
int n,m;
int max(int a,int b){return a>b?a:b;}
int main(){
    int T;
    int i,j,k;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        memset(a,0,sizeof(a));
        //memset(dp,0,sizeof(dp));
        for(i = 0;i<=n;i++){
            for(j = 0;j<=m;j++) dp[i][j] = -INF;
        }
        dp[1][0] = dp[0][1] = 0;   //这样做的目的是为了防止最左边的，还能够由它左边的相加来达到最大值。当然第一个除外。
        for(i = 1;i<=n;i++){
            for(j = 1;j<=m;j++){
                scanf("%d",&a[i][j]);
            }
        }
        for(i = 1;i<=n;i++){
            for(j = 1;j<=m;j++){
                dp[i][j] = max(dp[i-1][j],dp[i][j-1])+a[i][j]; //首先比较由左边，上边得到的较大值。
                for(k = 2;k<=j;k++){
                    if(j%k==0) dp[i][j] = max(dp[i][j],dp[i][j/k]+a[i][j]); // 看由前面跳跃到这一步是否可以得到更大的值。
                }
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}

记忆化：

#include<stdio.h>
#include<string.h>
int dp[25][1010],n,m;
int vis[25][1010],a[25][1010];
const int INF = 0x3fffffff;
int max(int a,int b){return a>b?a:b;}
int DP(int x,int y){
    if(x==n && y==m) return dp[x][y] = a[x][y];  //递归结束出口。
    if(vis[x][y]) return dp[x][y];  //这个结点已经知道从终点到这个点的最大值，可以直接返回。
    int ans = -INF;
    if(x<n) ans = max(ans,a[x][y]+DP(x+1,y));  //往右走一步。
    if(y<m) ans = max(ans,a[x][y]+DP(x,y+1)); //往左走一步。
    int k = 2;
    while(k*y<=m){
        ans = max(ans,a[x][y] + DP(x,k*y));  //跳倍数走。最终的结果保存在ans 里面。
        k++;
    }
    vis[x][y] = 1;  //这个点我们会知道结果了的。
    return dp[x][y] = ans;
}

int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        memset(dp,0,sizeof(dp));
        memset(vis,0,sizeof(vis));
        for(int i = 1;i<=n;i++){
            for(int j = 1;j<=m;j++) scanf("%d",&a[i][j]);
        }
        printf("%d\n",DP(1,1));
    }
    return 0;
}

另外一道非常类似的题。hdu2391_Filthy Rich   http://acm.hdu.edu.cn/showproblem.php?pid=2391

我也给出DP以及记忆化搜索的代码，具体的话就不多解释了，自己体会。

记忆化：

#include<cstdio>
#include<cstring>
int vis[1010][1010],dp[1010][1010],a[1010][1010];
int n,m;
int max(int a,int b){return a>b?a:b;}
int DP(int x,int y){
    if(x==n && y==m) {vis[x][y] = 1;return a[x][y];}
    if(vis[x][y]) return dp[x][y];
    int ans = 0;
    if(x+1<=n) ans = max(ans,DP(x+1,y)+a[x][y]);
    if(y+1<=m) ans = max(ans,DP(x,y+1)+a[x][y]);
    if(x+1<=n && y+1<=m) ans = max(ans,DP(x+1,y+1)+a[x][y]);
    vis[x][y] = 1;
    return dp[x][y] = ans;
}

int main(){
    int T;
    int i,j;
    int cnt = 1;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        memset(dp,0,sizeof(dp));
        memset(vis,0,sizeof(vis));
        for(i = 1;i<=n;i++){
            for(j = 1;j<=m;j++){
                scanf("%d",&a[i][j]);
            }
        }
        printf("Scenario #%d:\n%d\n\n",cnt++,DP(1,1));
    }
    return 0;
}

DP：

#include<cstdio>
#include<cstring>
int dp[1010][1010];
int a[1010][1010];
int n,m;
int max(int a,int b){return a>b?a:b;}
int main(){
    int T;
    int i,j;
    int cnt = 1;
    scanf("%d",&T);
    while(T--){
        memset(a,0,sizeof(a));
        memset(dp,0,sizeof(dp));
        scanf("%d%d",&n,&m);
        for(i = 1;i<=n;i++){
            for(j = 1;j<=m;j++) scanf("%d",&a[i][j]);
        }
        for(i = 1;i<=n;i++){
            for(j = 1;j<=m;j++){
                dp[i][j] = max(max(dp[i-1][j],dp[i][j-1]),dp[i-1][j-1])+a[i][j];
            }
        }
        printf("Scenario #%d:\n",cnt++);
        printf("%d\n",dp[n][m]);
        printf("\n");
    }
    return 0;
}

伟大的梦想成就伟大的人，从细节做好，从点点滴滴做好，从认真做好。