Задача про Разборчивую невесту на python

Нужно написать програму на python, которая решает стохастический вариант задачи про разборчивую невесту.

Условия задачи:

1. Существует фиксированное и известное число n претендентов на одну позицию, которые можно расположить в порядке качества.

2. Претенденты рассматриваются последовательно в случайном порядке.

3. Для каждого претендента j DM(decision maker) может узнать лишь относительный ранг претендента, то есть, насколько ценен он по сравнению с j - 1 ранее рассмотренными претендентами.

4. Однажды отвергнутого претендента нельзя вернуть. Если достигнут n-й претендент, его необходимо принять.

5'. DM получает вознаграждение pi(a) за выбор претендента с абсолютным рангом a, где pi(1) >=...>=pi(n).

Оптимальная стратегия заключается в том, что DM должна провести интервью с первыми t_1 – 1 претендентами и отбросить всех, затем между претендентом t_1 и претендентом t_2 – 1 она должна принимать только претендентов с относительным рангом 1 (тоесть лучшего среди встретившихся); между претендентом t_2 и претендентом t_3 – 1 она должна принимать претендентов с относительными рангами 1 или 2; и так далее. То есть мы задаем пределы t_1=<t_2=<...=<t_n<n для каждого этапа принятия решения

Особенность стохастического варианта заключается в том, что при построении модели присутствуют дополнительные параметры, такие как beta_r - дисперсия порогового распределения и mu_r - соответствует среднему значению порога для относительного ранга r. Это параметр, характеризующий центр распределения порогов.

В файле в пункте 3.2 описана модель и формулы которые используются для её построения. 

Нужно написать такой код, что бы результаты были максимально похожи на данные из таблицы 2 на стр. 10

Также нужны графики для визуализации похожие на графики с файла