TAOCP中的两个组合数算法

最近志哥让看《计算机程序设计艺术》中的组合数部分,目的是需要生成组合数。目前值看了一点点。虽然打了这么久的ACM但是不得不感叹。妙啊~~~秒~~。

这里先贴出使用c++写的L算法和T算法。

L算法和T算法:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
#include <iostream>
using namespace std;
const int M = 10000;
int c[M];
void visit(int t){
for(int i=1;i<=t;i++){
cout<<c[i];
}
cout<<endl;
}
void l(int n,int t){
int i = 0;
for(i=1;i<=t;i++){
c[i] = i-1;
}
c[t+1] = n;
c[t+2] = 0;
while(1){
visit(t);
i = 1;
while(c[i]+1 == c[i+1]){
c[i] = i-1;
i++;
}
if(i > t)
break;
c[i] = c[i]+1;
}
}
void t(int n,int t){
int i = 0;
for(i=1;i<=t;i++){
c[i] = i-1;
}
c[t+1] = n;
c[t+2] = 0;
int j = t;
int x = 0;
while(1){
T2:
visit(t);
if(j > 0){
x = j;
} else {
T3:
if(c[1]+1 < c[2]){
c[1] = c[1]+1;
continue;
} else {
j = 2;
while(1){
T4:
c[j-1] = j-2;
x = c[j] +1;
if(c[j+1] == x)
j = j+1;
else
break;
}
if(j > t){
break;
}
}
}
T6:
c[j] = x;
j = j-1;
}
}
int main(){
l(10,3);
t(10,3);
}