线性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
*/