- 80 名前：１３２人目の素数さん mailto:sage [2024/08/18(日) 06:09:54.38 ID:Qg46pTwQ.net]
- (* 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={1};
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]]]
