배경.

지난주에 델레스의 한정 3선이 가사에서 20만 명 정도를 폭사시키고 철수했다
스피커 대신 4장 정도 받았어요.
얻기 어려운 기회니까 확률이 얼마인지 시각화해 보려고요.
확실히 공짜 돌까지 합쳐서 750연발을 쳤어요.
며칠씩 나눠서 서서히 폭사했기 때문에 정확한 횟수를 기억하지 못한다.

컨디션

jupyter notebook on anaconda3-4.2.0먼저 다음을 수행하십시오.

import scipy.stats as st
import matplotlib.pyplot as plt

%matplotlib inline

우선 SSR 배출 확률

가봐.
SSR의 배출 확률은 1.5%입니다. 형상 분포의 누적 확률 분포를 살펴보겠습니다.

v = st.geom.cdf(range(0 , 500), 0.015)
plt.ylim(0, 1)
plt.grid()
plt.plot(v)

100연발은 20% 정도 할인해 드릴게요.
기대치는 1 / 0.015 = 66.6회다.
작년에 대포유에서 폭사한 P.
체감 5만 엔이면 1SSR 느낌.
감상을 서술하다.
9800엔=8400보석이기 때문에 50000엔으로 끝자리를 버리면 42000보석이다.
모두 10중대1를 넣으면 이쪽도 끝자리 16회(160중대)를 뺀다.

50000 // 9800 * 8400
# -> 42000
10 * (42000 // 2500)
# -> 160

그래프를 참조하면 160연발은 90% 정도인데 잘 어울릴 것 같아요.

SSR 배출 확률 제한

상응하는 와셔 한정 3선과 한정 스피커 2장을 정선했다
왜냐하면 배출 확률은 각각 0.3%일 거예요.
일종의 누적 확률 분포는 유한하다

v = st.geom.cdf(range(0 , 2000) , 0.003)
plt.ylim(0, 1)
plt.grid()
plt.plot(v)

모양이 같지만 가로축의 범위가 넓어지고 있으니 주의하세요.
500연발은 80퍼센트 정도 될 거예요.
기대치는 1 / 0.003 = 333.3회다.
직감적으로'500번 치면 환불되죠?'이런 느낌도 있지만.
실제로 5명 중 1명은'500번 돌려도 목표를 한정할 수 없다'.

최후

그럼 처음에 쓴 부분은 왜 스피커 4장을 줄였을까.
마이너스 2항식 분포2의 누적 확률 분포를 보면

v = st.nbinom.cdf(range(5000), 4, 0.003)
plt.ylim(0, 1)
plt.grid()
plt.plot(v)

에이
1800연발 정도면 8할의 확률이 있겠죠.
750연에서 4장을 뺀 것은 상당히 강한 운이라고 할 수 있는데, 확률은 10% 정도다.
그런데 그쪽이 아닌데..

덤

트위터에 교본을 붙이지 않은 게 생각나서 일단 여기에 붙여놨어요.
와셔로 찍는 습관이 없어 사후 일람화면이다.
아직 자습서가 끝나지 않았지만 앞으로도 잘 부탁드립니다.

실질 무료 하지만.  ↩
마이너스 2항식 분포에서 k=1은 기하학적 분포참고 자료↩이다.