一道非常好的DP，菜狗子多做这种题

题目描述

今年是国际数学联盟确定的“2000——世界数学年”，又恰逢我国著名数学家华罗庚先生诞辰90周年。在华罗庚先生的家乡江苏金坛，组织了一场别开生面的数学智力竞赛的活动，你的一个好朋友XZ也有幸得以参加。活动中，主持人给所有参加活动的选手出了这样一道题目:

设有一个长度N的数字串，要求选手使用K个乘号将它分成K+1个部分，找出一种分法，使得这K+1个部分的乘积能够为最大。

同时，为了帮助选手能够正确理解题意，主持人还举了如下的一个例子：

有一个数字串: 312，当N=3，K=1时会有以下两种分法：

1）3*12=36

2）31*2=62

这时，符合题目要求的结果是： 31*2=62

现在，请你帮助你的好朋友XZ设计一个程序，求得正确的答案。

输入格式

程序的输入共有两行:

第一行共有2个自然数N,K (6<=N<=40，1<=K<=30)

第二行是一个长度为N的数字串。

输出格式

相对于输入，应输出所求得的最大乘积（一个自然数）。

样例输入：

4 2

1231

样例输出：

6 

2

这道题和我之前做过的一道hihocoder非常类似，我在做这道题的时候想到的是搜索，说实话，我压根就没有想到用DP，这实在是一个悲伤的故事。持续更新我还是一个菜鸟，我的DP还是一如既往的菜。

这道题为什么可以用DP呢？因为这明显具有子问题，并且具备记忆化搜索的特征，至于DP的其他特征我是真没看出来。好了，那么这道题的关键就出来了，我们该如何DP？假设我们设dp[i][j]表示在1-i中插入j个乘号的最大值，那么dp[n][k]就是这道题的答案。看起来非常简单，但问题出来了，状态转移方程呢？dp[i][j]不再是普通的和dp[i-1][j]或者dp[i-1][j-1]有联系，他还和dp[i-k][j-1]也有联系(这里的k的范围要注意，我会在代码中给出解释)。

[cpp]  view plain  copy

 

1. #include<bits/stdc++.h>  
2. #define ll long long  
3. using namespace std;  
4. ll num[50][50];  
5. ll dp[50][50];  
6. string s;  
7. int N,K;  
8. int main()  
9. {  
10.     cin>>N>>K;  
11.     cin>>s;  
12.     int i,j,m;  
13.     memset(num,0,sizeof(num));  
14.       
15.     for(i=0;i<N;i++)  
16.     {  
17.         num[i][i]=s[i]-'0';               //预处理  
18.         for(j=i+1;j<N;j++)  
19.         {   
20.             num[i][j]=num[i][j-1]*10+(s[j]-'0');  //预处理   
21.         }  
22.     }  
23.       
24.     for(i=0;i<N;i++)  
25.     dp[i][0]=num[0][i];     //这里对于dp数组的初始化非常重要   
26.       
27.       
28.     for(i=0;i<N;i++)  
29.     {  
30.         for(j=0;j<=K && j<i+1 ;j++)   //这个j<i+1 的意思是，在前i+1个数中，你的乘号的数量最多是i个，这个不难理解   
31.         {  
32.             for(m=j;m<=i;m++)       //这一句才是这道题最最重要的地方！！！m的范围，这个你自己理解去吧，要学会自己思考！！！   
33.             {  
34.                 dp[i][j]=max(dp[i][j],dp[m-1][j-1]*num[m][i]);  
35.             }  
36.         }  
37.     }  
38.       
39.     printf("%lld\n",dp[N-1][K]);  
40.     return 0;  
41. }  

下面是网上一个人写的三维的DP，看起来比我这个要高大上，不过他的代码过不了最后一个测试点，我给他在输出的时候改了一下判断条件就过了。

#include <iostream>  
#include <string>  
#include <vector>  
using namespace std;  
  
int n,k;  
char a[50];    
long long num[50][50];  
long long dp[50][50][10];  
  
int main()  
{  
    cin >> n >> k;  
    cin >> a;  
    for(int i=0; i<n; i++)  
    {  
        num[i][i] = a[i]-'0';  
        dp[i][i][0] = num[i][i];  
        for(int j=i+1; j<n; j++)  
        {  
            num[i][j] = num[i][j-1]*10+(a[j]-'0');  
            dp[i][j][0] = num[i][j];  
        }  
    }  
    for(int i=0; i<n; i++)  
    {  
        for(int j=i+1; j<n; j++)  
        {  
            for(int kk=1; kk<=k; kk++)  
            {  
                long long maxV = 0;  
                for(int m=i; m<=j; m++)  
                {  
                    for(int left=0; left<=kk; left++)  
                    {  
                        dp[i][j][kk] =max(dp[i][j][kk],dp[i][m][left]*dp[m+1][j][kk-left-1]);
                    }  
                }  
                 
            }  
        }  
    }  
    if(dp[0][n-1][k]>99999999999999999999)
    cout<<0;
    else
    cout << dp[0][n-1][k];  
    return 0;  
}