- 84 名前：１３２人目の素数さん mailto:sage [2024/08/18(日) 13:29:06.79 ID:S1VoGAV5.net]
- >>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]
|
 |
