乘积最大(DP + 高精度)

题目链接 https://www.luogu.com.cn/problem/P1018

今年是国际数学联盟确定的“ 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 设计一个程序，求得正确的答案。（6≤N≤40,1≤K≤6）

由题目和数据范围容易想到本题应用DP方法，

dp[i][k[表示前i个数，分成k段时的最大值

所以我们需要枚举最后一个乘号是在哪里放的

转移时，直接计算一下最后一个乘号之前的数字分成k-1段时的最大值，再乘以最后一个乘号之后的数字就可以了。

注意：本题数字过大，应该用高精度计算，且因为我们要用到原数字中间的很多部分，所以可以先预处理原数字的所有子序列大小，最后的乘法即高精成高精

AC代码(细节部分有注释)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>

using namespace std;

typedef long long ll;
const int N = 50;
char s[N];
int a[N][N][N], f[N][10][50000];

void modify(int x, int s, int m)
{//f[m - 1][s - 1][0]因为从m处分所以前面为m-1，且转移前为s-1段 
	int ans[50000];//！！！原来p>=j对应此处m-1>=s-1 
	memset(ans, 0, sizeof(ans));
	for(int i = 1; i <= a[m][x][0]; i++)
		for(int j = 1; j <= f[m - 1][s - 1][0]; j++)//！！！！m-1注意可能会让数组越界为-1！！！ 
			ans[i + j - 1] += a[m][x][i] * f[m - 1][s - 1][j];//类似于乘法列竖式计算 
	for(int i = 1; i <= a[m][x][0] + f[m - 1][s - 1][0]; i++)
	{
		if(ans[i])
			ans[0] = i;
		ans[i + 1] += ans[i] / 10;//进位处理 
		ans[i] %= 10;
	}
	if(ans[0] < f[x][s][0])
		return;
	for(int i = ans[0]; i >= 1; i--)
	{
		if(ans[i] > f[x][s][i])
			break;
		if(ans[i] < f[x][s][i])
			return;
	}//判断是否更新 
	for(int i = 0; i <= ans[0]; i++)
		f[x][s][i] = ans[i];
	return;
}

int main()
{
	int n, k;
	scanf("%d%d", &n, &k);
	scanf("%s", s + 1);
	for(int i = 1; i <= n; i++)
		for(int j = i; j <= n; j++)
		{
			a[i][j][0] = j - i + 1;
			for(int p = j - i + 1; p >= 1; p--)
				a[i][j][j - i + 2 - p] = s[i + p - 1] - '0';
		}//字符串转高精的方法 
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= i && j <= k + 1; j++)//前i位最多k+1段或者i段 
		{
			if(j == 1)//j=1为边界，即不加乘号，只有一段 
			{
				f[i][j][0] = i;
				for(int p = 1; p <= i; p++)
					f[i][j][p] = a[1][i][p];
			}
			for(int p = i; p >= j; p--)//因为f[i][j]为前i位数分为几段，所以i最小也得为j 
				modify(i, j, p);//！！函数变量代表前i位分成j段，从p处转移 
		}
	for(int i = f[n][k + 1][0]; i >= 1; i--)
		printf("%d", f[n][k + 1][i]);
return 0;
}