[過去ﾛｸﾞ] 分からない問題はここに書いてね458 (1002ﾚｽ)
上下前次1-新
抽出解除 必死ﾁｪｯｶｰ(本家) (べ) 自ID ﾚｽ栞 あぼーん
このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ
638(1): 2020/03/15(日)04:25 ID:OTl1KJku(1/7) AAS
>>635
ｎ を 1〜２０でシミュレーションしてみた。
sim <- function(n,k=1e5){ # n:サイコロを振る回数 k:シミュレーション回数
sub <- function(n){
prod(sample(6,n,replace=TRUE))%%6==0 # n回の目の積の6で除算した剰余が0か？
}
mean(replicate(k,sub(n))) # 0となる割合を返す
省24
640: 2020/03/15(日)08:00 ID:OTl1KJku(2/7) AAS
>>638
1 : 1/6
2 : 5/12
3 : 133/216
4 : 325/432
5 : 6541/7776
6 : 4655/5184
省3
642(1): 2020/03/15(日)09:46 ID:OTl1KJku(3/7) AAS
>>635
1-(4^n+3^n-2^n)/6^n
643: 2020/03/15(日)09:52 ID:OTl1KJku(4/7) AAS
ｆ ＜− ｆｕｎｃｔｉｏｎ（ｎ）｛
ｌｉｂｒａｒｙ（ｇｍｐ）
ｎ＝ａｓ．ｂｉｇｑ（ｎ）
ｒ＝１−（４＾ｎ＋３＾ｎ−２＾ｎ）／６＾ｎ
ｒ２＝ｃａｐｔｕｒｅ．ｏｕｔｐｕｔ（ｒ）［２］
ｓｕｂｓｔｒ（ｒ２，５，ｎｃｈａｒ（ｒ２））
｝
省34
644: 2020/03/15(日)10:04 ID:OTl1KJku(5/7) AAS
分数表示 小数表示 シミュ値
1 1/6 0.1666666667 0.16695
2 5/12 0.4166666667 0.41535
3 133/216 0.6157407407 0.61633
4 325/432 0.7523148148 0.75147
5 6541/7776 0.8411779835 0.84065
6 4655/5184 0.8979552469 0.89742
省15
668: 2020/03/15(日)22:53 ID:OTl1KJku(6/7) AAS
>>661
いや、等しくなったけど。
> sim <- function(m){
+ n=2*m+1
+ i=1:(n-1)
+ a=sqrt(i*n)
+ b=floor(a)
省9
670: 2020/03/15(日)23:21 ID:OTl1KJku(7/7) AAS
10万＊２＋１までは成立することを確認。
> sim <- function(m,print=FALSE){
+ n=2*m+1
+ i=1:(n-1)
+ a=sqrt(i*n)
+ b=floor(a)
+ if(print) cat(sum(b%%2==0),sum(b%%2==1),'\n') # 偶数の数、奇数の数
省12
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.030s