洛谷 P1018 乘积最大（线性dp，高精度）

线性dp，高精度
本题要点：
1、 状态表示：
dp[i][j], 表示前i 个数，用 j个乘号得到的最大值
隐含条件，就是 数的个数 i - 1 >= j
2、 转移方程:
dp[i][j] = max{dp[k][j - 1], num[k + 1, i]}, j <= k < i
num[k + 1, i], 原来给出的字符串 s 从 k + 1 位 到 i位 所表示的数字。
3、因为数太大，超过 long long 范围，使用高精度。涉及的高精度运算有
大数 + 大数 bign add(bign a, bign b)
大数 * 大数 bign bign_multi(bign a, bign b)
大数 * 小数 bign multi(bign a, int b)

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int MaxN = 41;
char s[MaxN];
int n, K;

struct bign
{
	int d[1000];
	int len; 
	bign()
	:len(0)
	{
		memset(d, 0, sizeof d);
	}
};

bign dp[MaxN][7];	//  dp[i][j] 表示前i 个数，用 j个乘号得到的最大值 

bign str_to_bign(int L, int R)
{
	bign b;
	for(int i = R; i >= L; --i)
	{
		b.d[b.len++] = s[i] - '0';
	}
	// 去重前导 0， 不然后面比较大小会出错
	while(b.len > 1 && b.d[b.len - 1] == 0)
	{
		b.len--;
	}
	return b;
}

bign multi(bign a, int b)
{
	bign c;
	int carry = 0, tmp;
	for(int i = 0; i < a.len; ++i)
	{
		tmp = a.d[i] * b + carry;
		c.d[c.len++] = tmp % 10;
		carry = tmp / 10;
	}
	while(carry)
	{
		c.d[c.len++] = carry % 10;
		carry /= 10;
	}
	return c;
}

bign add(bign a, bign b)
{
	bign c;
	int carry = 0, tmp;
	for(int i = 0; i < a.len || i < b.len; ++i)
	{
		tmp = a.d[i] + b.d[i] + carry;
		c.d[c.len++] = tmp % 10;
		carry = tmp / 10;

	}
	while(carry)
	{
		c.d[c.len++] = carry % 10;
		carry = carry / 10;
	}
	return c;
}

bign bign_multi(bign a, bign b)
{
	bign c;
	for(int i = b.len - 1; i >= 0; --i)
	{
		c = multi(c, 10);
		c = add(c, multi(a, b.d[i]));
	}
	return c;
}

int cmp(bign a, bign b) // 0, -1 表示小于， 1 表示大于
{
	if(a.len != b.len)
	{
		if(a.len < b.len)	return -1;
		else return 1;
	}
	for(int i = a.len - 1; i >= 0; --i)
	{
		if(a.d[i] < b.d[i])	return -1;
		else if(a.d[i] > b.d[i])	return 1;
	}
	return 0;
}

int change(int L, int R)
{
	int ans  = 0;
	for(int i = L; i <= R; ++i)
		ans = ans * 10 + s[i] - '0';
	return ans;
}

void solve()
{
	bign tmp;
	for(int i = 1; i <= n; ++i)
	{
		dp[i][0] = str_to_bign(1, i);
	}
	for(int i = 2; i <= n; ++i)		// 前i个数
	{
		for(int j = 1; j < i && j <= K; ++j)	// j个乘号
		{
			for(int k = j; k < i; ++k)		// dp[i][j] = dp[k][j - 1] * change(k + 1, i)
			{
				tmp = bign_multi(dp[k][j - 1], str_to_bign(k + 1, i));
				if(cmp(dp[i][j], tmp) < 0)
					dp[i][j] = tmp;
			}
		}
	}
	for(int i = dp[n][K].len - 1; i >= 0; --i)
	{
		printf("%d", dp[n][K].d[i]);
	}
	printf("\n");
}

int main()
{
	scanf("%d%d", &n, &K);
	scanf("%s", s + 1);
	solve();
	return 0;
}

/*
4  2
1231
*/

/*
62
*/