Skip to Content

§4. Разорение игрока: блуждание между двумя стенками

Постановка

Два противника играют много независимых партий. Первый выигрывает партию с вероятностью pp, второй — с вероятностью q=1−pq = 1 - p (ничьих нет). Проигравший платит 1. Начальный капитал первого игрока kk, второго — M−kM - k. Игра кончается, когда один из них проиграет всё.

Ход игры — это траектория случайного блуждания Sn=Z1+⋯+ZnS_n = Z_1 + \dots + Z_n, поднятая на kk единиц вверх, где шаги ZiZ_i независимы, P(Zi=1)=p\mathbf P(Z_i = 1) = p, P(Zi=−1)=q\mathbf P(Z_i = -1) = q. Величина Sn+kS_n + k — капитал первого игрока после nn партий. Суммарный капитал MM не меняется, так что блуждание идёт по полосе между стенками 0 и MM. Найдём rkr_k — вероятность разорения первого игрока.

Смысл по-пацански: это та же картинка, что у Вальда в §3. Блуждание, две границы, остановка при первом касании. Только вместо «приняли гипотезу» — «пацан ушёл домой без денег».

Игра не длится вечно

Перенесём начало координат в точку (0,k)(0, k) и обозначим через ν\nu момент окончания игры. Так как MZ1=p−q\mathbf MZ_1 = p - q и DZ1=1−(p−q)2<∞\mathbf DZ_1 = 1 - (p - q)^2 < \infty, при p≠1/2p \ne 1/2 выполнены условия теоремы 3 из §3, и Mν<∞\mathbf M\nu < \infty.

Докажем, что и при p=1/2p = 1/2 (тогда MZ1=0\mathbf MZ_1 = 0) Mν\mathbf M\nu конечно. Возьмём чётное M=2mM = 2m и l=min⁡{k,M−k}l = \min\{k, M - k\}. С вероятностью 2−l≥2−m2^{-l} \ge 2^{-m} игра может закончиться за ll партий — достаточно ll выигрышей (или проигрышей) подряд. Разобьём блуждание на независимые отрезки длины mm. На каждом из них игра кончается с вероятностью не меньше 2−m2^{-m}, поэтому

P(ν>m)≤1−2−m,…,P(ν>jm)≤(1−2−m)j.\mathbf P(\nu > m) \le 1 - 2^{-m}, \quad \dots, \quad \mathbf P(\nu > jm) \le (1 - 2^{-m})^j.

Вероятности P(ν>i)\mathbf P(\nu > i) убывают с ростом ii, а при i=jmi = jm мажорируются членами геометрической прогрессии со знаменателем γ=1−2−m\gamma = 1 - 2^{-m}. Поэтому ряд ∑i=0∞P(ν>i)\sum_{i=0}^\infty\mathbf P(\nu > i) сходится, а по формуле (3) главы 1 его сумма равна Mν\mathbf M\nu. (Случай нечётного MM разберите самостоятельно.)

Ответ через тождество Вальда

Конечность Mν\mathbf M\nu позволяет применить тождество Вальда (теорема 5 из §3). При p=1/2p = 1/2 возьмём c0=−kc_0 = -k, c1=M−kc_1 = M - k, rk=P(Sν=c0)r_k = \mathbf P(S_\nu = c_0):

MSν=c0rk+c1(1−rk)=MZ1⋅Mν=0⋅Mν=0.\mathbf MS_\nu = c_0r_k + c_1(1 - r_k) = \mathbf MZ_1\cdot\mathbf M\nu = 0 \cdot \mathbf M\nu = 0.

Отсюда −krk+(M−k)(1−rk)=0-kr_k + (M - k)(1 - r_k) = 0 и

rk=1−kM.r_k = 1 - \frac kM.

Шансы на выигрыш прямо пропорциональны начальному капиталу. Перегнать в честной игре того, у кого денег в десять раз больше, можно лишь с вероятностью 1/111/11.

Ответ через первый шаг

Тот же ответ получим другим способом — разложением по первому шагу. По формуле полной вероятности после первой партии капитал станет k+1k + 1 (вероятность pp) или k−1k - 1 (вероятность qq), и дальше всё начинается заново:

rk=p rk+1+q rk−1,k=1,…,M−1,(18)r_k = p\,r_{k+1} + q\,r_{k-1}, \qquad k = 1, \dots, M - 1, \tag{18}

с граничными условиями r0=1r_0 = 1, rM=0r_M = 0. При p=1/2p = 1/2 соотношение (18) переписывается как rk+1−rk=rk−rk−1r_{k+1} - r_k = r_k - r_{k-1}: точки (k,rk)(k, r_k) лежат на одной прямой rk=A+Bkr_k = A + Bk. Из граничных условий A=1A = 1, B=−1/MB = -1/M.

Смысл по-пацански: «вероятность продуть отсюда» — среднее от «вероятностей продуть» в двух соседних клетках, куда ты попадёшь через партию. Если монетка честная, это просто середина отрезка — и вся картинка выпрямляется в линейку.

Аналогия с краевой задачей

Соотношение (18) похоже на краевую задачу для дифференциального уравнения второго порядка

y′′(x)=f(x,y),y(x0)=a,y(x1)=b.(19)y''(x) = f(x, y), \qquad y(x_0) = a, \quad y(x_1) = b. \tag{19}

При малом h>0h > 0 производную y′(x)y'(x) приближённо заменяют конечной разностью [y(x+h)−y(x)]/h[y(x + h) - y(x)]/h, а y′′(x)y''(x) — второй симметричной разностью

y(x+h)−y(x)h−y(x)−y(x−h)hh=y(x+h)−2y(x)+y(x−h)h2.(20)\frac{\frac{y(x+h) - y(x)}h - \frac{y(x) - y(x-h)}h}h = \frac{y(x + h) - 2y(x) + y(x - h)}{h^2}. \tag{20}

При p=1/2p = 1/2 соотношение (18) принимает вид rk+1−2rk+rk−1=0r_{k+1} - 2r_k + r_{k-1} = 0 — это правая часть (20) при x=kx = k, y(x)=rky(x) = r_k, h=1h = 1. А общее решение уравнения y′′=0y'' = 0 — та же прямая y=A+Bxy = A + Bx.

Метод стрельбы

Чтобы найти rkr_k при p≠1/2p \ne 1/2 (задача 6), воспользуемся идеей метода стрельбы для численного решения краевых задач. Допустим, мы умеем решать задачу Коши

y′′(x)=f(x,y),y(x0)=a,y′(x0)=c,(21)y''(x) = f(x, y), \qquad y(x_0) = a, \quad y'(x_0) = c, \tag{21}

например, методом Эйлера из §6 главы 2. Фиксируем aa и меняем cc — начальный наклон, «угол прицела». Пусть g(c)g(c) — значение решения (21) в точке x1x_1. Тогда краевая задача (19) сводится к поиску корня c∗c^* уравнения g(c)=bg(c) = b: бросаем мяч в кольцо, смотрим, перебросили или недобросили, и подправляем силу броска — делением отрезка пополам или методом Ньютона (§5 главы 9).

Для дискретной задачи (18) аналог y′y' — приращение Δk=rk+1−rk\Delta_k = r_{k+1} - r_k. Стреляем из r0=1r_0 = 1 с начальным наклоном c=Δ0c = \Delta_0 и подбираем cc так, чтобы попасть в rM=0r_M = 0:

Загрузка редактора…

Ответ задачи 6:

rk=λM−λkλM−1,λ=q/p.r_k = \frac{\lambda^M - \lambda^k}{\lambda^M - 1}, \qquad \lambda = q/p.

Устремим M→∞M \to \infty — второй игрок «бесконечно богат». Предел rkr_k равен 1, если p≤1/2p \le 1/2, и λk\lambda^k, если p>1/2p > 1/2. Искусный игрок может с положительной вероятностью 1−λk1 - \lambda^k никогда не проиграть даже бесконечно богатому противнику: при p=0,55p = 0{,}55 и k=10k = 10 это 1−(9/11)10≈0,871 - (9/11)^{10} \approx 0{,}87. Если он не проиграл сразу из-за случайности, то дальше его шансы резко растут благодаря «сносу» траектории вверх.

Смысл по-пацански: а вот обратная сторона. В рулетке ставка на красное выигрывает с вероятностью p=18/37≈0,486p = 18/37 \approx 0{,}486 — казино почти честное. Но если прийти с 50 фишками и играть до 100, разоришься с вероятностью 0,94, а не 0,5. Маленький перекос на длинной дистанции превращается в приговор.

Сколько длится игра

Найдём среднюю длительность игры lk=Mνl_k = \mathbf M\nu, k=1,…,M−1k = 1, \dots, M - 1. Те же рассуждения, что для (18), приводят к уравнению

lk=p lk+1+q lk−1+1,k=1,…,M−1,(22)l_k = p\,l_{k+1} + q\,l_{k-1} + 1, \qquad k = 1, \dots, M - 1, \tag{22}

с граничными условиями l0=0l_0 = 0, lM=0l_M = 0. Единица в правой части добавлена потому, что одна партия уже сыграна. Нетрудно проверить, что при p=1/2p = 1/2 решение lk=k(M−k)l_k = k(M - k). При больших MM и k≈aMk \approx aM среднее время игры равных по силе противников имеет порядок M2M^2.

Это становится менее удивительным, если вспомнить центральную предельную теорему: траектория симметричного блуждания SnS_n колеблется в среднем в пределах ±n\pm\sqrt n. Чтобы уйти на расстояние M/2M/2, нужно порядка M2/4M^2/4 шагов.

При p≠1/2p \ne 1/2 дело обстоит иначе: срабатывает снос (замечание 1 из §3), и игра идёт за время порядка MM (вопрос 4).

Разорение игрока

12 игр: капитал первого игрока после каждой партии. Красные кончились разорением (капитал 0), синие — выигрышем всех M = 20.Вероятность разорения r_k = 0.5000 (по 3000 играм — 0.5090). Средняя длительность игры l_k = 100.0 (по моделированию — 99.6).Нижний график — точная функция распределения длительности P(ν ≤ m). Медиана — 76 партий.

При k=10k = 10, M=20M = 20, p=1/2p = 1/2 средняя длительность l10=100l_{10} = 100. Но медиана — около 75 партий, а примерно в каждом пятом случае игра длится больше 150 партий: хвост у закона длинный. Сдвинь pp до 0,45 — разорение станет почти неизбежным (0,88), а игра в среднем короче (76 партий).

Проверь себя

Тест

Честная игра, p = 1/2. У тебя 10 рублей, у соперника 90. С какой вероятностью ты его разоришь?
Обновлено