高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
高校数学の質問スレ（医者・東大卒専用） Part438 http://rio2016.5ch.net/test/read.cgi/math/1723152147/
上
下
前
次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
484: １３２人目の素数さん [sage] 2024/12/29(日) 14:05:39.63 ID:bEHNzlJY 尿瓶ジジイ医者でも東大卒でもないのにこのスレばっか伸びてるのほんと図々しいなw http://rio2016.5ch.net/test/read.cgi/math/1723152147/484
485: １３２人目の素数さん [sage] 2024/12/29(日) 23:58:05.11 ID:/vJ/3Q0O 尿瓶ジジイまたダンマリかよ？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/485
486: １３２人目の素数さん [sage] 2024/12/30(月) 03:14:24.14 ID:DAgtZ3Zg # 遷移行列 P <- matrix(c( 0, 0.5, 0.5, 0, 0, # S0: 初戦 A vs B 0, 0, 0, 0.5, 0.5, # S1: A vs C 0, 0, 0, 0.5, 0.5, # S2: B vs C 0.5, 0, 0, 0, 0.5, # S3: 再試合 A vs B 0.5, 0, 0, 0.5, 0 # S4: 再試合 B vs A ), nrow = 5, byrow = TRUE) # 固有値と固有ベクトルの計算 eig <- eigen(t(P)) stationary <- Re(eig$vectors[, which(abs(Re(eig$values) - 1) < 1e-8)]) stationary <- stationary / sum(stationary)
# 各プレイヤーの勝率を計算 A_win <- stationary[2] + stationary[4] B_win <- stationary[3] + stationary[5] C_win <- stationary[1] cat("定常分布:\n") print(stationary) cat(sprintf("\nAの勝率: %.3f\n", A_win)) cat(sprintf("Bの勝率: %.3f\n", B_win)) cat(sprintf("Cの勝率: %.3f\n", C_win)) http://rio2016.5ch.net/test/read.cgi/math/1723152147/486
487: １３２人目の素数さん [sage] 2024/12/30(月) 03:17:30.79 ID:DAgtZ3Zg (* 遷移行列の定義 *) P = { {0, 0.5, 0.5, 0, 0}, (* S0: 初戦 A vs B *) {0, 0, 0, 0.5, 0.5}, (* S1: A vs C *) {0, 0, 0, 0.5, 0.5}, (* S2: B vs C *) {0.5, 0, 0, 0, 0.5}, (* S3: 再試合 A vs B *) {0.5, 0, 0, 0.5, 0} (* S4: 再試合 B vs A *) }; (* 遷移行列の転置 *) PT = Transpose[P]; (* 固有値と固有ベクトルの計算 *) {eigValues, eigVectors} = Eigensystem[PT]; (* 固有値1に対応する固有ベクトル *) stationaryVector = Select[eigVector
s, Norm[PT.# - #] & /@ eigValues == 1 ][[1]]; (* 正規化（確率の総和を1にする） *) stationaryVector = stationaryVector/Total[stationaryVector]; (* 各プレイヤーの勝率を計算 *) AWin = stationaryVector[[2]] + stationaryVector[[4]]; BWin = stationaryVector[[3]] + stationaryVector[[5]]; CWin = stationaryVector[[1]]; (* 結果表示 *) {stationaryVector, AWin, BWin, CWin} http://rio2016.5ch.net/test/read.cgi/math/1723152147/487
488: １３２人目の素数さん [sage] 2024/12/30(月) 03:31:36.28 ID:DAgtZ3Zg # 遷移行列 P <- matrix(c( 0, 0.5, 0.5, 0, 0, # S0: 初戦 A vs B 0, 0, 0, 0.5, 0.5, # S1: A vs C 0, 0, 0, 0.5, 0.5, # S2: B vs C 0.5, 0, 0, 0, 0.5, # S3: 再試合 A vs B 0.5, 0, 0, 0.5, 0 # S4: 再試合 B vs A ), nrow = 5, byrow = TRUE) # 固有値と固有ベクトルの計算 eig <- eigen(t(P)) stationary <- Re(eig$vectors[, which(abs(Re(eig$values) - 1) < 1e-8)]) stationary <- stationary / sum(stationary)
# 各プレイヤーの勝率を計算 A_win <- stationary[2] + stationary[4] B_win <- stationary[3] + stationary[5] C_win <- stationary[1] cat("定常分布:\n") print(stationary) cat(sprintf("\nAの勝率: %.3f\n", A_win)) cat(sprintf("Bの勝率: %.3f\n", B_win)) cat(sprintf("Cの勝率: %.3f\n", C_win)) http://rio2016.5ch.net/test/read.cgi/math/1723152147/488
489: １３２人目の素数さん [sage] 2024/12/30(月) 03:34:00.29 ID:DAgtZ3Zg (* 遷移行列の定義 *) P = { {0, 0.5, 0.5, 0, 0}, (* S0: 初戦 A vs B *) {0, 0, 0, 0.5, 0.5}, (* S1: A vs C *) {0, 0, 0, 0.5, 0.5}, (* S2: B vs C *) {0.5, 0, 0, 0, 0.5}, (* S3: 再試合 A vs B *) {0.5, 0, 0, 0.5, 0} (* S4: 再試合 B vs A *) }; (* 遷移行列の転置 *) PT = Transpose[P]; (* 固有値と固有ベクトルの計算 *) {eigValues, eigVectors} = Eigensystem[PT]; (* 固有値1に対応する固有ベクトル *) stationaryVector = Select[eigVector
s, Norm[PT.# - #] & /@ eigValues == 1 ][[1]]; (* 正規化（確率の総和を1にする） *) stationaryVector = stationaryVector/Total[stationaryVector]; (* 各プレイヤーの勝率を計算 *) AWin = stationaryVector[[2]] + stationaryVector[[4]]; BWin = stationaryVector[[3]] + stationaryVector[[5]]; CWin = stationaryVector[[1]]; (* 結果表示 *) {stationaryVector, AWin, BWin, CWin} http://rio2016.5ch.net/test/read.cgi/math/1723152147/489
490: １３２人目の素数さん [sage] 2024/12/30(月) 10:41:19.52 ID:dmM3LvbP >>482 自分の問題すら碌に解けないチンパンジーだったみたい http://rio2016.5ch.net/test/read.cgi/math/1723152147/490
491: １３２人目の素数さん [sage] 2024/12/30(月) 21:13:17.75 ID:DAgtZ3Zg >>490 類似問題のコードを弄るだけなのに クレクレ レス乞食かよ。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/491
492: １３２人目の素数さん [sage] 2024/12/30(月) 21:39:27.87 ID:dmM3LvbP >>491 クレクレレス乞食って自分で出した問題すら解けずに逃走するアンタのことだろ http://rio2016.5ch.net/test/read.cgi/math/1723152147/492
493: １３２人目の素数さん [sage] 2024/12/31(火) 07:15:10.51 ID:41bOMap+ (* 応用問題 600×nの三乗根が正整数になるｎの最小値を求めよ 演習問題 1234567890×nの五乗根が正整数になるｎの最小値を求めよ *) solve[x_,m_:3] := Times@@(FactorInteger[x][[All,1]]^Mod[(m-Mod[FactorInteger[x][[All,2]],m]),m]) solve[600,3] solve[1234567890,5] http://rio2016.5ch.net/test/read.cgi/math/1723152147/493
494: １３２人目の素数さん [sage] 2024/12/31(火) 07:15:43.90 ID:41bOMap+ >>492 レス乞食のために答を算出するRのコードを底辺シリツ医スレにアップしておいたから 実行してみ。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/494
495: １３２人目の素数さん [sage] 2024/12/31(火) 10:02:09.35 ID:6B6ig4bh >>494 これ答えじゃないよな？ レス乞食さん http://rio2016.5ch.net/test/read.cgi/math/1723152147/495
496: １３２人目の素数さん [sage] 2024/12/31(火) 14:55:01.89 ID:41bOMap+ EpiCycloid[R_,r_] := Module[ {circleA,circleB,redPoint,redPointTrace,frames,theta}, (* theta:Bの中心の偏角 *) circleA = Circle[{0, 0}, R]; circleB[theta_] := Circle[{(R+r) Cos[theta], (R+r) Sin[theta]}, r]; redPoint[theta_] := {(R+r) Cos[theta] + r Cos[(r+R)/r theta + Pi], (R+r) Sin[theta] + r Sin[(R+r)/r theta + Pi]}; (* 赤点の軌跡を保存するためのリスト *) redPointTrace = Table[redPoint[theta], {theta, 0, 2 Pi, 2 Pi/100}]; (* 円Bの中
心を円Aの円周に沿って移動させるためのアニメーション *) frames = Table[ Show[ Graphics[{EdgeForm[Black], FaceForm[None], circleA}], Graphics[{EdgeForm[None], FaceForm[Red], circleB[theta]}], Graphics[{Red, PointSize[Large], Point[redPoint[theta]]}], Graphics[{Red, Line[redPointTrace[[1 ;; Round[theta/ (2 Pi/100)]]]]}], PlotRange -> {{-2(R+r), 2(R+r)}, {-2(R+r), 2(R+r)}}, Axes -> True, AxesOrigin -> {0, 0} ], {theta, 0, 2 Pi, 2 Pi/100}]; (* アニメーション表示 *) ListAnimate[frames] ] EpiCycloid[3,1]
HypoCycloid[R_,r_] := Module[ {circleA,circleB,redPoint,theta,frames}, (* theta 円Bの中心の偏角 *) circleA = Circle[{0, 0}, R]; circleB[theta_] := Circle[{(R - r) Cos[theta], (R - r) Sin[theta]}, r]; redPoint[theta_] := {(R-r) Cos[theta] + r Cos[(1-R/r) theta], (R-r) Sin[theta] + r Sin[(1-R/r) theta]}; (* 円Bの中心を円Aの円周に沿って移動させるためのアニメーション *) frames = Table[ Show[ Graphics[{EdgeForm[Black], FaceForm[None], circleA}], Graphics[{EdgeForm[None], FaceForm[Red], circleB[theta]}],
Graphics[{Red, PointSize[Large], Point[redPoint[theta]]}], Graphics[{Red, Line[Table[redPoint[theta], {theta, 0, 2 Pi, 2 Pi/100}]]}], PlotRange -> {{-R, R}, {-R, R}}, Axes -> True, AxesOrigin -> {0, 0} ], {theta, 0, 2 Pi, 2 Pi/100} ]; (* アニメーションを表示 *) ListAnimate[frames] ] HypoCycloid[3,1] http://rio2016.5ch.net/test/read.cgi/math/1723152147/496
497: １３２人目の素数さん [sage] 2024/12/31(火) 14:56:07.64 ID:41bOMap+ >>495 実行すれば答がでるだろ。 証言をすべてみたす組み合わせのすべてで嘘つきになるのを列挙すればいい。 Rコード読めないの？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/497
498: １３２人目の素数さん [sage] 2024/12/31(火) 18:47:36.81 ID:B87FjDFd >>497 結局グダグダ言い訳してるが自分の問題すら答えられないアホってことね http://rio2016.5ch.net/test/read.cgi/math/1723152147/498
499: １３２人目の素数さん [sage] 2024/12/31(火) 23:12:32.74 ID:PjhLbT7v >>497 さすが答えの意味分かってないwww http://rio2016.5ch.net/test/read.cgi/math/1723152147/499
500: １３２人目の素数さん [sage] 2025/01/01(水) 08:56:34.48 ID:trf8HhUX ハイポサイクロイド https://i.imgur.com/Y699ttV.gif http://rio2016.5ch.net/test/read.cgi/math/1723152147/500
501: １３２人目の素数さん [sage] 2025/01/01(水) 11:12:05.21 ID:vYaE/f7n >>500 年明けても日本語読めないみたいだな http://rio2016.5ch.net/test/read.cgi/math/1723152147/501
502: １３２人目の素数さん [sage] 2025/01/03(金) 19:11:24.08 ID:KUn71WA0 https://i.imgur.com/PiWG2FN.gif http://rio2016.5ch.net/test/read.cgi/math/1723152147/502
503: １３２人目の素数さん [sage] 2025/01/03(金) 20:33:58.18 ID:gmyxpwxe 高校生にすら相手にされてないみたいだね 実に哀れ http://rio2016.5ch.net/test/read.cgi/math/1723152147/503
504: １３２人目の素数さん [sage] 2025/01/03(金) 22:47:45.54 ID:KUn71WA0 https://i.imgur.com/XfJfU12.gif http://rio2016.5ch.net/test/read.cgi/math/1723152147/504
505: １３２人目の素数さん [sage] 2025/01/03(金) 23:25:35.89 ID:rOKue8hC >>482 荒らすの止めてくれないか http://rio2016.5ch.net/test/read.cgi/math/1723152147/505
506: １３２人目の素数さん [sage] 2025/01/03(金) 23:25:35.96 ID:gmyxpwxe さっさと精神科行け http://rio2016.5ch.net/test/read.cgi/math/1723152147/506
507: １３２人目の素数さん [sage] 2025/01/05(日) 12:23:25.71 ID:tuC7O64g R=3; r=2; K=2; (* ロータ中心と頂点距離/ロータギア半径 *) apex[t_] := Module[ {p,c,d}, p = {(R-r) Cos[t] + R Cos[(t+Pi)/3], (R-r) Sin[t] + R Sin[(t+Pi)/3]}; c = {(R-r) Cos[t],(R-r) Sin[t]}; d = K(p-c)+c ] theta=Flatten@{Range[0,6 Pi, Pi/1000]}; ListPlot[apex[#]& /@ theta,AspectRatio ->1] x[t_]:=apex[t][[1]] y[t_]:=apex[t][[2]] length=Integrate[(D[x[t],t])^2+(D[y[t],t])^2,{t,0,6 Pi}] area=(1/2) Integrate[x[t] D[y[t],t] - y[t] D[x[t]
,t],{t,0, 6 Pi}] Clear[{R,r,K}] http://rio2016.5ch.net/test/read.cgi/math/1723152147/507
508: １３２人目の素数さん [sage] 2025/01/05(日) 12:47:54.51 ID:/2bPdu/v >>507 まだバカにされたりないのか尿瓶ジジイは しかし正規のスレじゃもうバカにすらしてもらえないみたいww http://rio2016.5ch.net/test/read.cgi/math/1723152147/508
509: １３２人目の素数さん [sage] 2025/01/07(火) 19:53:08.02 ID:rJfoFeSd 【問題】 固定円の半径を２、回転ローターの円の半径を３とするときに ルーローの三角の軌跡の内部でローターが移動できるには三角の頂点とロータ円の中心の距離が 回転ローターの円の半径の何倍以上であることが必要ですか？ > 2/3*(2+sqrt(3)) [1] 2.488034 http://rio2016.5ch.net/test/read.cgi/math/1723152147/509
510: １３２人目の素数さん [sage] 2025/01/07(火) 22:35:53.49 ID:zsZUtuIb >>509 高校数学で回転ローターとやらが定義されてるのかよwww 中卒以下の低脳でまともに数学用語使えないゴミなのはしょうがないけど、もう少し努力しろよwww http://rio2016.5ch.net/test/read.cgi/math/1723152147/510
511: １３２人目の素数さん [sage] 2025/01/08(水) 02:59:50.29 ID:0uIf2Djd ルーローって魯肉のことだな http://rio2016.5ch.net/test/read.cgi/math/1723152147/511
512: １３２人目の素数さん [sage] 2025/01/08(水) 04:11:37.72 ID:X/URcmTf rm(list=ls()) JugRiddle <- function( big=5, small=3, end=4){ # starting from the bigger jug movebig <- function(xy){ # start from c(big,0) x=xy[1] ; y=xy[2] # x==big if(x==big) re=c(big-(small-y),small) # x==0 if(x==0) re=c(big,y) # y==small if(y==small) re=c(x,0) # y==0 if(y==0 & x!=big){ if(x>=small) re=c(x-small,small) else re=c(0,x) } return(re) } STATUS=status=c(big,0) i=1 while(!identical(status,c(0,0))){# i=i+1 statu
s=movebig(status) STATUS=rbind(STATUS,status) } rownames(STATUS)=1:nrow(STATUS) colnames(STATUS)=c(paste0(big,'L'),paste0(small,'L')) (Bigger=STATUS) # starting from the smaller jug movesmall <- function(xy){ # start from c(0,small) x=xy[1] ; y=xy[2] if(y==small){ if(x<=(big-small)) re=c(x+small,0) else re=c(big, small-(big-x)) } if(y==0) re=c(x,small) if(x==big) re=c(0,y) if(x==0) re=c(y,0) return(re) } STATUS=status=c(0,small) i=1 while(!identical(status,c(0,0))){ # i=i+1 status=movesmall(status)
STATUS=rbind(STATUS,status) } rownames(STATUS)=1:nrow(STATUS) colnames(STATUS)=c(paste0(big,'L'),paste0(small,'L')) (Smaller=STATUS) if(all(end != c(Bigger,Smaller))) return(NA) min_Bigger=min(which(apply(Bigger,1,function(x) end %in% x))) min_Smaller=min(which(apply(Smaller,1,function(x) end %in% x))) list(Bigger=as.matrix(Bigger)[1:min_Bigger,], Smaller=as.matrix(Smaller)[1:min_Smaller,], min_Bigger=min_Bigger,min_Smaller=min_Smaller) } JugRiddle(5,3,4) JugRiddle(10,5,3) http://rio2016.5ch.net/test/re
ad.cgi/math/1723152147/512
513: １３２人目の素数さん [sage] 2025/01/08(水) 19:26:21.90 ID:SKdwPI7I コイントスの終了条件を Aくんは、 表が2回連続して出たら終了 Bくんは、 表裏の順に出たら終了（裏表の順では終了しない） することになった。 (1) Aくん、Bくんについて、コイントスの回数の期待値を求めよ。 (2) Aくんはの方が少ない回数で終了する確率を求めよ。小数解でよい。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/513
514: １３２人目の素数さん [sage] 2025/01/08(水) 19:29:44.12 ID:X/URcmTf >>513 （２）回数の差の分布 https://i.imgur.com/jDZubdl.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/514
515: １３２人目の素数さん [sage] 2025/01/08(水) 20:16:25.91 ID:X/URcmTf >>514 https://i.imgur.com/cmdTCLY.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/515
516: １３２人目の素数さん [sage] 2025/01/08(水) 20:27:25.85 ID:bjQ79ehj 日本語すら読めないチンパンジーが数学やった気になってるのか http://rio2016.5ch.net/test/read.cgi/math/1723152147/516
517: １３２人目の素数さん [sage] 2025/01/08(水) 20:27:57.47 ID:SKdwPI7I >>513 1000万回のシミュレーション https://i.imgur.com/JKpHsKb.png https://i.imgur.com/A1NiFF0.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/517
518: １３２人目の素数さん [sage] 2025/01/08(水) 20:28:49.92 ID:SKdwPI7I >>516 罵倒しか投稿できないのがFimoseくん。FランにちなんでPhimoseくんから昇格ｗ http://rio2016.5ch.net/test/read.cgi/math/1723152147/518
519: １３２人目の素数さん [sage] 2025/01/08(水) 20:30:20.34 ID:SKdwPI7I >>517 期待値の差よりも 回数の差の分布の方が直感に近いな。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/519
520: １３２人目の素数さん [sage] 2025/01/08(水) 20:33:00.33 ID:bjQ79ehj チンパンジーが図星で発狂中w http://rio2016.5ch.net/test/read.cgi/math/1723152147/520
521: １３２人目の素数さん [sage] 2025/01/09(木) 00:40:37.23 ID:bJYGsz4u >>518 他人がどうこうなんてお前が無能な事とは全く関係ないんだから、まずは国語の勉強してスレタイくらい読めるようになれよwwww 汚い言葉勉強してる場合じゃないぞwww http://rio2016.5ch.net/test/read.cgi/math/1723152147/521
522: １３２人目の素数さん [sage] 2025/01/09(木) 06:57:53.73 ID:0PPPIAFh やはり、下手な考え休むに似たりだな。 オッペンハイマーでの会話 Hahn and Strassmann in Germany. They split the uranium nucleus. How? Bombard it with neutrons. It's a nuclear fission. They did it, they split the atom. It's not possible. I'm gonna try to reproduce it. See? Can't be done. Very elegant. Quite clear. http://rio2016.5ch.net/test/read.cgi/math/1723152147/522
523: １３２人目の素数さん [sage] 2025/01/09(木) 07:00:06.18 ID:0PPPIAFh Fimoseくんはロータリーエンジンの動くシェーマの動画を作れないの？ 完成した動画をみてもハウジングが何かわからんようだし。 やっぱり、弄れるのはforeskinだけみたいだね。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/523
524: １３２人目の素数さん [sage] 2025/01/09(木) 08:14:29.03 ID:jj9pV73T >>522 アンタのことじゃん http://rio2016.5ch.net/test/read.cgi/math/1723152147/524
525: １３２人目の素数さん [sage] 2025/01/09(木) 08:14:47.34 ID:o9+qLk7l limit = 10000 max_count = 0 best_n = 0 for (n in 1:(limit - 1)) { count = 0 for (y in 1:(n - 1)) { if ((n*y) %% (n-y) ==0 & n!=y){ # x が整数になる条件 count = count + 1 } } if (count > max_count) { max_count = count best_n = n } } best_n max_count http://rio2016.5ch.net/test/read.cgi/math/1723152147/525
526: １３２人目の素数さん [sage] 2025/01/09(木) 08:16:10.89 ID:o9+qLk7l nを10000未満の正整数とする。 1/n+1/x=1/y を満たす{x,y}の正整数解の組み合わせの数が最大なnはいくつか？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/526
527: １３２人目の素数さん [sage] 2025/01/09(木) 12:34:42.59 ID:1fujxzuR こういう報道に接するとZ3が必須だな。 https://www.373news.com/_news/storyid/155006/ http://rio2016.5ch.net/test/read.cgi/math/1723152147/527
528: １３２人目の素数さん [sage] 2025/01/09(木) 12:35:57.90 ID:1fujxzuR >>524 映画みてないの？ 英語がわからないFラン？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/528
529: １３２人目の素数さん [sage] 2025/01/09(木) 17:28:07.32 ID:6xzpzcRA >>528 下手な考え休むに似たりって完全にアンタのことなんだけどまさか自覚ないんか？ww http://rio2016.5ch.net/test/read.cgi/math/1723152147/529
530: １３２人目の素数さん [sage] 2025/01/10(金) 02:04:06.88 ID:sfHDQvoL >>523 他人がどうこうとか自分の無能を誤魔化すばっかりwww 自分で定義分かってるとか思い込んでるから何にも理解出来ないんだよwww 思い込みじゃなくて自分できちんと日本語使って数学的な定義述べてみろよwww 出来ないんだろwww http://rio2016.5ch.net/test/read.cgi/math/1723152147/530
531: １３２人目の素数さん [sage] 2025/01/10(金) 02:06:10.84 ID:5hjYx106 Fimoseくんはロータリーエンジンの動くシェーマの動画を作れないの？ 完成した動画をみてもハウジングが何かわからんようだし。 やっぱり、弄れるのはforeskinだけ。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/531
532: １３２人目の素数さん [sage] 2025/01/10(金) 02:14:15.07 ID:H6m+2x8m 結局できなくて発狂かよ 自分で出したクソの処理すらできない模様 http://rio2016.5ch.net/test/read.cgi/math/1723152147/532
533: １３２人目の素数さん [sage] 2025/01/10(金) 02:18:01.93 ID:MdFFAMIP >>531 汚い言葉で誤魔化すばっかりかよwww 他人がどうこうなんて関係なくお前が無能な事を指摘して貰ってるんだから誤魔化さずに自分で努力してみろよwww http://rio2016.5ch.net/test/read.cgi/math/1723152147/533
534: １３２人目の素数さん [sage] 2025/01/10(金) 02:31:21.32 ID:H6m+2x8m 毎回毎回同じことしか言えないし知性のかけらもないのに自称医者東大卒とは笑わせるねww http://rio2016.5ch.net/test/read.cgi/math/1723152147/534
535: １３２人目の素数さん [sage] 2025/01/10(金) 05:36:36.58 ID:5hjYx106 >>534 俺は医科歯科卒 理１を蹴って入学。 あんたは動画すら作れないFランなんだろ。 このスレは出禁だぞ。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/535
536: １３２人目の素数さん [sage] 2025/01/10(金) 07:02:22.75 ID:nzxEBzqG >>535 医科歯科卒wwww 日本語すらまともに使えない中卒以下の低脳が何言ってんの？ Fラン受験出来るのが羨ましいのかよwww http://rio2016.5ch.net/test/read.cgi/math/1723152147/536
537: １３２人目の素数さん [sage] 2025/01/10(金) 08:29:22.85 ID://fsJZeZ 医科歯科ってスレタイ読めなくても卒業できるの？ ルールに則らないことを推奨してるの？ 医師国家試験でルールに則らないことを是とする回答をしたの？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/537
538: １３２人目の素数さん [sage] 2025/01/10(金) 08:30:21.85 ID://fsJZeZ 質問スレで出題ばかりするし質問には回答しないしで、 世の中のいなくなったほうが世界か良くなる人物の中の一人 http://rio2016.5ch.net/test/read.cgi/math/1723152147/538
539: １３２人目の素数さん [sage] 2025/01/10(金) 09:16:14.22 ID:zQblG/AR >>535 まあ匿名掲示板ならいくらでも経歴盛れるもんなそんな主張意味ないよ 本当に医科歯科卒ならそれ相応の能力を示せよ http://rio2016.5ch.net/test/read.cgi/math/1723152147/539
540: １３２人目の素数さん [sage] 2025/01/10(金) 11:01:30.43 ID:QR184n+C 今日は内視鏡が少なくて10:30に終了。予約の検査終了したら帰っていい契約なので( ・∀・)ｲｲ!! http://rio2016.5ch.net/test/read.cgi/math/1723152147/540
541: １３２人目の素数さん [sage] 2025/01/10(金) 11:16:40.37 ID:QR184n+C 熟練職人不足で新年から内視鏡週3回になった。 スタッフが可愛いので続けられる。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/541
542: １３２人目の素数さん [sage] 2025/01/10(金) 11:31:20.92 ID:zQblG/AR 本物の医者がこんな数学板の糞スレにいるわけないじゃん忙しいはずなのに妄想乙 http://rio2016.5ch.net/test/read.cgi/math/1723152147/542
543: １３２人目の素数さん [sage] 2025/01/10(金) 11:56:51.07 ID:02D377cs >>540 これが何かの証拠になるとか思ってんの？ 説得力皆無で今まで誰一人信じて貰ってないのが現実だろwwww 何の意味もない事繰り返し書き込んじゃうって知能低すぎん？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/543
544: １３２人目の素数さん [sage] 2025/01/10(金) 18:50:35.59 ID:nUmCuFV/ 尿瓶ジジイ急にダンマリかよ http://rio2016.5ch.net/test/read.cgi/math/1723152147/544
545: １３２人目の素数さん [sage] 2025/01/11(土) 07:06:37.50 ID:5YoQBcPi A君とB君がそれぞれ同じ規格立方体のサイコロを投げます。 終了条件は A君は、 同じ数の目が2回連続して出たら終了 B君は、 前の目の数より１多い数がでたら修了（但し、６→１でも終了するとする） 【問題】 (1) A君の投げた回数の期待値を求めよ (2) B君の投げた回数の期待値を求めよ (3) A君の投げた回数の方が多い確率をもとめよ (4) 二人の投げた回数が同数で終了する確率をもとめよ http://rio2016.5ch.net/test/read.cgi/math/1723152147/545
546: １３２人目の素数さん [sage] 2025/01/11(土) 08:13:55.14 ID:5YoQBcPi (* コイントスの終了条件を Aくんは、 表が2回連続して出たら終了 Bくんは、 表裏の順に出たら終了（裏表の順では終了しない） する。 *) pA[k_]:=Fibonacci[k-1]/2^k (* Aの終了回数ごとの確率 *) Sum[pA[k],{k,2,Infinity}] (* Σpmf==1 確認 *) Sum[k pA[k],{k,2,Infinity}](* 期待値 *) pB[k_]:=(k-1)/2^k (* Bの終了回数ごとの確率 *) Sum[pB[k],{k,2,Infinity}](* Σpmf==1 確認 *) Sum[k pB[k],{k,2,Infinity}](* 期待値 *) pABd[d_] := Sum[pA[k+d] pB[k],{k
,2,Infinity}] (* Aの方がd回多く終わる確率 *) pABd[0] (* 同じ回数で終わる確率 *) (* Aの方が多く終わる確率 *) Sum[pABd[d],{d,1,Infinity}] // Simplify pbA[k_]:=Sum[pB[j],{j,2,k-1}] (* A=ｋのときのB<Aの確率 *) Sum[pA[k] pbA[k],{k,2,Infinity}] http://rio2016.5ch.net/test/read.cgi/math/1723152147/546
547: １３２人目の素数さん [sage] 2025/01/11(土) 09:15:51.78 ID:EiIJz6LX >>545 結局誰にも信じられずに哀れだね http://rio2016.5ch.net/test/read.cgi/math/1723152147/547
548: １３２人目の素数さん [sage] 2025/01/11(土) 19:01:33.32 ID:WI+T/Kwt (* 1ヶ月間に馬に蹴られて死亡する兵士の数が母数λ=2のポアソン分布に従い、 1ヶ月間に補充される兵士の数はp=1/3の幾何分布に従うとする。 いずれの分布も定義域は非負整数。 1ヶ月後に兵士が増えている確率を算出し、シミュレーションとの合致を確認せよ。 *) p=1/3; lambda=2; pDiff[d_] := Sum[PDF[GeometricDistribution[p],k+d] PDF[PoissonDistribution[lambda],k],{k,0,Infinity}] Sum[pDiff[d],{d,1,Infinity}] % // N k=10^6; po=RandomVariate[PoissonDistri
bution[lambda],k]; ge=RandomVariate[GeometricDistribution[p],k]; d=ge-po; Boole[#>0& /@ d]//Mean // N http://rio2016.5ch.net/test/read.cgi/math/1723152147/548
549: １３２人目の素数さん [sage] 2025/01/11(土) 21:32:12.61 ID:EiIJz6LX >>548 マルチポストしてるのに誰にも相手にされてないみたい 実に哀れ http://rio2016.5ch.net/test/read.cgi/math/1723152147/549
550: １３２人目の素数さん [sage] 2025/01/11(土) 23:12:46.87 ID:WI+T/Kwt COVID-19の潜伏期はμ=1.6 、σ=0.5の対数正規分布(オレンジ色)に インフルエンザの潜伏期は形状母数k=2.0,尺度母数θ=1.0のガンマ分布(青色)に従うとする。 https://i.imgur.com/QmlK7yi.png 同時に感染したときにインフルエンザの方が先に発症する確率を求めなさい。 分布は独立でウイルス干渉はないものとする。 pdfz =\(z) integrate( \(x) dlnorm(x,1.6,0.5) * dgamma(x-z,shape=2,scale=1), 0,Inf,rel.tol = 1e-12)$value pdfz=Vectorize(pdfz) curve(pdfz(x),-10,
30) 1-integrate(pdfz,-Inf,0,rel.tol = 1e-12)$value k=1e7 Cov=rlnorm(k,1.6,0.5) Flu=rgamma(k,shape=2,scale=1) plotPost(Cov-Flu,compVal = 0,col=7) curve(pdfz(x),add=TRUE,lwd=2) mean(Cov>Flu) https://i.imgur.com/NKroLQc.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/550
551: １３２人目の素数さん [sage] 2025/01/11(土) 23:13:48.64 ID:WI+T/Kwt ここはRやWolframによる解法のコード置き場。 亀レスとの照合用。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/551
552: １３２人目の素数さん [sage] 2025/01/12(日) 08:08:57.37 ID:ht709xDC >>551 ちがいます >>5 http://rio2016.5ch.net/test/read.cgi/math/1723152147/552
553: １３２人目の素数さん [sage] 2025/01/12(日) 08:45:21.01 ID:jXkNnsH9 >>551 ここはお前専用の隔離スレだよ 医者・東大卒なんだろ？w もうここ以外に書き込むなよ http://rio2016.5ch.net/test/read.cgi/math/1723152147/553
554: １３２人目の素数さん [sage] 2025/01/12(日) 08:45:54.02 ID:Sb0hrOht >>5の日本語すら理解できないチンパンが紛れ込んでいる模様 http://rio2016.5ch.net/test/read.cgi/math/1723152147/554
555: １３２人目の素数さん [sage] 2025/01/12(日) 10:55:46.68 ID:qiFin7nw 【臨床応用問題】 下記のデータから適合する分布を選びそのパラメータを算出せよ。 https://i.imgur.com/N9lQ4Qu.png 算出例 AICやBICで判定すると非負量を定義域にする分布ではWeibull分布が最良だった。 そのパラメータは shape scale 6.404264 1.551445 インフルエンザの潜伏期は形状母数k=2.0,尺度母数θ=1.0のガンマ分布(青色)に従うとする。 https://i.imgur.com/YAI1lc2.png 【臨床応用問題】 発熱外来では呼吸器症状を呈するインフルエンザと消化器症状を呈
するノロウイルス感染に二分されたとする。 発熱外来で働く職員が同時に両方のウイルスに感染したときに同じ日に呼吸器症状と消化器症状が出現する確率を求めよ。 library(fitdistrplus) set.seed(2025) noro=c( runif(2,0,24), runif(2,25,28), runif(16,29,32), runif(16,33,36), runif(14,37,40), runif(2,41,44), runif(5,45,48) )/24 fit=fitdist(noro,'weibull') ; fit gofstat(fit) plot(fit) hist(noro,freq=F,main="",col=4,breaks = 'scott') fit$estimate curve(dweibull(x,shape=fit$estimate[1],scale=fit$estimate[2]),a
dd=TRUE,lwd=2) par=fit$estimate ; par Noro=rweibull(k,shape=par[1],scale=par[2]) plotPost(Noro,freq=FALSE,main='Norovirus',col=7) curve(dweibull(x,shape=fit$estimate[1],scale=fit$estimate[2]),add=TRUE) Flu=rgamma(k,shape=2,scale=1) par(mfrow=c(3,1)) plotPost(Flu,col=4) plotPost(Noro,col=2) plotPost(Flu-Noro,col=7,compVal = 0) mean(abs(Flu-Noro) < 1) http://rio2016.5ch.net/test/read.cgi/math/1723152147/555
556: １３２人目の素数さん [sage] 2025/01/13(月) 05:38:55.68 ID:cjaPfdCl https://i.imgur.com/0tg88Vd.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/556
557: １３２人目の素数さん [sage] 2025/01/13(月) 06:51:55.39 ID:oWSkSukX gr=Tuples[{"+","-","*","/",""},8]; join[x_] := ( {x1,x2,x3,x4,x5,x6,x7,x8}=x; "1"<>x1<>"2"<>x2<>"3"<>x3<>"4"<>x4<>"5"<>x5<>"6"<>x6<>"7"<>x7<>"8"<>x8<>"9" ) op[x_] := ToExpression@join@x re = op /@ gr; so
lve[n_] := join /@ Select[gr,op@#==n&] http://rio2016.5ch.net/test/read.cgi/math/1723152147/557
558: １３２人目の素数さん [sage] 2025/01/13(月) 07:07:22.30 ID:BfHaF1g6 下記のデータから潜伏時間の中央値の９５％信頼区間を求めよ。 https://i.imgur.com/N9lQ4Qu.png 算出方法にはいくつかの流儀があるが、好みの方法で 算出例 中央値を代表値としてブートストラップ法で算出した結果 BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS Based on 10000 bootstrap replicates CALL : boot.ci(boot.out = one.boot(dat, median, 10000, student = TRUE, M = 30)) Intervals : Level Normal Basic Studentized 9
5% (33.05, 35.86 ) (34.50, 34.50 ) (34.50, 34.50 ) Level Percentile BCa 95% (34.5, 34.5 ) (30.5, 30.5 ) Calculations and Intervals on Original Scale Warning : BCa Intervals used Extreme Quantiles Some BCa intervals may be unstable http://rio2016.5ch.net/test/read.cgi/math/1723152147/558
559: １３２人目の素数さん [sage] 2025/01/13(月) 08:16:08.96 ID:BfHaF1g6 Min[re] Max[re] solve[Min[re]] solve[Max[re]] Commonest[re] Length[solve[0]] http://rio2016.5ch.net/test/read.cgi/math/1723152147/559
560: １３２人目の素数さん [sage] 2025/01/13(月) 09:19:47.66 ID:GL4RAVeD https://i.imgur.com/uwF2H4N.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/560
561: １３２人目の素数さん [sage] 2025/01/13(月) 17:58:04.04 ID:kvwM3Liy gr=Tuples[{"+","-","*","/",""},8]; join[x_] := "1"<>x[[1]]<>"2"<>x[[2]]<>"3"<>x[[3]]<>"4"<>x[[4]]<>"5"<>x[[5]]<>"6"<>x[[6]]<>"7"<>x[[7]]<>"8"<>x[[8]]<>"9" op[x_] := ToExpression@join@x re = op /@ gr; sol
ve[n_] := join /@ Select[gr,op@#==n&] solve /@ Table[1111*m,{m,1,9}] {{1+23*45+6+78-9, 1*2*3*4*5/6*7*8-9, 1*23*45-6-7+89, 1*23/4*56+789, 1*234*5+6+7-8*9, > 1*234*5+6-7*8-9, 1*234*5-6*7-8-9, 1/2/3*4*5*6*7*8-9, 12+34*5*6+7+8*9, > 12*34-5+6+78*9}, {1+2+34*5/6*78+9}, {}, {}, {}, {12-3*45+6789}, {}, {}, {}} http://rio2016.5ch.net/test/read.cgi/math/1723152147/561
562: １３２人目の素数さん [sage] 2025/01/14(火) 06:30:12.29 ID:6B9y8ZmQ rm(list=ls()) library(fitdistrplus) solve=\(hours){ noro=c( runif(2,21,24), runif(2,25,28), runif(16,29,32), runif(16,33,36), runif(14,37,40), runif(2,41,44), runif(5,45,48) ) #hist(noro) distr=c("norm","lnorm","gamma","weibull") aic=NULL for(i in 1:length(distr)) aic=c(aic,gofstat(fitdist(noro,distr[i]))$aic) fitted=distr[which.min(aic)] fit=fitdist(noro,"fitted") # plot(fit) par=fit
$estimate #hist(noro,freq=F,main='Noro Virus incubation',col=4,axes=FALSE,xlab='hours',ylab='') ; axis(1) #curve(dgamma(x,par[1],par[2]),add=TRUE) #curve(dgamma(x,par[1],par[2]),0,5*24) pgamma(hours,par[1],par[2],lower.tail = FALSE) } solve(48) http://rio2016.5ch.net/test/read.cgi/math/1723152147/562
563: １３２人目の素数さん [sage] 2025/01/14(火) 18:27:07.83 ID:KQCqu9ZF " Fくんが２日前（４８時間とする）に喫食した食材でノロウイルスによる食中毒が発生したとする。 https://i.imgur.com/Aaz7iOi.png のデータから現在無症状のFくんが今後、発症する確率を算出しなさい。 算出に必要な条件は適宜設定してよい。 " rm(list=ls()) library(fitdistrplus) solve=\(hours,distr=c("lnorm","gamma","weibull"),verbose=FALSE){ noro=c( runif(2,21,24), runif(2,25,28), runif(16,29,32), runif(16,33
,36), runif(14,37,40), runif(2,41,44), runif(5,45,48) ) aic=NULL for(i in 1:length(distr)) aic=c(aic,gofstat(fitdist(noro,distr[i]))$aic) best_fit=distr[which.min(aic)] fit=fitdist(noro,best_fit) par=fit$estimate if(verbose){ hist(noro,freq=F,main='Noro Virus incubation',col=4,axes=FALSE,xlab='hours',ylab='') ; axis(1) curve(dgamma(x,par[1],par[2]),add=TRUE) } paste0('p',best_fit,'(hours,par[1],par[2],lower.tail = FALSE)') |> str2lang() |> eval() } solve(48,verbose=TRUE) k=1e4 p48=replicate(k,solve(
48,c("lnorm","gamma","weibull"))) round(HDInterval::hdi(p48),3)[1:2] median(p48) mean(p48) BEST::plotPost(p48,showCurve = F,col='lightgreen',xlab='p[hour>48]') http://rio2016.5ch.net/test/read.cgi/math/1723152147/563
564: １３２人目の素数さん [sage] 2025/01/14(火) 20:19:32.39 ID:KG5WXTP4 太郎くんと花子さんが一緒に食事をしてノロウイルスのいる生牡蠣を食べたとする。 https://i.imgur.com/lro8LGD.png のデータから太郎くんと花子さんが同じ日に発症する確率を求めなさい。 算出に必要な条件は適宜設定してよい、すなわち、好みの方法で計算してよい。 算出例：incubation periodがどんな分布に従うかを前提にしないで乱数発生させて算出。 R言語 calc=\(){ noro=c( runif(2,21,24), runif(2,25,28), runif(16,29,32), runif(16,33,36), runif(14,37,40),
runif(2,41,44), runif(5,45,48)) abs(diff(sample(noro,2,replace = TRUE))) < 24 } mean(replicate(1e6,calc())) Wolfram言語 calc[] := Module[{noro}, noro := Flatten@{ Table[Random[Real,{21,24}], 2], Table[Random[Real,{25,28}], 2], Table[Random[Real,{29,32}],16], Table[Random[Real,{33,36}],16], Table[Random[Real,{37,40}],14], Table[Random[Real,{41,24}], 2], Table[Random[Real,{45,48}], 5]}; (Abs@Differences@RandomChoice[noro,2])[[1]] < 24 // Boole ] Table[calc[],10^6] // Mean // N http://rio2016.5ch.net
/test/read.cgi/math/1723152147/564
ﾒﾓ帳
(0/65535文字)
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 427 ﾚｽあります
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
AAｻﾑﾈｲﾙ
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.021s