命运（动态规划）

Problem Description
穿过幽谷意味着离大魔王lemon已经无限接近了！
可谁能想到，yifenfei在斩杀了一些虾兵蟹将后，却再次面临命运大迷宫的考验，这是魔王lemon设下的又一个机关。要知道，不论何人，若在迷宫中被困1小时以上，则必死无疑！
可怜的yifenfei为了去救MM，义无返顾地跳进了迷宫。让我们一起帮帮执着的他吧！
命运大迷宫可以看成是一个两维的方格阵列，如下图所示：

yifenfei一开始在左上角，目的当然是到达右下角的大魔王所在地。迷宫的每一个格子都受到幸运女神眷恋或者痛苦魔王的诅咒，所以每个格子都对应一个值，走到那里便自动得到了对应的值。
现在规定yifenfei只能向右或者向下走，向下一次只能走一格。但是如果向右走，则每次可以走一格或者走到该行的列数是当前所在列数倍数的格子，即：如果当前格子是（x,y），下一步可以是（x+1,y），(x,y+1)或者(x,y*k) 其中k>1。
为了能够最大把握的消灭魔王lemon，yifenfei希望能够在这个命运大迷宫中得到最大的幸运值。

Input
输入数据首先是一个整数C，表示测试数据的组数。
每组测试数据的第一行是两个整数n,m，分别表示行数和列数(1<=n<=20,10<=m<=1000)；
接着是n行数据，每行包含m个整数，表示n行m列的格子对应的幸运值K ( |k|<100 )。

Output
请对应每组测试数据输出一个整数，表示yifenfei可以得到的最大幸运值。

Sample Input
1
3 8
9 10 10 10 10 -10 10 10
10 -11 -1 0 2 11 10 -20
-11 -11 10 11 2 10 -10 -10

Sample Output
52

思路：用动态规划得出到达每个位置获得的最大幸运值。要得到当前位置（i，j）的最大值，需要拿（i-1，j），（i，j-1），（i，j / k）（当j % k == 0时）的位置进行比较，取最大值，然后当前位置的最大值即为当前位置的值加上得到的最大值。其中 k 的取值范围为1 <=k<= j。
要注意的是二维数组 dp 的第0行第0列要初始化为小于100的数，因为第1行第1列的位置要从第0行或者第0列的值比较得到最大值，只有当第0行第0列的值足够小，在第一行或第一列的位置获得最大值时才不会干扰。而dp[1][1]要得到最大值，即要与dp[0][1]，dp[1][0]进行比较，再加上（1，1）位置的值，事实上dp[1][1] = mapp[1][1]，所以dp[0][1]，dp[1][0]要做处理等于0，这样dp[1][1] = dp[0][1]+mapp[1][1] = 0+mapp[1][1]，具体参见代码。

#include<stdio.h>
#include<string.h>
#define Max(x,y) x>y?x:y
int mapp[22][1002],dp[22][1002];

int main()
{
    int c,n,m,i,j,k;
    scanf("%d",&c);
    while(c--)
    {
        scanf("%d%d",&n,&m);
        for(i = 1; i<=n; i++)//记录数据是从第一行第一列开始的。
        {
            for(j = 1; j<=m; j++)
            {
                scanf("%d",&mapp[i][j]);
            }
        }
        /*memset(dp,0,sizeof(dp));//把dp[0][j]和dp[i][0]初始化为0，但是如果边界的数比0还小的时候那就会取0，所以出错
        比如mapp[1][2] = 5,dp[1][1] = -7,dp[0][2] = 0,dp[1][2]要得到最大值就会取dp[0][2]+mapp[1][2],而实际上dp[0][2]是不存在的，可以画一个二维数组的图理解
        把dp[i][0],dp[0][j]初始化为 <-100 的数。
        */
        for(i = 0; i<=n; i++)
        {
            dp[i][0] = -105;
        }
        for(j = 0; j<=m; j++)
        {
            dp[0][j] = -105;
        }
        dp[1][0] = 0,dp[0][1] = 0;//把起点处理一下，使起点dp[1][1]等于mapp[1][1]。
        for(i = 1; i<=n; i++)
        {
            for(j = 1; j<=m; j++)
            {
                dp[i][j] = Max(dp[i-1][j],dp[i][j-1]);
                for(k = 1; k<=j; k++)
                {
                    if(!(j%k))
                        dp[i][j] = Max(dp[i][j],dp[i][j/k]);//应该是取dp[i][k]，不是dp[i][j/k]？？？应该都可以
                        //k<=j时，取dp[i][j/k]; k<=j/2时，取dp[i][k] (取dp[i][j/k]，一些值取不到，例如 1 );
                }
                dp[i][j]+=mapp[i][j];
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}