【SSL】1007乘积最大

【SSL】1007乘积最大

Time Limit:1000MS

Memory Limit:65536K

Description

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

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

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

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

1) 312=36
　　2) 312=62

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

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

Input

程序的输入共有两行：
　　第一行共有2个自然数N，K（6≤N≤40，1≤K≤6）
　　第二行是一个长度为N的数字串。

Output

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

Sample Input

4 2
　　1231

Sample Output

62

思路

与加法最大类似。
引入记号：f(i,s)表示从1到i，s个乘号取得的最大值，g(i,j)表示从a[i]到a[j]的数字列，则有下列推导公式：
k=1时，f(n,1)=max{g(1,1)*g(2,n),g(1,2)*g(3,n),…,g(1,n-1)*g(n,n)}
k=2时，f(n,2)=max{f(n-1,1)*g(n,n),f(n-2,1)*g(n-1,n),…,f(2,1)*g(3,n)}
一般地,f(n,k)=max{f(n-1,k-1)*g(n,n),f(n-2,k-1)*g(n-1,n),…,f(k,k-1)*g(k+1,n)}
　　由此可见，为了计算f(n,k)必须算出f(n-1,k-1),f(n-2,k-1),…,f(k,k-1)。
f(n,k)=max{f(n-1,k-1)*g(n,n),f(n-2,k-1)*g(n-1,n),…,f(k,k-1)*g(k+1,n)}。
f[i][j]=max(f[l][j-1]*g[l+1][i]);1<=i<=n,1<=j<=k;j<=l<i;
预处理g(i,j)数字列。

代码

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<iomanip>
using namespace std;
char str[1010];
long long n,k,gs,f[1010][1010],cmpsz,s[1010],a[1010][1010]; 
void input()
{
	long long i,j;
	scanf("%lld%lld",&n,&k); 
	scanf("%s",&str);
	for(i=1;i<=n;i++)
	{
		s[i]=str[i-1]-'0';
		a[i][i]=s[i];
	}
	for(i=1;i<=n;i++)
		for(j=i;j<=n;j++)
			a[i][j]=a[i][j-1]*10+s[j];//预处理
	for(i=1;i<=n;i++)
		f[i][0]=a[1][i];
	return;
}
void DP()
{
	long long l,i,j;
	for(i=1;i<=n;i++)
		for(j=1;j<=k;j++)
			for(l=j;l<i;l++)
				f[i][j]=max(f[i][j],f[l][j-1]*a[l+1][i]);//状态转移方程
	return;
} 
int main()
{
	input();
	DP();
	printf("%lld",f[n][k]);
	return 0;
}