18/11/03 10:08:18.46 Ipwzsvmc.net
>>386
解析解は計算式1行で終わって味気ないのでシミュレーションで算出してみた
> sim <- function(b=2/10,c=3/10,d=4/10,B=1000,C=2000,D=3000){
+ a=1-(b+c+d)
+ br=c(0,B,C,D)
+ total=x=sample(br,1,prob=c(a,b,c,d))
+ while(x!=0){
+ x=sample(br,1,prob=c(a,b,c,d))
+ total=total+x
+ }
+ total
+ }
> re=replicate(1e3,mean(replicate(1e3,sim())))
> summary(re)
Min. 1st Qu. Median Mean 3rd Qu. Max.
17583 19539 20000 20006 20470 22049