代码片段_7215:简单的整数划分问题

总时间限制: 100ms 内存限制: 65536kB
描述
将正整数n 表示成一系列正整数之和，n=n1+n2+…+nk, 其中n1>=n2>=…>=nk>=1 ，k>=1 。
正整数n 的这种表示称为正整数n 的划分。正整数n 的不同的划分个数称为正整数n 的划分数。

输入
标准的输入包含若干组测试数据。每组测试数据是一个整数N(0 < N <= 50)。
输出
对于每组测试数据，输出N的划分数。
样例输入
5
样例输出
7
提示
5, 4+1, 3+2, 3+1+1, 2+2+1, 2+1+1+1, 1+1+1+1+1

注意对比题目：

（1）数的划分（把n分为k份，有点类似于放苹果，但放苹果这道题目里面盘子可以为空，它不允许为空。）

（2）复杂的整数划分

算法（一）递归法

把一个正整数n写成如下形式： n=m1+m2+...+mi; （其中mi为正整数，并且1 <= mi <= n），则m1,m2,...,mi为n的一个划分。如果m1,m2,...,mi中的最大值不超过m，即max(m1,m2,...,mi)<=m，则称它属于n的一个m划分。这里我们记n的m划分的个数为f(n,m);
例如当n=4时，他有5个划分，4,3,1,2,2,2,1,1,1,1,1,1;注意4=1+3 和 4=3+1被认为是同一个划分。

根据n和m的关系，考虑以下几种情况：
(1)当n=1时，不论m的值为多少（m>0)，只有一种划分即1;
(2) 当m=1时，不论n的值为多少，只有一种划分即n个1，1,1,1,...,1;
(3) 当n=m时，根据划分中是否包含n，可以分为两种情况：
　　(a). 划分中包含n的情况，只有一个即n；
　　(b). 划分中不包含n的情况，这时划分中最大的数字也一定比n小，即n的所有(n-1)划分。因此 f(n,n) =1 + f(n,n-1);
(4) 当n<m时，由于划分中不可能出现负数，因此就相当于f(n,n);
(5) 但n>m时，根据划分中是否包含最大值m，可以分为两种情况：
　　(a). 划分中包含m的情况，即m, x1,x2,...xi, 其中x1,x2,... xi 的和为n-m，可能再次出现m，因此是（n-m）的m划分，因此这种划分个数为f(n-m, m);
　　(b). 划分中不包含m的情况，则划分中所有值都比m小，即n的(m-1)划分，个数为f(n,m-1);因此 f(n, m) = f(n-m, m)+f(n,m-1);

综合以上情况，我们可以看出，上面的结论具有递归定义特征，其中（1）和（2）属于回归条件，（3）和（4）属于特殊情况，将会转换为情况（5）。而情况（5）为通用情况，属于递推的方法，其本质主要是通过减小m以达到回归条件，从而解决问题。其递推表达式如下：
f(n, m) =1; 　　　　　　　　　　 (n=1 or m=1)
=f(n, n); 　　　　　　　　 (n<m)
=1+ f(n, m-1); 　　　　 (n=m)
=f(n-m,m)+f(n,m-1);　　 (n>m)

 1 #include <stdio.h>
 2 long long GetPartitionCount(int n,int max)
 3 
 4     if(n==1||max==1) return 1;
 5     else if(n<max) return GetPartitionCount(n,n);
 6     else if(n==max) return 1+GetPartitionCount(n,max-1);
 7     else return GetPartitionCount(n,max-1)+GetPartitionCount(n-max,max);
 8 
 9 int main(int argc, char *argv[])
10 
11     int n;
12     while(scanf("%d",&n)!=EOF)
13     
14         printf("%lld\\n",GetPartitionCount(n,n));
15     
16     return 0;
17

算法（二）动态规划

这道题就是整数的划分，其实定义好了状态，就是简单的动态规划的递推。

我们定义dp[n][k]表示将n进行划分，最大的数不超过k的方案有多少种，那么我们可以得到如下的递推方案：

dp[n][k] = dp[n][k-1] + dp[n-k][k];

其中的dp[n][k-1]便是将n进行进行整数的划分，最大的数不超过k-1的方案数；dp[n-k][k]表示拿出一个k后，剩下的数被不超过k的数的表示的方案数。

其中当k>n的时候，则和dp[n][n]的值相同，下面通过递推的方式求出所有的解，然后对应的输入n，输出dp[n][n]就行了。

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdlib>
 4 #include<cstdio>
 5 using namespace std;
 6 
 7 const int MAX = 122;  
 8 int dp[MAX][MAX];
 9 
10 void Dynamic()
11 
12     for(int i=1; i<MAX; i++)
13     
14         dp[i][1] = dp[1][i] = dp[0][i] = 1;
15     
16     for(int i=2; i<MAX; i++)
17     
18         for(int j=2; j<MAX; j++)
19         
20             if(j<=i) dp[i][j] = dp[i][j-1] + dp[i-j][j];
21             else dp[i][j] = dp[i][i];
22         
23     
24   
25 int main()  
26   
27     int n;
28     Dynamic();
29     while(scanf("%d",&n)!=EOF)
30     
31         printf("%d\\n",dp[n][n]);
32     
33     return 0;  
34

这道题的动规思路，下面的论述比较清晰一些：

整数划分问题
数 n 的划分是将 n 表示成多个正整数之和的形式，划分可以分为两种情况：
第一种情况：划分的多个正整数中，正整数的数量是任意的。

这又可以分为划分的正整数中，正整数可以相同与不同两类

1. 划分的多个正整数可以相同, 递推方程可以表示为：
(1)   dp[n][m]= dp[n][m-1]+ dp[n-m][m]
分析dp[n][m]表示整数 n 的划分中，每个数不大于 m 的划分数。则划分数可以分为以下两种情况:
a. 划分中每个数都小于 m，相当于每个数不大于 m- 1, 故划分数为 dp[n][m-1].
           b. 划分中有至少一个数为 m. 那就在 n中减去 m , 剩下的就相当于把 n-m 进行划分，故划分数为 dp[n-m][m];
   (2)   dp[n][m]= dp[n][m+1]+ dp[n-m][m]
   其中dp[n][m]表示整数 n 的划分中，每个数不小于 m 的划分数。同理可证明该式。

2. 划分的多个正整数互不相同，递推方程可以表示为：
(1)    dp[n][m]= dp[n][m-1]+ dp[n-m][m-1]
   分析： dp[n][m]表示整数 n 的划分中，每个数不大于 m 的划分数。同样划分情况分为以下两种情况：
   a. 划分中每个数都小于 m, 相当于每个数不大于 m- 1,划分数为 dp[n][m-1].
   b. 划分中有一个数为 m. 在 n 中减去 m, 剩下相当对n- m 进行划分，并且每一个数不大于 m- 1，故划分数为 dp[n-m][m-1]
(2)    dp[n][m]= dp[n][m+1]+ dp[n-m][m]
   其中dp[n][m]表示整数 n 的划分中，每个数不小于 m 的划分数。

第二种情况：划分的多个正整数中，正整数的数量是固定的
把一个整数 n 无序划分成 k 份互不相同的正整数之和的方法总数。方程为： dp[n][k]= dp[n-k][k]+ dp[n-1][k-1];
   证明方法参考： http://www.mydrs.org/program/html/0369.htm
   另一种理解，总方法可以分为两类：
   第一类: n 份中不包含 1 的分法，为保证每份都 >= 2，可以先拿出 k 个 1 分
   到每一份，然后再把剩下的 n- k 分成 k 份即可，分法有: dp[n-k][k]
   第二类: n 份中至少有一份为 1 的分法，可以先那出一个 1 作为单独的1份，剩
   下的 n- 1 再分成 k- 1 份即可，分法有：dp[n-1][k-1]

类似问题：

M个小球装N个盒子，或者苹果装盘问题。比如：把M个球放到N个盒子,允许有空的盒子(不放球),有多少种放法?这些都属于典型的DP问题。

用F(m,n)表示有多少种放法:
如果m=0 或者 m=1 , F = 1
如果n=0 或者 n=1 , F =1
既F(0,0) = F(0,1) = F(1,0) = F(1,1) = 1
否则 F = F(m-n,n) + F(m,n-1)这就是DP的解空间递归解

关于整数的质因子和分解

【问题描述】
歌德巴赫猜想说任何一个不小于6的偶数都可以分解为两个奇素数之和。对此问题扩展，如果一个整数能够表示成两个或多个素数之和，则得到一个素数和分解式。对于一个给定的整数，输出所有这种素数和分解式。注意，对于同构的分解只输出一次（比如5只有一个分解2 + 3，而3 + 2是2 + 3的同构分解式）。

例如，对于整数8，可以作为如下三种分解：
(1) 8 = 2 + 2 + 2 + 2
(2) 8 = 2 + 3 + 3
(3) 8 = 3 + 5

【算法分析】
由于要将指定整数N分解为素数之和，则首先需要计算出该整数N内的所有素数，然后递归求解所有素数和分解即可。

原作者的C++代码：（感觉其实就是回溯的思路）

 1 #include <iostream>   
 2 #include <vector>   
 3 #include <iterator>   
 4 #include <cmath>   
 5 using namespace std;   
 6   
 7 // 计算num内的所有素数（不包括num）   
 8 void CalcPrimes(int num, vector<int> &primes)   
 9    
10     primes.clear();   
11     if (num <= 2)   
12         return;   
13        
14     primes.push_back(2);   
15     for (int i = 3; i < num; i += 2) 
16        
17         int root = int(sqrt(i));   
18         int j = 2;   
19         for (j = 2; j <= root; ++j)
20            
21             if (i % j == 0)   
22                 break;   
23            
24         if (j > root)   
25             primes.push_back(i);   
26        
27    
28   
29 // 输出所有素数组合(递归实现)   
30 int PrintCombinations(int num, const vector<int> &primes, int from, vector<int> &numbers)   
31    
32     if (num == 0) 
33        
34         cout << "Found: ";   
35         copy(numbers.begin(), numbers.end(), ostream_iterator<int>(cout, " "));   
36         cout << \'\\n\';   
37         return 1;   
38        
39        
40     int count = 0;   
41        
42     // 从第from个素数搜索，从而避免输出同构的多个组合   
43     int primesNum = primes.size();   
44     for (int i = from; i < primesNum; ++i)
45        
46         if (num < primes[i])   
47             break;   
48         numbers.push_back(primes[i]);   
49         count += PrintCombinations(num - primes[i], primes, i, numbers);   
50         numbers.pop_back();   
51        
52        
53     return count;   
54    
55   
56 // 计算num的所有素数和分解   
57 int ExpandedGoldbach(int num)   
58    
59     if (num <= 3)   
60         return 0;   
61        
62     vector<int> primes;   
63     CalcPrimes(num, primes);   
64        
65     vector<int> numbers;   
66     return PrintCombinations(num, primes, 0, numbers);   
67    
68   
69 int main()   
70    
71     for (int i = 1; i <= 20; ++i)
72        
73         cout << "When i = " << i << ":\\n";   
74         int count = ExpandedGoldbach(i);   
75         cout << "Total: " << count << "\\n\\n";   
76        
77

C语言代码如下：

 1 #include<stdio.h>
 2 #include<math.h>
 3 #include<string.h>
 4 
 5 #define maxNum 1000
 6 int primeNum[maxNum]=0,ansArr[maxNum]=0;
 7 int indexForPrimeNum,indexForAnsArr;
 8 
 9 int work(int num);//输出num的所有素数和分解的方案并返回其总方案数
10 int CalcPrimes(int num);//计算num内的所有素数（不包括num），结果保存在primeNum[]. 
11 
12 int PrintAns(int num,int from);
13 //回溯法的思想：从primeNum[]的第from个元素开始选择元素来累加构造num。
14 //构造结果放在ansArr[]中。寻找到一个构造方案后输出该方案.
15 //最终返回总的方案数
16 
17 int main(int argc, char *argv[])
18 
19     for (int i = 1; i <= 20; ++i)
20        
21         printf("When i = %d:\\n",i);   
22         int count = work(i);   
23         printf("Total: %d\\n\\n",count);   
24     
25     return 0;
26 
27 //计算num内的所有素数（不包括num），结果保存在primeNum[]. 返回数组元素个数 
28 int CalcPrimes(int num)
29 
30     int i,j,root,k=0;
31     
32     if(num<2) return k;
33     memset(primeNum,0,sizeof(primeNum));
34     
35     primeNum[k++]=2;
36     for(i=3;i<num;i++)
37     
38         root=sqrt(i);
39         for(j=2;j<=root;j++)
40         
41             if(i%j==0) break;
42         
43         if(j>root) primeNum[k++]=i;
44     
45     return k;
46 
47 //调用PrintAns()输出num的所有素数和分解的方案，返回其总方案数 
48 int work(int num)
49 
50     if(num<=3) return 0;
51     
52     indexForPrimeNum=CalcPrimes(num);
53     
54     memset(ansArr,0,sizeof(ansArr));
55     indexForAnsArr=0;
56     return PrintAns(num,0);
57 
58 
59 //回溯法的思想：从primeNum[]的第from个元素开始选择元素来累加构造num。
60 //构造结果放在ansArr[]中。寻找到一个构造方案后输出该方案.
61 //最终返回总的方案数
62 int PrintAns(int num,int from)
63 
64     int i;
65     if(num==0)
66     
67         printf("Found: ");
68         for(i=0;i<indexForAnsArr;i++) printf(" %d",ansArr[i]);
69         printf("\\n");
70         return 1;
71     
72     
73     int totalCount = 0;
74     //从第from个素数搜索，从而避免输出同构的多个组合
75     for(i=from;i<indexForPrimeNum;i++)
76     
77         if(num<primeNum[i]) break;
78         ansArr[indexForAnsArr++]=primeNum[i];
79         totalCount+=PrintAns(num-primeNum[i],i);
80         indexForAnsArr--;    
81     
82     return totalCount;
83

参考：

http://www.cppblog.com/superKiki/archive/2010/05/27/116506.html

http://blog.csdn.net/geniusluzh/article/details/8118683