题目描述
今年是国际数学联盟确定的“ 2000 ――世界数学年”,又恰逢我国著名数学家华罗庚先生诞辰 90 周年。在华罗庚先生的家乡江苏金坛,组织了一场别开生面的数学智力竞赛的活动,你的一个好朋友 XZ 也有幸得以参加。活动中,主持人给所有参加活动的选手出了这样一道题目:
设有一个长度为
N
N
N 的数字串,要求选手使用
K
K
K 个乘号将它分成
K
+
1
K+1
K+1 个部分,找出一种分法,使得这
K
+
1
K+1
K+1 个部分的乘积能够为最大。
同时,为了帮助选手能够正确理解题意,主持人还举了如下的一个例子:
有一个数字串:
312
312
312, 当
N
=
3
,
K
=
1
N=3,K=1
N=3,K=1 时会有以下两种分法:
-
3 × 12 = 36 3 \times 12=36 3×12=36
-
31 × 2 = 62 31 \times 2=62 31×2=62
这时,符合题目要求的结果是:
31
×
2
=
62
31 \times 2 = 62
31×2=62
现在,请你帮助你的好朋友 XZ 设计一个程序,求得正确的答案。
输入格式
程序的输入共有两行:
第一行共有
2
2
2 个自然数
N
,
K
N,K
N,K(
6
≤
N
≤
40
,
1
≤
K
≤
6
6≤N≤40,1≤K≤6
6≤N≤40,1≤K≤6)
第二行是一个长度为
N
N
N 的数字串。
输出格式
结果显示在屏幕上,相对于输入,应输出所求得的最大乘积(一个自然数)。
输入输出样例
输入样例 #1
4 2
1231
输出样例 #1
62
Link: Click Me
题解
这题有一定的难度,值得来仔细讨论。
关注到这题的数据范围,已经超过了 long long
,所以要使用高精度。笔者的代码没有使用高精度,但是其原理是一样的,供读者自行编写。(to do: 高精度)
这题是一道 区间DP 问题,基本上是 CSP
普及组难度的天花板了。我们把它按插入的乘号个数划分阶段,即
f
[
i
]
[
k
]
表
示
在
前
i
位
数
中
插
入
k
个
乘
号
所
得
的
最
大
值
f[i][k] 表示在前 i 位数中插入 k 个乘号所得的最大值
f[i][k]表示在前i位数中插入k个乘号所得的最大值
我们知道当该段获得最大值时,该段与该段后数字的乘积才能得到最大值(即满足 区间DP 的核心:从小区间推至大区间)。
我们还需要一个数组 a
用来存放第
i
i
i 位到第
j
j
j 位组成的自然数,例如:当输入字符串为
114514
114514
114514 时,a[2][5] = 1451(a 数组下标从 1 开始)
可得状态转移方程:
f
[
i
]
[
k
]
=
m
a
x
{
f
[
j
]
[
k
−
1
]
×
a
[
j
+
1
]
[
i
]
}
f[i][k] = max\{f[j][k-1]\times a[j+1][i]\}
f[i][k]=max{f[j][k−1]×a[j+1][i]}
其中
k
≤
j
≤
i
k\le j\le i
k≤j≤i。
C++ 代码(不带高精度):
#include <iostream>
typedef long long LL;
const int N = 49;
int n, k;
LL f[N][7], a[N][N];
// f[i][k] 表示在前 i 位数中插入 j 个乘号所得的最大值
char s[N];
int main() {
scanf("%d%d%s", &n, &k, s+1);
for (int i = 1; i <= n; i++)
for (int j = i; j<=n && j-i<19; j++) {
if (i == j) a[i][j] = s[j] - '0';
else a[i][j] = a[i][j-1]*10 + s[j] - '0';
}
for (int i = 1; i <= n; i++)
f[i][0] = a[1][i];
for (int i = 1; i <= n; i++)
for (int j = 1; j <= k; j++)
for (int t = j; t < i; t++)
f[i][j] = std::max(f[i][j], f[t][j-1]*a[t+1][i]);
printf("%lld", f[n][k]);
return 0;
}