-
Notifications
You must be signed in to change notification settings - Fork 0
/
CountByCars.py
55 lines (49 loc) · 1.72 KB
/
CountByCars.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
# CountByCars
# Daniel Mesham
# 2 February 2019
import numpy as np
import random
import sys
TOTAL_CARS_PASSED = 600
MAX_DIGITS = 6
FRACTION_NON_LOCAL_CARS = 0.2
NON_LOCAL_DIGITS = 3
NUM_TRIALS = 1000
def generate_licence(passed_local_licences = []):
"""Generate a sequence of numbers as they would appear in a licence number"""
random_frac = random.random()
if random_frac <= FRACTION_NON_LOCAL_CARS:
# Use a non-local licence
licence = random.randint(10 ** (NON_LOCAL_DIGITS-1), 10 ** NON_LOCAL_DIGITS)
return str(licence), False
# Use a local licence
licence = random.randint(1, 10 ** MAX_DIGITS)
while licence in passed_local_licences:
licence = random.randint(1, 10 ** MAX_DIGITS)
return str(licence), True
if __name__ == '__main__':
results = []
local_fractions = []
for t in range(0,NUM_TRIALS):
licences_with_length = {}
passed_local_licences = []
counted_licences = []
num_local_licences = 0
count = 1
for l in range(0, TOTAL_CARS_PASSED):
licence, is_local = generate_licence()
licences_with_length[len(licence)] = licences_with_length.get(len(licence), 0) + 1
if is_local:
num_local_licences += 1
passed_local_licences.append(licence)
if str(count) in licence:
count += 1
counted_licences.append(licence)
results.append(count-1)
local_fractions.append(num_local_licences/TOTAL_CARS_PASSED)
if '--list-all' in sys.argv:
print count-1
print np.mean(results)
if '--min-max' in sys.argv:
print "MIN = " + str(min(results))
print "MAX = " + str(max(results))