18/10/26 13:59:22.75 +rLvZj8m.net
100万回シミュレーションして確認
> # sim
> n=8
> k=3
> fo=readRDS('fast_order.rds')
> y=fo[sample(1:nrow(fo),1),]
> cat(LETTERS[y])
A D B E F G H C
> sim <- function(){
+ x=sample(1:n,k)
+ re=numeric(k)
+ for(i in 1:k){
+ re[i]=which(y==x[i])
+ }
+ max(re) - max(x)
+ }
> K=1e6
> re=replicate(K,sim())
> c(Earlier=mean(re<0),Later=mean(re>0),Equal=mean(re==0))
Earlier Later Equal
0.625315 0.267296 0.107389
319:卵の名無しさん
18/10/26 16:45:51.73 +rLvZj8m.net
東京医大、本来合格者入学許可へ 今年の受験生50人
2018年10月25日 02時06分
東京医科大=8月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生50人に対し、来年4月の入学を認める方針を固めたことが24日、関係者への取材で分かった。
昨年の本来合格者19人については、難しいとの意見が出ているもようだ。東京医大は50人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLリンク(www.nishinippon.co.jp)
URLリンク(www.t)
320:okyo-med.ac.jp/med/enrollment.htmlによると 学年 第1学年 第2学年 在学者数 133 113 昨年入学者の留年者や退学者が0として、 大学が公式認定した裏口入学者が少なくとも今年は133人中50人、昨年が113人中19人ということになる。 裏口入学率の期待値、最頻値、およびその95%信頼区間を求めよ。
321:卵の名無しさん
18/10/26 17:35:16.08 CYyweUPk.net
>>310
裏口入学者には裏金返して一旦、除籍、
来年の入試を受けさせて学力考査のみで合否判定、合格すればもとの学年に復学。
最低、これくらいはやるべき。
322:卵の名無しさん
18/10/27 00:52:44.70 lsV1Nek+.net
もう私立医大って裏口ばかりなんだろうな。
せめて慶応だけは裏口が無いと信じたい。
323:卵の名無しさん
18/10/27 10:45:04.64 mOp8hKBt.net
裏口私立医 でてこいや!
324:卵の名無しさん
18/10/27 14:29:59.93 NpJ0k/uw.net
六面体のサイコロでP君のサイコロは3面が1、Q君のサイコロは2面が1とする。
サイコロを降ってどちらか一方が1であればそちらが勝者。
どちらも1であるときはもしくはどちらも1でないならば どちらか一方だけが1が出て勝者が決まるまで繰り返す。
P君、Q君の勝つ確率を求めよ。
級数の和を使って
Σ[1,∞](1/4)^i = 1/4 *(1-1/4)= 1/3
Pr[P:win]=1/2+1/3*1/2=2/3
Pr[Q:win]1/4+1/3*1/4=1/3
で解ける(ド底辺シリツ医大卒の裏口バカには無理)が
数学板に書いたら速攻で級数計算なしでの正解が返ってきた。
世の中には頭のいい人がいるもんだ。
まあ、ド底辺シリツ医大に裏口入学するしかないバカもいるが。
325:卵の名無しさん
18/10/27 18:05:11.43 b2KaNqWj.net
興奮状態 大丈夫?
326:卵の名無しさん
18/10/27 18:45:07.04 NpJ0k/uw.net
>>315
>313はド底辺シリツのキーワードがないから別人。
327:卵の名無しさん
18/10/27 19:20:54.18 Kx89tomJ.net
理3 そこそこ頭いい
旧帝医 人間として普通
地方医、私立医 自分でも頭わるいこと自覚してるだろ?
328:卵の名無しさん
18/10/27 19:22:19.50 Kx89tomJ.net
>東京医大は50人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
うわひどい
これつまり年度をまたいで考えれば
本来受からない人間が(過去)受かって
受かるべき人間が(来年)落とされるってことじゃん
憲法違反じゃね?
329:卵の名無しさん
18/10/27 19:23:59.97 Kx89tomJ.net
私立医はバカだが地方国立医もたいがいだぞ
東京育ちの私立医は少なくとも文化資本はあるわな。
勉強しかできない(くせに勉強もたいしてできない)田吾作医よりも下とは言い切れない
330:卵の名無しさん
18/10/27 19:47:55.49 ZN2oZPvM.net
>>318
裏口入学の50人を除籍処分にしないかぎり、信頼の回復はないね。
つまり、いつまで経ってもシリツ医大卒=裏口バカと汚名は拭えない。
シリツ出身者こそ、裏口入学に厳しい処分を訴えるべきなのに、
僕も裏口入学だよ、と自慢するのが、シリツのクオリティー。
331:卵の名無しさん
18/10/27 19:51:57.10 ZN2oZPvM.net
>>316
ド底辺シリツ医大卒は自分が頭が悪いことすらわかってないぞ。
これが現実 1次方程式もできないド底辺シリツ医大卒の記録
URLリンク(imagizer.imageshack.com)
ド底辺シリツ医大卒の裏口バカは自分が馬鹿だという自覚がないから、与えられた数値から計算できない、問題が悪いと言い張るんだよなぁ。悪いのはそいつの頭なのに。
332:卵の名無しさん
18/10/27 19:52:52.05 ZN2oZPvM.net
おまけに自分が置かれている立場すら自覚していない。
東京医大の事件で点数加算の不正まで発覚しているのに、
>裏?
>そんなのあるかよばか
と言ってる椰子がいたな。
こいつも裏口入学したシリツ医だろうと推測されるね。
URLリンク(imagizer.imageshack.com)
333:卵の名無しさん
18/10/27 19:55:19.11 ZN2oZPvM.net
裏口入学を公言した高須院長を根拠もなく特待生認定してソースを求められたときの対応がこれな。
ド底辺シリツ医大卒の裏口バカ医って
患者から診断や治療の根拠を求められたときも
こんな対応をしているのかな?
高須院長は特待生だという裏口バカの根拠なき主張に対して
ソース出せと言われた時ののやりとり。
URLリンク(imagizer.imageshack.com)
URLリンク(imagizer.imageshack.com)
人間のクズだな。
334:卵の名無しさん
18/10/27 20:31:09.75 b2KaNqWj.net
いわゆる 勝ち組を自分の低いレベルまで引きずり落としたい人
335:卵の名無しさん
18/10/27 20:37:52.35 b2KaNqWj.net
女にモテて 経済力がある人に対するルサンチマン
336:卵の名無しさん
18/10/27 20:40:27.10 0SWbBA8x.net
>>324
ド底辺シリツ医大卒以下に堕ちるのは想像すらできん。
ナマポも避けたがるのがド底辺シリツ医大卒
実例↓
517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな?
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。
337:卵の名無しさん
18/10/27 20:41:25.17 0SWbBA8x.net
>>325
>>167
>>165
>>158
出身校も言えないような奴って気持ち悪いぞ。
んで、あんたどこ卒?
これに答えられない馬鹿だと患者も大変だぞ。
はよ、答えてみ。
66 卵の名無しさん sage 2018/10/08(月) 09:22:57.88 ID:mQHuWxdH
ド底辺聖マ卒 頭悪すぎ
という意見がでてるぞ。
338:卵の名無しさん
18/10/27 20:48:16.97 b2KaNqWj.net
ニーチェを読むことをお勧めします
339:卵の名無しさん
18/10/27 20:50:32.19 0SWbBA8x.net
>>328
>>278
>>259
あんた、どこ卒?
これであっているのか?
66 卵の名無しさん sage 2018/10/08(月) 09:22:57.88 ID:mQHuWxdH
ド底辺聖マ卒 頭悪すぎ
340:卵の名無しさん
18/10/27 20:55:11.15 b2KaNqWj.net
え? ニーチェ読めないの? 俺は帝京大学卒だがな
341:卵の名無しさん
18/10/27 20:57:09.93 0SWbBA8x.net
>>330
これを暗唱するといいぞ。
URLリンク(ishikisoku.com)
342:卵の名無しさん
18/10/27 20:58:34.00 0SWbBA8x.net
>>330
裏口バカが解けなかったこれの正解だしてみ
URLリンク(imagizer.imageshack.com)
343:卵の名無しさん
18/10/27 20:59:41.60 0SWbBA8x.net
>>330
これ人間のクズだと思うだろ?
高須院長は特待生だという裏口バカの根拠なき主張に対して
ソース出せと言われた時ののやりとり。
URLリンク(imagizer.imageshack.com)
URLリンク(imagizer.imageshack.com)
344:卵の名無しさん
18/10/27 21:02:21.59 kZ9du42U.net
高学歴スレに宣伝して来てやったぞ
ありがたく思えド底辺どもめ
スレリンク(newsplus板)
345:卵の名無しさん
18/10/27 21:13:42.31 b2KaNqWj.net
教養ないんだねw
346:卵の名無しさん
18/10/27 21:16:25.64 b2KaNqWj.net
キルケゴール ニーチェ カント ヘーゲル フロイト ラカン フーコー デリダ 読んだ?
347:卵の名無しさん
18/10/27 21:19:20.09 b2KaNqWj.net
それと フッサールは読んだ?レビ ストロース 読んだ?構造主義って知ってるかな?ソーカル事件て知ってる?柄谷行人て知ってる?
348:卵の名無しさん
18/10/27 21:27:35.20 b2KaNqWj.net
学問って広いんだね w
349:卵の名無しさん
18/10/27 21:34:33.98 b2KaNqWj.net
あ ハイデガーとアドルノ ホルクマイヤー を忘れてたわw あとはチョムスキーな?
350:卵の名無しさん
18/10/27 21:36:34.76 b2KaNqWj.net
ハーバーマスとレヴィナスは読むべきだね
351:卵の名無しさん
18/10/27 21:49:59.16 b2KaNqWj.net
構造主義はインパクトあったね フロイト マルクス ニーチェは当然読んでるよね 医学も含む学問の基本だよな?
352:卵の名無しさん
18/10/27 22:19:02.05 b2KaNqWj.net
こいつは教養ないな ちょこっといじれば知的レベルがバレるw
353:卵の名無しさん
18/10/28 05:21:48.30 kDnv47VV.net
外人の名前に権威を見出す俗物精神
勉強はできても教養はない
354:卵の名無しさん
18/10/28 05:25:27.95 kDnv47VV.net
問題が解けるとか
いろんなことを知ってるというだけでなく
それを知的な振る舞いに統合できて初めて教養と言うのだよね
テスト勉強とかポケモンの名前を全部覚えるとかなら子供でもできるわけで
355:卵の名無しさん
18/10/28 06:27:25.96 Q5x4cT59.net
知らないことは調べればいい、ただそれだけ。
こういうのは調べても出てこない。
まあ、計算スクリプトを俺がネットにあげたから調べれば出ては来るがな。
door to balloon time
356:が短縮するのに必要な計算。 診療所から病院に患者を救急搬送する。 病院から医師搭乗の救急車が診療所に向かっており10時到着予定と連絡が入った。 患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら 救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。 普通車から救急車への患者の乗り換えで10分余分に時間がかかる。 道路事情から救急車は病院から診療所への道は 平均時速60kmで、逆方向は平均時速45kmで定速走行する。診療所の普通車は信号待ちもあり平均時速30kmで定速走行する。 何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか?
357:卵の名無しさん
18/10/28 06:42:56.20 Q5x4cT59.net
東京医大、本来合格者入学許可へ 今年の受験生50人
2018年10月25日 02時06分
東京医科大=8月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生50人に対し、来年4月の入学を認める方針を固めたことが24日、関係者への取材で分かった。
昨年の本来合格者19人については、難しいとの意見が出ているもようだ。東京医大は50人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLリンク(www.nishinippon.co.jp)
URLリンク(www.tokyo-med.ac.jp)によると
学年 第1学年 第2学年
在学者数 133 113
昨年入学者の留年者や退学者が0として、
大学が公式認定した裏口入学者が少なくとも今年は133人中50人、昨年が113人中19人ということになる。
裏口入学率の期待値、最頻値、およびその95%信頼区間を求めよ。
358:卵の名無しさん
18/10/28 06:47:18.05 Q5x4cT59.net
医師国家試験って検索できれば小学生でも正答できるね。
医師国家試験
URLリンク(imagizer.imageshack.com)
こっちの方が難しい
問題
0,1,2だけを並べて整数を作っていきます。 ただし、左端には0はこないようにして、右端は必ず0にします。 また、となりあった数は異なるようにします。 このとき次の問に答えなさい。 (1)5けたの整数は何通り作れますか。
(2)8けたの整数は何通り作れますか。
359:卵の名無しさん
18/10/28 06:49:10.49 Q5x4cT59.net
一次方程式なしで解ける問題
マイクロソフトの入社試験問題
ボートに乗って、手荷物を船外に放り出すと、水位は上がるか、下がるか
あるブログで
>正解は「手荷物が浮くならば、全体の水位に変化はなく、手荷物が沈むならば水位は下がる。」
>意地悪な問題だな。
というのを読んだ。
答だけ覚えても楽しくないので荷物が沈むときどれだけ水位が低下するか計算してみる。
池の面積をA、水の比重は1、
手荷物の体積をV、比重とgとすると
水位はV(g-1)/A下がる。
中学入試程度の問題だが、医師国家試験よりずっと難しいね。
ド底辺シリツ医大卒の裏口バカ排除にこういう基礎学力を試す問題をだせばいいのに。
360:卵の名無しさん
18/10/28 06:50:49.78 Q5x4cT59.net
>>344
問題が解ける以前の低学力がシリツ医大卒。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こすわけなんだな。
URLリンク(i.imgur.com)
361:卵の名無しさん
18/10/28 07:53:26.49 88wXw622.net
>>330
たぶんそれ本を読んだと言うだけでニーチェの真意は理解してないと思う。
ましてや読んだのは翻訳した書物でしょ。
原文で読まないと著者の真意は分からないと思うよ。
帝京大卒くらいの基礎学力じゃなかなか理解するのは難しいと思う。
専門医と同じ。
マルチプルチョイスの専門医試験に受かったからといって、
みんながみんな実力があるとは限らない。
巷には実力の無い名前だけ専門医がウヨウヨいるのと同じw
362:卵の名無しさん
18/10/28 09:04:03.89 Q5x4cT59.net
>>350
帝京といえばこれ
URLリンク(ishikisoku.com)
363:卵の名無しさん
18/10/28 09:07:44.99 Q5x4cT59.net
普通は教育投資すればある程度の学力は身につく。
シリツ医大の高額学費が払える経済環境で教育投資したはず。
それでもド底辺シリツ医大にしか入れないのはそういうことだw.
364:卵の名無しさん
18/10/29 19:25:18.11 ykpN1OD3.net
東京裏口医大 とんでもねーことになったな~
365:卵の名無しさん
18/10/29 20:44:25.21 6vG4Elw2.net
世界まる見え!ハンパないSP★4
スレリンク(liventv板)
366:卵の名無しさん
18/10/29 21:28:26.20 72tLNU2F.net
シリツA医大、シリツB医大で試験をしたところ、裏口学生の平均点も正規学生の平均点もA医大の方が上だったのに
裏口正規合わせての平均点はB医大の方が上だったという。
このようなことがありえるのだろうか?
367:卵の名無しさん
18/10/30 00:37:45.25 Eer9XEMu.net
>>217
Haskellに移植してexeファイルにコンパイルしてコマンドラインから実行できるようにした。
import System.Environment
import Data.List
import Data.List.Split
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
main = do
argList <- getArgs -- m : 縦マス(短軸) n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
q = [0..m*n-1]
matQ = chunksOf n q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換
treasure = combinations k q -- 宝の組み合わせ
ip y = minimum $ map(\x -> elemIndices x p!!0) y -- 宝の、配列pでのindex列を求めて最小値を返す
iq y = minimum $ map(\x -> elemIndices x q!!0) y
idxp = map ip treasure -- 宝の組み合せで実行して
idxq = map iq treasure
p_q = zipWith (-) idxp idxq -- 差をとって大小判別
p1st = length $ filter (<0) p_q -- 短軸方向探索pが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw
368:卵の名無しさん
18/10/30 00:40:00.55 Eer9XEMu.net
きちんと動いていて( ・∀・)イイ!! これで安眠できる。
>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184
369:卵の名無しさん
18/10/30 08:08:15.89 Eer9XEMu.net
-- 全部の宝を先にみつけた方が勝者とするように書き換えてみた
import System.Environment
import Data.List
import Data.List.Split
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
main = do
argList <- getArgs -- m : 縦マス(短軸) n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
q = [0..m*n-1]
matQ = chunksOf n q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換
treasure = combinations k q -- 宝の組み合わせ
ip y = maximum $ map(\x -> elemIndices x p!!0) y -- 宝の、配列pでのindex列を求めて最小値を返す
iq y = maximum $ map(\x -> elemIndices x q!!0) y
idxp = map ip treasure -- 宝の組み合せで実行して
idxq = map iq treasure
p_q = zipWith (-) idxp idxq -- 差をとって大小判別
p1st = length $ filter (<0) p_q -- 短軸方向探索pが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw
370:卵の名無しさん
18/10/30 08:10:47.68 Eer9XEMu.net
>>358
バグの指摘歓迎
まあ、最大のバクは教育に金をかけてもド底辺シリツ医大にしか進学できなかった底辺頭脳だがなwwww
これね!!
URLリンク(ishikisoku.com)
371:卵の名無しさん
18/10/30 09:57:59.27 6vqzf6Lr.net
裏口私立医はバグ
372:卵の名無しさん
18/10/30 14:03:01.40 Eer9XEMu.net
>>360
1: It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。
というのがド底辺シリツ医大卒の第一法則なのだが、
不正が発覚した学生を除籍処分にしないシリツ医大もクズだよなぁ。
373:卵の名無しさん
18/10/30 21:03:47.41 kXLKssbE.net
>>356
数学板でc++に移植したコードが投稿された。
面白い問題をみると解決せずにはいられない人達に感謝。
shape(5x6), tcount=15 ---> pwin:25463979, qwin:25488051,eve
374:n:104165490 シミュレーションと比較してみる。 > p=25463979 > q=25488051 > e=104165490 > c(p,q,e)/(p+q+e) [1] 0.1641593 0.1643145 0.6715263 母集団1億5千万から10万のサンプリングなのでまずまずの近似だな。 > c(p1st_sim=mean(re15<0),q1st_sim=mean(re15>0),draw_sim=mean(re15==0)) p1st_sim q1st_sim draw_sim 0.16354 0.16515 0.67131
375:卵の名無しさん
18/10/30 21:15:23.20 kXLKssbE.net
いつも最後はcプログラマーのお世話になるなぁ。
お気楽言語で四苦八苦しているのを示すと
月光仮面のように現れてcのコードを投稿してくれる。
今回はc++なのでついていけそうもない。
376:卵の名無しさん
18/10/31 09:58:29.87 +mq0uM90.net
低脳私立医 でてこいや!
377:卵の名無しさん
18/10/31 10:29:16.05 q9Xt3MLv.net
>>355
B医大は裏口の割合が少なかった。
B医大の正規学生はA医大の正規学生に負けていたが、Aの裏口よりはずっと高かった
378:卵の名無しさん
18/10/31 11:53:23.56 bPxngJ5R.net
>>365
数字を伴った具体例が出せないとは、さてはシリツだな。
379:卵の名無しさん
18/10/31 15:20:14.36 qgJ05S6D.net
>>363
最適化された算出法のPythonでのコードが投稿されて
>PythonをHaskellにすればいい
>Haskellにもリスト内包表記があるんだから
とのアドバイスを頂いた。まさに移植の核心をつく助言だったな。サッカーでいうキラーパスだぜ。
import System.Environment
choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]
nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0
nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin n m c
main = do
argList <- getArgs -- m : 縦マス(短軸) n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)
380:卵の名無しさん
18/10/31 15:23:26.52 qgJ05S6D.net
コンパイルしてコマンドラインで動作するようにしたの汎用性があって( ・∀・)イイ!!
コンパイルできて不定長整数が扱えるHaskellならではだな。
こういうのも瞬時に計算してくれた、10×20部屋に宝箱100個
>takara 10 20 100
p1st = 15057759425309840160151925452579572328997602171271937639470,
q1st = 15057796557877993527038542474310161591275806044157319150135,
draw = 60432921540347294111327092128863840691952977587098698541050
381:卵の名無しさん
18/10/31 19:20:51.03 qgJ05S6D.net
>>367
{-
数学板で drawのtypo を指摘されてたのでデバッグ版を再度アップ。
数学板は細かいところまできちんと読んでいてくれる人がいて( ・∀・)イイ!!
-}
import System.Environment
choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]
nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0
nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin m n c
main = do
argList <- getArgs -- m : 縦マス(短軸) n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)
382:卵の名無しさん
18/10/31 19:22:13.69 qgJ05S6D.net
>>368
10×20部屋に宝箱100個の計算も修正
p1st = 15057759425309840160151925452579572328997602171271937639470
q1st = 15057796557877993527038542474310161591275806044157319150135
draw = 60432958672915447478213709150594429954231181459984080051715
383:卵の名無しさん
18/10/31 19:48:00.38 qgJ05S6D.net
14 x 15 宝105での算出値
Haskell
15050058965790721521451955450789741426719541629904231171092323
1505005898787002066617563667802976547834618559039773640553201
384:2 60392361586649265993221049300205393824371020946210111218855105 Python 15050058965790721521451955450789741426719541629904231171092323 15050058987870020666175636678029765478346185590397736405532012 60392361586649265993221049300205393824371020946210111218855105 C double参考値 15050058965790725315104920912133576716207522209102038099296256 15050058987870027325121443554283577687551916746692125429596160 60392361586649287873092121075176030813045740046471488224624640 HaskellとPythonは完全に一致
385:卵の名無しさん
18/10/31 19:50:07.01 qgJ05S6D.net
210から105個選ぶ組み合わせは
90492479540310008180848641429024900729436748166512078795479440
通り
386:卵の名無しさん
18/10/31 21:06:16.10 qgJ05S6D.net
{- 見づらいshellを起動するより、WinGHCiで対話型に算出しても速度は変わらないな。コンパイルしてコマンドラインから実行可能なようにもしておく-}
import System.Environment
choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]
nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0
nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin m n c
takara m n k = do
putStrLn $ "短軸p1st = " ++ show(mwin m n k)
putStrLn $ "長軸q1st = " ++ show(nwin m n k)
putStrLn $ "同等draw = " ++ show(draw m n k)
main = do
argList <- getArgs -- m : 縦マス(短軸) n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)
{-Prelude> takara 10 20 100
短軸p1st = 15057759425309840160151925452579572328997602171271937639470
長軸q1st = 15057796557877993527038542474310161591275806044157319150135
同等draw = 60432958672915447478213709150594429954231181459984080051715-}
387:卵の名無しさん
18/11/01 10:06:06.18 ZkYkKASz.net
数学板の質問スレに投稿していたら、勝率確率の算出プログラムをPythonからHaskellに移植せよとの課題が出された。
プログラム練習がてらに回答したら
引き分け計算式を間違えているとの指摘が返ってきた。
URLリンク(rio2016.2ch.net)
きちんと読んでもくれているんだなあと感心した。
388:卵の名無しさん
18/11/01 11:27:15.31 070VkLxq.net
今月国連WHOで解禁されようとしている医療大麻オイル
URLリンク(plaza.rakuten.co.jp)
医薬品としてはEpidiolexの名で承認されてます
国連で解禁されると日本の大麻取り締まり法がなくなるか、最低でも4条が改正され病院でも処方できるようになります
389:卵の名無しさん
18/11/01 22:07:12.63 8nUg0zIu.net
ほんと裏口私立医ばかり
順天とか昭和とか東京医とか馬鹿ばっかw
390:卵の名無しさん
18/11/02 00:49:56.65 CfCNBter.net
広島県の福山友愛病院で、
患者の病状と関係ない薬を大量に投与した。
しかも、意図的にです!
国は違うがドイツの場合
↓
裁判所で開かれた公判で、患者100人を殺害した罪を認めた。これでこの事件は、同国で戦後最悪級の連続殺人事件となった。
起訴されたのは、ドイツ北部デルメンホルストとオルデンブルクの病院で看護師をしていたニルス・ヘーゲル受刑者(41)。
当時勤務していたドイツ北部の2つの病院で2000~2005年にかけ、34~96歳の患者を殺害したことを認めた。
同受刑者は自分の蘇生措置の腕を同僚に見せびらかす目的や、退屈しのぎの目的で、
↓↓↓↓↓
患者に処方されていない薬を投与していたとされる。
↑↑↑↑↑
ドイツの場合。まあ日本は日本だが・・・
大口の病院は看護士の単独犯だったわけで捕まったが・・・
福山友愛病院は・・・・・
URLリンク(youtu.be)
391:卵の名無しさん
18/11/02 11:43:48.07 u7uevd67.net
内視鏡学会で今日は代診バイト、コストパフォーマンスが良くて( ・∀・)イイ!!。
ステーキハウスを予約した。
392:卵の名無しさん
18/11/02 12:32:58.23 k2/mPNL+.net
医学部入試不正はなぜ私立大学ばかり…背景に何が?
中原一歩2018.11.2 11:30AERA
東京医科大に慰謝料と入試成績の開示などを求めて会見した弁護団。
会見には顔と実名をさらして現役医大生も参加した (c)朝日新聞社
また補欠合格者の選定でも、臼井前理事長が学務課職員に指示し、
名簿で上位の5人を飛ばす形で、特定の受験生に電話連絡し、繰り上げ合格させていた。この5人はいずれも女性だった疑いがある。
驚くのは明らかな「性別」や「学歴」による「差別」を、
「調整」という言葉で正当化してきた大学側の意識だ。
一般常識とはかけ離れた「医学部ムラ」の論理は東京医大に限ったものではないと
社会学者・新雅史氏は指摘する。
また、こうした不祥事が私立大学で続いていることについて、国立大学に比べ、私立大学が圧倒的に経営のガバナンスが厳しく、
結果、理事長など経営陣の意向が入学試験にも反映されやすいからではないかと推測する。
「最先端の医療研究を主体とする国立大学に比べ、
大学経営の中心が臨床である私立大学では、理事会の意向を忖度し、
忠実に働いてくれる研修医が必要不可欠。
卒業後も大学病院で即戦力として長く働いてくれる人材を探している。
また手術を手がける外科を中心とする医局は、いまだに男の戦場という認識が強い。
老舗の医科大ほど手術室における女性の割合を気にする傾向がある」(新氏)
393:卵の名無しさん
18/11/02 18:51:13.22 3zC36uTy.net
広島県の福山友愛病院で、
患者の病状と関係ない薬を大量に投与した。
しかも、意図的にです!
国は違うがドイツの場合
↓
裁判所で開かれた公判で、患者100人を殺害した罪を認めた。これでこの事件は、同国で戦後最悪級の連続殺人事件となった。
起訴されたのは、ドイツ北部デルメンホルストとオルデンブルクの病院で看護師をしていたニルス・ヘーゲル受刑者(41)。
当時勤務していたドイツ北部の2つの病院で2000~2005年にかけ、34~96歳の患者を殺害したことを認めた。
同受刑者は自分の蘇生措置の腕を同僚に見せびらかす目的や、退屈しのぎの目的で、
↓↓↓↓↓
患者に処方されていない薬を投与していたとされる。
↑↑↑↑↑
ドイツの場合。まあ日本は日本だが・・・
大口の病院は看護士の単独犯だったわけで捕まったが・・・
福山友愛病院は・・・・・
URLリンク(youtu.be)
394:卵の名無しさん
18/11/03 00:17:51.85 Ipwzsvmc.net
>>373
P(短軸探索)が先、Q(長軸探索)、同時 の配置を表示するスクリプトを書いてみた。
数値を変えて実行できる。
m:短軸 n:長軸 k:宝の数
URLリンク(tpcg.io)
treasure_map(2,3,2)
P1st Q1st even
4 5 6
P (= long axis searcher) finds first.
■ ■ □
□ ■ ■
■ ■ ■
□ □ ■
■ ■ ■
□ ■ □
■ ■ ■
■ □ □
395:卵の名無しさん
18/11/03 00:18:49.30 Ipwzsvmc.net
Q (= short axis searcher) finds first.
■ □ □
■ ■ ■
■ □ ■
■ □ ■
■ □ ■
■ ■ □
■ ■ □
■ □ ■
■ ■ □
■ ■ □
Draw.
□ □ ■
■ ■ ■
□ ■ □
■ ■ ■
□ ■ ■
□ ■ ■
396:卵の名無しさん
18/11/03 00:26:24.89 Ipwzsvmc.net
真の論理力は、プログラミング技能で試されるといっても過言ではない。
論理的に完璧でなければプログラムは動かない代物だからだ。プログラムのバグが取り切れないうちは論理力などないに等しいと思った方がよい。
URLリンク(www.mag2.com) 👀
Rock54: Caution(BBR-MD5:7bff9ed63942b4cd01610d20b2c06e65)
397:卵の名無しさん
18/11/03 09:08:04.90 Ipwzsvmc.net
血液型でもいいんだが、このスレの趣旨に合わせた問題。
ド底辺シリツ医大では正規入学:縁故入学:裏口入学:ガチ裏入学の比率が4:3:2:1であったとする。
正規入学の学生を一人みいだすために必要な学生数の期待値はいくらか?
398:卵の名無しさん
18/11/03 09:42:42.39 Ipwzsvmc.net
ド底辺シリツ医大が寄付金の多寡で合格させているのを隠蔽するために裏口受験生用にクジを導入した。
クジは合格、小吉、中吉、大吉の4枚がある。
合格を引いたら寄付金は0、
小吉は寄付金を1000万を支払って合格がでるまでクジを引き続ける
中吉は寄付金を2000万を支払って合格がでるまでクジを引き続ける
大吉は寄付金を3000万を支払って合格がでるまでクジを引き続ける
途中で中止することはできない
裏口受験生が合格するのに必要な寄付金の期待値を求めよ。
399:卵の名無しさん
18/11/03 09:46:01.38 Ipwzsvmc.net
>>385
>385のクジを引く確率は合格10%小吉20%中吉30%大吉40%に設定されているとする。
裏口受験生が合格するのに必要な寄付金の期待値を求めよ。
400:卵の名無しさん
18/11/03 10:08:18.46 Ipwzsvmc.net
>>386
解析解は計算式1行で終わって味気ないのでシミュレーションで算出してみた
> sim <- function(b=2/10,c=3/10,d=4/10,B=1000,C=2000,D=3000){
+ a=1-(b+c+d)
+ br=c(0,B,C,D)
+ total=x=sample(br,1,prob=c(a,b,c,d))
+ while(x!=0){
+ x=sample(br,1,prob=c(a,b,c,d))
+ total=total+x
+ }
+ total
+ }
> re=replicate(1e3,mean(replicate(1e3,sim())))
> summary(re)
Min. 1st Qu. Median Mean 3rd Qu. Max.
17583 19539 20000 20006 20470 22049
401:卵の名無しさん
18/11/03 12:09:44.04 Ipwzsvmc.net
暇つぶしネタ
ド底辺シリツ医大受験生の親に裏口コンサルタントが訪れて裏金額に2つの決め方を提示した。
A: 定額で2000万円
B: サイコロを1の目がでるまでふったときの出た目を合計した値 × 100万円、 例 2,1と続けば300万、6,5,1なら1200万円
問題(1) AとBではどちらが有利か?
問題(2) Bを選択した場合5000万円以上必要になる確率はくらか?
402:卵の名無しさん
18/11/03 16:48:28.20 Ipwzsvmc.net
pythonを学ぼうとanacondaのインストール(容量が少なくなったのでSDカードのD:ドライブへインストール)を試みたら
これと全く同じ現象が出現
URLリンク(stackoverflow.com)
結局、cドライブにインストールにしたら問題なく終了できた。時間が無駄になったよ。
403:卵の名無しさん
18/11/03 18:22:06.20 mjfFp3DY.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
404:卵の名無しさん
18/11/03 19:09:03.11 pQVTtzkC.net
医者を騙るマジキチが暴れてるスレ
URLリンク(egg.2ch.sc)
405:卵の名無しさん
18/11/04 08:20:42.41 ONJCveSN.net
数学板でのレスに笑った。
宮崎大医学部志望なんですが
URLリンク(rio2016.2ch.net)
406:卵の名無しさん
18/11/05 01:21:37.36 NIiUrAnG.net
import Data.List
firstUnavailable x = let y = 0:x in head $([1..71] \\)$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
next x = [n:x|n<-[head x+1..firstUnavailable x]]
xss = iterate (\xs->concat [next x|x<-xs]) [[1]]
isGood x = let y = 0:x in (==70)$length $intersect [1..70]$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
main = do
print [x|x<-(xss !! 6),isGood x]
407:卵の名無しさん
18/11/05 01:32:39.42 NIiUrAnG.net
import Data.List
m = 69
sub x = do -- ans=[1,4,5,15,18,27,34]
let ijk = filter (<=70).nub $ sort [i+j+k| i<-x,j<-x,k<-x]
all (\y -> elem y ijk ) [0..70]
main = do
-- print $ [(1,4,5,e,f,g,h)| e<-[0..m],f<-[e..m],g<-[f..m],h<-[g..m],sub [0,1,4,5,e,f,g,h]] -- 動作確認用
print $ [(1,c,d,e,f,g,h)| c<-[1..m],d<-[c..m],e<-[d..m],f<-[e..m],g<-[f..m],h<-[g..m],sub [0,1,c,d,e,f,g,h]]
408:卵の名無しさん
18/11/05 07:48:20.82 GY2Jc0B6.net
高樹沙耶さんが大麻所持で逮捕された時に、
「医療大麻なんてありません!」と医者の権益保護に躍起になっていた西川史子の現在(笑)
URLリンク(goo.gl)
あははははは!ざまー(笑)
ガンかエイズか?
外国行って医療大麻吸った方が良いんじゃね?(笑)
409:卵の名無しさん
18/11/05 12:36:04.20 RNxpU/sa.net
救急車2台受けて新入院2人で諭吉2枚ゲット、
買いたい医学書もないなぁ。
ネット検索でたいていことが足りる。
こういう基礎学力が問われる問題は検索じゃ無理だな。
インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は99%と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
50人を無作為抽出してこの診断キットで診断したところ40人が陽性であった。
この母集団の有病率の期待値と95%CIはいくらか?
またこの診断キットの感度の期待値と95%CIはいくらか
410:卵の名無しさん
18/11/05 23:27:05.76 BePOdKHG.net
これ無茶苦茶難問だった。
サイコロを1000回ふったとき123456の順に並ぶ目がある確率は?
411:卵の名無しさん
18/11/05 23:31:25.29 BePOdKHG.net
>>397
シミュレーションはことの他の簡単。
diseq=function(x) grepl("123456",paste(sample(6,x,rep=T),collapse=''))
mean(replicate(1e7,diseq(1000)))
412:卵の名無しさん
18/11/06 00:06:05.69 QApCAMmZ.net
>>397
分数表示するとこうなるらしい。
29894670002765580717622018953762664878384906585883227072416001286367327613872462
36176982687935042063489589508768499096701653359496110507628202771697652461481898
15867083862983975849719649919755556602395081603445217810191450771872978378492822
20929434295896536747523490245609702844253652551469152963247478735449528154607738
94897010212340239653386134088594267247751133592603591616944486751154360658915673
06079634567518636788329151870647053752023096491715169208527344676777484044463919
37957902455021334206833507250332593780127396266576899688383930527654264882692836
52659851277089138926126233125349416055818289864095231536195189335702592341232170
57986389166894334399077869160455544669456052004648395254763440350805452857924718
777798944034868387931340048957755897293189482452979797088/6^1000
≒0.021102960211841870224...
413:卵の名無しさん
18/11/06 06:58:36.31 jkx5i2bQ.net
>>397
Haskellで瞬時に算出されてしまった。
import Data.List
import Data.Ratio
ps = map head $ iterate (\x->(tail x) ++ [x!!5 - 1%(6^6) * x!!0]) [1%1,1%1,1%1,1%1,1%1,1%1-1%(6^6)]
main = do
print $ fromRational $ 1%1 - (ps !! 999)
Prelude> :main
2.110296021184187e-2
414:卵の名無しさん
18/11/06 08:59:30.07 kuGN+vf1.net
>>400
コードを反芻したら漸化式をプログラムしただけと理解できたのでRに移植。
実数計算なので誤差がでる。
f = function(N){
p=numeric()
p[1]=p[2]=p[3]=p[4]=p[5]=p[6]=1
for(n in 7:(N+1)){
p[n]=p[n-1]-p[n-6]/(6^6)
}
1-p[N+1]
}
> f(1000)
[1] 0.02110296
415:卵の名無しさん
18/11/06 09:01:29.36 kuGN+vf1.net
表示桁を増やしてみた。
> options(digits=22)
> f(1000)
[1] 0.0211029602118424364221
416:卵の名無しさん
18/11/06 09:13:01.99 kuGN+vf1.net
サイコロを10000回ふったとき123456の順に並ぶめがある確率は?
> f(10000)
[1] 0.1928556782249950529717
と算出できて気持ち( ・∀・)イイ!!。
417:卵の名無しさん
18/11/06 15:52:59.42 jkx5i2bQ.net
コインを100万回投げた。連続して表がでる確率が最も高いのは何回連続するときか?
おい、ド底辺すべての場合は書き出して検証してくれ。
418:卵の名無しさん
18/11/06 15:53:52.55 jkx5i2bQ.net
seq_dice <- function(N=100,k=5,p=1/6){
P=numeric(N)
for(n in 1:(k-1)){
P[n]=0
}
P[k]=p^k
P[k+1]=p^k+(1-p)*p^k
for(n in (k+1):(N-1)){
P[n+1] = P[n] + (1-P[n-k])* p^(k+1)
}
return(P[N])
}
seq_dice()
seq_diceJ <- function(N=100,k=5,p=1/6){ # Just k sequence
seq_dice(N,k,
419:p)-seq_dice(N,k+1,p) } seq_diceJ() # vsdJ=Vectorize(seq_diceJ) NN=10^6 kk=1:30 p=0.5 y=vsdJ(NN,kk,0.5) which.max(y) plot(kk,y,bty='l',pch=19,xlab='sequence',ylab='probability')
420:卵の名無しさん
18/11/06 17:20:32.86 jkx5i2bQ.net
ド底辺シリツ医大卒の裏口バカのカウントが終わりそうもないので、
シミュレーションのコードも書いてみた。
# simulation
mhs = function(x){ # maximum head sequence
y=paste(x,collapse='')
str="1"
if(!grepl(str,y)) return(0)
else{
while(grepl(str,y)){
str=paste0(str,"1")
}
return(nchar(str)-1)
}
}
(x=sample(0:1,100,rep=T)) ; mhs(x)
sim <- function(n){
mhs(sample(0:1,n,rep=T))
}
re=replicate(1e3,sim(10000))
hist(re,col='lightblue',main='',ann=F,freq=F)
summary(re)
which.max(tabulate(re))
max(tabulate(re))/length(re)
421:卵の名無しさん
18/11/06 17:26:34.68 jkx5i2bQ.net
コインを100万回投げた。連続して表がでる確率が最も高いのは何回連続するときか?
おい、ド底辺すべての場合は書き出しは終わったか?
最大値をとるときの確率は 0.23685397402052999 になったが あってるか?
数字を0と1を書いて数を数えるだけだぞ。
ド底辺シリツ医大へ裏口入学でも数くらいは数えられるだと?
422:卵の名無しさん
18/11/06 18:32:53.38 RiP59Lv4.net
We hold these truths to be self-evident, that all "uraguchi" are created retards,
That they are endowed by their creator with certain unalienable traits, that among these are sloth, mythomania, and the pursuit of imbecility.
That to rectify these traits, Bottom Medical Schools (BMS) are instituted for retards, deriving their just powers from what has been referred as the arbitrary donation of the parents of the rectified,
That whenever any form of the retards becomes destructive of rectification,
it is the right of the BMS to suspend or expel them, and to impose additional tuition laying its foundation on such principles and organizing its powers in such form, as to them shall seem most likely to effect the profit of BMS .
423:卵の名無しさん
18/11/06 19:57:10.48 jkx5i2bQ.net
サイコロを1万回ふって1の目が5回以上続く確率の計算結果
3477607490736723/18014398509481984
= 0.19304599534124128
おい、ド底辺実験して実証してみろ。
目の出方、6の1万乗を列挙して数えるのでもいいぞ。
424:卵の名無しさん
18/11/06 22:22:43.93 zJMHEzfx.net
Prudence, indeed, will dictate that "uraguchi" long established cannot be changed for light and transient causes;
and accordingly all experience has shown, that BMS are more disposed to suffer, while evils are sufferable, than to right themselves by abolishing the forms from which they can pursue profit.
But when a long train of misbehavior and misconduct , showing invariably the same Imbecility evinces a design to reveal themselves as absolute Bona Fide Moron ,
it is their right, it is their duty, to expel such "uraguchi", and to provide new guards for their future security.
425:卵の名無しさん
18/11/06 22:23:52.75 zJMHEzfx.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
426:卵の名無しさん
18/11/06 22:36:16.20 QApCAMmZ.net
>>411
現在も通用する、不朽の名投稿
(quote)
私は昭和の時代に
427:大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。 常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。 本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8-9割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。当の本人には面と向かっては絶対にそんなことは言わないけどね。 (unquote)
428:卵の名無しさん
18/11/07 00:35:16.23 J7bMbWmD.net
>>401
有理数(分数表示)したいのでPythonに移植
from fractions import Fraction
def dice126(N):
P=list()
for n in range(6):
P.append(1)
P.append(1-1/(6**6))
for n in range(7,N+1):
P.append(P[n-1]-P[n-6]/(6**6))
return(1-P[N])
def dice123456(N):
print(Fraction(dice126(N)))
print(" = " + str(dice126(N)))
dice123456(1000)
429:卵の名無しさん
18/11/07 00:40:09.03 J7bMbWmD.net
>>405
# 有理数表示したかったのでPythonに移植
from fractions import Fraction
def seq_dice(N,k,p):
P=list()
for n in range(k-1):
P.append(0)
P.append(p**k)
P.append(p**k + (1-p)*p**k)
for n in range (k,N):
P.append(P[n]+(1-P[n-k])*p**(k+1))
return(P[N])
def seq_diceJ(N,k,p):
return(seq_dice(N,k,p) - seq_dice(N,k+1,p))
def dice(N,k,p):
print("Over " + str(k))
print(Fraction(seq_dice(N,k,p)))
print(" = " + str(seq_dice(N,k,p)))
print("Just " + str(k))
print(Fraction(seq_diceJ(N,k,p)))
print(" = " + str(seq_diceJ(N,k,p)))
dice(10000,5,1/6)
# URLリンク(tpcg.io)
430:卵の名無しさん
18/11/07 10:06:50.57 gLVWtWKy.net
# how to deal list in R
> a=NULL
> a[3]=3
> a
[1] NA NA 3
> a[1]
[1] NA
> for(i in 1:10) a[i]=as.character(i)
> a
[1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10"
> b[2]=2
Error in b[2] = 2 : object 'b' not found
>
431:卵の名無しさん
18/11/07 10:25:51.27 gLVWtWKy.net
# hiw to deal list in python
a=[]
for i in range(10):
a.insert(i,i)
print (a)
for i in range(10):
a[i]=i*i
print(a)
b=[0]*10
for i in range(10):
b[i]=str(i)
print (b)
del a[-3:-1]
print(a)
del b[0:2]
print(b)
# Error!
c=[]
for i in range(10):
c[i]=i
print (c)
432:卵の名無しさん
18/11/07 11:05:51.19 gLVWtWKy.net
pythonのlist 扱い
URLリンク(qiita.com)
433:卵の名無しさん
18/11/07 16:11:33.42 SGvkOphc.net
これ面白かった。
>半径3の円Pの外側に接している半径1の円QはPを一周するといくら回転しますか。
数値は違うけど100円玉2個使って2回転を実験しちゃったよ。
434:卵の名無しさん
18/11/07 18:26:33.76 1XbYddpG.net
ヒャッハー!
米国ミシガン州 住民投票
嗜好用大麻 合法化成る!
URLリンク(www.sankei.com)
435:卵の名無しさん
18/11/08 13:38:29.75 PGJy3ILP.net
ド底辺シリツ医大の裏口入学調査委員会が
裏口入学は高々10%と報告したとする。
その結果の検証に100人を調査したら少なくとも4人続けて裏口入学生であった、という。
この検証から裏口入学率が10%であるか否かを有意水準1%で検定せよ
436:。 漸化式から算出すると 0.008705348 https://egg.5ch.net/test/read.cgi/hosp/1540905566/61 検証のためにシミュレーションスクリプトを書いてみた。 起こる回数の少ない事象の確率を頻度で近似するにはシミュレーション回数が増えるなぁ。 m1s <- function (x) { # max one sequence c(1,1,0,1,1,1,0,0) => 3 if(sum(x)==0) return(0) chg <- c(TRUE, diff(x) != 0) chgidx <- c(which(chg), length(x)+1) run.length <- diff(chgidx) true.length <- run.length[x[chg] == 1] max(true.length) } prob=c(9/10,1/10) re1=replicate(1e6,m1s(sample(0:1,100,rep=T,prob=prob))) summary(re1) ; hist(re1,col='pink') > mean(re1>=4) [1] 0.008875 有効数字1桁で近似wまあ、1%未満であるのは一致したな。 理論値を分数で計算したいなぁ、Pythonを勉強するかなぁ。
437:卵の名無しさん
18/11/08 14:48:30.64 PGJy3ILP.net
ピタゴラス数の問題に嵌ってしまったぜい。
Haskellの練習課題dにいいな。
Prelude> m=50
Prelude> ps=[(a,b,c)|a<-[1..m],b<-[a..floor(m^2/2-1/2)],c<-[b..2*b],a^2+b^2==c^2]
Prelude> ps !! 99
438:卵の名無しさん
18/11/08 15:25:37.14 PGJy3ILP.net
開業医スレのシリツ三法則(試案、名投稿より作成)
1.私立医が予想通り糞だからしょうがない
>スレリンク(hosp板:101番)-102
>スレリンク(hosp板:844番)-853
ID:RFPm1BdQ
>スレリンク(hosp板:869番)
>スレリンク(hosp板:874番)-875
>スレリンク(hosp板:874番)-880
>スレリンク(hosp板:882番)
ID:liUvUPgn
2.馬鹿に馬鹿と言っちゃ嫌われるのは摂理
実例大杉!
3.リアルでは皆思ってるだけで口に出してないだけ
URLリンク(ishikisoku.com)
439:卵の名無しさん
18/11/08 16:39:33.18 PGJy3ILP.net
a^2+b^2=c^2を満たす3つの整数(a<b<c)
の組み合わせのうち(3,4,5)から数えて7番目は何になるかという問題がわかりません
鈍足のRでも計算できた(ド底辺シリツ医大卒には無理だろな)
> pita[7,]
a b c
9 40 41
> pita[77,]
a b c
42 56 70
> pita[777,]
a b c
216 288 360
440:卵の名無しさん
18/11/08 17:16:56.00 PGJy3ILP.net
>>423
こういう問題はシミュレーションで確信がもてないので嫌だね。
他の人が同じ答を出したら確信に繋がる。
でも、それがド底辺シリツ医大卒だと確信にはならんよなぁ。
薬剤を無駄にしないために非定型的な希釈をするときは国立卒の薬剤師に検算依頼しているね。
ド底辺シリツ医大卒なんぞは信頼できんよ。
441:卵の名無しさん
18/11/08 22:06:01.87 dwp/H1Qd.net
高樹沙耶さんが大麻所持で逮捕された時に、
「医療大麻なんてありません!」と医者の権益保護に躍起になっていた西川史子の現在(笑)
URLリンク(goo.gl)
あははははは!ざまー(笑)
ガンかエイズか?
外国行って医療大麻吸った方が良いんじゃね?(笑)
ヒャッハー!
米国ミシガン州 住民投票
嗜好用大麻 合法化成る!
URLリンク(www.sankei.com)
ヒャッハー!
大麻株、絶好調!
URLリンク(web.smartnews.com)
442:卵の名無しさん
18/11/09 07:47:17.21 guOlTCed.net
coin5 = function(x,n=5) any(diff(cumsum(x),n)==n)
mean(replicate(1e4,coin5(rbinom(100,1,0.5))))
数学板のRのスレにgrep使ったスクリプトをあげたら
別アルゴリズムで高速化してくれた。
頻度が小さい事象のシミュレーションは実行速度が問題になるのでありがたい。
ド底辺シリツ医大の裏口入学調査委員会が
裏口入学は高々10%と報告したとする。
その結果の検証に100人を調査したら4人以上続けて裏口入学生であった、という。
この検証から裏口入学率が10%であるか否かを有意水準5%で検定せよ。
実行速度比較
system.time(mean(replicate(1e5,any(diff(cumsum(rbinom(100,1,0.1)),4)==4))))
user system elapsed
19.210 0.050 21.897
>
> system.time(mean(replicate(1e5,with(rle(rbinom(100,1,0.1)), max(lengths[wh
<e(1e5,with(rle(rbinom(100,1,0.1)), max(lengths[which(values==1)])>=4))))
user system elapsed
42.360 0.020 43.736
443:卵の名無しさん
18/11/09 09:07:29.19 eBlXnuYD.net
東京医大 まさに裏口漬だったな
444:卵の名無しさん
18/11/09 13:09:42.77 gy6LZalS.net
ヒャッハー!
アメリカ中間選挙で大麻が大勝利!
全米で大麻肯定派知事が誕生!
2019年に医療大麻が全米で解禁される見込み!
URLリンク(forbesjapan.com)
445:卵の名無しさん
18/11/09 14:58:35.10 guOlTCed.net
>>420
import numpy as np
from fractions import Fraction
def seqNp(N = 10,K = 5,p = 0.5):
if N == K: return p**K
q = 1-p
a = [0]*(N+1)
for i in range(1,K+1):
a[i] = q/(p**(i))
for i in range(K,N):
a[i+1] = 0
for j in range(0,K):
a[i+1] = a[i+1]+a[i-j]
a[i+1] = q/p*a[i+1]
P0=[0]*(N+1)
for i in range(1,N+1):
P0[i] = a[i]*p**i
del P0[0]
MP = np.zeros([K,N])
MP[0] = P0
MP[1][0] = p
for i in range(K-3,K):
MP[i][0] = 0
for k in range(1,K):
for i in range(0,N-1):
MP[k][i+1] = p*MP[k-1][i]
re = 1-np.sum(MP,axis=0)
print(Fraction(re[N-1]))
print(re[N-1])
446:卵の名無しさん
18/11/09 14:59:21.85 guOlTCed.net
URLリンク(tpcg.io)
447:卵の名無しさん
18/11/09 15:17:45.55 gy6LZalS.net
ヒャッハー!
アメリカ中間選挙で大麻が大勝利!
全米で大麻肯定派知事が誕生!
2019年に医療大麻が全米で解禁される見込み!
URLリンク(forbesjapan.com)
448:卵の名無しさん
18/11/09 15:35:34.92 fYz4ML3G.net
裏口学生を除籍処分にしないのだから、シリツ=裏口の公式が成立だと思うね。
不正入学した学生を処分しろと積極的に発言すべきなのがシリツ卒だと思うのだが、「僕も裏口入学」と自慢するするのが追随してたよなぁ。
449:卵の名無しさん
18/11/09 15:37:00.05 fYz4ML3G.net
>>426
78410807303959/9007199254740992
0.008705348364830168
になるので有意水準1%でも棄却される。
450:卵の名無しさん
18/11/09 15:50:48.30 Ix1q0Sut.net
経営者が出てきてスレ荒らしまくってるよ!
【 きくち歯科クリニック 厚木 はヤバイ!】
スレリンク(body板)
【 横浜セントラルパーク歯科 はヤバイ!】
スレリンク(body板)
451:卵の名無しさん
18/11/09 22:40:48.73 fYz4ML3G.net
ド底辺シリツ医大の裏口入学調査委員会が100人を順に調査したら連続して裏口学生と判明した最大数が4人と判明した。
これだけのデータから裏口入学確率の95%信頼区間を推定せよ。
数値解なら算出可能だな。
seq2pCI <- function(N,K){
vp=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
curve(vp(x),bty='l') ; abline(h=0.05,lty=3)
peak=optimize(vp,c(0,1),maximum=TRUE)$maximum
lwr=uniroot(function(x,u0=0.05) vp(x)-u0,c(0.01,peak))$root
upr=uniroot(function(x,u0=0.05) vp(x)-u0,c(peak,0.99))$root
c(lower=lwr,upper=upr)
}
seq2pCI(100,4)
> seq2pCI(100,4)
lower upper
0.1662099 0.5685943
と計算できて気分が( ・∀・)イイ!!
452:卵の名無しさん
18/11/09 23:32:07.65 gy6LZalS.net
西洋医学はオワコンです!
これからはガン治療もエイズ治療も医療大麻です!
ヒャッハー!
アメリカ中間選挙で大麻が大勝利!
全米で大麻肯定派知事が誕生!
2019年に医療大麻が全米で解禁される見込み!
URLリンク(forbesjapan.com)
高樹沙耶さんが大麻所持で逮捕された時に、
「医療大麻なんてありません!」と医者の権益保護に躍起になっていた西川史子の現在(笑)
URLリンク(goo.gl)
あははははは!ざまー(笑)
ガンかエイズか?
外国行って医療大麻吸った方が良いんじゃね?(笑)
ヒャッハー!
米国ミシガン州 住民投票
嗜好用大麻 合法化成る!
URLリンク(www.sankei.com)
ヒャッハー!
大麻株、絶好調!
URLリンク(web.smartnews.com)
453:卵の名無しさん
18/11/10 07:39:14.10 Gl2x1zZV.net
開業医スレに答えられる頭脳はなさそう。
このスレにもいないだろうな。
ド底辺シリツ裏口調査団が100人を順次調査した。
裏口判明人数をそのまま公表はヤバすぎる結果であったため、
連続して裏口がみつかった最大数は4人であったとだけ公表した。裏口入学人数の95%信頼区間は?
454:卵の名無しさん
18/11/10 07:42:24.87 Gl2x1zZV.net
漸化式から確率計算、それを数値的に逆関数で求めたけど
もっとエレガントな方法があるんだろうな。
455:卵の名無しさん
18/11/10 08:30:24.32 BVAJndFJ.net
期待値(平均値�
456:j・最頻値なしで信頼区間でけ求めるのはカッコ悪いなぁ。 ド底辺シリツ裏口調査団が100人を順次調査した。 裏口判明人数をそのまま公表はヤバすぎる結果であったため、 連続して裏口がみつかった最大数は4人であったとだけ公表した。 公表結果が正しいとして裏口入学人数の期待値、最頻値、及び95%信頼区間を述べよ。 スクリプトに計算追加して期待値・最頻値も表示させることにした。 seq2pCI <- function(N,K,alpha=0.05){ vp=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p)) curve(vp(x),lwd=2,bty='l') ; abline(h=0.05,lty=3) peak=optimize(vp,c(0,1),maximum=TRUE)$maximum mean=integrate(function(x)x*vp(x),0,1)$value/integrate(function(x)vp(x),0,1)$value lwr=uniroot(function(x,u0=alpha) vp(x)-u0,c(0.01,peak))$root upr=uniroot(function(x,u0=alpha) vp(x)-u0,c(peak,0.99))$root c(lower=lwr,mean=mean,mode=peak,upper=upr) } seq2pCI(100,4)
457:卵の名無しさん
18/11/10 08:35:58.90 BVAJndFJ.net
seq2pCI <- function(N,K,alpha=0.05){
vp=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
curve(vp(x),lwd=2,bty='l') ; abline(h=0.05,lty=3)
peak=optimize(vp,c(0,1),maximum=TRUE)$maximum
mean=integrate(function(x)x*vp(x),0,1)$value/integrate(function(x)vp(x),0,1)$value
lwr=uniroot(function(x,u0=alpha) vp(x)-u0,c(0.01,peak))$root
upr=uniroot(function(x,u0=alpha) vp(x)-u0,c(peak,0.99))$root
c(lower=lwr,mean=mean,mode=peak,upper=upr)
}
seq2pCI(100,4)
# URLリンク(tpcg.io)
458:卵の名無しさん
18/11/10 08:56:45.04 BVAJndFJ.net
>>440
これをPythonに移植するの敷居が高い。
Rの固有関数を使いまくり。
まあ、分数表示する必要がないし、処理速度も問題ないからモーティベーションもわかないし。
459:卵の名無しさん
18/11/10 09:12:19.51 BVAJndFJ.net
seq2pCI <- function(N,K,alpha=0.05){
vp=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
curve(vp(x),lwd=2,bty='l',xlab='Pr[head]',ylab=paste('Pr[max',K,'-head repetition]'))
abline(h=alpha,lty=3)
peak=optimize(vp,c(0,1),maximum=TRUE)$maximum
mean=integrate(function(x)x*vp(x),0,1)$value/integrate(function(x)vp(x),0,1)$value
lwr=uniroot(function(x,u0=alpha) vp(x)-u0,c(0.01,peak))$root
upr=uniroot(function(x,u0=alpha) vp(x)-u0,c(peak,0.99))$root
c(lower=lwr,mean=mean,mode=peak,upper=upr)
}
seq2pCI(100,5)
460:卵の名無しさん
18/11/10 13:47:51.57 eRQ4an/O.net
>>439
連続して裏口がみつかった最大数を変化させて、裏口入学者数の期待値、最頻値、95%信頼区間をグラフにしてみた。
URLリンク(i.imgur.com)
461:卵の名無しさん
18/11/10 15:31:51.93 eRQ4an/O.net
>>435
確率密度関数が左右対象でないから、
QuontileでなくHDI(Highest Density Interval)での95%信頼区間をpdfからcdfの逆関数を作って算出してみる。
# pdfからcdfの逆関数を作ってHDIを表示
pdf2hdi <- function(pdf,xMIN=0,xMAX=1,cred=0.95){
nxx=1001
xx=seq(xMIN,xMAX,length=nxx)
xx=xx[-nxx]
xx=xx[-1]
xmin=xx[1]
xmax=xx[nxx-2]
AUC=integrate(pdf,xmin,xmax)$value
PDF=function(x)pdf(x)/AUC
cdf <- function(x) integrate(PDF,xmin,x)$value
ICDF <- function(x) uniroot(function(y) cdf(y)-x,c(xmin,xmax))$root
hdi=HDInterval::hdi(ICDF,credMass=cred)
c(hdi[1],hdi[2])
}
462:卵の名無しさん
18/11/10 15:32:12.62 eRQ4an/O.net
seqNp <- function(N=100,K=5,p=0.5){
if(N==K) return(p^K)
q=1-p
a=numeric(N) # a(n)=P0(n)/p^n , P0(n)=a(n)*p^n
for(i in 1:K) a[i]=q/p^i # P0(i)=q
for(i in K:(N-1)){ # recursive formula
a[i+1]=0
for(j in 0:(K-1)){
a[i+1]=(a[i+1]+a[i-j])
}
a[i+1]=q/p*a[i+1]
}
P0=numeric(N)
for(i in 1:N) P0[i]=a[i]*p^i # P0(n)=a(n)*p^n
MP=matrix(rep(NA,N*K),ncol=K)
colnames(MP)=paste0('P',0:(K-1))
MP[,'P0']=P0
MP[1,'P1']=p
for(i in (K-2):K) MP[1,i]=0
for(k in 2:K){
for(i in 1:(N-1)) MP[i+1,k]=p*MP[i,k-1]
} # Pk(n+1)=p*P(k-1)(n)
ret=1-apply(MP,1,sum)
ret[N]
}
463:卵の名無しさん
18/11/10 15:38:44.97 eRQ4an/O.net
2つのサブルーチンを定義してから、
# N試行で最大K回連続成功→成功確率pの期待値、最頻値と95%HDI
# max K out of N-trial to probability & CI
mKoN2pCI <- function(N=100 , K=4 , conf.level=0.95){
pmf=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
mode=optimize(pmf,c(0,1),maximum=TRUE)$maximum
auc=integrate(pmf,0,1)$value
pdf=function(x) pmf(x)/auc
mean=integrate(function(x)x*pdf(x),0,1)$value
curve(pdf(x),lwd=3,bty='l',xlab='Pr[head]',ylab='density')
lu=pdf2hdi(pdf,cred=conf.level)
curve(pdf(x),lu[1],lu[2],type='h',col='lightblue',add=T)
re=c(lu[1],mean=mean,mode=mode,lu[2])
print(re,digits=4)
invisible(re)
}
> mKoN2pCI(100,4)
lower mean mode upper
0.1747676 0.3692810 0.3728936 0.5604309
> seq2pCI(100,4,0.05,T)
lower mean mode upper
0.1662099 0.3692810 0.3728936 0.5685943
はやり、HDIの方が信頼区間幅が小さいのがみてとれる。
まあ、四捨五入で人数計算すれば、推定裏口入学者数の数はほぼ同じ。
464:卵の名無しさん
18/11/10 17:41:54.09 eRQ4an/O.net
確率密度関数の面積からのクォンタイルでの信頼区間の計算
mKoN2pCIq <- function(N=100 , K=4 , alpha=0.05){
pmf=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
mode=optimize(pmf,c(0,1),maximum=TRUE)$maximum
auc=integrate(pmf,0,1)$value
pdf=function(x) pmf(x)/auc
curve(pdf(x),bty='l')
mean=integrate(function(x)x*pdf(x),0,1)$value
cdf=function(x) MASS::area(pdf,0,x)
vcdf=Vectorize(cdf)
lwr=uniroot(function(x)vcdf(x)-alpha/2,c(0,mode))$root
upr=uniroot(function(x)vcdf(x)-(1-alpha/2),c(mode,1))$root
c(lower=lwr,mean=mean,mode=mode,upper=upr)
}
mKoN2pCIq(100,4)
> mKoN2pCIq(100,4)
lower mean mode upper
0.1748351 0.3692810 0.3728936 0.5605172
> seq2pCI(100,4)
lower mean mode upper
0.1662099 0.3692810 0.3728936 0.5685943
465:卵の名無しさん
18/11/10 17:55:20.02 eRQ4an/O.net
> mKoN2pCI(100,4)
lower mean mode upper
0.1747676 0.3692810 0.3728936 0.5604309
> mKoN2pCIq(100,4)
lower mean mode upper
0.1748351 0.3692810 0.3728936 0.5605172
> mKoN2pCI(100,5)
lower mean mode upper
0.2445409 0.4469764 0.4589692 0.6418295
> mKoN2pCIq(100,5)
lower mean mode upper
0.2390957 0.4469764 0.4589692 0.6368895
466:卵の名無しさん
18/11/10 21:05:30.69 raIjTO8G.net
底辺に行くのは、二浪三浪して仕方なく、なんだろうなやはり
467:卵の名無しさん
18/11/10 21:16:17.54 eRQ4an/O.net
>>439
裏口入学確率pのときに100人中最大連続が4人になる確率をグラフにする。
これを標準化(グラフ下方の面積で割るだけ)して確率密度関数化して原点周りの一次モーメントで平均値を算出、
平均値周りの2次モーメントで分散を計算。
> (mean=integrate(function(x)x*pdf(x),0,1)$value)
[1] 0.369281
> (var=integrate(function(x)(x-mean)^2*pdf(x),0,1)$value)
[1] 0.009878284
から
> (sd=sqrt(var))
[1] 0.09938955
正規分布として下2.5%、上2.5%のとなるpの値を求めると
> qnorm(0.025,mean,sd)
[1] 0.174481
> qnorm(0.975,mean,sd)
[1] 0.5640809
累積密度関数の逆関数を数値的に作って求めたHighest Densiti Intervalは
> pdf2hdi(pdf,Print=F)
lower upper
0.17477 0.56043
とほどんど同じになった。
最大4人連続で続いたときの裏口入学者の期待値は37人
95%信頼区間は17人から56人という結果が得られた。
468:卵の名無しさん
18/11/11 07:13:25.44 A3mMnPP+.net
>>449
最初から推薦でシリツ専願というのがあるらしいぞ。
国立卒には無縁なのでよく知らん。
そういえば底辺シリツにも地域枠ってあるらしいね。
額に地底と刺青してほしいねw
469:卵の名無しさん
18/11/11 07:14:21.60 A3mMnPP+.net
ド底辺シリツ医大卒より恥ずかしいことはないよ。
東大卒のお言葉です。
URLリンク(ishikisoku.com)
470:卵の名無しさん
18/11/11 07:15:11.37 A3mMnPP+.net
ナマポも避けたがるのがド底辺シリツ医大卒
ナマポの症例報告です。
実例↓
517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな?
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。
471:卵の名無しさん
18/11/11 07:15:43.10 A3mMnPP+.net
裏口入学が都市伝説だとか過去の話というのは
嘘だと全国民が知りつくしたのが東京医大の事件だね。
縁故や寄付金じゃなくて試験データの改竄。
シリツ=裏口バカが確定。
まあ、当事者はずっと前に認識していたようだ。
ガチの裏と呼んでその存在を認めている。
URLリンク(imagizer.imageshack.com)
ガチの裏wwwwwww
472:卵の名無しさん
18/11/11 12:00:36.04 RT24eksu.net
>>453
ここが重要ポイントな。
>初診時に言ってた
実際に診察を受けなくても不利益を被ると看破されている聡明なナマポですなぁ。
473:卵の名無しさん
18/11/11 14:17:50.63 RT24eksu.net
あるド底辺シリツ医大の裏口入学者には背中に1から通し番号がふられている。
一学年の定員は100人である。
裏口入学5人の番号のうち最大の番号が60であった。
この学年の裏口入学者数の期待値を分数で答えよ。
お~い、ド底辺。
計算結果は5030556272103101/70368744177664 になったが、お前の母校で実験検証してくれ。
分数表示はともかく、この計算できない国立医学部卒がいたのにはびっくりした。
二項分布の足し算もできない学力でも臨床医がやれるんだよなぁ。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こすわけなんだな。
URLリンク(i.imgur.com)
474:卵の名無しさん
18/11/11 19:26:54.80 RT24eksu.net
>>450
Cで書いてみたけど10万個のコイン投げまでは計算してくれたけど
100万個にすると暴走してしまった。
数字の宣言が実に面倒。小数を整数で冪乗しようとするとエラーが出てくるのでいちいちcastしなくちゃならん。
インタープリターに慣れちゃうとこの辺りがとても面倒に感じる。
475:卵の名無しさん
18/11/11 19:28:27.77 RT24eksu.net
/*
C:\pleiades\workspace>coin5 100000 15
Over 15 heads at 100000 flips : 0.782609
Just 15 heads at 100000 flips : 0.248904
*/
#include <stdio.h>
#include<stdlib.h>
#include<math.h>
#define p 0.5
double flip(long N,double k){
double P[N];
long i;
for(i = 0; i < k-1; i++) {
P[i] = 0;
}
P[(int)k-1] = pow(p,k);
P[(int)k] = pow(p,k) + (1-p)*pow(p,k);
for(i = (int)k; i < N; i++){
P[i+1] = P[i] + (1-P[i-(int)k])*pow(p,(double)(k+1));
}
return P[N];
}
int main(int argc,char *argv[]){
long N = atol(argv[1]);
double k = atof(argv[2]);
double PN =flip(N,k);
double PNJ = flip(N,k) - flip(N,k+1);
printf("Over %d heads at %d flips : %f\n", (int)k ,N ,PN);
printf("Just %d heads at %d flips : %f\n", (int)k ,N ,PNJ);
return 0;
}
476:卵の名無しさん
18/11/11 20:36:53.59 VtTy0gP3.net
>>456
追加の設問を忘れてた。
この学年には何人の裏口入学者数の95%信頼区間を整数で述べよ。
477:卵の名無しさん
18/11/11 20:42:06.75 VtTy0gP3.net
>>459
>>456
日本語修正w
この学年には何人の裏口入学者がいると推論されるか。
その人数の95%信頼区間を整数で述べよ。
478:卵の名無しさん
18/11/11 22:25:33.05 RT24eksu.net
>>458
doubleをlong doubleにしても100万コイントスは計算できなかった。
pythonでは時間がかかったけど、1億コイントスの確率を分数計算してくれた。
479:卵の名無しさん
18/11/11 22:44:07.45 UCMkX/If.net
Last but not least, three laws of Do-Teihen(lowest-tier) Medical School, currently called Gachi'Ura by its graduates.
It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater they had gone through.
The Do-Teihen graduates are so ashamed of having bought their way into the exclusively lowest-tier medical school
that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
480:卵の名無しさん
18/11/11 23:15:47.97 RT24eksu.net
>>462
お~い、ド底辺。
計算結果は5030556272103101/70368744177664 になったが、お前の母校で実験検証してくれ。
あるド底辺シリツ医大の裏口入学者には背中に1から通し番号がふられている。
一学年の定員は100人である。
裏口入学5人の番号のうち最大の番号が60であった。
この学年の裏口入学者数の期待値を分数で答えよ。
この学年には何人の裏口入学者がいると推論されるか。
その人数の95%信頼区間を整数で述べよ。
481:卵の名無しさん
18/11/11 23:41:15.05 RT24eksu.net
(1)ド底辺シリツ裏口調査団が100人を順次調査した。
裏口判明人数をそのまま公表はヤバすぎる結果であったため、
連続して裏口がみつかった最大数は4人であったとだけ公表した。
(2)あるド底辺シリツ医大の裏口入学者には背中に1から通し番号がふられている。
一学年の定員は100人である。
偶々みつかった裏口入学5人の番号のうち最大の番号が60であった。
この二つの情報から、裏口人数の期待値と95%信頼区間を求めてみよう。
離散量なので単に足し算すれば求まる。
vsJ=Vectorize(function(p) seqNpJ(N=100,K=4,p))
par(mfrow=c(2
482:,2)) curve(vsJ,bty='l') n=60:100 plot(n/100,vsJ(n/100),bty='l',type='h') pmf=choose(60-1,5-1)/choose(n,5) plot(n,pmf,bty='l',type='h') plot(n,vsJ(pp)*pmf,type='h',bty='l') pdf=vsJ(pp)*pmf/sum(vsJ(pp)*pmf) > n[HDIofGrid(pdf)$`indices`] [1] 60 61 62 63 64 65 66 (1)(2)を前提にすると裏口入学の学生数の95%信頼区間は60~66人と算出できた。
483:卵の名無しさん
18/11/11 23:46:55.27 RT24eksu.net
人数の期待値は
> sum(n*pdf)
[1] 61.87202
484:卵の名無しさん
18/11/12 00:46:14.52 G7wSOjAG.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
485:卵の名無しさん
18/11/12 17:06:21.06 gzeh1x/+.net
>>463
f <- function(Mmax,M=5,Nmax=100){
n=Mmax:Nmax
pmf=choose(Mmax-1,M-1)/choose(n,M)
pdf=pmf/sum (pmf)
mean=sum(n*pdf)
upr=n[which(cumsum(pdf)>0.95)[1]]
lwr=Mmax
c(lwr,mean,upr)
}
x=5:100
y=sapply(x,function(n) f(n)[2])
z=sapply(x,function(n) f(n)[3])
plot(x,y,type='l',lwd=2,bty='l')
segments(x,x,x,z,col='gray')
cbind(x,y,z)
486:卵の名無しさん
18/11/12 22:53:02.12 3/93pSbh.net
Last but not least, three laws of Do-Teihen(lowest-tier) Medical School, currently called Gachi'Ura by its graduates.
It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater they had gone through.
The Do-Teihen graduates are so ashamed of having bought their way into the exclusively lowest-tier medical school
that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
487:卵の名無しさん
18/11/13 09:02:16.20 UTrDDuDZ.net
本日、東京医科大学が記者会見を行い、2017年、2018年度の入試の際に不当に不合格と判定された101人に意向を確認した上で、追加合格者を認定することを発表しました。
101人の選考方法
2017年度と2018年度の入試において、繰り上げ合格とされた者よりも、上位の成績だったにもかかわらず不合格とされた受験生を調べたところ、2017年度は32人(男16人、女16人)、2018年度は69人(男18人、女53人)だったとのこと。
医学部の定員は120名
今回の問題とは全く関係なく、東京医科大学の医学部の1学年の定員は120名となっている。
URLリンク(www.savingdoctor.com)
すると2年間で101/240が裏口入学と認めたことになるのではないか?
問 : 1学年あたり裏口入学者数95%信頼区間を求めよ、また、99%信頼区間も求めよ。
488:卵の名無しさん
18/11/13 09:03:43.61 UTrDDuDZ.net
>>468
ド底辺スレなのでちゃんと原文の日本語もつけろ。
Last but not least, three laws of Do-Teihen(lowest-tier) Medical School, currently called Gachi'Ura by its graduates.
最後にド底辺医大の三法則を掲げましょう。
1: It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。
2: The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater they had gone through.
ド底辺シリツ医大卒は恥ずかしくて、学校名を皆さま言いません。
3: The Do-Teihen graduates are so ashamed of having bought their way into the exclusively lowest-tier medical school
that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
ド底辺特殊シリツ医大卒は裏口入学の負い目から裏口馬鹿を暴く人間を偽医者扱いしたがる。
489:卵の名無しさん
18/11/13 09:05:12.90 UTrDDuDZ.net
本日、東京医科大学が記者会見を行い、2017年、2018年度の入試の際に不当に不合格と判定された101人に意向を確認した上で、追加合格者を認定することを発表しました。
101人の選考方法
2017年度と2018年度の入試において、繰り上げ合格とされた者よりも、上位の成績だったにもかかわらず不合格とされた受験生を調べたところ、2017年度は32人(男16人、女16人)、2018年度は69人(男18人、女53人)だったとのこと。
医学部の定員は120名
今回の問題とは全く関係なく、東京医科大学の医学部の1学年の定員は120名となっている。
URLリンク(www.savingdoctor.com)
すると2年間で101/240が裏口入学と認めたことになるのではないか?
問 : 1学年あたり裏口入学者数95%信頼区間を求めよ、また、99%信頼区間も求めよ。
99%信頼区間を計算すると、上限は50%を超えちゃうね。
シリツ医大=裏口 でいいよな?
490:卵の名無しさん
18/11/13 14:15:43.70 UTrDDuDZ.net
ド底辺シリツ医大にある年の裏口受験票があることが判明した。
枚数は1000枚以下であることは判明しているが、その数を推定したい。
裏口入学調査員が無作為に10枚選んで番号を記録して元に戻した。
別の裏口入学調査員が無作為に20枚選んで番号を記録した。
二人の調査員の記録した番号を照合すると3枚の番号が一致していた。
この情報からこの年の裏口受験証の数の最頻値、期待値とその95%信頼区間を求めよ。
491:卵の名無しさん
18/11/13 14:58:16.38 EZEaUTrv.net
Last but not least, three laws of Do-Teihen(lowest-tier) Medical School, currently called Gachi'Ura by its graduates.
It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater they had gone through.
The Do-Teihen graduates are so ashamed of having bought their way into the exclusively lowest-tier medical school
that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
492:卵の名無しさん
18/11/13 15:34:13.41 UTrDDuDZ.net
>>472
速攻で数学板で正解されてしまった。
信頼区間を求めよを追加すると、レスがつきにくいので外したのがよかったな。
rm(list=ls())
g=10 # b=N
493:-g s=20 x=3 Nmax=1000 N2p= function(N){ if(N < g + s - x) return(0) b = N-g dhyper(x,g,b,s) # choose(g,x)*choose(b,s-x)/choose(N,s) } pmf=Vectorize(N2p) # HDI source('hdi.R') N=1:Nmax y=pmf(N) plot(N,y,bty='l',type='l',ylab='probability') z=y/sum(y) # pdf 相当 plot(N,z,type='h',col='gray',ylab='density',bty='l') sum(N*z) min(HDIofGrid(z)[[1]]) ; max(HDIofGrid(z)[[1]]) # quantile zc=cumsum(z) min(which(zc>0.025)) ; min(which(zc>0.975))
494:卵の名無しさん
18/11/13 15:39:05.37 UTrDDuDZ.net
シミュレーション結果と一致したので答には自信があったのだが、数学板でのレスと一致したので確信がもてた。
Rのスクリプトはこれ
Nmin=g+s-x
sub <- function(){
n=sample(Nmin:Nmax,1)
a=sample(n,g)
b=sample(n,s)
if(sum(a %in% b)==x) return(n)
}
sim <- function(){
r=sub()
while(is.null(r)) r=sub()
return(r)
}
re=replicate(1e5,sim())
summary(re) ; hist(re,freq=F,col="lightblue",main='',xlab='N')
495:卵の名無しさん
18/11/13 15:44:29.15 UTrDDuDZ.net
>>475
> summary(re) ; hist(re,freq=F,col="lightblue")
Min. 1st Qu. Median Mean 3rd Qu. Max.
28.00000 72.00000 107.00000 149.49666 174.00000 1000.00000
> MODE <- function(x) {
+ dens <- density(x)
+ mode_i <- which.max(dens$y) # densityの頂点となるindex
+ mode_x <- dens$x[mode_i]
+ mode_y <- dens$y[mode_i]
+ mode_x
+ }
> MODE(re)
[1] 68.90811110786575
496:卵の名無しさん
18/11/13 16:51:34.16 UTrDDuDZ.net
東京医大、本来合格者入学許可へ 今年の受験生50人
2018年10月25日 02時06分
東京医科大=8月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生50人に対し、来年4月の入学を認める方針を固めたことが24日、関係者への取材で分かった。
昨年の本来合格者19人については、難しいとの意見が出ているもようだ。東京医大は50人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLリンク(www.nishinippon.co.jp)
本日、東京医科大学が記者会見を行い、2017年、2018年度の入試の際に不当に不合格と判定された101人に意向を確認した上で、追加合格者を認定することを発表しました。
101人の選考方法
2017年度と2018年度の入試において、繰り上げ合格とされた者よりも、上位の成績だったにもかかわらず不合格とされた受験生を調べたところ、2017年度は32人(男16人、女16人)、2018年度は69人(男18人、女53人)だったとのこと。
URLリンク(www.savingdoctor.com)
19人50人と言ってたのが、32人69人に増えてんじゃん。
497:卵の名無しさん
18/11/13 17:23:07.50 VaMoMifM.net
>>477
どっちにしてもそんなの氷山の一角
過去からの累積入学者なんて膨大な人数なんじゃねえの?
と言うか、元々高額な授業料を払えない受験生を排除した入学試験を
実施している時点で裏口みたいなものだから、元々大学側はそういうこと
に対しては鈍感というか、あまり気にしてなかったんじゃねえの?
国が大学の補助金が云々言い出したので慌てふためいているだけで。
498:卵の名無しさん
18/11/13 18:55:32.34 UTrDDuDZ.net
>>478
ド底辺シリツ医大卒の裏口バカの言い分をご覧ください。
東京医大の事件で点数加算の不正まで発覚しているのに、
>裏?
>そんなのあるかよばか
だそうです。
URLリンク(imagizer.imageshack.com)
499:卵の名無しさん
18/11/13 18:58:14.37 UTrDDuDZ.net
>>478
元は文科省の汚職なのだが、いつの間にやら男女比がどうやらこうやら。
まあ、裏口入学学生を除籍処分にしないのだから、
シリツ=裏口バカの公式は健在だよな。
500:卵の名無しさん
18/11/13 19:46:47.63 WZBh4CIg.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
501:卵の名無しさん
18/11/14 00:17:40.02 lTkyjX2F.net
ニュー速+のこれが正義だと思う。
251 名無しさん@1周年 2018/11/09(金) 03:01:30.51 ID:T61ryGbd0
不正合格者退学にして空けろよ
252 名無しさん@1周年 2018/11/09(金) 03:27:59.34 ID:Wfw1oWQ60
裏口で入ったのをやめさせろよ。
502:卵の名無しさん
18/11/14 07:58:51.21 1BnjGtIH.net
>>478
縁故や「寄付金」での加点も裏口と呼ぶならば裏口でない学生の方が少ないかもね。
503:卵の名無しさん
18/11/14 08:01:55.34 1BnjGtIH.net
>>481
単発idのこいつもド底辺シリツに裏口入学なんだろうなぁ。
ちゃんと原文もつけろよ。
現在も通用する、不朽の名投稿
(quote)
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8-9割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
(unquote)
504:卵の名無しさん
18/11/14 09:23:05.13 8Fqyq1WG.net
コピペだがこれに納得。
私立大学では、誰を入学させるかは大学の自由などと言ってる人
「大学設置基準」(学校教育法に基づく文部科学省令)の明確な違反(法令違反)だ。
2条の2 入学者の選抜は、公正かつ妥当な方法により、適切な体制を整えて行うものとする。
ちなみに、これに基づく入試の公正確保のための事務次官通達もでている。
また、学校法人は、私立学校法により、公の認可を受けた公益法人。税制面の優遇も受けている。
営利企業とは全く異なる存在。さらに、東京医大を例にとれば、平成28年で20億以上の補助金が配布されている。
一人当たり年間300~400万円、6年で2000万円前後の税金が使われている。
こうした前提を無視し、適当な言説をばらまくのは控えるべき。
(参考資料)
私立大学における入学者選抜の公正確保等について
URLリンク(www.mext.go.jp)
505:卵の名無しさん
18/11/15 17:51:53.99 u639IIxL.net
ABC
ACB
BAC
BCA
+CAB
--------
3123
A,B,Cは?
Haskell
[(a,b,c)|a<-[1..9],b<-[1..9],c<-[1..9], (a+a+b+b+c)*100+(b+c+a+c+a)*10+(c+b+c+a+b)==3123]
Python
n = range(1,10)
for A in n:
for B in n:
for C in n:
if (A+A+B+B+C)*100+(B+C+A+C+A)*10+(C+B+C+A+B)==3123:
print( [A,B,C] )
R
n=1:9
N=c(100,10,1)
for(A in n){
for(B in n){
for(C in n){
S=sum(c(A,B,C)*N+c(A,C,B)*N+c(B,A,C)*N+c(B,C,A)*N+c(C,A,B)*N)
if(S==3123) ans=c(A,B,C)
}
}
}
ans
506:卵の名無しさん
18/11/15 18:36:42.99 u639IIxL.net
C
#include <stdio.h>
int A,B,C,i;
int main() {
for(A=1;A<10;A++){
for(B=1;B<10;B++){
for(C=1;C<10;C++){
i=(A+A+B+B+C)*100+(B+C+A+C+A)*10+(C+B+C+A+B);
if(i==3123) printf("%d %d %d\n",A,B,C);
}
}
}
}
507:卵の名無しさん
18/11/15 20:59:32.76 ziKbMPPq.net
Last but not least, three laws of Do-Teihen(lowest-tier) Medical School, currently called Gachi'Ura by its graduates.
It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater they had gone through.
The Do-Teihen graduates are so ashamed of having bought their way into the exclusively lowest-tier medical school
that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
508:卵の名無しさん
18/11/16 06:15:26.91 DH6oVlEk.net
でも私立医大卒の医者って良いよな~。
だって親がお金持ちじゃなかったら本当は
医者じゃなかったかもしれない連中なんだぜw
509:卵の名無しさん
18/11/16 07:19:25.14 CpQxI30j.net
[(a,b,c,d,e,f,g)|a<-[1..9],b<-[1..9],c<-[1..9],d<-[1..9],e<-[1..9],f<-[1..9],g<-[1..9],a/b*c/d=(10*e+b)/g]
510:卵の名無しさん
18/11/16 08:32:47.11 OhMUxuM/.net
1問目は1から9を多くて1回づつ使って等式を完成させる
URLリンク(i.imgur.com)
□/□ * □/□ = □□/□
(a/b)*(c/d) = (10*e+f) /g
左辺の分数は互換なのでa>cとして
Prelude> map f r
[[7,2,3,6,1,4,8],[7,6,3,2,1,4,8],[8,2,3,6,1,4,7],[8,2,7,4,6,3,9],[8,2,7,6,1,4,3],[8,4,7,2,6,3,9],[8,4,7,6,2,1,9],[8,6,3,2,1,4,7],[9,2,7,4,6,3,8],
[9,2,8,4,6,3,7],[9,4,7,2,6,3,8],[9,4,7,6,2,1,8],[9,4,8,2,6,3,7],[9,4,8,6,2,1,7],[9,6,7,4,2,1,8],[9,6,8,4,2,1,7]]
お~い、ド底辺、他に答がないか検算しておけ。
511:卵の名無しさん
18/11/16 12:28:38.79 OhMUxuM/.net
>>491
Rでのスパテゲティーコード
実用的な時間で算出できるかなぁ?
# □/□ * □/□ = □□/□
# a/b * c/d == ef/g (c>a)
F=function(a,b,c,d,e,f,g){
a!=b & a!=c & a!=d & a!=e & a!=f & a!=g &
b!=c & b!=d & b!=e & b!=f & b!=g & c!=d &
c!=e & c!=f & c!=g & d!=e & d!=f & d!=g & e!=f & e!=g & f!=g & (a/b)*(c/d)==(10*e+f)/g
}
n=1:9
ans=NULL
for(a in n){
for(b in n){
for(c in a:9){
for(d in n){
for(e in n){
for(f in n){
for(g in n){
if(F(a,b,c,d,e,f,g)){
ans=rbind(ans,c(a,b,c,d,e,f,g))}}}}}}}}
ans
512:卵の名無しさん
18/11/16 12:44:40.76 OhMUxuM/.net
分母に1を許さないことにすると
> ans
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 3 2 7 6 1 4 8
[2,] 3 2 8 6 1 4 7
[3,] 3 6 7 2 1 4 8
[4,] 3 6 8 2 1 4 7
[5,] 7 2 8 4 6 3 9
[6,] 7 2 9 4 6 3 8
[7,] 7 4 8 2 6 3 9
[8,] 7 4 9 2 6 3 8
[9,] 7 4 9 6 2 1 8
[10,] 7 6 8 2 1 4 3
[11,] 7 6 8 4 2 1 9
[12,] 7 6 9 4 2 1 8
[13,] 8 2 9 4 6 3 7
[14,] 8 4 9 2 6 3 7
[15,] 8 4 9 6 2 1 7
[16,] 8 6 9 4 2 1 7
513:卵の名無しさん
18/11/16 13:01:06.28 OhMUxuM/.net
分母に整数を許すと34個あるな。(左辺の分数は区別しない)
514:卵の名無しさん
18/11/16 14:41:49.22 OhMUxuM/.net
>>494
-- □/□ * □/□ = □□/□
r0 = [[a,b,c,d,e,f,g]|a<-[1..9],b<-[2..9],c<-[(a+1)..9],d<-[2..9],e<-[1..9],f<-[1..9],g<-[2..9],(a/b)*(c/d)==(10*e+f)/g,a/=b,a/=c,a/=d,a/=e,a/=f,a/=g,b/=c,b/=d,b/=e,b/=f,b/=g,c/=d,c/=e,c/=f,c/=g,d/=e,d/=f,d/=g,e/=f,e/=g,f/=g]
r1 = [[a,b,c,d,e,f,g]|a<-[1..9],b<-[1..9],c<-[(a+1)..9],d<-[1..9],e<-[1..9],f<-[1..9],g<-[2..9],(a/b)*(c/d)==(10*e+f)/g,a/=b,a/=c,a/=d,a/=e,a/=f,a/=g,b/=c,b/=d,b/=e,b/=f,b/=g,c/=d,c/=e,c/=f,c/=g,d/=e,d/=f,d/=g,e/=f,e/=g,f/=g]
f x = map floor x
ans0 = map f r0
ans1 = map f r1
f0 = do
putStrLn "One is NOT permitted as denominator"
putStrLn $ show(length(ans0)) ++ " combinations"
print ans0
f1 = do
putStrLn "permissive of 1 as denominator"
putStrLn $ show(length(ans1)) ++ " combinations"
print ans0
main = do
f0
f1