%%time

%matplotlib inline
import math
import random
import matplotlib.pyplot as plt
import numpy as np
import bisect

stats = np.loadtxt('us_births_69_88.csv', skiprows = 1, delimiter=',', dtype = int)

maxi = max(stats[:,2])
print(maxi)

def funkcja(v, tab):
    tab = [0 for i in range(0, v + 1)]
    k = 0
    while True:
        p = random.randint(0, v - 1)
        q = random.randint(0, maxi)
        if stats[p][2] >= q:
            #return p
            i = p;
            tab[i] = tab[i] + 1
            if tab[i] > 1:
                return k
            k = k + 1
    return -1

v = len(stats)
tab = [0 for i in range(0, v + 1)]
tab2 = [funkcja(v, tab) for i in range(100000)]

plt.hist(tab2, bins=30)
plt.show()

# ilość kroków aż będzie dana data
# samplować ale od dni a nie kiedy pierwsze wystąpienie