楼主:
Serenede (Serenede)
2014-12-21 00:27:29sav = 0; iter = 10000000;
For[i = 1, i <= iter, i++,
If[Random[]^2 < 1/2, sav = sav + 1;]
] // Timing
2 sav/iter // N
这样居然也算了23秒,觉得有点太久,请问要怎样写才会快呢?
AbsoluteTiming@Total@NestList[If[Random[]^2<0.5,1,0]&,0,iter]