Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- """
- Car park problem: http://mathworld.wolfram.com/RenyisParkingConstants.html
- theoretical solution: 0.7475979202...
- """
- import random
- uni = random.uniform
- THEORY = 0.7475979202 # for large L
- def park(L):
- """ How many cars can you park on a street with length L,
- if drivers park their car at random spots
- L : length of parking space - length of car=1
- Recursive solution
- """
- if L < 1:
- return 0
- r = uni(0, L - 1)
- return 1 + park(r) + park(L - r - 1)
- def main():
- L = 1000000.0
- cars = park(L)
- sol = cars / L
- print L, sol, "Theory:", THEORY, "Diff %:", (sol - THEORY) / THEORY * 100
- main()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement