import random
import time
import sys
m=0
n=0
p=int(sys.argv[1])
while True:
    r=random.randint(1,100)
    if r <= p:
        i = 1
    else:
        i = 0

    a=random.randint(1,100)
    if a <= 80:
        x = i
    else:
        x = 1 - i

    b=random.randint(1,100)
    if b <= 80:
        y = i
    else:
        y = 1 - i

    if x == y and x == 1:
        if i == 1:
            n = n + 1
        m = m + 1
    if m > 0:
        print("hit " , n , " in " , m , " ", n/m)
    time.sleep(0.001)
【 在 dawei78 的大作中提到: 】
: 程序模拟，xyp/(xyp+(1-p)(1-x)(1-y))居然是对的！
:
:
--
FROM 223.71.29.*