一个组合数的求解方法
来源:文书网 2.6W
有从10个不同的数中取出7个,如何求出所有组合?初看这个问题,确实不太好下手,虽然我们理解这个问题很容易,但要让计算机“理解”可得花点功夫了。
先分析:首先想到,如果由10个元素中的7个组成一个序列,并且这7个元素互不相等,这就比较接近于正解了;然后考虑,组合中7元素是不分先后的,我们如何剔除多余的数据呢(如四选二中(1,3)与(3,1)是相同解)?
我们可以在编程时作一种限定,7个元素的排列顺序满足我们的一个规定,这样,就可以依据相同位置的值不能相同来排除不正确的解。
这样我们的第一个解呼之欲出了:可以利用数组来进行选取,不同的数组下标顺序就代表了不同的解,而且我们约定这个下标序列必须是升序,这就利于我们排除冗余值。
在本文中,约定这10个数是如下形式的数组,并且已经赋值。计算的结果在一个TMemo控件中显示。
var
Value:array[1..10]ofinteger;
解法一、
按钮1的`点击事件处理。
on1Click(Sender:TObject);
var
idx1,idx2,idx3,idx4,idx5,idx6,idx7:integer;
tmpStr:string;
begin
r;
foridx1:=1to4do
foridx2:=idx1+1to5do
foridx3:=idx2+1to6do
foridx4:=idx3+1to7do
foridx5:=idx4+1to8do
foridx6:=idx5+1to9do
foridx7:=idx6+1to10do
begin
tmpStr:=IntToStr(idx1)+''+IntToStr(idx2)+''+IntToStr(idx3)
+''+IntToStr(idx4)+''+IntToStr(idx5)+''+IntToStr(idx6)
+''+IntToStr(idx7);
(tmpStr);
end;
end;
解中只显示了下标的组合,实际应用把下标改为Value数组即可:tmpStr:=IntToStr(Value[idx1])+''+IntToStr(Value[idx2])+...;。
这个解的一个亮点就是每一个循环变量的初值都是它前一个变量加1;这就保证后一个下标一定不等于前一个下标,请体会一下为什么循环控制变量的终值为4~10。
解法二、
中学的数学教材中就讲到C(10,7)=C(10,3),这个很好理解,从10中选7,剩下3个,所有选七的组合完成,也就是所有选三的组合完成,反之亦然。
按钮2的点击事件处理:
on2Click(Sender:TObject);
var
idx1,idx2,idx3,idx4:integer;
tmpStr:string;
begin
r;
foridx1:=1to8do
foridx2:=idx1+1to9do
foridx3:=idx2+1to10do
begin
tmpStr:='';
foridx4:=1to10do
if(idx4<>idx1)and(idx4<>idx2)and(idx4<>idx3)//加入一个元素
thentmpStr:=tmpStr+IntToStr(Value[idx4])+'';
(tmpStr);
end;
end;
这个解是十选三,然后显示剩下的七个,是不是有点意思呢?
以上加入元素的判断可以改写为:
if((idx1+100)*(idx2+100)*(idx3+100)mod(idx4+100)>0)
请体会一下同余理论的运用,这里利用了101到110这十个数中任一个数都不被其它三个数的乘积整除的特性
解法三、
以上两解也可用集合来实现。
再想想,以上两解都是用了多重循环,必须定义多个循环变量,通用性似乎差了点。
观察这些循环(尤其是解一),形式是何等的一致!
也许诸位已经想到我要用递归调用来解决这个问题了。
要设计递归,首先要确定哪些变量是必须向这个函数传递的,容易知道,调用时要知道当前元素的起始下标与终止下标,当然还有两个必不可少的参数就是我这个函数计算的是从多少选多少的,这样就确定了四个参数。
还有一个要注意的地方,是递归何时显示数据?显然应当在求最后一个元素时。
(关于递归调用的一些详细的论述,可以参看乔林的《参透Delphi/Kylix》一书或其它教材)
请看实现:
全局变量
var
GIDX:array[1..20]ofinteger;//记录下标的数组。
procedureMyRecur(Total,SelCnt,BLoc,ELoc:integer);
var
idx,jjj:integer;
tmpStr:string;
begin
if(ELoc=Total)//终止状态
thenbegin
foridx:=BLoctoELocdo
begin
GIDX[SelCnt+ELoc-Total]:=idx;//注意此处应与下面8888处一致
tmpStr:='';
forjjj:=1toSelCntdotmpStr:=tmpStr+IntToStr(Value[GIDX[jjj]])+'';
(tmpStr);
end;
end
elsebegin
foridx:=BLoctoELocdo
begin
GIDX[SelCnt+ELoc-Total]:=idx;//8888
MyRecur(Total,SelCnt,idx+1,ELoc+1);
end;
end;
end;
procedureSelNumber(Total,SelCnt:integer);