高校数学の質問スレ（医者・東大卒専用） Part438 (899ﾚｽ)
上下前次1-新
抽出解除 必死ﾁｪｯｶｰ(本家) (べ) 自ID ﾚｽ栞 あぼーん
ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡
108(2): 2024/08/21(水)08:28 ID:trPnwZW4(1/6) AAS
>>106
で、３Dグラフの動画はまだできんのだが？
110: 2024/08/21(水)10:42 ID:trPnwZW4(2/6) AAS
>>107
レスありがとうございます。
コードを解読しました。
(* a 人でジャンケンして参加者が b 人になる確率,a==b ならアイコ *)
p[a_,b_]:=p[a,b]=If[a==b,1-(2^a-2)/3^(a-1),Binomial[a,b]/3^(a-1)]
m=11;
(* 参加者数の変遷を列挙 *)
省7
111: 2024/08/21(水)11:14 ID:trPnwZW4(3/6) AAS
デバッグして、数値が合致しました。
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= (* jan[n_,m_] := n 人でジャンケンして参加者が m 人になる確率 *)
In[2]:= 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]
省8
112: 2024/08/21(水)11:27 ID:trPnwZW4(4/6) AAS
(* 各回ごとの確率 *)
pn={jan[n,1]};
Total@Flatten@Table[jan[n,x1]jan[x1,1],{x1,2,n}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,1],{x1,2,n},{x2,2,x1}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,1],{x1,2,n},{x2,2,x1},{x3,2,x2}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4}]//AppendTo[pn,#]&;
省5
113(1): 2024/08/21(水)11:44 ID:trPnwZW4(5/6) AAS
シミュレーションによる検証
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= j[n_] :=( (* n人でジャンケンして勝者が決まるまでの回数と勝者の数*)
count=0;
Until[Length@Union@a==2,a=RandomChoice[Range[3],n];count++];
b=Sort@Union@a;
省25
117(4): 2024/08/21(水)17:57 ID:trPnwZW4(6/6) AAS
>>115
Flattenの使い方とか熟練者のコード解析は勉強になる。
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.029s