P(M,N)  從 M 個內取 N 個組合。    您可以使用下列函示 GetN(N:integer);      for i:=1 to N do
  GetN(i);    就可以取得所有子集合。 :)        // 在 M 個元件內取 N 個並加以排列，可重複選取    //SetAry[0..M-1]={1,2,3,4...M}; 取 N 個，可能組合有 M^N 種 = round(Exp(Ln(M)*N))。  CombAry[M^N];    curID:integer;    curID:=0;    procedure GetN(N:integer);
var
  i,j,k,subCombination:byte;
begin
  if (N>1) then
  begin
    subCombination:=1;
    for j:=1 to N-1 do subCombination:=subCombination*M; //其下階 所有 for loop 的組合總數
    
    for i:=0 To M-1 do
    begin        
      for k:=0 to SubCombination-1 do
      begin
        CombAry[curID k]=SetAry[i];
      end;          GetN(N-1); //遞減
    end;
  end
  else
  begin
    for i:=0 To M-1 do
    begin
      CombAry[curID]=SetAry[i];
      inc(curID);
    end;
    exit;
  end;
end;