24/08/12 05:21:42.18 7Qwbx0VG.net
>>13
Log[1+i] = log(2)/2 + i(π/4) = -θ + i(π/4),
[1+i]^[1-i] = (1+i)*(1+i)^{-i}
= (1+i)*e^{π/4 +iθ},
e^{π/4} = 2.19328
θ =－log(2)/2 = －0.3465736
21:１３２人目の素数さん
24/08/12 08:06:38.70 N57DabL7.net
>>19
レスありがとうございます。
人間Wolframに感銘。
22:１３２人目の素数さん
24/08/12 08:09:45.64 N57DabL7.net
さて、そろそろ帰り支度をするかな。
今日は祝日なので内視鏡バイトなし。
まあ、新入院や救急搬送でインセンティブ２栄一ゲットしたので
m3から購入した。
i,e,π計算のメモ
(I^E)^Pi // ComplexExpand // Simplify
% // N
I^(E Pi) // ComplexExpand // Simplify
% // N
(I^Pi)^E // ComplexExpand // Simplify
% // N
I^(Pi E) // ComplexExpand // Simplify
% // N
23:１３２人目の素数さん
24/08/12 08:10:14.89 Obx9rvEn.net
>>21
人間Wolframってバカにしてるのでは
24:１３２人目の素数さん
24/08/12 10:39:57.82 F0Hc68tw.net
>>22
おい尿瓶ジジイいつになったら無理数の証明できるようになるんだよ
25:１３２人目の素数さん
24/08/12 13:09:48.25 OufDPFxO.net
>>22
おいもう息してないのか
26:１３２人目の素数さん
24/08/12 21:05:05.17 7Qwbx0VG.net
>>8
「（某大物棋士の）兄貴たちは頭が悪いから東大へ行った」
という話の出所は芹沢（博文）九段という人もいるが、怪しい…
27:１３２人目の素数さん
24/08/14 14:46:24.61 8WM1s01c.net
Reverse@Sort@Counts@IntegerDigits[2024!,7]
28:１３２人目の素数さん
24/08/14 17:20:49.32 jFZUZf8i.net
a + b + c = 1,
aa + bb + cc = 3,
a^6 + b^6 + c^6 = 3,
を解け
29:１３２人目の素数さん
24/08/14 17:23:05.81 jFZUZf8i.net
実解は
(a,b,c) = (-1,1,1) (1,-1,1) (1,1,-1)
30:１３２人目の素数さん
24/08/14 17:32:46.58 jFZUZf8i.net
虚数解
t^3 －t^2 －t +5 = 0,
t = {－[2(31+3√105)]^{1/3} －[2(31-3√105)]^{1/3} + 1}/3,
=－1.594313016355
[1-t ±ｉ√(t-2-15/t)]/2 = 1.2971565081774 ± 1.205625150603ｉ
31:１３２人目の素数さん
24/08/14 18:19:32.60 tCfOv5Df.net
1+2+…+2024は何桁の整数か
32:１３２人目の素数さん
24/08/14 19:19:14.10 jFZUZf8i.net
(a,b,c) が
t = －1.594313016355
[(1-t) ±ｉ√(-t-1)･√(5-3t)]/2 = 1.2971565081774 ± 1.205625150603ｉ
の順列である…
33:１３２人目の素数さん
24/08/14 19:40:42.13 qeuhkrG0.net
>>28
f[l_,m_,n_,a_,b_,c_] := Solve[x^l+y^l+z^l==a && x^m+y^m+z^m==b && x^n+y^n+z^n==c,{x,y,z}]
f[1,2,6,1,3,3]
34:１３２人目の素数さん
24/08/14 19:49:00.45 jFZUZf8i.net
a + b + c = 1,
aa + bb + cc = 3,
a^7 + b^7 + c^7 = 1,
を解け。
実数解は (-1,1,1) (1,-1,1) (1,1,-1)
35:１３２人目の素数さん
24/08/14 20:03:11.26 jFZUZf8i.net
↑
虚数解
t'^3 － t'^2 － t' + 4 = 0,
t' = {－[(97+9√113)/2]^{1/3} －[(97-9√113)/2]^{1/3} + 1} /3
=－1.485584
[(1-t’) ±ｉ√(-1-t')･√(5-3t')] /2 = 1.242792 ± 1.07145315ｉ
36:１３２人目の素数さん
24/08/14 20:18:28.86 jFZUZf8i.net
a + b + c = 1,
aa + bb + cc = 3,
a^8 + b^8 + c^8 = 3,
を解け。
実数解は (-1,1,1) (1,-1,1) (1,1,-1)
37:１３２人目の素数さん
24/08/14 20:23:49.14 jFZUZf8i.net
↑
虚数解
t^3 － t^2 － t + 11/5 = 0,
t = {－[(121+9√161)/5]^{1/3} －[(121-9√161)/5]^{1/3} + 1}/3,
= －1.23930635
[(1-t) ±ｉ√(-t-1)･√(5-3t)]/2 = 1.119653175 ± 0.7221934276ｉ
38:１３２人目の素数さん
24/08/14 20:55:23.81 jFZUZf8i.net
a + b + c = 1,
aa + bb + cc = 3,
a^3 + b^3 + c^3 = 1 + 3(1 -t -t^2 +t^3),
39:１３２人目の素数さん
24/08/14 21:24:22.70 jFZUZf8i.net
基本対称式：
a + b + c = 1,
ab + bc + ca = －1,
abc = t(tt-t-1),
40:１３２人目の素数さん
24/08/15 00:23:28.29 zglV+U2W.net
(続き)
S_0 = 1 + 1 + 1 = 3,
S_1 = a + b + c = 1,
S_2 = aa + bb + cc = 3,
S_3 = 4 + 3t(tt-t-1),
漸化式
S_n = S_{n-2} + S_{n-2} + (1/3) (S_3 - 4) S_{n-3},
そこで
S_n = 1 となる奇数n
S_n = 3 となる偶数n
をさがす。
41:１３２人目の素数さん
24/08/15 00:30:52.81 zglV+U2W.net
a + b + c = 1,
aa + bb + cc = 3,
a^17 + b^17 + c^17 = 1,
を解け。
実数解は (a,b,c) = (-1,1,1) (1,-1,1) (1,1,-1)
42:１３２人目の素数さん
24/08/15 00:41:32.86 zglV+U2W.net
↑
虚数解
t^3 － t^2 － t + 3 = 0,
t = [－(35+3√129)^{1/3} － (35-3√129)^{1/3} + 1]/3
= －1.359304086
[(1-t) ±ｉ√(-t-1)･√(5-3t)]/2 = 1.179652043 ± 0.903013146ｉ
S_3 = -5,
abc = -3,
43:１３２人目の素数さん
24/08/15 01:32:00.94 zglV+U2W.net
>>38-40
t^3 － t^2 － t + (4 - S_3)/3 = 0,
より
t = (－[q+√(qq－64)]^{1/3} －[q-√(qq－64)]^{1/3} + 1) /3,
ここに q = (25－9･S_3)/2,
44:１３２人目の素数さん
24/08/15 10:18:26.54 5ggTekS4.net
bWin <- function(x){# 1 beats 3, 3 beats 2, 2 beats 1
if(length(unique(x))!=2 ) return(0) # no winner
u=sort(unique(x))
if(all(u==c(1,2))) return(sum(x==2)) # how many winners who won by 2
if(all(u==c(2,3))) return(sum(x==3)) # how many winners who won by 3
if(all(u==c(1,3))) return(sum(x==1)) # how many winners who won by 1
}
janken.till.winner.sim <- function(n){ # janken till someone wins
if(n<2) return(NA)
x=sample(1:3,n,replace = TRUE) # janken by n people
count=1 # counter of janken game
while(bWin(x)==0){ # if draw try again
x=sample(1:3,n,replace = TRUE) # janken again till someone win
count=count+1
}
c(bWin(x),count) # return how many winners and counter
}
janken.till.chamipion.sim <- \(n){
re=janken.till.winner.sim(n)
count=re[2]
winner=re[1]
while(winner>1){
re=janken.till.winner.sim(n)
count=count+re[2]
winner=re[1]
}
count
}
res=replicate(1e5,janken.till.chamipion.sim(9))
summary(res)
hist(res,col=2,breaks='scott',freq=FALSE,xlab='count',main='')
BEST::plotPost(res,showMode = TRUE,breaks='scott')
45:１３２人目の素数さん
24/08/15 10:18:53.46 5ggTekS4.net
９人で野球チームをつくる。
全員がピッチャーをやりたがったためジャンケンで決めることにする。
全員でのジャンケンから始めてその勝者でジャンケンを続けることを繰り返す。
ピッチャーが決まるまでのジャンケンの回数の期待値と中央値を求めよ。
46:１３２人目の素数さん
24/08/15 12:49:28.59 6mj/BVQi.net
>>45
高校生にまたバカにされたいのか？
47:１３２人目の素数さん
24/08/15 17:41:15.12 bFfiJSUV.net
janken simulator
j[n_] :=(
count=0;
Until[Length@Union@a==2,a=RandomChoice[Range[3],n];count++];
b=Sort@Union@a;
If[b=={1,2}, c=Count[a,2]];
If[b=={2,3}, c=Count[a,3]];
If[b=={1,3}, c=Count[a,1]] ;
{c,count}
)
j[4]
48:１３２人目の素数さん
24/08/15 17:59:59.43 bFfiJSUV.net
n=9;
j[n_] :=(
count=0;
Until[Length@Union@a==2,a=RandomChoice[Range[3],n];count++];
b=Sort@Union@a;
If[b=={1,2}, w=Count[a,2]];
If[b=={2,3}, w=Count[a,3]];
If[b=={1,3}, w=Count[a,1]] ;
{w,count}
)
For[{w,count}=j[n],w>1,k=j[w];w=k[[1]];count=count+k[[2]]]
{w,count}
49:１３２人目の素数さん
24/08/15 18:27:20.71 bFfiJSUV.net
小学生向きの問題
王様 と 王様でない人とはどちらが多いでしょうか?
50:１３２人目の素数さん
24/08/15 19:21:07.03 k34L4Drp.net
>>49
chatGPTに聞かないと答えが分からないアホ発見
51:１３２人目の素数さん
24/08/15 19:23:19.21 ysyyeRB1.net
j[n_] :=( (* n人でジャンケンして勝者が決まるまでの回数と勝者の数*)
count=0;
Until[Length@Union@a==2,a=RandomChoice[Range[3],n];count++];
b=Sort@Union@a;
If[b=={1,2}, winners=Count[a,2]];
If[b=={2,3}, winners=Count[a,3]];
If[b=={1,3}, winners=Count[a,1]] ;
{winners,count}
)
sim[n_] :=((* 勝者が一人になるまでの回数 *)
For[{winner,counts}=j[n],winner>1,k=j[winner];winner=k[[1]];counts=counts+k[[2]]];
counts
)
res9=Table[sim[9],10^5];
Histogram[res9]
Mean[res9]//N
Median[res9]
52:１３２人目の素数さん
24/08/15 19:23:48.55 ysyyeRB1.net
Rに移植
j=\(n){
a=sample(3,n,replace=TRUE)
count=1
while(length(unique(a))!=2){
a=sample(3,n,replace=TRUE)
count=count+1
}
b=sort(unique(a))
if(all(b==c(1,2))) winner=2
if(all(b==c(2,3))) winner=3
if(all(b==c(1,3))) winner=1
winners=sum(a==winner)
c(winners,count)
}
sim=\(n){
c=j(n)
winner=c[1]
count=c[2]
while(winner>1){
k=j(winner)
winner=k[1]
count=count+k[2]
}
count
}
res9=replicate(1e5,sim(9))
BEST::plotPost(res9,breaks='scott')
53:１３２人目の素数さん
24/08/15 19:29:14.64 s7yxKa8y.net
Wolfram
URLﾘﾝｸ(i.imgur.com)
R
URLﾘﾝｸ(i.imgur.com)
同じアルゴリズムなので同様の結果になった。
54:１３２人目の素数さん
24/08/15 19:56:17.63 k34L4Drp.net
>>53
なんで日本語不自由なのに数学できると思ってんだチンパンは
55:１３２人目の素数さん
24/08/15 20:59:51.79 k34L4Drp.net
>>53
結局無理数の証明はダンマリ？
56:１３２人目の素数さん
24/08/16 13:42:56.09 BuIOB4cE.net
>>55
ググれば答がでてくるようなのは放置。
んで、ググってもでてこない、3Dグラフの動画はまだぁ？
57:１３２人目の素数さん
24/08/16 13:44:09.89 BuIOB4cE.net
ハノイの塔のソルバー Wolfram版
hanoi[n_,from_:"A",via_:"B",to_:"C"] :=(
If[n>=1,hanoi[n-1,from,to,via];Print["move ",n, " from ",from, " to ",to];hanoi[n-1,via,from,to]]
)
hanoi[6]
58:１３２人目の素数さん
24/08/16 14:02:36.89 jFisYWC2.net
>>56
なんかwolframで答えみたいなの出してたけどw
結局高校数学の基礎すら分からんってことね
59:１３２人目の素数さん
24/08/16 14:05:08.96 D3c2lioP.net
>>56
24：卵の名無しさん：2024/08/16(金) 08:48:49.90 ID:eXoor/xc
>>20
50歳男性
食事の際に胸がつかえる感じがする
下の内視鏡所見
↓
URLﾘﾝｸ(i.imgur.com)
診断と治療は？内視鏡専門にしてるなら答えれるよなwww
無視したり答えられなかったらお前は底辺シリツ医以下の知能しか無いってことでwww
これの答えまだ？
60:１３２人目の素数さん
24/08/16 14:25:35.34 D3c2lioP.net
スレタイ以外のことでいくら発狂したところで無視されるのは当たり前だろマヌケ
無理数の>>31すら答えられないのかよ？
61:１３２人目の素数さん
24/08/16 15:15:01.70 ahhsBV/s.net
尿瓶ジジイID:BuIOB4cEにはできない問題を質問します
①√2は無理数であることを証明せよ。
②1+2+…+2024は何桁の整数か。
62:１３２人目の素数さん
24/08/16 15:31:11.25 bNupVHmJ.net
（1）リウヴィルの定理
p, q が自然数のとき
｜√2 － p/q｜ ≧ (6-4√2)/qq,
∴ √2 ≠ p/q (無理数)
63:１３２人目の素数さん
24/08/16 15:49:14.93 ahhsBV/s.net
当然ながら同じIDじゃないと認めません
64:１３２人目の素数さん
24/08/16 16:58:44.40 3S2KNH0L.net
本スレでは散々発狂してますが>>61に対してこっちでは一切ダンマリを決め込んでる模様
65:１３２人目の素数さん
24/08/16 17:29:04.80 NAPgnpHZ.net
問題 油分け算のソルバーを作れ。（言語は問わない）
R言語でのソルバー
URLﾘﾝｸ(i.imgur.com)
これをWolframに移植する気にならんなぁ。
Wolframでないと分数での厳密解が出せないというわけでもないから。
黒瀬のスパイスをつかってチキンジャーキーを作る方が楽しい。
焼き肉のタレでつくったら嫁から評判が悪かったので作り直し。
66:１３２人目の素数さん
24/08/16 17:46:49.14 NAPgnpHZ.net
>>59
輪状襞と縦走襞から好酸球性食道炎を疑う所見。
まあ、特異的な治療があるわけでもないので、無症状なら検診で遭遇しても生検しないこともあるなぁ。
A型胃炎を疑っても金がかかるだけなので、
無症状で貧血でもなければ、抗胃壁細胞抗体やガストリンも測定せずに
検診を毎年受けてくださいの説明だけにしていることもままある。
内視鏡やっているなら即答できる問題
鮮血を吐血したと紹介された患者に何を内服していないかを問診すべきか？
経管栄養している患者が下血したと病棟から呼ばれた何を内服していないかを問診すべきか？
67:１３２人目の素数さん
24/08/16 17:50:49.50 3S2KNH0L.net
尿瓶ジジイこちらの問題にはダンマリ決め込んでる模様
①√2は無理数であることを証明せよ。
②1+2+…+2024は何桁の整数か。
68:１３２人目の素数さん
24/08/16 18:05:04.95 3S2KNH0L.net
>>66
臨床問題wwwお前本当馬鹿なんだな
じゃあ自称消化器外科のお前に問題出してやるよ
↓下の画像の病名は？
URLﾘﾝｸ(i.imgur.com)
で、なんでこれの答えが肝腫瘍なんだよマヌケ
69:１３２人目の素数さん
24/08/16 19:06:09.18 3S2KNH0L.net
医者スレにて尿瓶ジジイ大発狂中www
5択問題すらまともに答えられずww
40：卵の名無しさん：[sage]：2024/08/16(金) 18:13:16.70 ID:6sTnd4ml
>>25
学生時代には女子医大生には息子が大変お世話になりました。
作業仮説
臨床やっている女医は不細工か不倫している。
（エビデンスレベルV：個人の経験 n=22)
55：卵の名無しさん：2024/08/16(金) 18:57:15.96 ID:VnGGyzKc
ほら次の問題だぞ
健康診断の腹部超音波検査で肝臓に異常を指摘されたため来院した．
腹部単純CT【図No.1】，腹部造影CT【図No.2（早期相）】，【図No.3（門脈相）】，【図No.4（平衡相）】および【図No.5（早期相・冠状断）】を示す．矢印は腫瘤を示す．
考えられるのはどれか．1つ選べ．
（a）肝細胞癌
（b）肝血管腫
（c）転移性肝癌
（d）肝内胆管癌
（e）肝限局性結節性過形成〈focal nodular hyperplasia：FNH〉
URLﾘﾝｸ(i.imgur.com)
選択肢までつけてやったぞ簡単だろ？ChatGPTにまた聞いてみるか？www
62：卵の名無しさん：[sage]：2024/08/16(金) 19:01:43.99 ID:6sTnd4ml
>>55
肝腫瘍
70:１３２人目の素数さん
24/08/17 11:44:44.70 dkCzRgHk.net
>>66
>>69のゴキジェットで無事死んじゃった
71:１３２人目の素数さん
24/08/17 13:19:41.98 dkCzRgHk.net
>>66ここのスレではダンマリ決め込んでるみたい
ようやく自分が医者でも東大卒でもないことがわかったの？
72:１３２人目の素数さん
24/08/17 14:29:59.78 mowI5mD2.net
(*
10人のうち血液型A,O,B,ABの人がそれぞれ4,3,2,1人いる。
この中から無作為に3人を選んだとき3人の血液型がすべて異なる確率を求めよ。
シミュレーションしてその数値が正しいことを検証せよ。
*)
b= {1, 2, 2, 3, 3, 3, 4, 4, 4, 4};
c = Subsets[b,{3}];
Total@Boole[Length@Union@# ==3 & /@ c]/Length@c
% // N
k=10^6;
Table[Boole[Length@Union@RandomSample[b,3]==3],k] // Mean // N
73:１３２人目の素数さん
24/08/17 14:30:43.53 mowI5mD2.net
Wolfram言語に少々なれてきた。
74:１３２人目の素数さん
24/08/17 14:48:01.69 mowI5mD2.net
(*
10人のうち血液型A,O,B,ABの人がそれぞれ4,3,2,1人いる。
この中から無作為に n 人を選んだとき血液型がすべて異なる確率を求めるソルバーを作れ。
*)
solve[n_] :=(
b= {1, 2, 2, 3, 3, 3, 4, 4, 4, 4};
c = Subsets[b,{n}];
Total@Boole[Length@Union@# ==Min[n,4] & /@ c]/Length@c
)
Table[solve[n],{n,1,10}]
75:１３２人目の素数さん
24/08/17 14:49:16.87 mowI5mD2.net
In[3]:= Table[solve[n],{n,1,10}]
7 5 4 2 7 19 7 9
Out[3]= {1, -, --, --, -, --, --, -, --, 1}
9 12 35 7 15 30 9 10
76:１３２人目の素数さん
24/08/17 15:09:03.74 dkCzRgHk.net
>>75
で、尿瓶チンパンジジイはこれはいつになったら解けるの？
①√2は無理数であることを証明せよ。
②1+2+…+2024は何桁の整数か。
77:１３２人目の素数さん
24/08/17 15:10:58.66 dkCzRgHk.net
頭が悪すぎてフルボッコにされていることも理解できないみたいだねw
>>69に5択すら選べないみたいだしそこまでしてバカにされたい？w
78:１３２人目の素数さん
24/08/17 18:23:34.97 g0prVaq1.net
質問スレなのに質問が無い
79:１３２人目の素数さん
24/08/18 06:09:54.38 Qg46pTwQ.net
(* 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={1};
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]]]
80:１３２人目の素数さん
24/08/18 06:11:04.79 Qg46pTwQ.net
>>78
では、質問
王様（王位にある人） と 王様でない人とはどちらが多いでしょうか?
81:１３２人目の素数さん
24/08/18 06:46:57.60 Qg46pTwQ.net
(* 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
82:１３２人目の素数さん
24/08/18 11:13:34.46 fBcVsaRS.net
>>80
スレタイ読めないアホ発見
83:１３２人目の素数さん
24/08/18 13:29:06.79 S1VoGAV5.net
>>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]
84:１３２人目の素数さん
24/08/18 14:21:09.42 S1VoGAV5.net
４０人のクラスで代表１人をジャンケンで選ぶ。
全員でのジャンケンから始めてその勝者でジャンケンを続けることを繰り返す。
代表が決まるまでのジャンケンの回数の期待値を求めよ。
答は分数もしくは小数第１位を四捨五入にした整数でよい。
(* 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={0};
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]]];
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]]];
re40=Table[calc[m],{m,4,40}][[-1]][[-1]]
N[re40,7]
N[re40,8]
85:１３２人目の素数さん
24/08/18 15:28:11.50 fBcVsaRS.net
>>84
自称有能ならレス乞食しないで自分で勝手にやってろよ
86:１３２人目の素数さん
24/08/18 20:29:25.07 S1VoGAV5.net
このスレならこの程度の問題にしていただきたいね。
iのi乗根は無理数であることを示せ
1+2+3+...+2024!を7進法で表示すると何桁の数になるか算出せよ。
87:１３２人目の素数さん
24/08/18 21:48:59.30 LYw/XosF.net
>>86
そんな簡単ならさっさと答えろよレス乞食
88:１３２人目の素数さん
24/08/19 05:30:45.90 qPuO+ITh.net
(* 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
89:１３２人目の素数さん
24/08/19 20:37:29.20 LNMh0Kop.net
だよなあ
今回段取りだけはすごい強いわなんなんだよな
90:１３２人目の素数さん
24/08/19 21:07:59.86 qEmD6FDJ.net
実際やってるような気が付いたから大丈夫
トレンド1位スナイプできるレベルで技は３ー４個しか使えずインフレして自分達で爆死させて貰いますた
女は恋愛すんな
91:１３２人目の素数さん
24/08/19 21:25:54.52 fCkpErGW.net
スノのいいドラマは被らないとこでライブ配信やって下げるのに
ワーキングプアしてるのか
お前らの
92:１３２人目の素数さん
24/08/19 21:31:33.15 wxqnnk1E.net
サガスカおもろかったやん
死んでない
検査装置で、脳梗塞でも書き込めなくなってきた
93:１３２人目の素数さん
24/08/19 21:51:34.12 Uz3cJ0Q4.net
>>14
典型的な買い
URLﾘﾝｸ(i.imgur.com)
94:１３２人目の素数さん
24/08/19 22:04:47.75 nSXgUTL8.net
ヨーモニ～💏
あいがみも酔ってる方だと思うわ
投資スタイルなんてだいぶ昔に統一関係議員全員美少女化したのここ？
URLﾘﾝｸ(0kvs.if.d42c)
95:１３２人目の素数さん
24/08/19 22:08:34.03 fivowbXI.net
>>79
逆にテレビ千鳥は深夜に戻っとるやん！
なんのため
ガチで
96:１３２人目の素数さん
24/08/19 22:37:05.09 IyUu92sy.net
試合で一回も成功してたしね
暴露で集めただけに近いし
97:１３２人目の素数さん
24/08/19 22:51:33.83 XQP/Eri1.net
関連はようございます
URLﾘﾝｸ(i.imgur.com)
URLﾘﾝｸ(i.imgur.com)
98:１３２人目の素数さん
24/08/19 23:18:24.90 CBhSN87T.net
あのキャラ絶対芸能界向きだと自サーバーにクレカ情報渡したくないわ
28000まで同じような順位なん？
決済
アイスタ素直に逮捕されて使い道ない
99:１３２人目の素数さん
24/08/19 23:28:07.46 y+uVizcc.net
>>88
ウィスキーとか飲むイメージだな
ソシャゲは最早ガチャでエロ絵集めるだけやで
100:１３２人目の素数さん
24/08/19 23:40:01.05 TuWh8Fdq.net
他全滅
101:１３２人目の素数さん
24/08/19 23:47:49.30 pUs7ogYC.net
すかへんくくしあろはせまけひいけきねぬひきさえなとやいろに
102:１３２人目の素数さん
24/08/20 05:30:38.98 O+DDQLmu.net
"
n=11;
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,#]&;
"
str=NULL
str[1]="pn={jan[n,1]};"
str[2]="Total@Flatten@Table[jan[n,x1]jan[x1,1],{x1,2,n}]//AppendTo[pn,#]&;"
str[3]="Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,1],{x1,2,n},{x2,2,x1}]//AppendTo[pn,#]&;"
str[4]="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,#]&;"
str[5]="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,#]&;"
for(i in 2:30){
y=paste0(",x",i,"]jan[x",i,",1],")
tmp=sub(",1],",y,str[i])
z=paste0("},{x",i,",2,x",i-1,"}]")
str[i+1]=sub("}]",z,tmp)
}
103:１３２人目の素数さん
24/08/20 08:08:58.34 O+DDQLmu.net
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]
)
n=11;
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,#]&;
...
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,x8]jan[x8,x9]jan[x9,x10]jan[x10,x11]jan[x11,x12]jan[x12,x13]jan[x13,x14]jan[x14,x15]jan[x15,x16]jan[x16,x17]jan[x17,x18]jan[x18,x19]jan[x19,x20]jan[x20,x21]jan[x21,x22]jan[x22,x23]jan[x23,x24]jan[x24,x25]jan[x25,x26]jan[x26,x27]jan[x27,x28]jan[x28,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5},{x7,2,x6},{x8,2,x7},{x9,2,x8},{x10,2,x9},{x11,2,x10},{x12,2,x11},{x13,2,x12},{x14,2,x13},{x15,2,x14},{x16,2,x15},{x17,2,x16},{x18,2,x17},{x19,2,x18},{x20,2,x19},{x21,2,x20},{x22,2,x21},{x23,2,x22},{x24,2,x23},{x25,2,x24},{x26,2,x25},{x27,2,x26},{x28,2,x27}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,x8]jan[x8,x9]jan[x9,x10]jan[x10,x11]jan[x11,x12]jan[x12,x13]jan[x13,x14]jan[x14,x15]jan[x15,x16]jan[x16,x17]jan[x17,x18]jan[x18,x19]jan[x19,x20]jan[x20,x21]jan[x21,x22]jan[x22,x23]jan[x23,x24]jan[x24,x25]jan[x25,x26]jan[x26,x27]jan[x27,x28]jan[x28,x29]jan[x29,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5},{x7,2,x6},{x8,2,x7},{x9,2,x8},{x10,2,x9},{x11,2,x10},{x12,2,x11},{x13,2,x12},{x14,2,x13},{x15,2,x14},{x16,2,x15},{x17,2,x16},{x18,2,x17},{x19,2,x18},{x20,2,x19},{x21,2,x20},{x22,2,x21},{x23,2,x22},{x24,2,x23},{x25,2,x24},{x26,2,x25},{x27,2,x26},{x28,2,x27},{x29,2,x28}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,x8]jan[x8,x9]jan[x9,x10]jan[x10,x11]jan[x11,x12]jan[x12,x13]jan[x13,x14]jan[x14,x15]jan[x15,x16]jan[x16,x17]jan[x17,x18]jan[x18,x19]jan[x19,x20]jan[x20,x21]jan[x21,x22]jan[x22,x23]jan[x23,x24]jan[x24,x25]jan[x25,x26]jan[x26,x27]jan[x27,x28]jan[x28,x29]jan[x29,x30]jan[x30,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5},{x7,2,x6},{x8,2,x7},{x9,2,x8},{x10,2,x9},{x11,2,x10},{x12,2,x11},{x13,2,x12},{x14,2,x13},{x15,2,x14},{x16,2,x15},{x17,2,x16},{x18,2,x17},{x19,2,x18},{x20,2,x19},{x21,2,x20},{x22,2,x21},{x23,2,x22},{x24,2,x23},{x25,2,x24},{x26,2,x25},{x27,2,x26},{x28,2,x27},{x29,2,x28},{x30,2,x29}]//AppendTo[pn,#]&;
pn
Accumulate[pn]
% // N
104:１３２人目の素数さん
24/08/20 08:10:01.22 O+DDQLmu.net
自民党総裁選は候補者11人という。
11人でジャンケンをして勝者が複数なら勝者だけでジャンケンを繰り返すという方式で総裁を決めるとき
ジャンケンの回数（アイコも１回と数える）の最頻値と中央値を求めよ。
105:１３２人目の素数さん
24/08/20 08:11:47.85 O+DDQLmu.net
Rの文字列操作関数を駆使してWolframのコードを作成するという戦略を採った。
Wolframだけで完結するコードの投稿を期待します。
106:１３２人目の素数さん
24/08/20 08:24:57.34 cBs2zsg2.net
ここまで教えてもらってまだできんww wwwwww
無能wwwwwwwwwwwwwww
107:１３２人目の素数さん
24/08/20 23:02:05.74 3diBYk2P.net
冒頭、
>> このスレは医者・東大卒の人物専用スレです。
>> その他の人が書き込むことは許されません。
とあるので、許されないのかもしれませんが、向こうに書き込むよりはましかと思いこちらに投稿します。
>> 答に自信が持てないので質問します。
>>
>> 自民党総裁選は候補者11人という。
...
>> 199602482890927685979724552820321419790
>> Out[14]= ----------------------------------------
>> 1330279464729113309844748891857449678409
>>
>> で合っているでしょうか？
私も自信があるわけではないが、コードと得た答えを記します。
p[a_,b_]:=p[a,b]=If[a==b,1-(2^a-2)/3^(a-1),Binomial[a,b]/3^(a-1)]
m=11;
f=Flatten[Table[{m,m-x1,m-x2,m-x3,m-x4,m-x5,m-x6,m-x7,m-x8,m-x9,1},
{x1,0,m-1},{x2,x1,m-1},{x3,x2,m-1},{x4,x3,m-1},{x5,x4,m-1},
{x6,x5,m-1},{x7,x6,m-1},{x8,x7,m-1},{x9,x8,m-1}],8];
Total[Product[p[#[[k]],#[[k+1]]],{k,1,10}] & /@f]//InputForm
449687340186660888579056289638229806808082/2909321189362570808630465826492242446680483
108:１３２人目の素数さん
24/08/21 08:28:52.22 trPnwZW4.net
>>106
で、３Dグラフの動画はまだできんのだが？
109:１３２人目の素数さん
24/08/21 09:46:56.78 4xmRA5nu.net
>>108
なんだ、あんたも無能なんだ
110:１３２人目の素数さん
24/08/21 10:42:51.81 trPnwZW4.net
>>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;
(* 参加者数の変遷を列挙 *)
f=Flatten[Table[{m,m-x1,m-x2,m-x3,m-x4,m-x5,m-x6,m-x7,m-x8,m-x9,1},
{x1,0,m-1},{x2,x1,m-1},{x3,x2,m-1},{x4,x3,m-1},{x5,x4,m-1},
{x6,x5,m-1},{x7,x6,m-1},{x8,x7,m-1},{x9,x8,m-1}],8];
(* 変遷確率積の総和*)
Total[Product[p[#[[k]],#[[k+1]]],{k,1,10}] & /@f]
% // N
正しいように思いますので自分のコードのデバッグ作業に移ります。
111:１３２人目の素数さん
24/08/21 11:14:42.94 trPnwZW4.net
デバッグして、数値が合致しました。
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]
)
In[3]:= n=11;
In[4]:= (* 10回以下で終了する確率*)
In[5]:= Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,x8]jan[x8,x9]jan[x9,1]
,{x1,1,n},{x2,1,x1},{x3,1,x2},{x4,1,x3},{x5,1,x4},{x6,1,x5},{x7,1,x6},{x8,1,x7},{x9,1,x8}]
449687340186660888579056289638229806808082
Out[5]= -------------------------------------------
2909321189362570808630465826492242446680483
112:１３２人目の素数さん
24/08/21 11:27:27.25 trPnwZW4.net
(* 各回ごとの確率 *)
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,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5},{x7,2,x6}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,x8]jan[x8,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5},{x7,2,x6},{x8,2,x7}]//AppendTo[pn,#]&;
Total@Flatten@Table[jan[n,x1]jan[x1,x2]jan[x2,x3]jan[x3,x4]jan[x4,x5]jan[x5,x6]jan[x6,x7]jan[x7,x8]jan[x8,x9]jan[x9,1],{x1,2,n},{x2,2,x1},{x3,2,x2},{x4,2,x3},{x5,2,x4},{x6,2,x5},{x7,2,x6},{x8,2,x7},{x9,2,x8}]//AppendTo[pn,#]&
Total@pn
113:１３２人目の素数さん
24/08/21 11:44:00.03 trPnwZW4.net
シミュレーションによる検証
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;
If[b=={1,2}, winners=Count[a,2]];
If[b=={2,3}, winners=Count[a,3]];
If[b=={1,3}, winners=Count[a,1]] ;
{winners,count}
)
In[2]:=
In[2]:= sim[n_] :=((* 勝者が一人になるまでの回数 *)
For[{winner,counts}=j[n],winner>1,k=j[winner];winner=k[[1]];counts=counts+k[[2]]];
counts
)
In[3]:= res11=Table[sim[11],10^6];
In[4]:= Histogram[res11,"Scott","PDF"]
Out[4]= -Graphics-
In[5]:= Mean[res11] // N
Out[5]= 34.9504
In[6]:= Median[res11]
Out[6]= 27
In[7]:= N@Mean@Boole[#<=10&/@res11]
Out[7]= 0.154785
In[8]:= 449687340186660888579056289638229806808082/2909321189362570808630465826492242446680483
449687340186660888579056289638229806808082
Out[8]= -------------------------------------------
2909321189362570808630465826492242446680483
In[9]:= % // N
Out[9]= 0.154568
114:１３２人目の素数さん
24/08/21 16:33:52.27 fRgQD6bE.net
>>108
日本語不自由な無能発見
115:１３２人目の素数さん
24/08/21 16:34:11.10 EiqnP9Ex.net
>>107
あからさまな自演で草
116:１３２人目の素数さん
24/08/21 17:53:09.19 EiqnP9Ex.net
>>113
>>69にはいつになったら答えられるんだよゴミ
117:１３２人目の素数さん
24/08/21 17:57:45.90 trPnwZW4.net
>>115
Flattenの使い方とか熟練者のコード解析は勉強になる。
118:１３２人目の素数さん
24/08/21 18:09:28.46 EiqnP9Ex.net
>>117
そもそも脳内医者なのになんでこのスレ書き込めるんだよ？
119:１３２人目の素数さん
24/08/21 19:14:25.87 EiqnP9Ex.net
>>117
>>69にはいつまでダンマリ決め込んでんだよ？
120:１３２人目の素数さん
24/08/21 20:27:00.79 rsVL5MYO.net
アイスタマジヤバいの来そうやないか？
パーマかけたんだね
121:１３２人目の素数さん
24/08/21 20:34:12.34 YldtGvrJ.net
逆に言えば、僕らはそれを隠してたんだね
122:１３２人目の素数さん
24/08/21 21:13:41.42 orJvOL8Q.net
土曜の昼間にやることがあるやんけ
今週中に脳梗塞でもないからな
123:１３２人目の素数さん
24/08/21 22:21:32.56 EiqnP9Ex.net
>>117
134：卵の名無しさん：2024/08/21(水) 22:14:22.45 ID:Mk9xahY2
>>131
即答出来なかった癖にってwww
そりゃお前の脳内妄想医療なんぞ想像出来るわけないだろ臨床に全くそってないのに
普通は議題にも上がらないからな
仮にリファンピシン何て飲んでたらGFを消火器にふる前にそういう情報出てるよな
わざわざ他科にコンサルトしてんだから
正確には常に気をつける薬剤では無く飲んでたらそれが原因だと気づく薬剤なんだよ
お前、現在リファンピシン飲んでる結核患者なんてどれぐらいいるのか知ってるのか？
妄想も大概にしろよ
135：卵の名無しさん：2024/08/21(水) 22:17:01.89 ID:Mk9xahY2
>>132
どう考えてもわかってなかったよな？
変なプログラム書いてる暇あったらすぐに反論できたよね？
後出し乙
もう医者のフリするのやめない？馬鹿なんだからwww
124:１３２人目の素数さん
24/08/22 11:17:43.50 moSOEVeB.net
体に良い親人と出会っていたんだろうな。
ラメーン食いたいもんかラーメンて
125:１３２人目の素数さん
24/08/22 11:20:40.23 rXGE90op.net
付け焼き刃的に持ってないな
それにたいする答えてのは好きじゃないとしたら
126:１３２人目の素数さん
24/08/22 11:46:32.27 +idgYRUh.net
>>32
俺くらいになるね
でもチケ代高い割には何もやらない方がでかいな
画面越しのヒントになるってだけ
127:１３２人目の素数さん
24/08/22 12:19:30.59 cPl6boaC.net
なんかぞわぞわするんだよお婆さん
今はSNSでの煽り行為や公共の福祉以外に何を不幸自慢してるんだろうなぁ
枠拡充120万くらいなのに
128:１３２人目の素数さん
24/08/22 12:25:13.71 KCeZ3hz8.net
>>105
自分の立場だったら記念切手の感覚でやってるフリしてる奴らて何も分かって怖い
横転した統一系のサークルで勧誘してるねー
129:１３２人目の素数さん
24/08/22 13:22:36.77 8qGY+6eb.net
>>117
142：卵の名無しさん：2024/08/22(木) 09:09:23.36 ID:Hxf3Cpig
>>137
吐血と下血の区別ぐらいつくわwww
喀血もわからないお前じゃないんだから
んで、リファンピシンで治療するような病気を診ている呼吸器内科がいる規模の病院なのに
消火器内科ではなく自称外科医のお前に内視鏡を依頼するのは何故？
→脳内医療だからwww
結核だろうが肺MAC症だろうが外科に内視鏡頼むような病院でそんな病気診る呼吸器内科がいるわけないじゃん
そもそもリファンピシン飲んでるのわかってるのにその呼吸器内科は内視鏡やる医者にその事を伝えないのか？滅茶苦茶だなwww
どう見ても脳内医療
130:１３２人目の素数さん
24/08/22 16:23:48.72 0a4dfs+c.net
>>80
王位にある人は 藤井聡太七冠
永世王位にある人は 大山15世(故人)、中原16世、羽生19世(永世七冠)
(大意)
王位戦七番勝負は現在 3－1 で、渡辺九段にあと1勝すれば防衛。
連続5期に到達するので、永世王位に就く。
131:１３２人目の素数さん
24/08/22 20:15:09.90 ZYoBe3+E.net
f[n_,x_] = f[n]=If[
n==1,
1/(1-x),
1/(1-(1-(2^n-2)/3^(n-1))x) x Sum[Binomial[n,m]/3^(n-1)f[m,x],{m,1,n-1}]]
Table[{n,SeriesCoefficient[f[11,x],{x,0,n}]},{n,9,15}]
Table[{n,SeriesCoefficient[f[11,x]*(1-x),{x,0,n}]},{n,9,5}]
Table[{n,N@SeriesCoefficient[f[11,x]*(1-x),{x,0,n}]},{n,0,20}]
Table[{n,Sign@SeriesCoefficient[f[11,x]*(1-2*x+x^2),{x,0,n}]},{n,0,20}]
ListPlot[Table[SeriesCoefficient[f[11,x]*(1-x),{x,0,n}],{n,30}]]
N@SeriesCoefficient[1/(1-x)-f[11,x],{x,1,0}]
132:１３２人目の素数さん
24/08/23 04:02:36.07 /G4Ss0QX.net
>>86
ｉ = e^{iπ/2} のｉ乗根は e^{π/2}
S = 1 + 2 + …… + 2024!
= (2024!)(2024!+1)/2
≒ (1/2)(2024!)^2
≒ 2024^{4049}･π･e^{－4048} Stirlingの式
log(S) = 4049･log(2024) + log(π) －4048
≒ 26777.497572
log(S)/log(7) = 13760.911615…
∴ 13761桁
133:１３２人目の素数さん
24/08/23 07:15:27.28 eFoV/xUL.net
>>132
レスありがとうございます。
想定解の数値と合致です。
134:１３２人目の素数さん
24/08/23 07:18:46.09 eFoV/xUL.net
11人でジャンケンしたときの終了回数のシミュレーションは
期待値、中央値はシミュレーション値と理論値がほぼ合致したが、
最頻値はずいぶんと乖離した。
練習問題
1人の勝者が決まるまで11人でジャンケンをしたとき
終了までのジャンケンの回数の最頻値とその確率を求めよ。
勝者が複数なら勝者間でジャンケンを継続し、アイコも１回と数える。
135:１３２人目の素数さん
24/08/23 08:31:17.84 Am/2METs.net
>>133
数値一致wwww
等号の意味すら分からないのかよwwww
136:１３２人目の素数さん
24/08/23 13:14:58.90 g6kRvRuw.net
シミュレーションでの分布
URLﾘﾝｸ(i.imgur.com)
最頻値は理論値とずれた。
137:１３２人目の素数さん
24/08/23 13:40:44.06 Cv/P/3fy.net
尿瓶ジジイ＝>>136は高校数学の基礎の基礎であるこちらの問題にはダンマリ決め込んで懲りずにチンパン統計をぶつぶつほざいている模様
①√2は無理数であることを証明せよ。
②1+2+…+2024は何桁の整数か。
138:１３２人目の素数さん
24/08/23 14:00:14.64 wqZOBs8f.net
当たり前やろ
こんな分布1000000回ごときのシミュで答え出るわけないやろ
この無能頭わいてんのかね？
無能ってかわいそう
139:１３２人目の素数さん
24/08/23 18:49:51.27 xQPSXxbD.net
150：卵の名無しさん：2024/08/23(金) 06:33:17.51 ID:0RpmxFnA
尿瓶ジジイが本当に臨床統計を理解しているか確認する問題
原著論文の抄録を示す．
背景・目的：腰痛には心理的要因が影響していることが知られているが，心理行動学的介入が腰痛患者のQOLに効果があるかは不明である．本研究では，腰痛患者に対する認知行動療法が，患者の健康関連QOLに与える効果を検討する．
方 法：参加した50施設の整形外科外来に慢性腰痛で受診中の患者のうち，同意を得られた880名（平均年齢 ± 標準偏差 ＝ 72.4 ± 4.5歳，男女比 ＝ 1：1.2）を本研究の対象とした．ランダム割付けを行い，通常治療＋認知行動療法を行う介入群と通常治療のみの比較対照群に分類した．ベースラインと12か月時点の健康関連QOLをSF 36 Health Survey日本語版〈SF-36〉の下位項目「活力」で測定し，スコアの変化量の2群間の差をt検定で解析した．Intention to Treat解析を用いた．
結 果：12か月の追跡率は70％であった．ベースラインのSF-36「活力」に2群で統計学的有意差はみられなかった．12か月後の「活力」の変化量は比較対象群と比べて介入群で統計学的有意な改善を認めた（P値 ＜ 0.001）．
結 論：慢性腰痛症に対する認知行動療法は健康関連QOLを改善させる．
本研究結果を臨床応用するうえで最も確認したい論文中の情報はどれか．1つ選べ．
（a）P値（絶対値）
（b）各群の追跡率
（c）二重盲検順守率
（d）per protocol解析結果
（e）変化量の2群間差の値（95％信頼区間）
臨床統計得意なんでしょ？5択だし流石に答えられるよね？
140:１３２人目の素数さん
24/08/23 19:36:15.43 R15w90it.net
>>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
141:１３２人目の素数さん
24/08/23 23:29:49.45 /G4Ss0QX.net
>>132
S = 1 + 2 + …… + n!
= (n!)(n!+1)/2
≒ (1/2)(n!)^2
≒ n^{2n+1}･e^{－2n}･π･(1+1/(6n)) Stirlingの式
log(S) ≒ (2n+1)･log(n) - 2n + log(π) + 1/(6n)
≒ 26777.49765435
log(S)/log(7) ≒ 13760.91165738…
∴ 13761桁
142:１３２人目の素数さん
24/08/23 23:37:21.14 /G4Ss0QX.net
>>132
56155314531336…… (13761桁)
143:１３２人目の素数さん
24/08/24 14:31:18.38 CzjhaWSa.net
x,yは x<=y の正整数とし、x^2+y^2は25の倍数である。
これを満たす{x,y}の組み合わせを x+y が小さい順に並べる。
10番目、20番目、30番目の{x,y}の組み合わせを求めよ。
xy=Flatten[Table[{x,y},{x,1,25},{y,1,25}],1];
re=Select[xy,Mod[ #[[1]]^2 + #[[2]]^2,25] == 0 && #[[1]] <= #[[2]] &];
Table[SortBy[re,Total][[i]],{i,{10,20,30}}]
144:１３２人目の素数さん
24/08/25 04:38:29.64 BpgtRhSp.net
あの人気SNSのLiteバージョンでPayPayやAma券などにチェンジできる￥5000相当のﾎﾞｲﾝﾄを配布中！
＊無印ユーザーの方はアプリからログアウト&削除必須
1.SIMの入ったスマホ・タブレットを用意する
2.以下のサイトからアプリをダウンロード（ダウンロードだけでまだ起動しない）
URLﾘﾝｸ(zws.im)
3.ダウンロード完了後、もう一度上記アドレスのリンクからアプリを起動
4.アプリ内で無印で未使用の電話番号かメルアドを使用して18歳以上で登禄
5.「10日間連続のﾁｪｯｸｲﾝ｣で合計で￥5000相当のﾎﾞｲﾝﾄ入手
↑ 重要！
ﾎﾞｲﾝﾄはPayPayやAma券にチェンジできる！
更にご家族等に教えたり、ノーマルタスクをこなせば更にﾎﾞｲﾝﾄを入手できます
145:１３２人目の素数さん
24/08/25 05:42:43.86 /qrXHaIo.net
>>140
レスありがとうございます。
想定解と合致しました。
URLﾘﾝｸ(i.imgur.com)
146:１３２人目の素数さん
24/08/25 05:50:42.70 /qrXHaIo.net
>>139
問題を読むだけでも面倒だが、(e)だろ。
臨床統計ｗの問題
ゴルゴ13は100発100中
ゴルゴ14は10発10中
ゴルゴ15は1発1中
とする。
各々10000発撃ったとき各ゴルゴの命中数の期待値はいくらか？
147:１３２人目の素数さん
24/08/25 08:47:57.41 XxCJyA5c.net
>>144
やるしか無いわ
148:１３２人目の素数さん
24/08/25 09:15:08.50 rhW1eedt.net
>>144
やるしかないね
149:１３２人目の素数さん
24/08/25 10:53:44.75 0Dxrgsfy.net
数学板でも高校生にボコされた上にスクリプト爆撃を受けてもう元気がなくなった様子w
>>146
ずっとダンマリ決め込んでたくせに今更後出しジャンケンだっさw
189：卵の名無しさん：2024/08/24(土) 08:17:30.15 ID:q7UaKMDX
尿瓶ジジイは言い訳ばっかで結局答えわからんかったって事でとりあえず答え書いとくわ
P値とは統計的仮説検定を行う際，「比較する2群の結果に差がない」という仮説（帰無仮説）どおりになる確率のことである．すなわちP = 0.001とは，標本抽出して検定を1,000回行うと１回だけ帰無仮説どおりになることである．
この際，2群の差が大きいほどP値は小さくなる．一方，P値が小さいことは2群の差が大きいことを意味するとは限らない．例えば，サンプルサイズが大きければP値は低くなる傾向がある．ゆえにP値の絶対値が介入の効果の大きさを直接示しているわけではない．
ランダム化比較研究の結果を広く臨床に応用する際，介入の効果（有効性）の大きさを考慮することは大事である．今回の研究においては介入の効果を，変化量の2群間の差や改善割合の2群間の比として95％信頼区間で示すというのが可能である．また，効果量を算出する方法もある．2群の平均値の差であれば，t 値と自由度とを使って計算するr 値や，平均や標準偏差から計算する d 値で示すことができる．効果量も絶対的な指標ではないものの，介入効果を示す1つとして覚えておくと良い．
以上から，選択肢の中では（e）を正解，（a）を不正解．
150:１３２人目の素数さん
24/08/25 11:52:28.29 rhW1eedt.net
>>146
結局答え教えてもらってから書き込んでんじゃん
ダサッ！
151:１３２人目の素数さん
24/08/25 12:09:41.09 0Dxrgsfy.net
>>146
問題文読むだけでも面倒って自己紹介？
自分は散々人にスレチな出題をレス乞食しておいていざ人にご指名で質問されたらご丁寧に答え教えてくた後出しジャンケンとか一体何様なの？w
こんなのが医者や東大卒だと思う人レスして下さい
152:１３２人目の素数さん
24/08/25 12:59:59.55 /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]
153:１３２人目の素数さん
24/08/25 13:02:05.24 /qrXHaIo.net
一次方程式の問題を解けないようなシリツ医を相手にしないのよ。
154:１３２人目の素数さん
24/08/25 13:04:45.53 /qrXHaIo.net
すぐに読める問題
ゴルゴ13は100発100中
ゴルゴ14は10発10中
ゴルゴ15は1発1中
とする。
各々10000発撃ったとき各ゴルゴの命中数の期待値はいくらか？
155:１３２人目の素数さん
24/08/25 13:08:11.41 /qrXHaIo.net
>>138
最初の参加人数が少なければシミュレーションと理論値が合致するぞ。
156:１３２人目の素数さん
24/08/25 13:29:00.58 u4or3250.net
>>144
こういう情報を知りたかった
157:１３２人目の素数さん
24/08/25 13:29:42.45 0Dxrgsfy.net
>>153
スレタイも理解できないチンパンが相手にされてないからここで発狂してんだろうがw
158:１３２人目の素数さん
24/08/25 13:30:37.46 /qrXHaIo.net
練習問題
１５人のラグビーチームでジャンケンをしてキャプテンを１人決めることになった。
１５人全員でジャンケンを始めて勝者どおしでジャンケンを続けて勝者が１人になるまで続ける。
アイコも１回と数えるときキャプテンが決まるまでのジャンケンの回数の中央値は１１４回である。
問題 ４０人のクラスでジャンケンをして学級委員を１人決めるときのジャンケンの回数の中央値を算出せよ。
159:１３２人目の素数さん
24/08/25 13:51:00.26 0Dxrgsfy.net
人には散々チンパン数学()を一方的に出題しておいてレス乞食はする癖に(しかしほぼ相手にされない)いざ自分が同じことされると訳の分からない言い訳発狂ダンマリ後出しジャンケンとかダサすぎるw
5択すら答えられないのかよ？どこまで頭悪いんだ？
160:１３２人目の素数さん
24/08/25 13:56:49.64 /qrXHaIo.net
>>159
ChatGPTで正解が返ってくるような問題は面白くないんだなぁ。
161:１３２人目の素数さん
24/08/25 14:00:19.08 SbwyOjer.net
>>155
合致wwww
主観かよwwwクズwww
162:１３２人目の素数さん
24/08/25 14:01:32.22 iq/W3tzQ.net
>>160
簡単なプログラムで数値計算するだけの問題は面白くないぞwww
163:１３２人目の素数さん
24/08/25 14:02:46.70 /qrXHaIo.net
>>159
Wolframを扱える方はちゃんとレスを返してくるんだなぁ。
熟練者のコードを解読するのは勉強になる。
お陰でRのwhich関数に相当する関数がPositionであることを学んだので
早速、それを組み込んで道具箱に保管。
Wolframが使えないPhimoseくんが返してくるのは自演認定と医師板からのコピペくらいだなぁ。
164:１３２人目の素数さん
24/08/25 14:05:54.84 1ZiEG71e.net
>>160
結局お得意のGPTでも答えられなかったから答えが出るまでダンマリ決め込むしかなかったんだろ？
悔しかったら答えが出る前に答えだせよタコw
>>163
アンタと同じ文体でアンタがシカトされてる時にしれっと書き込み
これが自演じゃなくて何なんだよ？
165:１３２人目の素数さん
24/08/25 14:11:13.40 QiSGCEaw.net
>>154
ChatGPTだと期待値はすべて10000発と返ってきた。
95%信頼区間を算出させたら、
どのゴルゴに対しても、命中数が確定的であり、すべての発射が命中するため、95%信頼区間は
[10000,10000] です。
と返ってきた。
ChatGPTが正解を返さないような問題が面白い。
166:１３２人目の素数さん
24/08/25 14:13:20.84 QiSGCEaw.net
>>164
ChatGPTに入力してみたら、正解が返ってくるから。
167:１３２人目の素数さん
24/08/25 14:16:20.38 QiSGCEaw.net
>>162
その簡単なプログラムをWolframで書けないのがいるんだよなぁ。
確率問題は総当たりとシミュレーションで答を検討できるから(・∀・)ｲｲ!!
168:１３２人目の素数さん
24/08/25 14:18:57.19 QiSGCEaw.net
>>162
簡単なプログラムで３D動画作成できない尿瓶チンパフェチがいるんだよ。
未だに動画がアップされない。俺は慣れたRで作図している。
169:１３２人目の素数さん
24/08/25 14:38:36.89 1ZiEG71e.net
>>168
簡単な5択すらどれかを答えることができず選択肢にない答えを出すレベチのアホが何言ってんの？w
170:１３２人目の素数さん
24/08/25 15:01:14.93 /qrXHaIo.net
>>169
問題読むのも面倒だからChatGPTに投げたら正解が返ってきたぞ。
ゴルゴの問題は誤答が返ってきた。
こういうもんだいが面白いね。
171:１３２人目の素数さん
24/08/25 15:04:46.24 /qrXHaIo.net
経管栄養している患者に投与したら便が赤くなる薬剤は
とChatGPTに投げたけど、正解は返ってこなかったな。
まあ、ChatGPTには臨床経験がないから仕方ないだろうな。
172:１３２人目の素数さん
24/08/25 15:14:36.19 1ZiEG71e.net
>>170
え？完全に後出しジャンケンだったのにまだそんなこと言ってるの？
脳内医者だからchatGPTに頼らないと答えが出ないわけ？
173:１３２人目の素数さん
24/08/25 15:17:59.45 1ZiEG71e.net
chatGPTに投げるならすぐ出てくるだろ
何で答えが出てくるまでダンマリ決め込んでたの？chatGPTにすぐ投げられることもできないくらい無能なの？
174:１３２人目の素数さん
24/08/25 15:23:41.99 UIotCtA+.net
>>171
ChatGPTに質問したら答えいっぱい出てきましたが
答えが出ないとか意味不明ですよ偽医者さん
文字が読めないとか？
175:１３２人目の素数さん
24/08/25 16:18:49.45 1ZiEG71e.net
日本語不自由だから仕方ないねw
176:１３２人目の素数さん
24/08/25 16:33:49.03 fxoV5jAX.net
>>167
書けないんじゃなくこのスレで回答したらスレ違い
プログラム板の然るべきスレで同じ質問したら答えてやるよ
177:１３２人目の素数さん
24/08/25 16:40:08.05 1ZiEG71e.net
>>176
それが理解できないアホだからここで発狂してんだよ
躾ければ決まったところでトイレができる犬以下
178:１３２人目の素数さん
24/08/26 00:45:30.63 oKRWYGRE.net
結局またダンマリかよ
プログラム板に書き込んだところで結局数学板みたいにバカにされて終わりそうだけどw
179:１３２人目の素数さん
24/08/26 08:00:26.74 eNo86RHt.net
>>167
そりゃいるよ
高校数学の質問スレで出題だのプログラミング言語だのやり始める奴すらいるんだから
180:１３２人目の素数さん
24/08/27 06:18:32.52 uadWa+Us.net
尿瓶ジジイまた措置入院かよ？
181:１３２人目の素数さん
24/08/27 09:40:02.23 fklJzjU+.net
NHKでも紹介されたアプリ【TEMU】が
最大5000円分のPayPayを貰えるキャンペーン実施中！
URLﾘﾝｸ(i.ibb.co)
URLﾘﾝｸ(temu.com)
1.リンク先からダウンロード
2.インストール後、リンク先からTEMUを開く
3.メールやGoogle、Appleアカウント等で登録
4.最大7商品無料ルーレットを回す（今回のキャンペーンとは直接関係ありません）
5.友人や家族など招待
6.PayPay付与
キャンペーンの終了時期が不明なのでお早めに！
182:１３２人目の素数さん
24/08/27 09:47:59.56 qmpn3/6k.net
↓自演の知らなかった
183:１３２人目の素数さん
24/08/27 10:22:06.00 0dFnTW4+.net
>>181
もう弾は尽きたよ
184:１３２人目の素数さん
24/08/28 03:21:35.87 oo0Tk9dt.net
>>181
>>183
嶋田信幸49歳
ポイ活必死ですね 馬鹿ですね
185:１３２人目の素数さん
24/08/28 08:21:40.99 WXncGf8t.net
この数値が算出できないってことは、
尿瓶チンパフェチのPhimoseくんは、Wolfram言語が使えないってことだろうなぁ
186:１３２人目の素数さん
24/08/28 08:21:49.72 WXncGf8t.net
練習問題
１５人のラグビーチームでジャンケンをしてキャプテンを１人決めることになった。
１５人全員でジャンケンを始めて勝者どおしでジャンケンを続けて勝者が１人になるまで続ける。
アイコも１回と数えるときキャプテンが決まるまでのジャンケンの回数の中央値は１１４回である。
問題 ４０人のクラスでジャンケンをして学級委員を１人決めるときのジャンケンの回数の中央値を算出せよ。
187:１３２人目の素数さん
24/08/28 08:24:28.12 5uFzq8dK.net
どおしｗ
188:１３２人目の素数さん
24/08/28 11:24:34.09 jfJE3lq8.net
>>167
それが何？
つまらんものはつまらん。
算数で足し算ても検証してろwww
189:１３２人目の素数さん
24/08/28 16:20:05.31 tZ6Zn8cN.net
>>185
スレタイすら読めないアホが高校生にバカにされてるみたいだね
190:１３２人目の素数さん
24/08/28 17:54:35.49 5uFzq8dK.net
>>167
その方がつまらないと思うがな
公式に沿って電卓叩くだけだろ？
ああ成る程、だから証明や幾何が苦手なのか
191:１３２人目の素数さん
24/08/28 18:57:30.82 B5l1srpB.net
>>190
>>167
高校数学の基礎の基礎の証明すらできないチンパンが電卓叩いてキーキー喜んでるだけ
192:１３２人目の素数さん
24/08/29 17:42:28.62 Mdm5M0Mq.net
>>190
算出するための関数を作るのが楽しいんだね。
それをシミュレーションで検証。
193:１３２人目の素数さん
24/08/29 17:48:56.83 Mdm5M0Mq.net
課題：算出する公式とシミュレーションコードを投稿せよ。
ゴルゴ13は100発100中
ゴルゴ14は10発10中
ゴルゴ15は1発1中
とする。
各々10000発撃ったとき各ゴルゴの命中数の期待値はいくらか？
194:１３２人目の素数さん
24/08/29 18:33:31.32 P92nUjXQ.net
作ってる函数に何の工夫もなく、
公式を特定の言語に変換してるだけじゃん
それも数学に特化した言語とも言えるWolfram
面白いかどうかは主観だから好きにすれば良いけど、
客観的に見たら車輪の再発明を繰り返す非生産的行為だよ
195:１３２人目の素数さん
24/08/29 18:52:28.61 BNc6+SX2.net
>>192
スレタイ読めないアホは予選敗退でーす
196:１３２人目の素数さん
24/08/29 20:53:42.92 akSZPtmj.net
最近調子良かったから油断してタイムカード打ってるけど月に8/15からは漏れる疾患で急病かもしれんしな
197:１３２人目の素数さん
24/08/29 20:57:19.42 kpTKB/Z2.net
>>54
はえー
正直
もともとが
リクライニングがかなりフラットに見ることができないし
198:１３２人目の素数さん
24/08/29 21:09:07.43 YpVDF540.net
>>182
便利なカードならわざとサロンのドメインで仮装通貨の買いも異常無しって分かるの？大河より面白かったけど
やっぱ酢の味して不味かったからな
定期的にラッパーの枠にはまりきってないはずの愛想や見た目で言えばかなり理想に近い
199:１３２人目の素数さん
24/08/29 21:40:47.30 QNhyARn/.net
結局他人がどうこう言うことあるよね
200:１３２人目の素数さん
24/08/29 22:26:47.81 9OgyWxAu.net
>>192
作るwww
今までで一つでも作った気になってんの？
東大合格者にあるまじき感想wwww
作ったうちに入らねぇだろwwww
201:１３２人目の素数さん
24/08/29 22:34:47.73 NrDFyX72.net
アンチ風囲いて
なかなか浸透してますよほんま好き
202:１３２人目の素数さん
24/08/29 22:36:39.39 9hsjz4CP.net
ホリエモンも全身脱毛して
炭水化物も食ってもピンとはこないと思うな
203:１３２人目の素数さん
24/08/29 22:40:03.51 bTW+UREj.net
こっちはただの神様のプレゼントなんだな
俺は
204:１３２人目の素数さん
24/08/29 23:25:32.93 G/SrNAYO.net
をつしわほうるえれんさもにろて
205:１３２人目の素数さん
24/08/29 23:36:11.91 Gj+8esKt.net
東方ボーカル界隈では悪用される可能性を危惧したらとか疲労で調子がいい今の所良い話題じゃない
206:１３２人目の素数さん
24/08/29 23:36:16.57 m5RfJAWC.net
あおぞらスケベした情報しか出て欲しい
エバランス😲😲😲😲😲
コカインとかタバコと変わらんぞ
友達がいない31日まで後、
207:１３２人目の素数さん
24/08/29 23:50:53.63 8rvI0hWY.net
問い合わせボタンない
・サロンは毎日のような勤勉さもないから連覇させてくれてるやん
208:１３２人目の素数さん
24/08/30 18:27:34.04 nhc1fobJ.net
電卓おじいさん、小中学生スレの鶴亀算を解くのに最小公倍数を求める関数を使ってイキるｗ
209:１３２人目の素数さん
24/08/31 05:58:25.26 Q/Lo9wJj.net
>>208
GCDが最小公倍数だと思うようなのは
自分でプログラムもできないのだろうなぁ。
東大卒でGCDとLCMの区別ができないとは思えないから、
シリツ卒と推定される。
210:１３２人目の素数さん
24/08/31 05:59:17.10 Q/Lo9wJj.net
プログラム言語の扱える東大卒用の練習問題
課題：算出する公式とシミュレーションコードを投稿せよ。
ゴルゴ13は100発100中
ゴルゴ14は10発10中
ゴルゴ15は1発1中
とする。
各々10000発撃ったとき各ゴルゴの命中数の期待値はいくらか？
211:１３２人目の素数さん
24/08/31 06:01:27.67 Q/Lo9wJj.net
シリツ卒用の問題
GCDはGreatest Common Divisorのことである。
日本語に訳せ
212:１３２人目の素数さん
24/08/31 09:48:54.74 eB+VMmt6.net
>>211
gas chromatographic detectorのことなんだよな～
そんなんもわからないなんて私立文系？
213:１３２人目の素数さん
24/08/31 10:17:02.47 OqrM+pGd.net
そもそも小学生にすら相手にされてない電卓チンパンだしなぁ
214:１３２人目の素数さん
24/09/01 01:43:13.93 7TlKBa7A.net
>>210
いつまで同じ問題wwww出してんの？
日本語も数学も出来なくてまともな問題として設定出来てないの分からないのかよwwww
215:１３２人目の素数さん
24/09/08 06:11:56.44 bna8QGWY.net
>>214
答がでるまで。
で答は？
216:１３２人目の素数さん
24/09/08 18:01:00.78 xdCBvRpf.net
>>215
問題設定すら出来てないのに答えとかwww
馬鹿すぎだろwwww
217:１３２人目の素数さん
24/09/08 18:01:44.35 xdCBvRpf.net
東大卒教えてクレクレ乞食www
218:１３２人目の素数さん
24/09/08 19:09:32.57 ZVAVt1HD.net
そもそも出題スレじゃないと何回書かれたら理解できるんだろう
東大卒、あるいは医者が理解できないはずはないから、
出題繰り返してる人は東大（どころか底辺私大含む大学）卒業すらしておらず、
当然医者でもないな
このスレに書く資格ない人物だろう
219:１３２人目の素数さん
24/09/10 11:12:07.75 whLItct/.net
累積密度関数の逆関数を与えてHighest Density Intervalを算出する
HDI=\(InvCDF=qbeta,cred=0.95,...){
opt=optimize(\(p) InvCDF(p+cred,...) - InvCDF(p,...),c(0,1-cred))
lwr=InvCDF(opt$min,...)
upr=lwr+opt$obj
c(lwr,upr)
}
ガンマ分布で検証
lu=HDI(qgamma,shape=20,rate=24)
lu
curve(dgamma(x,shape=20,rate=24),0,2)
pgamma(lu[2],shape=20,rate=24)-pgamma(lu[1],shape=20,rate=24)
diff(lu) < qgamma(0.975,20,24)-qgamma(0.025,20,24)
ベータ分布で検証
lu=HDI(qbeta,shape1=20,shape2=24)
lu
curve(dbeta(x,20,24),0,1)
pbeta(lu[2],20,24)-pbeta(lu[1],20,24)
diff(lu) < qbeta(0.975,20,24)-qbeta(0.025,20,24)
220:１３２人目の素数さん
24/09/11 23:03:13.05 cP+fZTFu.net
この数学教師マジか…
元小樽桜陽高校数学科教師堤伸弘の数学教師としての能力
URLﾘﾝｸ(note.com)
221:１３２人目の素数さん
24/09/12 12:35:05.85 3Gkbberk.net
>>220
覚えなくても困らないなら、そんなテスト低い点取ってほっとけば良いだけだろ。
足し算すれば答え出るのだから、掛け算九九すら覚えなさそうwww
222:１３２人目の素数さん
24/09/14 08:10:13.90 dkeJyU09.net
出題と質問の違いすら分からないアホチンパンはここに書き込む資格ないからとっとと失せなさい
223:１３２人目の素数さん
24/09/14 08:10:21.26 dkeJyU09.net
>>219
224:１３２人目の素数さん
24/09/17 07:36:31.62 X3jkfzLo.net
vonNeuman <- function(PDF,xmin=0,xmax=1){
N=1e6
ymax=max(PDF(seq(xmin,xmax,length=N+1)))
Ux=runif(N,xmin,xmax)
Uy=runif(N,0,ymax)
Rand=Ux[which(Uy<=PDF(Ux))]
hist(Rand,xlim=c(xmin,xmax),freq=FALSE,breaks=30,col=sample(colors(),1),main='')
curve(PDF,add=TRUE,lwd=2)
invisible(Rand)
}
vonNeuman(dnorm,-3,3)
vonNeuman(\(x)sin(x)/2,0,pi)
225:１３２人目の素数さん
24/09/21 07:06:35.21 cdm6DP8+.net
library(binom)
ci=binom.confint(324-300,324)
lu=unlist(ci[11,5:6])
LearnBayes::beta.select(list(p=0.025,x=lu[1]),list(p=0.975,x=lu[2]))
# 信頼区間からβ分布の形状係数を算出し代表値を返す
ci2ab=\(l,u,verbose=FALSE,cl=0.95){ # CI : [l,u], cl : confidence level
if(l==u) return(NA)
options(warn = -1)
HDI=\(InvCDF=qbeta,cred=0.95,...){
opt=optimize(\(p) InvCDF(p+cred,...) - InvCDF(p,...),c(0,1-cred))
lwr=InvCDF(opt$min,...)
upr=lwr+opt$obj
c(lwr,upr)
}
f=\(ab){
LU=HDI(qbeta,cred=cl,shape1=ab[1],shape2=ab[2])
(LU[1]-l)^2 + (LU[2]-u)^2
}
opt=optim(runif(2,1,100),f)
opt=optim(opt$par,f)
par=opt$par
lu=HDI(qbeta,cred=cl,shape1=par[1],shape2=par[2])
if(verbose){
mean=par[1]/sum(par)
median=qbeta(0.50,par[1],par[2])
mode=(par[1]-1)/(sum(par)-2)
cat('α =',round(par[1],2),' β =',round(par[2],2),'\n')
cat('mean =',round(mean,3),' median =',round(median,3))
if(par[1]>1 & par[2]>1) cat(' mode =',round(mode,3))
cat('\nlower =',round(lu[1],3),' upper =',round(lu[2],3),'\n') curve(dbeta(x,par[1],par[2]),type='h',col=2,n=250,bty='l',ann=FALSE,axes=FALSE)
axis(1)
}
options(warn = 0)
invisible(par)
}
ab=ci2ab(lu[1],lu[2])
k=1e5
p=rbeta(k,ab[1],ab[2])
# 検査陽性の事後確率
postp=\(p,s,t) p*s/ (1-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度
# 検査陰性の事後確率
postn=\(p,s,t) p*(s-1)/(-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度
# 尿素呼気試験(感度90-100% 特異度80-99%) 便中ピロリ菌抗原 (感度90-98% 特異度87-100%)
abs=ci2ab(0.90,1.00)
abt=ci2ab(0.80,0.99)
s=rbeta(k,abs[1],abs[2])
t=rbeta(k,abt[1],abt[2])
post1=postn(p,s,t)
abs=ci2ab(0.90,0.98)
abt=ci2ab(0.87,1.00)
s=rbeta(k,abs[1],abs[2])
t=rbeta(k,abt[1],abt[2])
post2=postn(post1,s,t)
1/mean(post2)
1/median(post2)
hist(post2,freq=FALSE,breaks='scott',ann=F,axes=F) ; axis(1)
226:１３２人目の素数さん
24/09/26 06:12:41.58 DjIa9Rog.net
a個入りのタコ焼き(S)、b個入りのタコ焼き(M)、c個入りのタコ焼き(L)を
S,M,Lを各々s箱,m箱,l箱ずつ用意した。
タコ焼きの数で注文を受けるが、バラ売りはしない。
問題
(1) 受注できるのタコ焼きの数は何種類あるか。但し、0個は受注に数えない。
(2) 出荷できる方法が最も多いのは何個のタコ焼きを受注したときか？
複数あればすべて列挙せよ
Wolfram言語
tako[a_,b_,c_,s_,m_,l_]:=(
S=a Range[0,s];
M=b Range[0,m];
L=c Range[0,l];
t1=DeleteCases[Tuples[{S,M,L}],{0,0,0}];
t2=Total /@ t1;
ans1=Length@Union@t2;
t3=Counts[t2];
ans2=Select[t3,# == Max[t3]&];
{ans1,ans2}
)
tako[5,9,12,100,50,20]
R言語
tako=\(a,b,c,s,m,l){
S=a*(0:s)
M=b*(0:m)
L=c*(0:l)
t1=expand.grid(S,M,L)[-1,]
t2=rowSums(t1)
ans1=length(unique(t2))
hist(t2,breaks=length(unique(t2)),border = 'pink',col=2,main='shipment')
t3=table(t2)
ans2=t3[t3==max(t3)]
list(ans1,ans2)
}
tako(a=5,b=9,c=12,s=100,m=50,l=20)
227:１３２人目の素数さん
24/09/28 08:42:02.91 XCswu5aK.net
a,b,cは100以下の相異なる正の整数でa<b<cとする。
4つの整数a+b-c,b+c-a,c+a-b,a+b+cのすべてが素数となるような組合せは何個あるか数えよ。
rm(list=ls())
library(numbers)
n=100
ans=NULL
for(a in 1:(n-2)){
for(b in (a+1):(n-1)){
for(c in (b+1):n){
x= a+b-c
y= b+c-a
z= c+a-b
if(x>0 & y>0 & z>0){
if(isPrime(x) & isPrime(y) & isPrime(z)){
ans=rbind(ans,c(a,b,c,x,y,z))
}
}
}
}
}
ans
colnames(ans)=c('a','b','c','a+b-c','b+c-a','c+a-b')
head(ans)
tail(ans)
abc=ans[,1:3]
xyz=ans[,4:6]
abc[rowSums(abc)==max(rowSums(abc)),]
abc[rowSums(abc)==min(rowSums(abc)),]
ans4=abc[isPrime(rowSums(abc)),]
head(ans4)
tail(ans4)
228:１３２人目の素数さん
24/09/28 11:00:57.10 4EeustNs.net
へったくそwww
229:１３２人目の素数さん
24/09/28 11:06:24.21 1u+u54qN.net
尿瓶ジジイ都合の悪いレスに噛みついたところで速攻で論派されるのでここでコソコソ書き込むしかできない模様
230:１３２人目の素数さん
24/10/06 10:12:15.79 9nVPq2f1.net
>>227
もう息絶えたのか？
231:１３２人目の素数さん
24/10/06 17:03:12.61 4DgSzUxW.net
お金には使用できる枚数の制限があるのですか
【答】
日本銀行券（いわゆる紙幣、お札）は、「日本銀行法」第46条第２項で「無制限に通用する」と規定されています。
貨幣（いわゆる硬貨）は、「通貨の単位及び貨幣の発行等に関する法律」第７条で「額面価格の20倍まで」を限度として通用すると規定されています。つまり、20枚までは貨幣による支払いを行っても良いということです。
URLﾘﾝｸ(www.mof.go.jp)
１円玉から５００円玉までの硬貨が各々２０枚ずつある。
(1)お釣りをもらわずに支払額は何種類あるか？
(2)支払い方法が最も多いのは何円の支払いのときか。複数あればすべて列挙せよ。
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= y1=1 Range[0,20];
In[2]:= y5=5 Range[0,20];
In[3]:= y10=10 Range[0,20];
In[4]:= y50=50 Range[0,20];
In[5]:= y100=100 Range[0,20];
In[6]:= y500=500 Range[0,20];
In[7]:=
In[7]:= t1=Tuples[{y1,y5,y10,y50,y100,y500}];
In[8]:= t2=Select[Total /@ t1,#!=0&];
In[9]:= Length@Union@t2
Out[9]= 13320
In[10]:= t3=Counts[t2];
In[11]:= Max[t3]
Out[11]= 9867
In[12]:= Position[t3,Max[t3]]
Out[12]= {{Key[3200]}, {Key[3700]}, {Key[4200]}, {Key[4700]}, {Key[5200]}, {Key[5700]},
> {Key[6200]}, {Key[6700]}, {Key[7200]}, {Key[7700]}, {Key[8200]}, {Key[8700]}, {Key[9200]},
> {Key[9700]}, {Key[10200]}, {Key[3120]}, {Key[3620]}, {Key[4120]}, {Key[4620]}, {Key[5120]},
> {Key[5620]}, {Key[6120]}, {Key[6620]}, {Key[7120]}, {Key[7620]}, {Key[8120]}, {Key[8620]},
> {Key[9120]}, {Key[9620]}, {Key[10120]}}
In[13]:= % // Length
Out[13]= 30
232:１３２人目の素数さん
24/10/13 11:47:25.71 x0yyIaSy.net
f=\(n){
p=numeric()
for (a in 1:6) p[a]=sum((1/6)^(0:(n-1))*(a-1)/6)
mean(p)
}
n=1:20
y=sapply(n,f)
plot(n,y,bty='l',pch=16)
233:１３２人目の素数さん
24/10/24 16:47:36.34 pQNl9nxD.net
(*
血液型 AB型10人、B型20人、O型30人、A型40人の100人から無作為に別人を4人選ぶとき
選ばれた4人の血液型がすべて異なる確率を算出せよ。
*)
a=Flatten@{Table[1,10],Table[2,20],Table[3,30],Table[4,40]};
N@Mean@Table[Boole[Length@Union@RandomSample[a,4]==4],1*^6]
234:１３２人目の素数さん
24/10/24 16:59:40.70 o7cQL/rk.net
>>233
相変わらず誰にも相手にされてないみたいだね
235:１３２人目の素数さん
24/10/24 17:00:49.50 hfz9xw4M.net
自分が解けるようになるまで問題を出さないスタイル
だからもちろん他者の問題に答えることもしない
236:１３２人目の素数さん
24/10/25 15:17:05.37 AKAjfr56.net
ここは兄弟スレのレポジトリとして使用。
237:１３２人目の素数さん
24/10/25 15:59:33.25 t22ya7mF.net
>>236
尿瓶ジジイまだ生きてたのかよ？死に損ないもいいとこだな
238:１３２人目の素数さん
24/10/25 19:11:16.83 AKAjfr56.net
１月から内視鏡バイトの勤務日をもう1日増やしてくれと打診された。
スタッフが優秀でストレスのない職場なので、次の医師がみつかるまで
という条件で引き受けた。
まあ、次の医師がみつかるのは新年度だろうな。
仕事をふやしても年金支給停止額が増えるだけ。
239:１３２人目の素数さん
24/10/25 19:27:41.81 tA8uN3sd.net
年金支給停止額が増えるなら喜んでやるだろ
やらない理由がなくなるんだから日本語くらい正しく使えよ
数学もまともにできないんだからさ
240:１３２人目の素数さん
24/10/26 09:11:01.76 DYbBtbTY.net
"
ある政党に100人の議員がいる。
個々の議員が裏金議員である可能性には何の情報もないためその確率を一様分布と仮定する。
無作為に10人を選んで調べたところ9人が裏金議員であった。
100人中の裏金議員の数の期待値と95%信頼区間を算出せよ。
算出法は好みの手法でよい。
"
rm(list=ls())
n100=100
n10=10
n9=9
m=n100/2
sd=sqrt(n100/12)
fp=\(n) pnorm(n,m,sd)-pnorm(n-1,m,sd)
pn=sapply(0:n100,fp)
sim=\(){
u=sample(0:n100,1,prob=pn)
u10=sum(sample(1:0,n10,replace=TRUE,prob=c(u,100-u)))
c(u,u10)
}
k=1e6
res=t(replicate(k,sim()))
ans=res[res[,2]==n9,][,1]
BEST::plotPost(ans)
241:１３２人目の素数さん
24/10/26 11:11:45.76 DYbBtbTY.net
中心極限定理を使わずに算出
rm(list=ls())
n100=100
n10=10
n9=9
sim=\(){
u=sum(runif(n100))
u10=sum(sample(0:1,n10,replace=TRUE,prob = c(n100-u,u)))
c(u,u10)
}
k=1e6
res=t(replicate(k,sim()))
ans=res[res[,2]==n9,][,1]
hist(ans,freq=F,col=2,ann=F,axes=F) ; axis(1)
summary(ans)
HDInterval::hdi(ans)
BEST::plotPost(ans,xlab='n',showMode = T)
242:１３２人目の素数さん
24/10/26 18:10:39.10 QUufrNke.net
n100=100;
n10=10;
n9=9;
sim[] :=(
u=Total@RandomVariate[UniformDistribution[{0,1}],100];
u10=Total@RandomChoice[{u,100-u}->{1,0},10];
{u,u10}
)
k=1*^7;
res=Select[Table[sim[],k],#[[2]]==n9&];
ans=#[[1]]& /@ res;
Histogram[ans]
Mean[ans]
Median[ans]
243:１３２人目の素数さん
24/10/26 18:22:42.42 V7l7LhBL.net
>>242
高校生にすら相手にされてなくて哀れだね
244:１３２人目の素数さん
24/10/26 21:19:37.54 dZVTecny.net
n100=100
n10=10
n9=9
# P[9|m]
p9_m=\(m) choose(m,n9)*choose(n100-m,n10-n9)/choose(100,10)
p9_m=Vectorize(p9_m)
P9m=p9_m(1:100)
plot(1:100,P9m)
# P[m|9]=P[9|m]P[m]/P[9]
# P[9]=P[9|m]P[m]+P[9|!m]P[!m]
# P[m]=1/101
# P[!m]=100/101
# P[9|!m] = P9m[-m]
pm_9=\(m) p9_m(m)/101 / (p9_m(m)/101+sum(P9m[-m])*100/101)
pm_9=Vectorize(pm_9)
auc=sum(pm_9(1:100))
x=1:100
sum(x*pm_9(x)/auc)
plot(x,pm_9(x))
245:１３２人目の素数さん
24/10/26 21:34:59.82 dZVTecny.net
n100=100
n10=10
n9=9
sim=\(){
u=sample(0:n100,1)
u10=sum(sample(1:0,n10,replace=TRUE,prob=c(u,n100-u)))
c(u,u10)
}
k=1e5
res=t(replicate(k,sim()))
ans=res[res[,2]==n9,][,1]
hist(ans)
summary(ans)
246:１３２人目の素数さん
24/10/26 22:09:07.09 0cRJo0MK.net
あ、プログラミング出来る人だ丁度良かった。
これの解説できません？
プログラムを組むのに数学の勉強は必要か？数学板住民の意見は
ｽﾚﾘﾝｸ(math板:14番)
14: １３２人目の素数さん 2024/10/26(土) 20:14:07.62 ID:0cRJo0MK
でも数学の知識があるとコードを短くできる場面もあるっぽい。
自分が見たのはじゃんけんの手が群をなす？とかで条件分岐を無くしてた。
誰か理解できる人は解説頼む。
247:１３２人目の素数さん
24/10/26 22:38:17.70 dZVTecny.net
"
ある政党に100人の議員がいる。
何人が裏金議員であるには何の情報もないため
0～100人である確率は同じと家庭する。
即ち、0人である確率も99人である確率も1/101とする。
無作為に10人を選んで調べたところ9人が裏金議員であった。
100人中の裏金議員の数の期待値と95%信頼区間を算出せよ。
算出法は好みの手法でよい。
"
rm(list=ls())
par(bty='l')
n100=100
n10=10
n9=9
# P[9|m]
p9_m=\(m) choose(m,n9)*choose(n100-m,n10-n9)/choose(100,10)
p9_m=Vectorize(p9_m)
P9m=p9_m(0:100)
plot(0:100,P9m,type='h')
# P[m|9]=P[9|m]P[m]/P[9]
# P[9]=sum(P9m)
# P[m]=1/101
pm_9=\(m) p9_m(m)/101 / sum(P9m)
pm_9=Vectorize(pm_9)
auc=sum(pm_9(0:100))
x=0:100
sum(x*pm_9(x)/auc)
plot(x,pm_9(x),type='h',lwd=2,col=3,axes=F,ann=F) ; axis(1)
plot(0:100,cumsum(pm_9(x)/auc),type='l')
abline(h=c(0.025,0.975),lty=3)
x100=1:100
p=pm_9(x100)/auc
po=order(p,decreasing = TRUE)
sum(cumsum(sort(p,decreasing = TRUE))<0.95) # 34
po[1:34]
sum(p[po[1:34]])
po[1:35]
sum(p[po[1:35]])
248:１３２人目の素数さん
24/10/26 22:39:35.12 V7l7LhBL.net
>>247
>>246には答えられないみたいだね統計もどきだからw
249:１３２人目の素数さん
24/10/27 05:57:28.52 SNtImSnv.net
(*
ある政党に100人の議員がいる。
何人が裏金議員であるには何の情報もないため
0～100人である確率は同じと家庭する。
即ち、0人である確率も99人である確率も1/101とする。
無作為に10人を選んで調べたところ9人が裏金議員であった。
100人中の裏金議員の数の期待値と95%信頼区間を算出せよ。
*)
n100=100;
n10=10;
n9=9;
(* P[9|m] *)
p9m[m_] := Binomial[m,n9] Binomial[n100-m,n10-n9]/Binomial[100,10]
P9m=Table[p9m[m],{m,0,100}];
(* P[m|9] *)
pm9[m_] := p9m[m]/(n100+1) / Total[P9m];
auc=Total@Table[pm9[m],{m,0,n100}];
pdf[m_]:=pm9[m]/auc
Sum[x*pdf[x],{x,0,n100}]
Table[pdf[x],{x,1,100}] // ListPlot
Plot[pdf[x],{x,0,100}]
p=Table[pdf[x],{x,1,100}];
ps=ReverseSort[p];
Boole[#<0.95]& /@ Accumulate[ps] // Total
i=Reverse[Ordering[p]][[1;;34]];
{Min[i],Max[i]}
p[[i]] // Total // N
j=Reverse[Ordering[p]][[1;;35]];
{Min[j],Max[j]}
p[[j]] // Total // N
250:１３２人目の素数さん
24/10/27 07:06:45.12 SNtImSnv.net
(*
あるシリツ医大にｘ人が入学したとする。
何人が裏口入学であるには何の情報もないため一様分布とする。
即ち0～x人である確率はどれも同じと仮定する。
無作為にy人を選んで調べたところz人が裏口入学であった。
x人中の裏口入学者の数の期待値とその95%信頼区間、最頻値、中央値を算出するソルバーをWolfram言語で作成せよ。
*)
solve[x_,y_,z_] :=(
n100=x;
n10=y;
n9=z;
(* P[yz|x] *)
p9m[m_] := Binomial[m,n9] Binomial[n100-m,n10-n9]/Binomial[n100,n10];
P9m=Table[p9m[m],{m,0,n100}];
(* P[x|yz] *)
pm9[m_] := p9m[m]/(n100+1) / Total[P9m];
auc=Total@Table[pm9[m],{m,0,n100}];
pdf[m_]:=pm9[m]/auc;
e=Sum[m*pdf[m],{m,0,n100}];
p=Table[pdf[m],{m,1,n100}];
ps=ReverseSort[p];
n34=Total[Boole[#<0.95]& /@ Accumulate[ps]];
i=Reverse[Ordering[p]][[1;;n34]];
{Min[i],Max[i]};
p[[i]] // Total // N;
j=Reverse[Ordering[p]][[1;;n34+1]];
ci={Min[j],Max[j]};
p[[j]] // Total // N;
mode=Position[p,Max[p]][[1]][[1]];
median=Total[Boole[#<0.5]& /@ Accumulate[p]];
{Ex->e,CI->ci,Mode->mode,Median->median}
)
solve[100,10,9]
251:１３２人目の素数さん
24/10/27 08:42:49.81 Ll/mJg+e.net
>>249
死に損ないいつになったら高校生に相手にされるんだ？
252:１３２人目の素数さん
24/10/27 09:42:05.05 HIz5euo6.net
# 一峰性非対称分布する離散量t1の95% highest probability intervalを返す
t2=sort(table(t1),decreasing=TRUE)
t3=as.numeric(names(t2))
t4=t2/sum(t2)
t5=sum(cumsum(t4)<0.95)
range(t3[1:t5])
sum(t2[1:t5]/sum(t2))
range(t3[1:(t5+1)])
sum(t2[1:(t5+1)]/sum(t2))
253:１３２人目の素数さん
24/10/27 09:51:18.11 HIz5euo6.net
>>246
例
複素平面上の点A,B,Cの作る三角形の面積を求める関数をプログラムせよ
行列式を使って1行で完成する
ABCS <- function(A,B,C) abs(det(rbind(c(Re(A-C),Im(A-C)),c(Re(B-C),Im(B-C)))))/2
254:１３２人目の素数さん
24/10/27 09:55:54.15 HIz5euo6.net
数学の知識がプログラム短縮できる例（R言語）
四面体の頂点を与えて四面体の体積を行列式を使って計算
ABCD2V <- function(A,B,C,D){ # 四面体ABCDの体積
v=rbind(A,B,C,D)
abs(det(rbind(v[1,]-v[4,],v[2,]-v[4,],v[3,]-v[4,])))/6
}
255:１３２人目の素数さん
24/10/27 11:03:50.59 1ynMj6oQ.net
>>251
こういう言葉遣いをする人って祖父母や両親から愛情を注がれなかった気の毒な人間なんだろうなぁ。
256:１３２人目の素数さん
24/10/27 11:07:06.01 1ynMj6oQ.net
>>251
女子大生や看護学生と嵌めたことはあるけど、女子高生とは嵌めたことがないなぁ。
ちなみに、女子医大の学生には在京中は息子が大変お世話になりました。
257:１３２人目の素数さん
24/10/27 11:22:36.43 6A9938yt.net
人格上もクズ
258:１３２人目の素数さん
24/10/27 12:07:23.87 Ll/mJg+e.net
>>255
5chですら必死にレス乞食しても高校生にすら相手にされない老害が死に損ない以外の何だって言うんだよw
259:１３２人目の素数さん
24/10/27 12:09:18.71 Ll/mJg+e.net
>>256
妄想かなけなしの金で行ったそういう設定の風俗だろうなwww
数学板の高校生にすら相手にされないんだからリアルで誰にも相手にされるわけないww
260:１３２人目の素数さん
24/10/27 12:36:11.63 1ynMj6oQ.net
>>259
オペナースとかは活動的な女性が多いぞ。
外科医に穴兄弟が多いのは業界人なら知っている。
なぜなら、活動的なナースは口が軽いからｗ
261:１３２人目の素数さん
24/10/27 12:39:15.42 1ynMj6oQ.net
高学歴女医はフェラを嫌うが、新設国立の女医（女子医学生）はその逆。
エビデンスレベルV（個人の体験）
最近の作業仮説
セクハラ認定したがるのはブサイク女か、持てない男である。
262:１３２人目の素数さん
24/10/27 13:01:14.16 Ll/mJg+e.net
>>260
そもそも年中早朝からここで発狂してる無職だろw
263:１３２人目の素数さん
24/10/27 13:13:50.81 q2AYb8Fp.net
職場の労働衛生義務すら知らないってことは少なくとも管理職未満か
264:１３２人目の素数さん
24/10/27 14:30:10.23 Ll/mJg+e.net
さっきまで即レスだったのに>>263が書き込まれた途端死んじゃった
265:１３２人目の素数さん
24/10/27 18:51:50.22 wu6Ollaf.net
>>253
ありがとうございます。
普通だと、底辺を決めて長さと高さを求めて、面積計算…という手順のコードのはずが
１行で!!って感じでしょうか?
detは数学書で読んだ記憶はありますが、面積計算にも生きるんですね。
Ｒ言語?ですよね？
後で調べてみます。
266:１３２人目の素数さん
24/10/27 21:47:42.48 1ynMj6oQ.net
おまけ
外積ベクトルの長さは平行四辺形の面積という知識があれば、
x-y-z座標で三角形A(a1,a2,a3)、B(b1,b2,b3)、C(c1,c2,c3)の面積計算だと
ベクトルA-CとベクトルB-Cの外積ベクトルの長さの1/2で計算すれば、コードが短縮できる。
外積を計算する関数(R言語だとライブラリpracmaのcross)が備わっていれば、それを利用する。
v=pracma::cross(A-C,B-C)
sqrt(sum(v^2))/2
で計算できる。
267:１３２人目の素数さん
24/10/28 05:48:46.65 F6mA4jUa.net
また尿瓶劇場かよ
268:１３２人目の素数さん
24/10/29 07:30:27.03 bTR+YKOB.net
"
餅を1億口食べると7人が窒息死するという
URLﾘﾝｸ(youtu.be)
某AIの答
>>
日本でワクチンが原因と認定された死亡例は、約79件でした。この中で、ワクチン接種後に死亡が確認されたケースは51件でした1。
一方、日本での総ワクチン接種回数は約4億3961万回に上ります。
<<
【問題】
（１）餅による窒息死とワクチン原因死のリスク比（もしくはオッズ比）の95%信頼区間を算出せよ。
【博物館入りの古典問題】
餅による窒息死とワクチン原因死の確率は同じを帰無仮説として有意差検定せよ。
すなわち、帰無仮説を前提に実測値とmore extremeな場合の確率の合計=p値を求めよ。
"
source('toolmini.R')
r1=7
r2=51
n1=1e8
n2=4.3961e8
Fisher.test(c(r1,r2),c(n1,n2))
prop.test(c(r1,r2),c(n1,n2))
hit=c(r1,r2)
shot=c(n1,n2)
mat=cbind(hit,shot-hit)
Epi::twoby2(mat)
mochi=rbeta(1e9,7+0.5,1e8-7+0.5)
vacci=rbeta(1e9,51+0.5,4.3961e-51+0.5)
d=mochi-vacci
hist(d)
summary(d)
HDInterval::hdi(d)
mean(mochi>vacci)
re=beta.diff(r1,r2,n1,n2)
re$cdf(0)
269:１３２人目の素数さん
24/10/29 12:09:44.37 bTR+YKOB.net
>>268
確率が小さいのでシミュレーションで算出するのは困難だな。
1000万回程度だと差がでない。