高校数学の質問スレ（医者・東大卒専用） Part438

(* jan[n_,m_] := n 人でジャンケンして参加者が m 人になる確率 *)
jan[n_,m_] := (
If[m==0 || m>n, Return[0]];
If[m==n,1 - 3*(2^n-2)/3^n,3*Binomial[n,m]/3^n]
)
(* je : 勝者が1人になるまでのジャンケン回数の期待値リスト *)
solve[n_]:=(
je={0};
calc[m_] := AppendTo[je,x /. Solve[x == jan[m,m](x+1)+jan[m,1] + Sum[jan[m,i](1+je[[i]]),{i,2,m-1}],x][[1]]];
Table[calc[m],{m,2,n}][[-1]];
je
)
solve[10][[-1]]
% // N