Top to random simul.py

From Werner KRAUTH

import random
NCards = 3
data = {}
for iter1 in range(1000000):
    L= [k for k in range(NCards)]
    for iter2 in range(10):
        a = L.pop(0)
        L.insert(random.randint(0, len(L)), a)
    L = tuple(L)
    data[L] = data.get(L, 0) + 1
for k in data:
    print(k, data[k])