运行很久，能出来，但是不知道为啥要这么长时间。
N        组合数        时间s
100        4950        0.028
200        19900        0.761
300        44850        3.806
400        79800        5.9
500        124750        35.257
600        179700        52.437
700        244650        154.352
800        319600        355.408
900        404550        635.694
1000        499500        1177.313
1100        604450        1455.569
1200        719400        2257.302
1300        844350        3201.859
1400        979300        4494.663
1500        1124250        5625.264

自己重新写了个组合的程序
proc mycom(N);
local x,y,z,z1,z2,rest1,rest2,rest;
x=ones(N,N);
y=seqa(1,1,N);
z=x.*y;
z1=z[1:N-1,2:N];
z2=z[2:N,1:N-1];
rest1=vech(z1');
rest2=vech(z2);
retp(rest1~rest2);
endp;

自己重新写了个组合的程序
proc mycom(N);
local x,y,z,z1,z2,rest1,rest2,rest;
x=ones(N,N);
y=seqa(1,1,N);
z=x.*y;
z1=z[1:N-1,2:N];
z2=z[2:N,1:N-1];
rest1=vech(z1');
rest2=vech(z2);
retp(rest1~rest2);
endp;

结果一下就出来了：

N	组合数	时间s
100	4950	0
200	19900	0
300	44850	0.015
400	79800	0.016
500	124750	0.031
600	179700	0.031
700	244650	0.031
800	319600	0.063
900	404550	0.062
1000	499500	0.078
1100	604450	0.094
1200	719400	0.125
1300	844350	0.124
1400	979300	0.188
1500	1124250	0.202

鼓励自己动手修编程序！