洛谷P1249 最大乘积
题目描述
一个正整数一般可以分为几个互不相同的自然数的和,如 3 = 1 + 2 3=1+2 3=1+2, 4 = 1 + 3 4=1+3 4=1+3, 5 = 1 + 4 = 2 + 3 5=1+4=2+3 5=1+4=2+3, 6 = 1 + 5 = 2 + 4 6=1+5=2+4 6=1+5=2+4。现在你的任务是将指定的正整数 n n n分解成若干个互不相同的自然数的和,且使这些自然数的乘积最大。
输入格式
只一个正整数 n , ( 3 ≤ n ≤ 10000 ) n,(3\leq n\leq 10000) n,(3≤n≤10000)。
输出格式
第一行是分解方案,相邻的数之间用一个空格分开,并且按由小到大的顺序。
第二行是最大的乘积。
输入#1:
10
输出#1:
2 3 5
30
开始没看到是不同的正整数。如果正整数随意,我们知道可以尽量多地分成3来处理。但是如果要正整数各不相同,而且个数也不定,似乎就比较棘手。 n n n的范围会上万,直接搜索的话肯定是不行的,所以去思考是否存在某种取数的策略。
首先从小学我们就知道,一个正整数拆成两个正整数之和,这两个数只要不是特别小,积普遍是大于和的,所以直觉上应该要尽量地拆这个数。然后应该尽量拆成比较小的数,因为大的数分成两个数的积之后比这个数更大。然后就产生一种贪心策略:把 n n n拆成 2 + 3 + ⋯ + k 2+3+\cdots+k 2+3+⋯+k的形式,这样每个数都尽可能小,积也会更大。但是 n n n不一定能刚好表示成这个形式,所以还需要处理一下边界情况。
那么用 i i i表示当前加到的数,从 2 2 2开始增大,计算当前的和。如果刚好等于 n n n,那么ok很完美,直接去计算积了。如果大于 n n n,也就是说我们得到 ∑ k = 2 i − 1 k < n \sum\limits_{k=2}^{i-1}k<n k=2∑i−1k<n并且 ∑ k = 2 i k > n \sum\limits_{k=2}^ik>n k=2∑ik>n。那么就存在一个取舍的问题,需要用一个大的数替换掉 2 2 2到 i − 1 i-1 i−1之间的一个数。显然应该用 i i i替换掉 ∑ k = 2 i k − n \sum\limits_{k=2}^ik-n k=2∑ik−n。因为,考虑把 i + 1 i+1 i+1放进去替换,假设之前替换掉 t t t,现在为了保证和为 n n n就要替换掉 t + 1 t+1 t+1,比较前后两种替换不同的数字之积存在大小关系 i ( t + 1 ) > ( i + 1 ) t i(t+1)>(i+1)t i(t+1)>(i+1)t,所以前者更大。同理可证用更大的数替换积会更小。
那么贪心策略就证明了正确性,虽然不是很严谨,但是基本没问题。可以举几个例子验证一下,一举就发现了问题,假设 n o w s u m nowsum nowsum表示当前的和,如果 n o w s u m = = n + 1 nowsum==n+1 nowsum==n+1,那需要替换1,但是之前我们不选1,所以这种情况需要特判一下。显然此时为了替换掉2,不得不放弃用 i i i而使用 i + 1 i+1 i+1。
最后的最后,还没有结束,由于 n n n可以达到 10000 10000 10000,乘下来的结果绝对爆了。所以要用string做一下大数乘法,代码如下:
#include <iostream>
#include <stdio.h>
#include <string>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <string.h>
#include <vector>
using namespace std;
int n;
string mul(string a,string b){
int m=a.size(),n=b.size();
int c[m+n+1];
memset(c,0,sizeof(c));
for(int i=m-1;i>=0;i--){
for(int j=n-1;j>=0;j--){
c[(m-1-i)+(n-1-j)]+=(a[i]-'0')*(b[j]-'0');
}
}
for(int i=0;i<m+n;i++){
c[i+1]+=c[i]/10;
c[i]%=10;
}
int flag=0;
string res;
for(int i=m+n;i>=0;i--){
if(c[i]||i==0){
flag=1;
}
if(flag){
res.push_back(c[i]+'0');
}
}
return res;
}
int main() {
cin>>n;
vector<int> v;
int nowsum=0,i=2,ignore=-1;
while(1){
nowsum+=i;
v.push_back(i);
if(nowsum==n){
break;
}
else if(nowsum==n+1){
ignore=2;
v.pop_back();
v.push_back(i+1);
break;
}
else if(nowsum>n){
ignore=nowsum-n;
break;
}
i++;
}
string res="1";
for(auto x:v){
if(x!=ignore){
cout<<x<<(x==v.back()?"":" ");
res=mul(res,to_string(x));
}
}
cout<<endl<<res<<endl;
return 0;
}