Advertisement
Guest User

Untitled

a guest
Jun 23rd, 2017
50
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.40 KB | None | 0 0
  1. import random
  2. import math
  3.  
  4. def prob(n, p):
  5. return math.pow(n, -1/(p))*math.log(10)
  6.  
  7. def cwp(p):
  8. while True:
  9. x = random.random()
  10. n = random.choice([1, 2, 3, 4, 5, 6, 7, 8, 9, 10])
  11. pp = prob(n, p)
  12. if pp > x:
  13. return n%10
  14.  
  15.  
  16. def rn():
  17. ndig = cwp(1)+1
  18. w = []
  19. for x in xrange(ndig):
  20. w.append(str(cwp(x+1)))
  21. return int(''.join(w))
  22.  
  23. for x in xrange(100):
  24. print rn()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement