监狱长给犯人们出的难题

ggggwhw

1. (*总人数*)
   
2. renshu = 100;
   
3. (*实验次数:*)
   
4. cishu = 10000;
   
5. 
   
6. chenggong = 0;
   
7. For[kk = 1, kk <= cishu, kk++,
   
8. qiutu = {};
   
9. hezi = {};
   
10. haoma = Table[i, {i, 1, renshu}];
    
11. For[ii = 1, ii <= renshu, ii++,
    
12. fanghaoma = RandomInteger[{1, Length[haoma]}];
    
13. AppendTo[hezi, haoma[[fanghaoma]]];
    
14. haoma = Delete[haoma, fanghaoma];
    
15. ];
    
16. For[ii = 1, ii <= renshu, ii++,
    
17. For[jj = 1, jj <= renshu/2, jj++;,
    
18. If[Length[qiutu] < ii,
    
19. AppendTo[qiutu, {hezi[[ii]]}],
    
20. AppendTo[qiutu[[ii]], hezi[[Last[qiutu[[ii]]]]]];
    
21. ];
    
22. If[Last[qiutu[[ii]]] == ii, Break[]];
    
23. ];
    
24. If[jj > renshu/2, Break[]];
    
25. ];
    
26. If[Length[qiutu] == renshu, chenggong += 1(*; Print[qiutu]*);];
    
27. ];
    
28. Print["成功次数" <> ToString[chenggong]];
    
29. Print["实验总次数" <> ToString[cishu]];
    
30. Print["成功的概率" <> ToString[chenggong/cishu // N]];

复制代码

虽然不是很清楚所谓的数字环的理论,但是我觉得这种方法还是不错的.上面是我按照该思想做的Mathematica 程序.当然未必是最短的程序,但是实验次数不是很多时,还是可以出结果的.
成功次数3108
实验总次数10000
成功的概率0.3108
经历的时间是174秒,我还是可以接受这么长时间的.
程序中的(*; Print[qiutu]*)改成Print[qiutu]后可以输出成功时的抽取方法.但是如果实验次数很大时,输出会很长的.如果实验10次成功2~4次那么输出一下可以形象地看出抽取的过程.

[ 本帖最后由 ggggwhw 于 2009-10-21 09:27 编辑 ]