【题解】NOIP2000 提高组 乘积最大

题目描述

今年是国际数学联盟确定的“ 2000 ――世界数学年”，又恰逢我国著名数学家华罗庚先生诞辰 90 周年。在华罗庚先生的家乡江苏金坛，组织了一场别开生面的数学智力竞赛的活动，你的一个好朋友 XZ 也有幸得以参加。活动中，主持人给所有参加活动的选手出了这样一道题目：
设有一个长度为 $N$ 的数字串，要求选手使用 $K$ 个乘号将它分成 $K+1$ 个部分，找出一种分法，使得这 $K+1$ 个部分的乘积能够为最大。
同时，为了帮助选手能够正确理解题意，主持人还举了如下的一个例子：
有一个数字串：$312$， 当 $N=3,K=1$ 时会有以下两种分法：

1. $3\times 12=36$

2. $31\times 2=62$

这时，符合题目要求的结果是: $31\times 2=62$
现在，请你帮助你的好朋友 XZ 设计一个程序，求得正确的答案。

输入格式

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

输出格式

结果显示在屏幕上，相对于输入，应输出所求得的最大乘积（一个自然数）。

输入输出样例

输入样例 #1

4  2
1231

输出样例 #1

62

Link: Click Me

题解

这题有一定的难度，值得来仔细讨论。

关注到这题的数据范围，已经超过了 long long，所以要使用高精度。笔者的代码没有使用高精度，但是其原理是一样的，供读者自行编写。(to do: 高精度)

这题是一道 区间DP 问题，基本上是 CSP 普及组难度的天花板了。我们把它按插入的乘号个数划分阶段，即
$$f[i][k]表示在前i位数中插入k个乘号所得的最大值$$
我们知道当该段获得最大值时，该段与该段后数字的乘积才能得到最大值（即满足 区间DP 的核心：从小区间推至大区间）。

我们还需要一个数组 a 用来存放第 $i$ 位到第 $j$ 位组成的自然数，例如：当输入字符串为 $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\le j\le 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;
}