.整数的拆解组合问题c语言,整数拆分的两种解法（已完成）

4.1 集合的分划和第二类Stirling数

定义1  (集合的划分)

设A 是有限集.A 的一族子集{Ai}称为是集A 的一个划分

,如果满足:

(1)每个子集Ai 都非空;

(2)这些子集两两不交;

(3)它们的并为A.

每个Ai称为划分的一个块,有k个块的划分称为k- 划分.

定义2 (第二类Stirling 数)

一个n元集的所有k- 划分的个数,用S(n,m) 表示,称为第二类Stirling 数.

注:将一个n元集划分成m块,就相当于将n个有区别的球放到m个相同的盒子中使得无一空盒.

例:求3元集和4元集的所有划分数.

由S(n,m) 的定义易知:

S(n,1)=1;

S(n,n)=1;

S(n,k)=0 (k ＞n)

S(n,k)对任何正整数n和k都有定义.

4.2第二类Stirling数的性质

定理1

第二类Stirling数S(n,k)满足递推关系S(n+1,k)=S(n,k-1)+kS(n,k).

证明:易验证k=n+1成立;k＞n+1 时也成立.当1≤k≤n时取n+1 元集A 的某一元素a, 将A 的k- 划分分成两类:一类是a 作为单独一块的,一类是a 不是单独一块的.第一类的划分数为S(n,k-1), 第二类的划分分成两步来实现,一步为将A 中除了元素a 划分成k块,第二步将a 放入某一块中.由乘法原理得第二类的划分数为kS(n,k).最后由加法原理,定理得证.

将第二类Stirling数列成表,它与杨辉三角类似.

定理2

第二类Stirling数S(n,k)满足下列性质

(1) S(n,2)=2的n-1次方-1;

(2) S(n,n-1)=C(n,2).

证明1:可由定理1递推得到.

证明2:由组合意义:

(1)S(n,2)是n元集A={a1,a2,…,an}的2-划分数.先取a1,则对于a2,…,an各有两种选择,或者和a1在同一块,或者不在同一块.但不能都和a1在同一块.

(2)S(n,n-1)是n元集A={a1,a2,…,an}的n-1-划分数.这样的划分必定一块有2个元素,其余的块都只有1个元素.只要确定了2个元素的块就确定了这个划分.

定理3

第二类Stirling数S(n,k)满足S(n+1,k)=C(n,k-1)S(k-1,k-1)+C(n,k)S(k,k-1)+…+C(n,n)S(n,k-1).

证明：S(n+1,k)是集合A={a1,a2,..,an，an+1}的k分划数。对于A的一个k分划，设包含an+1的那块是B,则其余的k-1块构成了A\B的一个k-1分划。反过来，给定A的一个含an+1 的子集B，若A\B含有至少k-1个元素，则A\B的一个k-1分划加上B就构成A的一个k分划。令C=A\B，易看出结论成立。

4. 3 正整数的有序分拆

定义3  (有序分拆)正整数n的一个k分拆把n表示成k个正整数的和n=n1+n2+…+nk(k ≥1)若分拆不仅与ni的值有关,也与它们的顺序有关.则称为有序k分拆.

例如4=2+1+1

=1+2+1

=1+1+2

是4的所有3个有序3分拆.注:正整数n的一个有序k-分拆,就相当于将n个无区别的球放到k个不同的盒

子中使得无一空盒.

定理4

正整数n的有序k分拆的个数为C(n-1,k-1).证明: 正整数n的有序k分拆n=n1+n2+…+nk等价于在n个竖杠的n-1个空隙中插入k-1个杠子.这样的方法数为C(n-1,k-1).

4. 4 正整数的无序分拆

定义4 (无序分拆)

正整数n的一个k分拆把n表示成k个正整数的和n=n1+n2+…+nk(k≥1)若分拆只与ni的值有关,与它们的顺序无关.则称为无序k分拆.我们用B(n,k)表示n的无序k分拆的个数,用B(n)表示n的所有分拆的个数.

注:正整数n的一个无序k- 分拆,就相当于将n个无区别的球放到k个相同的盒子中使得无一空盒.

例如: 4=4

=3+1=2+2

=2+1+1

=1+1+1+1

显然有

(1)B(n,k)=0 (k＞n);

(2)B(n,1)=1;

(3)B(n,n)=1;

(4)B(n)=B(n,1)+B(n,2)+…+B(n,n).

定理5

正整数n的无序k分拆的个数B(n,k)满足递推关系B(n+k,k)=B(n,1)+B(n,2)+…+B(n,k).

证明:我们考虑所有n的分成至多k个分部的分拆,这样的分拆总数为B(n,1)+B(n,2)+…+B(n,k).n的每个分成至多k个分部的分拆可表示为n=n1+n2+…+nm+0+…+0,这里n1≥n2≥…≥nm1≤m≤k

这个和式包含k项.它与n+k的下述k分拆一一对应.n+k=(n1+1)+(n2+1)+…+(nm +1)+1+…+1,这里n1≥n2≥…≥nm1≤m≤k

分类： 算法

2012-09-30 17:04

263人阅读

评论(0)

前几天在算法书上看到一个整数拆分的题目，觉得挺有意思，记录如下：

题目：给定一个整数n，输出这个整数拆分的可能总数

例如：n==6有

6

5+1

4+2    4+1+1

3+3    3+2+1     3+1+1+1

2+2+2     2+2+1+1        2+1+1+1+1

1+1+1+1+1+1

共11种分解方法，所以输出应该为11。

分析一

拆分按照因子从大到小排列，每一次拆分都可视为问题规模的减少，所以可以使用递归解决。

设q(n,m)为整数n使用不大于m的整数进行拆分的所有情况总数，因此有

1)当n==m时

可以分为两种情况，一个是使用n本身，只有一种情况。二个是使用不大于n-1的整数进行拆分。

所以此时q(n,m)=1+q(n,n-1)；

2)当n

使用比n大的数进行拆分没有意义，所以此时q(n,m)=q(n,n)

3)当n>m时

这时候有两种情况，一个是使用m对n进行拆分，一个是使用小于m的数对n进行拆分

对于第一个，使用m对n进行拆分，所以拆分出来的情况最大的数是m，剩下的所有数加起来为n-m，问题变成使用不大于m的整数对n-m进行拆分，所以此时为q(n-m,m)

对于第二个，使用小于m的数对n进行拆分，表示为q(n,m-1)

4)当n==1或者m==1时 只有一种情况，返回1。注意，这里可能会漏掉m==1的情况，实际上这种情况是存在的。

所以可以使用递归函数编程如下：

#include 

#include 

usingnamespacestd;

intq(intn,intm){

if(n==1||m==1){

return1;

}

if(n

if(n==m){

returnq(n,m-1)+1;

}

if(n>m)

returnq(n,m-1)+q(n-m,m);

}

voidmain(){

intn;

while(scanf("%d",&n)!=EOF){

printf("%d\n",q(n,n));

}

}

分析二

使用母函数法

整数分解用母函数可以这样理解，分别用任意个1，2，3，4，……，n可以加起来可以表示成n的种数。又因为当使用整数m对n进行分解时，所使用的次数不能多于n/m次，所以可以写下母函数如下：

G(x)＝(1+x^1+x^2+……+x^n)*(1+x^2+x^4+……+x^((n/2)*2))*……

＊(1+x^m+x^(2*m)+x^(3*m)+……+x^((n/m)*m))*……*(1+x^n)

在程序中使用set数组表示每一轮乘法后得到系数，c数组表示到现在为止乘法得到的系数总和。最后算出结果后x^n对应的系数则为可分解的可能数。代码如下：

#include

constintnum=1000;

intset[num];

intc[num];

voidinit(){

for(inti=0; i

c[i]=1;

set[i]=0;

}

}

intmain(){

intn;

intsum;

init();

while(scanf("%d",&n)!=EOF){

sum=0;

init();

for(inti=2; i<=n; i++){

for(intj=0; j<=n; j+=i){

for(intk=0; j+k<=n; k++){

set[k+j]+=c[k];

}

}

for(intx=0; x<=n; x++){

c[x]=set[x];

set[x]=0;

}

}

printf("%d\n",c[n]);

}

return0;

}

题目二：给定一个整数n，输出这个整数拆分的可能形式(即输出全部情况)

使用递归情况(输出实在麻烦。。弄了很久(>﹏

整个输出类似于一颗树，以分解6为例，过程如下图

代码如下：

#include 

//使用一个数组记录在递归过程中产生的前面需要重复输出的值

intset[100];

//用于在递归过程中判断是否递归到最深处，输出回车

intk;

//此函数表示使用不大于m的整数对n进行拆分的情况，i用于表示set数组已经存在的记录数长度

voidq(intn,intm,inti){

if(n==k&&n!=m){

//此时递归栈已经退回到某一分支的最上层，输出回车

//并重置计数器i为0

printf("\n");

i=0;

}

if(n==1){

//当n为1，意味者着只能表示1

printf("1 ");

return;

}

elseif(m==1){

//当m为1，意味着要输出n个m相加

for(inti=0; i

printf("1+");

printf("1 ");

return;

}

if(n

q(n,n,i);

}

if(n==m){

//当n等于m时，到达本次递归求和的一个叶子，此时需要输出多一个空格，表示下一次输出为另一个叶子

printf("%d ",n);

//在递归输出另一个叶子之前，将之前记录的在叶子之上的数一并输出，如上图示过程1

for(intj=0; j

printf("%d+",set[j]);

q(n,m-1,i);

}

if(n>m){

//如果n大于m，使用m作为分解，则要首先输出m+的形式

printf("%d+",m);

//记录下作为树干节点m的值并使i自增

set[i++]=m;

//递归输出m+以后的分解

q(n-m,m,i);

//递归完毕后需要将数组记录后退一个，回到上一个节点，如上图示过程2

i--;

//执行另一个分支，在下一次递归之前输出记录的数据，如上图示过程3

for(intj=0; j

printf("%d+",set[j]);

//递归输出另一分支情况

q(n,m-1,i);

}

}

voidmain(){

intn;

while(scanf("%d",&n)!=EOF){

if(n<=0){

printf("Please input a positive interger.\n\n");

continue;

}

k=n;

q(n,n,0);

printf("\n\n");

}

}