=分治法求整数划分
把一个正整数写成n=m1+m2+…+mi,其中m1>=m2>=m3>=…>=mi,1<=mi<=n,i>=1,则{m1,m2…,mi}叫做n的一个划分。
算法思想
设正整数n的不同划分个数称为正整数的划分数,记作p(n)。
如果{m1,m2,…mi}中最大的值不超过m,称他为n的一个m划分,记作q(n,m)。
例如n=4有5个划分:{4};{3+1};{2+2};{2+1+1};{1+1+1+1}。
该问题是求出n的所有划分个数,即f(n,n)。
算法分析
1、n==1:m取任意值,只有一种划分{1}
2、m==1:n取任意值,只有一种划分{1+1+1+…+1}
3、n==m:根据划分中是否包含n,分为:
(1)划分中包含n,只有一种{n}
(2)划分中不包含n,这时划分中最大的数字也一定比n小,即n所有的n-1划分
(3)因此q(n,n)=1+q(n,n-1)
4、n<m:实际上不允许出现n<m情况,这时即为q(n,n)
5、n>m:根据划分中是否包含最大值m,可以分为
(1)包含m的情况,即{m,{x1,x2,….xi}并且x1+x2+…+xi=n-m,因此这种情况下为q(n-m,m)。
(2)划分中不包含m且划分中所有值都比m小,即n的所有(m-1)划分,为q(n,m-1)
6、由以上分析得:
q(n,m)=1; n==1|m==1
q(n,m)=q(n,n); n<m
q(n,n)=1+q(n,m-1); n==m
q(n,m)=q(n-m,m)+q(n,m-1); m<n
源代码如下:
#include<iostream>
using namespace std;
int f(int n1,int m1)
{
if(n1==1||m1==1)
return 1;
else if(n1==m1)
{
return 1+f(n1,n1-1);
}
else if(n1>m1)
{
return f(n1-m1,m1)+f(n1,m1-1);
}
else if(n1<m1)
{
return f(n1,n1);
}
}
int main()
{
int p;
cin>>p;
cout<<f(p,p)<<endl;
return 0;
}
输入与输出:
5
7
6
11
4
5