题目链接:http://bailian.openjudge.cn/practice/4117/
总时间限制: 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 ",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 ",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++代码:(感觉其实就是回溯的思路)
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
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 << ' '; 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 << ": "; 74 int count = ExpandedGoldbach(i); 75 cout << "Total: " << count << " "; 76 } 77 }
C语言代码如下:
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
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: ",i); 22 int count = work(i); 23 printf("Total: %d ",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(" "); 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