20/03/22 23:00:24.46 liILqu/N.net
>>871
n=10のときで総当たりで計算してみた
> rm(list=ls())
> n=10
> dec2nw <- function(num, N, digit = n){
+ r=num%%N
+ q=num%/%N
+ while(q > 0 | digit > 1){
+ r=append(q%%N,r)
+ q=q%/%N
+ digit=digit-1
+ }
+ return(r)
+ }
> d=t(sapply(0:(2^n-1),function(num) dec2nw(num,2,n))) # 10枚のコインの裏表の順列（０を取り除く表とする）
> head(d) ; tail(d)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 0 0 0 0 0 0 0 0 0 0
[2,] 0 0 0 0 0 0 0 0 0 1
[3,] 0 0 0 0 0 0 0 0 1 0
[4,] 0 0 0 0 0 0 0 0 1 1
[5,] 0 0 0 0 0 0 0 1 0 0
[6,] 0 0 0 0 0 0 0 1 0 1
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1019,] 1 1 1 1 1 1 1 0 1 0
[1020,] 1 1 1 1 1 1 1 0 1 1
[1021,] 1 1 1 1 1 1 1 1 0 0
[1022,] 1 1 1 1 1 1 1 1 0 1
[1023,] 1 1 1 1 1 1 1 1 1 0
[1024,] 1 1 1 1 1 1 1 1 1 1
> d1=apply(d,1,sum) # 各行ごとの合計＝裏がでた枚数
> sum((1/2)^d1)/2^n # 各行ごとに (1/2) ^ 裏の枚数を計算して合算する
[1] 0.056314
> (1-1/4)^n
[1] 0.056314
合致した。