Python语言:
Codee#24529
01
# -*- coding: cp936 -*-
02 #monta cola蒙特卡洛方法算PI
03 #younger.x.shen@gmail.com
04 #2011/12/08
05
06 import random
07
08 def getThePi ( randsum ):
09
10 inner = 0.0
11 outer = 0.0
12
13 for i in range ( randsum ) :
14 temp1 = random . random ()
15 temp2 = random . random ()
16
17 if temp1 ** 2 + temp2 ** 2 <= 1 :
18 inner += 1
19
20 perc = inner / randsum
21 print ( inner )
22 circle = 4 * perc
23 # s = PI * R * R
24 # PI = s / (R*R)
25
26 return circle
02 #monta cola蒙特卡洛方法算PI
03 #younger.x.shen@gmail.com
04 #2011/12/08
05
06 import random
07
08 def getThePi ( randsum ):
09
10 inner = 0.0
11 outer = 0.0
12
13 for i in range ( randsum ) :
14 temp1 = random . random ()
15 temp2 = random . random ()
16
17 if temp1 ** 2 + temp2 ** 2 <= 1 :
18 inner += 1
19
20 perc = inner / randsum
21 print ( inner )
22 circle = 4 * perc
23 # s = PI * R * R
24 # PI = s / (R*R)
25
26 return circle
不知道为何高亮这么不好用,插件又装不上去,很郁闷,写多了java ee ,这点代码我也是敲了10多分钟的,算法功力大退,虽然是差不多第一次使用python吧,感觉用python描述算法确实是比c好使,不过还是很慢,之后的时间,我会像云风那样把代码优化到极限,希望可以算得更加精确一点,不过蒙特卡洛方法的先天特性,还是很慢,没有办法。