P4157 [SCOI2006]整数划分（数学推导+大数）

最新推荐文章于 2021-11-26 16:05:46 发布

原创最新推荐文章于 2021-11-26 16:05:46 发布 · 415 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#ACM

数论专栏收录该内容

23 篇文章

订阅专栏

本文介绍了如何利用数学技巧，如均值不等式和对数分析，解决一道关于将正整数拆分为多个部分以最大化乘积的问题。通过讨论n对3的余数，作者给出了一个简单的算法，优先选择3和2的组合，展示了如何在有限时间内找到最优解。

原题题面

读入一个正整数 $n$ （ $10 \leq n \leq 31000$ ）。要求将 $n$ 写成若干个正整数之和，并且使这些正整数的乘积最大。

例如， $n = 13$ ，则当n表示为 $4 + 3 + 3 + 3$ （或 $2 + 2 + 3 + 3 + 3$ ）时，乘积= $108$ 为最大。

输入格式

只有一个正整数 $n$ （ $10 \leq n \leq 31000$ ）

输出格式

第1行输出一个整数，为最大乘积的位数。
第2行输出最大乘积的前100位，如果不足100位，则按实际位数输出最大乘积。
（提示：在给定的范围内，最大乘积的位数不超过5000位）。

输入样例

13

输出样例

3
108

题面分析

简单概括一下就是:已知 $∑xi=N\sum{x_i}=N$ ，求 $max{∏xi}max\{\prod{x_i}\}$ .
看到 $∑xi\sum{x_i}$ 和 $∏xi\prod{x_i}$ ~~很容易~~ 就能想到均值不等式。
假设 $x_i$ 有 $n$ 项，我们有
$(∑xin)n≥∏xi(\frac{\sum{x_i}}{n})^n\geq\prod{x_i}$
因为取到等号时， $x_i$ 两两相等，因此 $x_i$ 要取的尽量平均，于是问题变成了求出
$(∑xin)n=(Nn)n(\frac{\sum{x_i}}{n})^n=(\frac{N}{n})^n$ 的最大值
我们可以用一些高数知识对它进行处理
记 $y=(Nn)ny=(\frac{N}{n})^n$ ,两边套 $I n$ ，得到 $I n y = n I n N - n I n n$
两边再求导，得到 $y′y=InN−Inn−1\frac{y'}{y}=InN-Inn-1$ (这里不再赘述为什么左边不是 $1y\frac{1}{y}$ )
移项得到 $y′=(InN−Inn−1)∗(Nn)ny'=(InN-Inn-1)*(\frac{N}{n})^n$
令 $y^{'} = 0$ ，即 $(InN−Inn−1)∗(Nn)n=0(InN-Inn-1)*(\frac{N}{n})^n=0$
因为 $(Nn)n(\frac{N}{n})^n$ 必定大于0，故原式等价于求 $(I n N - I n n - 1) = 0$ ，易知 $n=Nen=\frac{N}{e}$
即把每个数尽可能地分成e。
且 $y$ 在 $(0,Ne](0,\frac{N}{e}]$ 上单调递增, $[Ne,∞)[\frac{N}{e},\infty)$ 上单调递减。
但考虑到在整数范围， $Ne\frac{N}{e}$ 取不到，因此候选项有 $N2,N3\frac{N}{2},\frac{N}{3}$ 两个数，于是问题就转化为哪个优先取的问题。
代入方程，得到 $2N2<3N32^{\frac{N}{2}}<3^{\frac{N}{3}}$
如何证明？
两边取对数得到 $N2In2和N3In3\frac{N}{2}In2和\frac{N}{3}In3$
注意到 $In22<In33\frac{In2}{2}<\frac{In3}{3}$
故 $2N2<3N32^{\frac{N}{2}}<3^{\frac{N}{3}}$ ,
故优先取3，然后取2。 证毕。
于是这就变成了憨憨题，我们根据 $3n\ mod\ 3$ 的不同结果来讨论取的方法。
$3=0n\ mod\ 3=0$ 时，无脑拿3,答案就是 $3N33^\frac{N}{3}$
$3=1n\ mod\ 3=1$ 时，拿出一个3变两个2,答案就是 $3N−13−1∗223^{\frac{N-1}{3}-1}*2^2$
$3=2n\ mod\ 3=2$ 时，拿一个2之后全变3,答案就是 $3N−23∗23^{\frac{N-2}{3}}*2$

AC代码(70ms)

import java.math.BigInteger;
import java.util.Scanner;

import static java.lang.Math.min;
/**
 * @author DrGilbert
 */
public class Main
{
    public static void main (String[] args)
    {
        Scanner cin = new Scanner (System.in);
        int n;
        while(cin.hasNext ())
        {
            n=cin.nextInt ();
            BigInteger ans;
            if(n%3==0)
            {
                ans=BigInteger.valueOf (3).pow ((n/3));
            }
            else if(n%3==1)
            {
                ans=BigInteger.valueOf (3).pow ((n/3-1)).multiply (BigInteger.valueOf (4));
            }
            else
            {
                ans=BigInteger.valueOf (3).pow ((n/3)).multiply (BigInteger.valueOf (2));
            }
            String ans1=ans.toString ();
            System.out.println (ans1.length ());
            for(int i=0;i<min(100,ans1.length ());i++)
            {
                System.out.print (ans1.charAt (i));
            }
            System.out.println ();
        }
    }
}