hdu2571 简单dp的两种解决思路

命运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 Input13 89 10 10 10 10 -10 10 1010 -11 -1 0 2 11 10 -20-11 -11 10 11 2 10 -10 -10Sample Output

52

这是第一种从左上角往右下角打表的方法。

#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>

using namespace std;
#define maxx 100000;

const int ma1=50;
const int ma2=2500;
int data[ma1][ma2];
int dp[ma1][ma2];

void init()
{
    memset(data,0,sizeof(data));
   memset(dp,0,sizeof(dp));
}

int main()
{
   int K;
   int m,n;
   scanf("%d",&K);
   while(K--)
   {
       scanf("%d%d",&n,&m);
       init();
       for(int i=1;i<=n;i++)
          for(int j=1;j<=m;j++)
             {
                 dp[i][j]=-maxx;//这一步初始化很关键
                 scanf("%d",&data[i][j]);
             }
        dp[1][1]=data[1][1];
        for(int i=2;i<=n;i++)
              dp[i][1]=max(dp[i][1],dp[i-1][1]+data[i][1]);//处理第一列，只从上面来
        for(int j=2;j<=m;j++)//处理第一行,看清处理方法，特别是整除部分,
        {
            dp[1][j]=max(dp[1][j],dp[1][j-1]+data[1][j]);//从左边来

          for(int k=1;k<=j/2;k++)//下面计算部分与这条语句思想相同
                if(j%k==0)
                   dp[1][j]=max(dp[1][j],dp[1][k]+data[1][j]);
        }
        //计算开始
        for(int i=2;i<=n;i++)
            for(int j=2;j<=m;j++)
              {
                  dp[i][j]=max(dp[i][j],dp[i-1][j]+data[i][j]);
                  dp[i][j]=max(dp[i][j],dp[i][j-1]+data[i][j]);
                  for(int k=1;j<=j/2;k++)
                       if(j%k==0)
                          dp[i][j]=max(dp[i][j],dp[i][k]+data[i][j]);
              }
        cout<<dp[n][m]<<endl;
              }
    return 0;
}

这是第二种从右下向左上角打表的方法

#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>

using namespace std;
#define maxx 100000;

const int ma1=50;
const int ma2=2500;
int data[ma1][ma2];
int dp[ma1][ma2];

void init()
{
    memset(data,0,sizeof(data));
   memset(dp,0,sizeof(dp));
}

int main()
{
   int K;
   int m,n;
   scanf("%d",&K);
   while(K--)
   {
       scanf("%d%d",&n,&m);
       init();
       for(int i=1;i<=n;i++)
          for(int j=1;j<=m;j++)
                 scanf("%d",&data[i][j]);
        dp[n][m]=data[n][m];
        for(int i=n-1;i>=1;i--)
              dp[i][m]=dp[i+1][m]+data[i][m];//处理最后一列
        for(int j=m-1;j>=1;j--)//处理最后一行
        {
            dp[n][j]=dp[n][j+1]+data[n][j];
            for(int k=j*2;k<=m;k++)
                  if(k%j==0)
                    dp[n][j]=max(dp[n][j],dp[n][k]+data[n][j]);
        }
        //计算开始
        for(int i=n-1;i>=1;i--)
            for(int j=m-1;j>=1;j--)
              {
                  dp[i][j]=max(dp[i+1][j]+data[i][j],dp[i][j+1]+data[i][j]);
                  for(int k=j*2;k<=m;k++)
                       if(k%j==0)
                          dp[i][j]=max(dp[i][j],dp[i][k]+data[i][j]);
              }
        cout<<dp[1][1]<<endl;
              }
    return 0;
}

与该题相同每到DP的题都会有两种解法，一种是从前往后推的递推，还有一种是从后往前推的递归，两种算法的时间复杂度是相同的，对部分题两种方法可能会有代码长短的区别，但耗时上没太大区别，因此大多数情况下我们两种方法都是可行的。DP的核心思想是防止重复计算，因此记录下每次的计算结果是关键，这样可以大大缩短运算时间，因此这是一种高效的算法，同样，更深层次的DP不光会缩短时间复杂度，同样也会考虑空间复杂度的缩短，因为每次计算中都会出现很多多余的数据，每次计算都可以覆盖多余的数据以此来缩短空间复杂度，简单的零一背包问题就是这样，读者可以自己理解。

对于本题，当然也可以缩短空间复杂度，将从左上到右下角的代码改为：

#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

int data[22][1500];


void init(int n,int m)
{
    for(int i=0;i<=n;i++)
      for(int j=0;j<=m;j++)
         data[i][j]=-500000;
    data[0][1]=data[1][0]=0;

}

int main()
{
 int K;
 int m,n;
 scanf("%d",&K);
 while(K--)
  {
      scanf("%d%d",&n,&m);
      init(n,m);
      for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&data[i][j]);
              int temp=data[i][j];
             data[i][j]+=max(data[i-1][j],data[i][j-1]);
             for(int k=1;k<=j/2;k++)
                if(j%k==0)
                   data[i][j]=max(data[i][j],data[i][k]+temp);
        }
        cout<<data[n][m]<<endl;

  }
}

读者可自行测试一下，与前面的代码长度和memery上有明显的缩短，因为去掉了打表用的dp数组，空间上肯定节省。