首页 > 其他 > 详细

P1018 乘积最大

时间:2020-04-25 21:48:18      阅读:66      评论:0      收藏:0      [点我收藏+]

p1018 乘积最大

题目描述

今年是国际数学联盟确定的“2000――世界数学年”,又恰逢我国著名数学家华罗庚先生诞辰90周年。在华罗庚先生的家乡江苏金坛,组织了一场别开生面的数学智力竞赛的活动,你的一个好朋友XZ也有幸得以参加。活动中,主持人给所有参加活动的选手出了这样一道题目:

设有一个长度为N的数字串,要求选手使用K个乘号将它分成K+1个部分,找出一种分法,使得这K+1个部分的乘积能够为最大。

同时,为了帮助选手能够正确理解题意,主持人还举了如下的一个例子:

有一个数字串:312, 当N=3,K=1时会有以下两种分法:

3×12=36

31×2=62

这时,符合题目要求的结果是:31×2=62

现在,请你帮助你的好朋友XZ设计一个程序,求得正确的答案。

输入格式

程序的输入共有两行:

第一行共有2个自然数N,K(6≤N≤40,1≤K≤6)

第二行是一个长度为N的数字串。

输出格式

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

输入输出样例

输入

4  2
1231

输出

62

题解

在n个数中插入k个乘号.

首先看到位数范围40位, 就能想到, 乘法操作后, 可能会爆long long, 保险起见, 上高精.

对于这个题的设定, 第i个乘号确定后将数分成的两部分, 乘号数量一定(左i-1个和右k-i个), 乘积互不影响. 所以只要保证两边单独的乘积最大, 就能保证最后的乘积最大.

这样就设计出了状态:

f(i,j)表示第j个乘号放在第i个数字之后, 前i个数字的乘积最大值.

设A[l,r]表示用A的第l位到第r位形成的高精数

那么就可以写出状态转移方程:

\[f[i][j]=max(f[j-1][j-1]*A[j,i],f[j][j-1]*A[j+1,i],...,f[i-1][j-1]*A[i,i]) \]

最后答案也可以通过f数组推出.

\[ans=max(f[k][k]*A[k+1,n],...,f[n-1][k]*A[n,n]) \]

代码

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int n, K, A[50], C[1005], b[50][10][1005], d[1005], ji[1005], ans[1005];
char a[50];
bool flg[1005];
void take(int l, int r) {//从A中取[l,r]区间的数字组成一个数, 存入高精数d.
	d[0] = r - l + 1;
	for (int ti = l; ti <= r; ti++) {
		d[ti - l + 1] = A[ti];
	}
	return;
}
void print(int x[]) {//打印高精数x
	for (int pi = x[0]; pi >= 1; pi--) {
		cout << x[pi];
	}
	cout << endl;
	return;
}
void times(int x[], int y[]) {//计算高精数x,y的乘积, 并且将积存入高精数ji
	memset(ji, 0, sizeof(ji));
	int xx = 0;
	for (int ti = 1; ti <= x[0]; ti++)//枚举因数A的每一位 
	{
		xx = 0;
		for (int tj = 1; tj <= y[0]; tj++) {//枚举因数B的每一位 
			ji[ti + tj - 1] += xx + (x[ti] * y[tj]);
			xx = ji[ti + tj - 1] / 10;
			ji[ti + tj - 1] %= 10;
		}
		ji[ti + y[0]] += xx;
	}
	ji[0] = x[0] + y[0];
	while (ji[ji[0]] == 0 && ji[0] > 1) {
		ji[0]--;
	}
	return;
}
bool cmp(int x[],int y[]) {//比较高精数x,y, 若x>y则返回true, 反之返回false
	if (x[0] < y[0]) {
		return false;
	}
	if (x[0] > y[0]) {
		return true;
	}
	for (int ci = x[0]; ci >= 1; ci--) {
		if (x[ci] < y[ci]) {
			return false;
		}
		if (x[ci] > y[ci]) {
			return true;
		}
	}
	return false;
}
int main() {
	cin >> n >> K;
	for (int i = 1; i <= n; i++) {
		cin >> a[i];
		A[n - i + 1] = a[i] - 48;//倒着存, 高精常规操作, 不解释
	}
	A[0] = n;//高精数第0位存长度
	for (int i = 1; i < n; i++) {//枚举
		take(1, i);//这是取b[i][1]的值
		for (int ij = 0; ij <= d[0]; ij++) {//只放一个乘号, 结果固定
			b[i][1][ij] = d[ij];
		}
		for (int j = 2; j <= K; j++) {//枚举乘号
			for (int k = j; k < i; k++) {//枚举区间右界
				take(k + 1, i);//取A数的第k+1位到第i位作为一个乘数
				//print(b[k][j - 1]);//注释语句均为调试所用, 保留原因见最后一行
				times(b[k][j - 1], d);//代入第二个乘数,计算乘法
				//print(b[k][j - 1]);
				//cout << i << " " << j << endl;
				if (cmp(ji, b[i][j])) {//比较并更新(转移)
					for (int ci = 0; ci <= ji[0]; ci++) {
						b[i][j][ci] = ji[ci];
					}
				}
			}
		}
	}
	for (int i = K; i <= n - 1; i++) {//枚举放完所有的乘号后, 所得的答案
		take(i + 1, n);
		times(b[i][K], d);//乘法
		//cout << i << " " << K<< endl;
		//print(d);
		//print(b[i][K]);
		//print(ji);
		if (cmp(ji, ans)) {//比较并更新答案
			for (int j = ji[0]; j >= 0; j--) {
				ans[j] = ji[j];
			}
		}
	}
	print(ans);//输出答案
	return 0;
}//保留调试语句刚好100行

P1018 乘积最大

原文:https://www.cnblogs.com/Wild-Donkey/p/12775410.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!