Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import random
- import math
- def prob(n, p):
- return math.pow(n, -1/(p))*math.log(10)
- def cwp(p):
- while True:
- x = random.random()
- n = random.choice([1, 2, 3, 4, 5, 6, 7, 8, 9, 10])
- pp = prob(n, p)
- if pp > x:
- return n%10
- def rn():
- ndig = cwp(1)+1
- w = []
- for x in xrange(ndig):
- w.append(str(cwp(x+1)))
- return int(''.join(w))
- for x in xrange(100):
- print rn()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement