今年是国际数学联盟确定的“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的范围要注意,我会在代码中给出解释)。
下面是网上一个人写的三维的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;
}