[過去ﾛｸﾞ] 分からない問題はここに書いてね458 (1002ﾚｽ)
上下前次1-新
抽出解除 必死ﾁｪｯｶｰ(本家) (べ) 自ID ﾚｽ栞 あぼーん

---

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ

---

493: １３２人目の素数さん [sage] 2020/03/06(金) 23:04:18.50 ID:EFqGY3yx(1/2) AAS
>>484
> N=6^10
> A0=5^10
> B0=5^10
> B1=10*5^9
> A0B0=4^10
> A0B1=10*4^9
> 1-(A0+B0+B1-(A0B0+A0B1))/N
[1] 0.41467302314603127

---

494: １３２人目の素数さん [sage] 2020/03/06(金) 23:42:01.10 ID:EFqGY3yx(2/2) AAS
"サイコロを10回振るとき、1が1回以上かつ2が2回以上でる確率"
rm(list=ls())
options(digits=22)

# ¬(1が1回以上　∧　2が2回以上でる) == 1が０回　∨　(２が0回 ∨ 2が１回）
N=6^10 # すべての順列
A0=5^10 # 1が0回の場合
B0=5^10 # 2が０回の場合
B1=10*5^9 # 2が１回の場合

A0B0=4^10 # 1が0回の場合∧2が０回の場合
A0B1=10*4^9 # 1が0回の場合∧2が1回の場合

1-(A0+B0+B1-(A0B0+A0B1))/N

N - (A0+B0+B1-(A0B0+A0B1))
N
25073692/60466176
6268423/15116544

# シラミ潰しに数え上げる
library(gtools)
pm=permutations(6,10,rep=T)
f <- function(x){
sum(x==1)>0 & sum(x==2)>1
}
r=sum(apply(pm,1,f))　# 25073692

sim <- function(){
f(sample(6,10,rep=TRUE))
}
mean(replicate(1e6,sim()))
１００万回のシミュレーションでの割合
> mean(replicate(1e6,sim()))
[1] 0.41464000000000001

---

上下前次1-新書関写板覧索設栞歴

---

ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ

---

ぬこの手 ぬこTOP 0.042s