高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞

---

ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡

---

81(1): １３２人目の素数さん [sage] 2024/08/18(日)06:46 ID:Qg46pTwQ(3/3)
(* n 人でジャンケンして勝者がm人になる確率 *)
ja[n_,m_] := (
If[m>=n,Return[0]];
If[m==0,1 - 3*(2^n-2)/3^n,3*Binomial[n,m]/3^n]
)
(* 勝者が1人になるまでのジャンケン回数の期待値 *)
je={0};
AppendTo[je,x /. Solve[x == ja[2,0]x +1 ,x][[1]][[1]]]
AppendTo[je,x /. Solve[x == ja[3,0]x + ja[3,2]je[[2]]+1][[1]][[1]] ]
AppendTo[je,x /.Solve[x == ja[4,0]x + ja[4,3]je[[3]] + ja[4,2]je[[2]] + 1,x][[1]][[1]]]
AppendTo[je,x /.Solve[x == ja[5,0]x + Sum[ja[5,i],{i,2,4}] + 1,x][[1]][[1]] ]
AppendTo[je,x /.Solve[x == ja[6,0]x + Sum[ja[6,i],{i,2,6-1}] + 1,x][[1]][[1]] ]
calc[m_] := AppendTo[je,x /.Solve[x == ja[m,0]x + Sum[ja[m,i],{i,2,m-1}] + 1,x][[1]][[1]] ]
calc[7]
calc[8]
calc[9]
calc[10]
calc[11]
% // N

---

83: １３２人目の素数さん [sage] 2024/08/18(日)13:29 ID:S1VoGAV5(1/3)
>>81
アルゴリズムのバグ発見したのでデバッグ

(* ja[n_,m_] := n 人でジャンケンして勝者がm人になる確率 *)
ja[n_,m_] := (
If[m>=n,Return[0]];
If[m==0,1 - 3*(2^n-2)/3^n,3*Binomial[n,m]/3^n]
)
(* je : 勝者が1人になるまでのジャンケン回数の期待値リスト *)
je={1};
AppendTo[je,x /. Solve[x == ja[2,0](x +1)+ja[2,1] ,x][[1]]]
AppendTo[je,x /. Solve[x == ja[3,0](x+1)+ja[3,1]+ja[3,2](1+je[[2]]),x][[1]]]
AppendTo[je,x /. Solve[x == ja[4,0](x+1)+ja[4,1]+ja[4,2](1+je[[2]])+ja[4,3](1+je[[3]]),x][[1]]]
AppendTo[je,x /. Solve[x == ja[5,0](x+1)+ja[5,1]+ja[5,2](1+je[[2]])+ja[5,3](1+je[[3]])+ja[5,4](1+je[[4]]),x][[1]]]
calc[m_] := AppendTo[je,x /. Solve[x == ja[m,0](x+1)+ja[m,1] + Sum[ja[m,i](1+je[[i]]),{i,2,m-1}],x][[1]]]
calc[6]
calc[7]
calc[8]
calc[9]
calc[10]
calc[11]

---

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

---

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

---

ぬこの手 ぬこTOP 0.048s