一、 题目
最优分解问题
二、 问题描述
设n是一个正整数。现在要求将n分解为若干互不相同的自然数的和,且使这些自然数的乘积最大。
算法设计:对于给定的正整数n,计算最优分解方案。
数据输入:由文件 input.txt 提供输入数据。文件的第一行是正整数n。
结果输出:将计算出的最大乘积输出到文件output.txt。
输入文件示例 | 输出文件示例 |
input.txt | output.txt |
10 | 30 |
三、 算法基本思想
先对整数分解分析可以发现如下结论:
若 a + b = const,则 |a - b| 越小,a·b越大。
根据原问题的描述,需要将正整数n分解为若干互不相同的自然数的和,同时又要使自然数的乘积最大。当n<4 时对n的分解的乘积是小于n的;当n大于或等于4时,n = 1 + (n-1)因子的乘积也是小于n的,所以n = a + (n-a), 2≤a≤n-2,可以保证乘积大于n,即越分解乘积越大。
因此基于之前发现的结论和上面的分析,可以采用如下贪心策略:将n分成从2开始的连续自然数的和,如果最后剩下一个数,将此数在后项优先的方式下均匀地分给前面各项。
该贪心策略首先保证了正整数所分解出的因子之差的绝对值最小,即|a - b|最小;同时又可以将其分解成尽可能多的因子,且因子的值较大,确保最终所分解的自然数的乘积可以取得最大值。
四、 算法实现(C语言)
1 /*************************************************
2 File Name: optimal_decomposition.cpp
3 Functions: 最优分解问题求解
4 将正整数n分解为若干个自然数的和,
5 使这些自然数的乘积最大。
6 Author: Jeccey
7 Created: Jeccey
8 Last Change: 2013-07-03
9 **************************************************/
10
11 #include <stdio.h>
12 #include <stdlib.h>
13 #include <string.h>
14
15 const int MAX = 51;
16 int a[MAX - 1]; //记录所分解的自然数
17
18 int PIntdecomp(int n);
19
20 int main()
21 {
22 int n;
23 char ch;
24
25
26 FILE *fp1, *fp2;
27 if( (fp1=fopen("input.txt","r")) == NULL){
28 printf("cannot open file!\n");
29 exit(0);
30 }
31 fscanf(fp1,"%d",&n);
32
33 memset(a, 0, MAX*sizeof(int));
34
35 if( (fp2=fopen("output.txt","w")) == NULL){
36 printf("cannot open file!\n");
37 exit(0);
38 }
39 fprintf(fp2,"%d\n",PIntdecomp(n));
40
41 return 0;
42 }
43
44 /****************************************************
45 Function name: PIntedecomp
46 Description: 将正整数n分解为若干个自然数的和,
47 保证所分解自然数的乘积最大并返回。
48 Arguments: n (int类型,且n<=437)
49
50 Returns: mresult, int 类型
51 n所分解的自然数的乘积,返回0表示错误
52 *****************************************************/
53 int PIntdecomp(int n)
54 {
55 int k,j,mresult;
56
57 if (n < 1 ) //非正整数返回0
58 return 0;
59 if (n < 5) //若n<5,结果是其本身
60 return n;
61 else{ //若 n>= 5
62 k = 0;
63 a[k] = 2;
64 n -= 2;
65 /* 贪心策略:先从2开始分成连续自然数的和 */
66 for (; n > a[k]; ){
67 a[++k] = a[k-1] + 1;
68 n -= a[k];
69 }
70 /* 如果剩下一个数,将其按后项优先的方式
71 均匀分给前面各项 */
72 if (n == a[k]){
73 a[k]++;
74 n--;
75 }
76 for (j = 0; j < n; j++){
77 a[k-j]++;
78 }
79 /* 计算连乘积 */
80 for(mresult = 1, j = 0; j <= k; j++){
81 mresult *= a[j];
82 }
83
84 return mresult;
85 }
86 }
View Code
五、关于自然数因子可重复情况
n<4 —— 乘积<n
n≥4—— 1. 1和(n-1)→乘积<n,无效;
2. a和(n-a),且 2≤a≤n-2,分解使乘积增加(可以证明)。
若因子大于4则继续分解,直至a和(n-a)都小于4。
因此最终分解结果:全是2或3;又3×3>2×2×2,所以有三个2换成两个3。
最终结果:“M个3”或“M个3和一个2”或“M个三和两个2”
即 R=3^M 或 R=(3^M)×2或 R=(3^M)×4