高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
高校数学の質問スレ（医者・東大卒専用） Part438 http://rio2016.5ch.net/test/read.cgi/math/1723152147/
上
下
前次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
抽出解除
ﾚｽ栞
ﾘﾛｰﾄﾞ規制
です｡10分ほどで解除するので､
他のﾌﾞﾗｳｻﾞ
へ避難してください｡
81: １３２人目の素数さん [sage] 2024/08/18(日) 06:46:57.60 ID:Qg46pTwQ (* 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 http://rio2016.5ch.net/test/read.cgi/math/1723152147/81
83: １３２人目の素数さん [sage] 2024/08/18(日) 13:29:06.79 ID:S1VoGAV5 >>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] http://rio2016.5ch.net/test/read.cgi/math/1723152147/83
ﾒﾓ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.046s