【洛谷】P1249 最大乘积

洛谷P1249 最大乘积

题目描述

一个正整数一般可以分为几个互不相同的自然数的和，如 $3=1+2$，$4=1+3$，$5＝1+4=2+3$，$6=1+5＝2+4$。现在你的任务是将指定的正整数 $n$分解成若干个互不相同的自然数的和，且使这些自然数的乘积最大。

输入格式

只一个正整数 $n,(3\le n\le 10000)$。

输出格式

第一行是分解方案，相邻的数之间用一个空格分开，并且按由小到大的顺序。

第二行是最大的乘积。

输入#1：

10

输出#1:

2 3 5
30

       开始没看到是不同的正整数。如果正整数随意，我们知道可以尽量多地分成3来处理。但是如果要正整数各不相同，而且个数也不定，似乎就比较棘手。$n$的范围会上万，直接搜索的话肯定是不行的，所以去思考是否存在某种取数的策略。

       首先从小学我们就知道，一个正整数拆成两个正整数之和，这两个数只要不是特别小，积普遍是大于和的，所以直觉上应该要尽量地拆这个数。然后应该尽量拆成比较小的数，因为大的数分成两个数的积之后比这个数更大。然后就产生一种贪心策略：把$n$拆成$2+3+\cdots +k$的形式，这样每个数都尽可能小，积也会更大。但是$n$不一定能刚好表示成这个形式，所以还需要处理一下边界情况。

       那么用$i$表示当前加到的数，从$2$开始增大，计算当前的和。如果刚好等于$n$，那么ok很完美，直接去计算积了。如果大于$n$，也就是说我们得到$\sum _{k=2}^{i-1}k<n$并且$\sum _{k=2}^{i}k>n$。那么就存在一个取舍的问题，需要用一个大的数替换掉$2$到$i-1$之间的一个数。显然应该用$i$替换掉$\sum _{k=2}^{i}k-n$。因为，考虑把$i+1$放进去替换，假设之前替换掉$t$，现在为了保证和为$n$就要替换掉$t+1$，比较前后两种替换不同的数字之积存在大小关系$i(t+1)>(i+1)t$，所以前者更大。同理可证用更大的数替换积会更小。

       那么贪心策略就证明了正确性，虽然不是很严谨，但是基本没问题。可以举几个例子验证一下，一举就发现了问题，假设$nowsum$表示当前的和，如果$nowsum==n+1$，那需要替换1，但是之前我们不选1，所以这种情况需要特判一下。显然此时为了替换掉2，不得不放弃用$i$而使用$i+1$。

       最后的最后，还没有结束，由于$n$可以达到$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;
}