rm(list=ls())

g=rep(0:3,c(2700,100,100,100))
"
確率の問題で分からないことがあるので計算方法を教えて下さい

あるガチャガチャがあります。
中身は
　アタリA　3.3%
　アタリB　3.3%
　アタリC　3.3%
　ハズレ　90.0%
となっておりアタリやハズレを引いても無限に補充され続けます（常にこの確率です）。
また10回毎にアタリ確定がありABCのどれかが1/3で出てきます。

このガチャガチャでアタリ3種を5個づつ手に入れるには何回回せばいいでしょうか？

＞アタリやハズレを引いても無限に補充され続けます（常にこの確率です）。
という設定が現実的ではないので
開始時に
アタリA　100
アタリB　100
アタリC　100
ハズレ　 2700
があり、補充はされないとする。
"

f=\(g){
i=sample(length(g),10)
list(g[i],g[-i])
}
sim=\(n=5){
count=10
re=f(g)
p=re[[1]]
g=re[[2]]
while((sum(p==1)<n | sum(p==2)<n | sum(p==3)<n)){
re=f(g)
p=c(p,re[[1]])
g=re[[2]]
count=count+10
}
count
}

count=replicate(10^6,sim())
BEST::plotPost(count,breaks='scott')