combinate函数对N是否有限制？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

电脑计算combinate(500,2)还可以，但是combinate(1000,2)就一直出不来，理论上后者的行数约是前者的4倍，时间上应该4倍左右也能出来，但是等了10多分钟，还没有出来，求原因！
谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：ATE COM bin MBI 函数

运行很久，能出来，但是不知道为啥要这么长时间。
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

spy1889 发表于 2013-5-14 09:10
运行很久，能出来，但是不知道为啥要这么长时间。
N        组合数        时间s
100        4950        0.028

自己重新写了个组合的程序
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;

spy1889 发表于 2013-5-14 09:10
运行很久，能出来，但是不知道为啥要这么长时间。
N        组合数        时间s
100        4950        0.028

自己重新写了个组合的程序
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

鼓励自己动手修编程序！