負けた場合に3連敗以上している確率を計算していたことに気づいた。

> p=1/6
> n=407
> sim=\(){
+ rbinom(n,1,p) |> rle() -> x
+ cbind(x$values,x$lengths) -> d
+ any((d[d[,1]==1,])[,2]>=3)
+ }
> k=1e6
> replicate(k,sim()) |> mean()
[1] 0.79583

シミュレーションでも8割弱の結果が得られた。