sav = 0; iter = 10000000;
For[i = 1, i <= iter, i++,
If[Random[]^2 < 1/2, sav = sav + 1;]
] // Timing
2 sav/iter // N
這樣居然也算了23秒，覺得有點太久，請問要怎樣寫才會快呢?

Nest[UnitStep[(1/2-Random[]^2)]+#&,0,10000000]

RandomVariate[NormalDistribution[1/Sqrt[2], 1/Sqrt[12]//Sqrt[10000000]]]

AbsoluteTiming@Total@NestList[If[Random[]^2<0.5,1,0]&,0,iter]