Python：每日一题 137

一只小章鱼

趋势如果成正太分布就对了，但我不会怎么模拟出这个场景

冰封雪舞

import random
k=int(input ('请输入需要测试球的总数'))
num = [0]*10
for i in range(k):
    j=0
    for t in range(9):
        j=random.randint(j,j+1)
    num[j] += 1
print(num)

====================== RESTART: D:\Python34\linshi4.py ======================
请输入需要测试球的总数100000
[237, 1774, 7044, 16355, 24581, 24547, 16538, 6947, 1796, 181]
>>>

shigure_takimi

久疤K 发表于 2018-1-13 22:54
显然这不是模拟出来的，这是算概率后计算出来的，而且最终的球数目不对。

1. #  是的，当时不知道怎么模拟，只好计算春概率了。
   
2. #  因为四舍五入，最终球的数目对不上。
   
3. #  看了别人的代码，明白了怎么搞随机掉落。
   
4. #  试着写了下，random.random()<0.5掉在左边，否则掉在右边，结果符合预期了。
   
5. 
   
6. import random
   
7. def drop():
   
8.     a = [[0 for i in range(10)] for j in range(10)]
   
9.     for i in range(10):
   
10.         if i == 0:
    
11.             index = 0
    
12.             a[i][index] = 1
    
13.         elif i > 0:
    
14.             index = a[i-1].index(1) if random.random() < 0.5 else a[i-1].index(1)+1
    
15.             a[i][index] = 1
    
16.     return a[9].index(1)
    
17. 
    
18. result = {}
    
19. for i in range(10):
    
20.     result[i] = 0
    
21. 
    
22. n = 100000
    
23. for i in range(n):
    
24.     result[drop()] += 1
    
25. 
    
26. for k in result:
    
27.     print(k, '-->', result[k])
    
28. 
    
29. ##    0 --> 204
    
30. ##    1 --> 1680
    
31. ##    2 --> 7027
    
32. ##    3 --> 16255
    
33. ##    4 --> 24737
    
34. ##    5 --> 24709
    
35. ##    6 --> 16488
    
36. ##    7 --> 7008
    
37. ##    8 --> 1680
    
38. ##    9 --> 212

复制代码

h20060304

来凑下热闹：

可以实现 任意多行，任意多次试验。

直接上代码：

1. import random
   
2. import time
   
3. 
   
4. 
   
5. def zero_or_one():
   
6.     if random.random() < 0.5:
   
7.         t = 0
   
8.     else:
   
9.         t = 1
   
10.     yield t
    
11. 
    
12. 
    
13. # 最终位置: 往左不变, 往右加 1
    
14. def get_str(num_of_sphere=10, r=10):
    
15.     while num_of_sphere:
    
16.         final_position = 0
    
17.         for i in range(r - 1):
    
18.             direction = next(zero_or_one())
    
19.             if direction:
    
20.                 final_position += 1
    
21.         num_of_sphere -= 1
    
22.         yield final_position
    
23. 
    
24. 
    
25. total = 1000
    
26. row = 10
    
27. print('试验总球数为: {}'.format(total))
    
28. t_list = list(get_str(total, row))
    
29. starttime = time.time()
    
30. for num in range(row):
    
31.     print('落在第 {} 个位置的球个数为: {}'.format(num, t_list.count(num)))
    
32. print('花费的时间：{}'.format(time.time() - starttime))
    

复制代码

Elastcio

1. import random
   
2. result = []
   
3. for i in range(10):
   
4.     result.insert(i,0)
   
5. t=0
   
6. while t<100000:
   
7.     ball = 0
   
8.     for i in range(9):
   
9.         ball += random.randint(0,1)
   
10.     result[ball] += 1
    
11.     t += 1
    
12. print(result)
    

复制代码

塔利班

1. import random as r
   
2. n=100000
   
3. L=[0]*10
   
4. def step():
   
5.     return r.randint(0,1)-0.5
   
6. while n:
   
7.     start=4.5
   
8.     t=9
   
9.     while t:
   
10.         start += step()
    
11.         t -= 1
    
12.     L[int(start)]+=1
    
13.     n -= 1
    
14. print(L)

复制代码

咕咕鸡鸽鸽

算法稍微改了一下，不知道是不是这样算

1. import random
   
2. 
   
3. def fun137():
   
4.     res_list = []
   
5.     for test in range(100000):
   
6.         #先算一个球的趋向
   
7.         str1 = ""
   
8.         for i in range(10):            
   
9.             str1 += random.choice(["+","-"])
   
10.         res = eval("10" + "1".join(str1)[:-1])
    
11.         res_list.append(res//2)
    
12. 
    
13.     count_dict = {}
    
14.     for i in range(10):
    
15.         count = res_list.count(i)
    
16.         count_dict[i] = count
    
17.         
    
18. 
    
19.     return count_dict
    
20. 
    
21. print(fun137())

复制代码