P1018 [NOIP2000 提高组] 乘积最大

题目

题目链接

题解

状态定义：dp[i][j]表示前i个数分成j段（即需要j+1个*）的最大乘积。
状态转移：dp[i][j] = max(dp[k-1][j-1]*a[k][i], dp[i][j])，表示在第k-1和第k个数之间加上一个*得到的最大值，其中前k-1个数采用的是dp值，即被分为了j-1段后的最大值，而后面的数（从k到i）则视为一整个数，与前面的dp值相乘。遍历全部k的位置，取最大值即为 dp[i][j]。

---

首先，暴力是不行的，暴力枚举每个*的位置，O(40^6)=O((2^12)*(10^6))=O(4,096,000,000)。

再者，由于要保存和计算的数最高为40位，因此还需要采用高精度。（之前写的高精度，有空重新整理一下）

---

整体思路就是先写出来非高精度版，再将其中的乘、加改成高精度的函数即可。

高精度函数：

1. max：计算两个高精度数的较大者，并返回string类型的高精度数。
2. add：计算一个高精度数和一个int类型数的和，因为用不到两个高精度数相加，所以就写了这一个。
3. mul：重载成两个函数了。第一个是计算两个高精度数的乘积，返回string类型的高精度数；第二个是计算一个高精度数与10相乘的结果。

---

另：如果在提交时出现RE，或者测试我给出的样例时出现了奇怪错误，要么是数组开小了，要么是高精度乘法最后位数的处理出现了问题。（经验之谈）

代码

未采用高精度的60%代码：

#include<bits/stdc++.h>
using namespace std;
const int N = 50, M = 10;

int n, m;
long long dp[N][M], a[N][N];
string num;

int main()
{
	cin >> n >> m;
	cin >> num;
	
	for (int i = 0;i < n;i ++)
		for (int j = i;j < n;j ++)
			a[i+1][j+1] = a[i+1][j]*10 + (num[j]-'0');
	
	for (int i = 1;i <= n;i ++)
		dp[i][1] = a[1][i];
	
	for (int i = 1;i <= n;i ++)
		for (int j = 2;j <= m+1 && j <= i;j ++) 
			for (int k = j;k <= i;k ++)
				dp[i][j] = max(dp[k-1][j-1] * a[k][i], dp[i][j]);
		
	cout << dp[n][m+1] << endl;
	
	return 0;
}

采用高精度的AC代码：

#include<bits/stdc++.h>
using namespace std;
const int N = 100, M = 20;

int n, m;
long long dp[N][M], a[N][N];
string num, f[N][M], b[N][N];

string max (string s1, string s2) {
	// 先比较长度，长者大；否则比较高位数，大者大；否则一样。
	int n1 = s1.size(), n2 = s2.size();
	if (n1 != n2) return n1>n2?s1:s2;
	for (int i = 0;i < n1;i ++)
		if (s1[i] != s2[i]) 
			return s1[i] > s2[i] ? s1 : s2;
	return s1; 
}

string add (string s, int c) {
	int nn = s.size();
	for (int i = nn-1;i >= 0;i --) {
		int t = s[i] - '0';
		s[i] = (t + c) % 10 + '0';
		c = (t + c) / 10;
	}
	return (c == 0 ? s : "1" + s);
}

string mul (string s1, string s2) {
	// 比较标准的高精度加法模板
	int len1 = s1.size(), len2 = s2.size();
	int c[len1+len2+1], a[len1+1], b[len2+1];
	s1 = '.'+s1; s2 = '.'+s2;
	memset(c, 0, sizeof c);
	
	for(int i=1;i<=len1;i++)
        a[len1-i+1]=s1[i]-'0';
    for(int i=1;i<=len2;i++)
        b[len2-i+1]=s2[i]-'0';
 
    for(int i=1;i<=len1;i++)   
        for(int j=1;j<=len2;j++){
            c[i+j-1]+=a[i]*b[j];
            c[i+j]+=c[i+j-1]/10;
            c[i+j-1]%=10;
        }
 
    int len=len1+len2-1;
    if(c[len+1]!=0) len++;
 
 	string res = "";
    for(int i=len;i>=1;i--) {
    	res += to_string(c[i]); // to_string将int转为string，string只能拼接string，不能拼接char
	}
    return res;

	/*
	// 这种方式也可以，不标准
	reverse(s1.begin(), s1.end());
	reverse(s2.begin(), s2.end());
	int n1 = s1.size(), n2 = s2.size(), c = 0;
	int res[n1+n2+1], len = n1 + n2 - 1;
	memset(res, 0, sizeof res);
	
	for (int i = 0;i < n1;i ++)
		for (int j = 0;j < n2;j ++) {
			int t = (s1[i]-'0') * (s2[j]-'0');
			res[i+j] += t;
			res[i+j+1] += res[i+j]/10;
			res[i+j] %= 10;
		}
	
	while(res[len]) len ++; // 将len初始化为n1+n2-1是因为两个数乘积最少位数为两个数位数之和-1
	string ans(len, '0'); // 必须这样初始化，否则直接将字符赋值到string上进行输出是空的
	for (int i = 0;i < len;i ++) 
		ans[len-i-1] = res[i] + '0';
	return ans;
	*/ 
}

string mul (string s) {
	// 默认乘以10
	// 乘以10，直接加个0就行 
	return s + "0" ;
}

// main部分的代码几乎没变，就是加了几个函数、改了几个变量名
int main()
{
	cin >> n >> m;
	cin >> num;
	
	for (int i = 0;i < n;i ++)
		for (int j = i;j < n;j ++)
			b[i+1][j+1] = add(mul(b[i+1][j]), num[j]-'0');
			
	for (int i = 1;i <= n;i ++)
		f[i][1] = b[1][i];
	
	for (int i = 1;i <= n;i ++)
		for (int j = 2;j <= m+1 && j <= i;j ++) 
			for (int k = j;k <= i;k ++) {
				f[i][j] = max(mul(f[k-1][j-1], b[k][i]), f[i][j]);		
			}

		
	cout << f[n][m+1] << endl;
	
	return 0;
}

/*
几组样例
 
15 5
234567891987127
93730005427968

40 5
2323456789194567892345678919871271987127
1799210433424844361866616035388793782468

40 5
2323456789194567892345678919871271900000
1658633212909120759464865210877632800000
*/