高校数学の質問スレ(医者・東大卒専用) Part438 (979レス)
高校数学の質問スレ(医者・東大卒専用) Part438 http://rio2016.5ch.net/test/read.cgi/math/1723152147/
上
下
前次
1-
新
通常表示
512バイト分割
レス栞
抽出解除
レス栞
リロード規制
です。10分ほどで解除するので、
他のブラウザ
へ避難してください。
140: 132人目の素数さん [sage] 2024/08/23(金) 19:36:15.43 ID:R15w90it >>134 M=Table[If[n==m,1-(2^m-2)/3^(m-1),Binomial[m,n]/3^(m-1)],{m,1,11},{n,1,11}]; u=Table[Boole[i==11],{i,1,11}]; v=Table[Boole[i== 1],{i,1,11}]; p=Differences@Table[u.MatrixPower[M,i].v,{i,0,20}]; max=Max[p] 12597340638043496499971512842764330462294598371/515377520732011331036461129765621272702107522001 Position[p,max] 11 http://rio2016.5ch.net/test/read.cgi/math/1723152147/140
145: 132人目の素数さん [sage] 2024/08/25(日) 05:42:43.86 ID:/qrXHaIo >>140 レスありがとうございます。 想定解と合致しました。 https://i.imgur.com/b226zjT.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/145
152: 132人目の素数さん [sage] 2024/08/25(日) 12:59:59.55 ID:/qrXHaIo >>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] http://rio2016.5ch.net/test/read.cgi/math/1723152147/152
メモ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.031s