- 154 名前：１３２人目の素数さん mailto:sage [2024/08/25(日) 12:59:59.55 ID:/qrXHaIo.net]
- >>140
行列で計算させると算出時間が爆速なのにびっくりしました。
達人のスクリプトを改造して道具箱に保存しておきます。
(* j 人でジャンケンをしたときの終了までの回数の最頻値とその確率を返す *)
calc[j_]:=(
M=Table[If[n==m,1-(2^m-2)/3^(m-1),Binomial[m,n]/3^(m-1)],{m,1,j},{n,1,j}];
p=Differences@Table[MatrixPower[M,i][[j,1]],{i,0,10j}];
max=Max@p;
Flatten@{Position[p,max],max,N[max]}
)
In[3]:= calc[15]
Out[3]= {30, 65101358743766874914341259145354001254712997240185483777481087387163094008455949207206\
> 93198405902336999941273392239247993407703628004425130823658476122231689832344483758966236574\
> 375695 / 11947838420050013668726696739307151046843799152024135169583095938840977078626722578\
> 97327618239887790786549346048626664496721871548575328400043101228717425477619608889629973635\
> 327326175449, 0.0054488}
calc[15]
|
 |
