誕生日のパラドックスを分かりやすく解説（Pythonプログラム付き） - Qiita

k = np.arange(1,1000+1) p = np.ones(len(k))*364/365 p = 1 - np.multiply.accumulate(p) fig,ax = plt.subplots(facecolor='white') ax.plot(k,p) ax.set_yticks([0.25,0.5,0.75,0.9,1]) ax.set_xlim(1,k[-1]) ax.set_ylim(0,1) ax.tick_params(direction='in') ax.set_xlabel('集団の構成人数（人）') ax.set_ylabel('集団のなかに自分と同じ誕生日の人が\n少なくとも1人は存在する確率（%）') for pp in [0.5,0.75,0.9]: idx = np.where(p>=pp)[0][0] ax.scatter(k[idx],