18/10/21 21:20:01.39 Z8EDA9Rx.net
# 縦3マス、横4マスの12マスのうちランダムに選ばれた2マスにそれぞれ宝が眠っている。
# AEIBFJ…の順で縦に宝を探していく方法をとるP君と、
# ABCDEFGH…の順で横に宝を探していく方法をとるQ君が、同時に地点Aから探索を開始した。
# どっちの方が有利？
#
# ABCD
# EFGH
# IJKL
y=LETTERS[1:12]
z=matrix(y,ncol=4,byrow=T)
as.vector(z) # P
as.vector(t(z)) # Q(=y)
x=c(1,1,rep(0,10))
is.P1st <- function(){
Q=sample(x)
z=matrix(Q,ncol=4,byrow=T)
P=as.vector(z)
which.max(P) < which.max(Q)
}
k=1e3
re=replicate(k,mean(replicate(k,is.P1st())))
summary(re)
802:卵の名無しさん
18/10/21 21:45:36.09 Z8EDA9Rx.net
> x=c(1,1,rep(0,10))
> PQ <- function(){
+ Q=sample(x)
+ z=matrix(Q,ncol=4,byrow=T)
+ P=as.vector(z)
+ c( even=which.max(P) == which.max(Q),
+ p1st=which.max(P) < which.max(Q),
+ q1st=which.max(P) > which.max(Q))
+
+ }
> k=1e6
> re=replicate(k,PQ())
> mean(re['even',]) ; 13/(26+27+13)
[1] 0.197025
[1] 0.1969697
> mean(re['p1st',]) ; 26/(26+27+13)
[1] 0.393803
[1] 0.3939394
> mean(re['q1st',]) ; 27/(26+27+13)
[1] 0.409172
[1] 0.4090909
803:卵の名無しさん
18/10/22 02:17:14.19 YJqFFfdN.net
>>748
一般解
treasure <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
c(P1st=sum(re<0),Q1st=sum(re>0),even=sum(re==0))
}
treasure(3,4,2)
treasure(5,10,3)
804:卵の名無しさん
18/10/22 22:25:36.36 JjE8cyU0.net
|￣￣￣￣￣￣| ..--―-- __
| 国 | .:.:.:.:.:.:.:.:..:.:.:.:.:.:.:.:.:.｀.
| 試 |.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:＼
| 浪 |:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:ヽ
| 人 |.:.:./:.:.:.:.:.:.:.:.:.:.:ヽ:.:.:.:.:.:.:..:..:.:ﾊ
| ！？ |.:.:.|:.:.:.:.:.:.::.:.:.:.:.:.|:.:.:.:.:.:.:.:.:.:.:.:.}
|＿＿＿＿＿＿|:.:.:l.:.:l:.:.:.:.:l.:.:.:.:.:ｌ:.:.:.:.:.:ｌ:.:.:.:.:.:j
l:.:.:.:.i:.:.:.:ヽl∧ﾊ＼:.:.l:.:/∨ハ:.:/l:.:.:.:ｌ:.:.ｌ
|:.:.:.:.l:.:.:.:l斗―ミ､ ヽ!/ ,. -―ﾐ､ !:.:./:.:.:ｌ
iﾄ､:.:.|:.:.:.|（ ・）` '´（・ ）` l:.ノ:.:.:.:.:l
ﾉ:.:ﾊ:.:.:l ` 一' , ｰ一'´ﾉ:.:.:.:.|.:.:人
／:.:.:.У:人ι ﾑｲ:.:.:l:.:.:.:./ヽ
(/i:.:{::.ｲ:.:.:ﾊ.丶 -―- ／:.:l:.:.:.:l.:.:.:八
. ﾚ'´}:.:.:./:.:.}ﾘ}＞ ＿_ .イ人:.:ﾊ:.:.:|:.:ﾘ}
__)ノ ﾄソ／ ＼_ _／ ＼ ﾘ 乂へ
イ /##｛ ＿/∨＼＿ }### i ＼
/ /###レイ. ヽ人／ ヽ/###.ｌ ヽ
. ｌ /####l / ヽ /####ｉ }
805:卵の名無しさん
18/10/22 22:27:47.24 JjE8cyU0.net
国試浪人の事務員については専用スレがあるのでそちらへ
ｽﾚﾘﾝｸ(hosp板:29番)
806:卵の名無しさん
18/10/22 22:28:47.06 JjE8cyU0.net
妬みマンチクリン がどういう意味なんだろうと思ってググってみたら
こんなスレが出てきました
自称医科歯科卒が専門医をねたむスレ [転載禁止]©2ch.net
ｽﾚﾘﾝｸ(hosp板)
このスレを辿ってみたら白黒コピーの医師免許の一部が
大部分を隠してアップされてましたが
ｽﾚﾘﾝｸ(hosp板:437番)
これって病院で事務用に保管しているコピーですよね
807:卵の名無しさん
18/10/23 03:47:27.97 0jwkOV+G.net
あれ？
朝鮮殺戮殺人学会が警察にテロ工作員を送り込み、
なんと
警察から犯罪ライセンスを与えられ、
監禁罪、薬物大量投与テロ、
などやり、
テロ工作拠点だとバレた
福山友愛病院については？
福山友愛病院 事件 でググれば？
URLﾘﾝｸ(sp.nicovideo.jp)
808:卵の名無しさん
18/10/23 07:05:26.69 RnIgsz+6.net
>>751
スクリプトの最適化とか示せばいいのに
どれでもいいから掲載してあるRのコードをCで高速化してみ。
809:卵の名無しさん
18/10/23 07:11:17.37 RnIgsz+6.net
>>751
専門医スレで誰も正解出せなかった、
これを解くコードを書いて仮想国試浪人程度の頭脳があることを示してくれ。
door to balloon timeが短縮するのに必要な計算。
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？
810:卵の名無しさん
18/10/23 07:42:35.45 RnIgsz+6.net
国立卒ならこれくらい答えてシリツとの差をみせてほしいね。
インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は９９％と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
５０人を無作為抽出してこの診断キットで診断したところ４０人が陽性であった。
この母集団の有病率の期待値と９５％ＣＩはいくらか？
またこの診断キットの感度の期待値と９５％ＣＩはいくらか
811:卵の名無しさん
18/10/23 11:55:53.89 7flOAmF/.net
digi = function(x){ # 1000 -> 4 , 999 -> 3
n=ceiling(log10(x))
ifelse(10^n==x,n+1,n)
}
n2a <- function(num){ # nmu to array 122 -> c(1,2,2)
N=10
r=num%%N
q=num%/%N
while(q>0){
r=append(q%%N,r)
q=q%/%N
}
return(r)
}
one2n <- function(x){ # 121 -> 13
a=n2a(x)
k=digi(x)
p=2^((k-1):0)
sum(a*p)
}
x=1212121212
one2n(x)
812:卵の名無しさん
18/10/23 13:18:37.27 cl3+ifwu.net
このスレは事務員が日がな一日妄想を垂れ流し
見物人たちがそれを見てフルボッコにするスレである
事務員 とは？
同一者の別名として、薄汚いジジイ、国試浪人、統計野郎、自称医科歯科、がある
自称医科歯科卒、実際は九州の駅弁国立出身で、卒業後は実家の東京に帰り国試浪人となり23年間経過、家庭教師バイトなどを経て現在は事務員、とされている
本人は、勤務医でたまに開業医の手伝いと内視鏡バイト、専門医なし、と主張しているが、連日連夜の異常な書き込み数の多さからは、勤務医とは考え難く、彼の職業がたとえ何であったとしても、人としてド底辺なことは間違いがない
自ら事務員であると告白したこともある
ｽﾚﾘﾝｸ(hosp板:108番)-109
自作自演も頻繁に行なっている
事務仕事の際に手にした他人の給与明細や同窓会費振込票を盗撮し、自分のものとしてアップロードしたこともある、犯罪なんじゃねーの？
統計絶対主義でマウントを取ろうとするが、大学卒業後の人生の大部分においてそれが何の役に立つのかはなはだ疑問である
統計に基づいた重病の治療成績を患者に説明しても、患者の希望する治療は患者の事情によって異なることは一定年数以上の臨床経験があるものなら誰でも知っていることだ
それが理解できない彼は、彼自身が素性を明らかにしているようなものである
コンプ薬屋同様に、私立医への執拗な誹謗中傷を数年余にわたって続けている
裏口批判をしているが、実際にそれが可能な資産家の子弟がいるであろう右派スレには近寄らない、一度書き込んだ形跡があったが指摘したら別人だと言って逃げた
弱いものが多そうなスレしか狙わないのもコンプ薬屋と共通する異常者の特徴である
劣等コンプレックスの塊なんだろうな
自分の主張に対して都合の悪い突込みが入ると、相手を私立医と決めつけて、さてはシリツだな、裏口バカ、というセリフでワンパターンに返すことが多い
国試問題を挙げて簡単だとほざいてますが、本番で通らなければお話になりません
いくらネットでわめいていても、医師国試の本番で通らなければお話になりません
大事なことなので２回言いました
813:卵の名無しさん
18/10/23 14:23:32.78 7flOAmF/.net
-- 不定長整数が扱えて遅延IOなのでメモリ不足になりにくいな --
main = do
putStr $ "Type number: "
str <- getLine
let xs = concat $ iterate (\x->[1:n| n<-x] ++ [2:n|n<-x]) [[1],[2]]
n = read str
putStrLn $ concat $ map show (xs !! (n-1))
814:卵の名無しさん
18/10/23 14:54:10.02 7flOAmF/.net
>>760
コンパイルして実行したけど、１億めはメモリ不足に陥った。
815:卵の名無しさん
18/10/23 15:44:29.51 7flOAmF/.net
{1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0}
URLﾘﾝｸ(www.wolframalpha.com)
816:卵の名無しさん
18/10/23 16:49:32.31 RnIgsz+6.net
>>761
> digit12(10^8) # 1億め
12222212122221111211111112
> digit12(10^12) # 1兆め
221211122121211212112121112111111111112
Rのコードはここ
URLﾘﾝｸ(tpcg.io)
817:卵の名無しさん
18/10/23 18:48:55.68 cl3+ifwu.net
>>763 は国試に２３回も落ちた挙句
婚活にも失敗してる素人童貞で
くるくるぱーの裏口バカに
なっちゃってるのらぁあぁぁ
Ｆラン事務員の濃ゆぅぅい生ガキ汁
ド底辺の臭いが落ちないよぉ
んほおぉぉぉおぉぉ
818:卵の名無しさん
18/10/23 21:27:28.15 RnIgsz+6.net
>>763
１０の６８乗を無量大数というらしい
無量大数＋１を２進数表示できるかやってみた。
Prelude> :main
Input integer : 100000000000000000000000000000000000000000000000000000000000000000001
1110110101100011101000100011000111010100110001001111101100100111010011001010011110101010101010000110001111101110010010111101110101001000010101101100010111000100000000000000000000000000000000000000000000000000000000000000000001
さすが不定長整数を扱えるHaskell。
819:コンブ薬屋
18/10/23 23:11:54.77 aUOG+wre.net
>>765
んで、あんた、どこ卒？
さては国試浪人の事務員だな
820:卵の名無しさん
18/10/24 07:56:03.66 cJkK7XDW.net
>>766
俺、医科歯科卒。二期校時代に入学。
国立卒なのでこういう問題にも答が出せる。
インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は９９％と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
５０人を無作為抽出してこの診断キットで診断したところ４０人が陽性であった。
この母集団の有病率の期待値と９５％信用区間はいくらか？
またこの診断キットの感度の期待値と９５％信用区間はいくらか
821:卵の名無しさん
18/10/24 13:25:08.39 R65m5krb.net
>>767は国試に２３回も落ちた挙句
婚活にも失敗してる素人童貞で
くるくるぱーの裏口バカに
なっちゃってるのらぁあぁぁ
Ｆラン事務員の濃ゆぅぅい生ガキ汁
ド底辺の臭いが落ちないよぉ
んほおぉぉぉおぉぉ
822:卵の名無しさん
18/10/24 14:57:07.89 IciU1D6F.net
無職ではまともな女は相手しないだろうw
823:卵の名無しさん
18/10/24 18:00:20
824:.55 ID:wCLV17Kd.net
825:卵の名無しさん
18/10/24 18:00:40.97 cJkK7XDW.net
>>768
>国立卒なのでこういう問題にも答が出せる。
を受けて正解を投稿すれば、
ド底辺シリツ医大卒でも国立卒レベルの数学ができることを示す機会なのに
馬鹿を晒す機会にするとはド底辺シリツ医大卒の裏口ガイジらしいぜｗ
826:卵の名無しさん
18/10/24 18:11:39.40 cJkK7XDW.net
>>769
無職のナマポでさえ避けたがるのがド底辺シリツ医大卒
ナマポの症例報告です。
実例↓
517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな？
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。
827:卵の名無しさん
18/10/25 07:04:57.74 /nor6Fod.net
裁判所事務官採用試験（大阪高裁管轄） の結果
URLﾘﾝｸ(www.courts.go.jp)
男880 1次通過者224 合格者43
女736 1次通過者159 合格者87
調整の出来ない1次試験の通過率 男25.5% 女21.6％
調整し放題の面接と称する試験合格率 男19.2% 女54.7%
f <- function(ratio=1.2,man=800,woman=200,pass=100){
total=man+woman
i=0:pass
sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > ratio*(pass-i)/woman))
}
f((87/736)/(43/880),736,880,43+87)
f((87/159)/(43/224),159,224,43+87)
prop.test(c(43,87),c(880,736))
prop.test(c(224,159),c(880,736))
prop.test(c(43,87),c(224,159))
828:卵の名無しさん
18/10/25 07:38:38.53 /nor6Fod.net
> binom.test(43,43+87,224/(224+159),alt="less")
829:卵の名無しさん
18/10/25 14:30:04.77 VIrGA62+.net
さあシリツのじかんだ
URLﾘﾝｸ(cdn-ak.f.st-hatena.com)
830:卵の名無しさん
18/10/25 18:32:09.06 /nor6Fod.net
この方がすっきり。
男880 1次通過者224 合格者43
女736 1次通過者159 合格者87
x=c(880,736)
y=c(224,159)
z=c(43,87)
fisher.test(cbind(x-y,y))
fisher.test(cbind(y-z,z))
fisher.test(cbind(x-z,z))
prop.test(y,x)
prop.test(z,y)
prop.test(z,x)
831:卵の名無しさん
18/10/25 18:35:26.32 /nor6Fod.net
>>772
>初診時に言ってた
シリツ医大卒＝裏口バカと認識しているわけだな。
>>
東京医大、本来合格者入学許可へ 今年の受験生５０人
2018年10月25日 02時06分
東京医科大＝８月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
<<
URLﾘﾝｸ(www.nishinippon.co.jp)
入学者数って100人程度のだよな？
半分が裏口ってこと？
832:卵の名無しさん
18/10/25 19:31:40.10 MruIcscN.net
treasure <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
print(c(P1st=sum(re<0),Q1st=sum(re>0),even=sum(re==0)))
invisible(list(P1st=tre[,which(re<0)
833:],Q1st=tre[,which(re>0)],even=tre[,which(re==0)])) } t342=treasure(3,4,2) print(matrix(LETTERS[t342$P1st],nrow=2),quote=F) print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F) print(matrix(LETTERS[t342$even],nrow=2),quote=F) t452=treasure(4,5,2) print(matrix(LETTERS[t452$P1st],nrow=2),quote=F) print(matrix(LETTERS[t452$Q1st],nrow=2),quote=F) print(matrix(LETTERS[t452$even],nrow=2),quote=F)
834:卵の名無しさん
18/10/25 19:32:05.84 MruIcscN.net
縦nマス、横n+1マスのn(n+1)マスのうちランダムに選ばれた2マスにそれぞれ宝が眠っている。
縦1列を探し終えたらすぐ右の1列に移って宝を探していく方法をとるP君と、横1行を探し終えたらすぐ下の1行に移って宝を探していく方法をとるQ君が、同時に左上の地点から探索を開始した。
例えば、n=3の時はP君はAEIBFJCGKDHLの順で探す。Q君はABCDEFGHIJKの順で探すことになる。
ABCD
EFGH
I JK L
1つの地点を捜索するのにかかる時間は同じで、相手が1度探し終えた地点を重複して調べることも当然ある。
相手より先に宝を見つけた方を勝者とする。同時の場合は引き分けとする。
どちらの方が有利になるだろうか？
835:卵の名無しさん
18/10/25 19:48:03.86 MruIcscN.net
>>779
> t(sapply(1:20,treasure1))
P1st Q1st even
[1,] 0 0 1
[2,] 4 5 6
[3,] 26 27 13
[4,] 84 83 23
[5,] 203 197 35
[6,] 413 398 50
[7,] 751 722 67
[8,] 1259 1210 87
[9,] 1986 1910 109
[10,] 2986 2875 134
[11,] 4320 4165 161
[12,] 6054 5845 191
[13,] 8261 7987 223
[14,] 11019 10668 258
[15,] 14413 13972 295
[16,] 18533 17988 335
[17,] 23476 22812 377
[18,] 29344 28545 422
[19,] 36246 35295 469
[20,] 44296 43175 519
836:卵の名無しさん
18/10/25 20:43:08.89 VIrGA62+.net
どれ、たまには算数の問題を投稿してやろう
URLﾘﾝｸ(seiga.nicovideo.jp)
837:卵の名無しさん
18/10/25 21:18:10.86 /nor6Fod.net
>>780
4以上だと短軸探索が有利になるようにみえるな。
838:卵の名無しさん
18/10/25 22:33:55.27 MruIcscN.net
>>778
n=2
ABC
DEF
の場合
短軸方向探索Pが先に宝を発見する埋め方：４通り
> print(matrix(LETTERS[t232$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4]
[1,] C D D E
[2,] D E F F
長軸方向探索Qが先に宝を発見する埋め方：５通り
> print(matrix(LETTERS[t232$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5]
[1,] B B B C C
[2,] C E F E F
同時に宝を発見する埋め方：６通り
> print(matrix(LETTERS[t232$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] A A A A A B
[2,] B C D E F D
839:卵の名無しさん
18/10/25 23:25:11.32 0YTgsLE2.net
算数の問題です
URLﾘﾝｸ(sp.seiga.nicovideo.jp)
840:卵の名無しさん
18/10/25 23:37:44.75 0YTgsLE2.net
算数の問題です
URLﾘﾝｸ(sp.seiga.nicovideo.jp)
841:卵の名無しさん
18/10/26 00:02:56.11 +rLvZj8m.net
>>779
数学板で２つの宝を先に発見した方を勝者にしたらどうなるかというレスがあったのでプログラムを改変してみた。
treasure2 <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
max(p)-max(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
print(c(P1st=sum(re<0),Q1st=sum(re>0),even=sum(re==0)))
invisible(list(P1st=tre[,which(re<0)],Q1st=tre[,which(re>0)],even=tre[,which(re==0)]))
}
t342=treasure2(3,4,2)
#短軸方向探索Pが先に２つの宝を発見する埋め方
print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
#長軸方向探索Qが先に２つの宝を発見する埋め方
print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
#同時に２つめの宝を発見する埋め方
print(matrix(LETTERS[t342$even],nrow=2),quote=F)
842:卵の名無しさん
18/10/26 00:04:20.43 +rLvZj8m.net
> t342=treasure2(3,4,2)
P1st Q1st even
27 26 13
> #短軸方向探索Pが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B B C C C D E E E E F
[2,] E F I J K E F I J K I J K K F I J K I
[,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27]
[1,] F F G G G I I J
[2,] J K I J K J K K
> #長軸方向探索Qが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B C C C C C D D D D D
[2,] B C D G H C D G H D E F G H E F G H I
[,20] [,21] [,22] [,23] [,24] [,25] [,26]
[1,] E E F F G H H
[2,] G H G H H I J
> #同時に２つめの宝を発見する埋め方
> print(matrix(LETTERS[t342$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,] A B C D D E F G H H I J K
[2,] L L L J L L L L K L L L L
やはり、直観と違ってイーブンにはならないな。
843:卵の名無しさん
18/10/26 00:14:13.73 +rLvZj8m.net
縦５マス、横１０マスで宝が３マスに埋まっているときに
全部の宝を発見した方が勝者とすると、縦方向探索Pと横方向探索Qとして
勝者となる埋め方の場合の数は
> treasure2(5,10,3)
P1st Q1st even
9142 8832 1626
と算出できた。
部屋割りカウントできない底辺頭脳でもこれくらいは数えられるかな？
844:卵の名無しさん
18/10/26 08:17:44.20 R30QBCXi.net
算数の問題です
URLﾘﾝｸ(seiga.nicovideo.jp)
中卒国試浪人の事務員には解けないだろうな
845:卵の名無しさん
18/10/26 08:21:49.38 46gTk/Hn.net
昼休みの薬の説明会はフィブラストだった。
熱傷診療ガイドライン(www.jsbi-burn.org/members/guideline/pdf/guideline2.pdf)に引用された論文としてこんなのが商用パンフに掲載されていた。
URLﾘﾝｸ(imagizer.imageshack.com)
早期群: n1=17 m1=36.2 sd1=21.5
後期群: n2=12 m2=56.1 sd2=30.1
で早期群の方が上皮化期間が有意(t検定でp<0.05)に短いと主張している。
生データがないので標本数・平均・標準偏差から検定してみる。
ttest2 = function(n1,m1,s1, n2,m2,s2){ # ni:標本数 mi:平均 si:標準偏差
n=n1+n2-2
u=((n1-1)*s1^2+(n2-1)*s2^2)/n
t=(m1-m2)/sqrt(u/n1+u/n2)
pe=2*pt(-abs(t),n)
t=(m1-m2)/sqrt(s1^2/n1 + s2^2/n2)
n=(s1^2/n1+s2^2/n2)^2 / ( (s1^2/n1)^2/(n1-1) + (s2^2/n2)^2/(n2-1))
pu=2*pt(-abs(t),n)
p.values=c(T.test=pe,Welch.test=pu)
return(p.values)
}
> ttest2(n1,m1,sd1, n2,m2,sd2)
T.test Welch.test
0.04701090 0.06461679
等分散を仮定したｔ検定でやっと0.05未満になるが、等分散を仮定しないWelchの方法だと有意差がでない。
846:卵の名無しさん
18/10/26 08:22:44.14 46gTk/Hn.net
URLﾘﾝｸ(imagizer.imageshack.com)
のグラフをみると後期群の最大値は外れ値に見える。
後期群が n2=12 m2=56.1 sd2=30.1
を満たすようなシミュレーションデータを作成して
そのデータの最大値を除いた標本群を修正後期群として
早期群と検定（等分散を仮定するｔ検定と仮定しないWelchの検定）してｐ値を出す。
これを1万回繰り返してｐ値がどんな値になるかやってみた。
その結果は、
等分散ｔ検定でｐ値は平均 0.105(95%CI 0.0916 0.116)で有意とは言えない。
Welochの検定でｐ値は平均 0.126(95%CI 0.114 0.129)で有意とは言えない。
シミュレーションデータによる検定のRスクリプトはここに書いた。
ｽﾚﾘﾝｸ(hosp板:82番)
ガイドラインに引用されている論文すら、あんまり信用できないと思った。
“Statistics are like bikinis. What they reveal is suggestive, but what they conceal is vital.”
847:卵の名無しさん
18/10/26 08:39:18.67 46gTk/Hn.net
>>789
中卒で受けられる国試って何だろね？
その一文だけで頭が悪いのがわかる。
848:卵の名無しさん
18/10/26 11:29:32.95 46gTk/Hn.net
こういうのの実行結果を数学板に投稿すると
レスが返って来て勉強になるなぁ。
# ABCDEFGに対してなら
# BCDEFGAが一番勝率高い気がする
library(gtools)
n=7
k=2
perm=permutations(n,n)
Q=perm[1,]
np=nrow(perm)
p1st=numeric(np)
for(i in 1:np){
P=perm[i,]
tre=combn(n,k)
nt=ncol(tre)
re=numeric()
for(j in 1:nt){
re[j]=min(which(tre[1,j]==P),which(tre[2,j]==P))-
min(which(tre[1,j]==Q),which(tre[2,j]==Q))
}
p1st[i]=sum(re<0)
}
plot(p1st)
p1st[which.max(p1st)]
(p.max=which(p1st==15))
print(matrix(LETTERS[perm[p.max,]],ncol=7),quote=F)
#
849:卵の名無しさん
18/10/26 11:48:48.32 46gTk/Hn.net
宝を２個先に～い見つけた方が勝者になるとき
library(gtools)
n=8
k=2
perm=permutations(n,n)
Q=perm[1,]
np=nrow(perm)
p1st=numeric(np)
for(i in 1:np){
P=perm[i,]
tre=combn(n,k)
nt=ncol(tre)
re=numeric()
for(j in 1:nt){
re[j]=max(which(tre[1,j]==P),which(tre[2,j]==P))-
max(which(tre[1,j]==Q),which(tre[2,j]==Q))
}
p1st[i]=sum(re<0)
}
plot(p1st)
mv=p1st[which.max(p1st)]
(p.max=which(p1st==mv))
print(matrix(LETTERS[perm[p.max,]],ncol=8),quote=F)
850:卵の名無しさん
18/10/26 12:35:57.72 +rLvZj8m.net
数学板の、分からない問題はここに書いてね スレで話題のテーマで問題を作ってみた。
ド底辺シリツ医大生が８人いてそのうち３人は裏口入学であることが判明している。
８人をA,B,C,D,E,F,G,Hとする。
A,B,C,D,E,F,G,Hの順に調査して３人みつかれば調査は終了とする。
A,B,C,D,E,F,G,Hの順の調査と比べて早くみつかる可能性の探索順を列挙せよ。
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,] A B D E F G H C
[2,] A C D E F G H B
[3,] A D B E F G H C
[4,] A D C E F G H B
[5,] B A D E F G H C
[6,] B C D E F G H A
[7,] B D A E F G H C
[8,] B D C E F G H A
[9,] C A D E F G H B
[10,] C B D E F G H A
[11,] C D A E F G H B
[12,] C D B E F G H A
[13,] D A B E F G H C
[14,] D A C E F G H B
[15,] D B A E F G H C
[16,] D B C E F G H A
[17,] D C A E F G H B
[18,] D C B E F G H A
おい、ド底辺頭脳よ、これであっているか検証してみろ。
851:卵の名無しさん
18/10/26 12:56:44.55 +rLvZj8m.net
treasure <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
return(c(短軸有利=sum(re<0),長軸有利=sum(re>0),同等=sum(re==0)))
}
sapply(1:12,function(k) treasure(3,4,k))
sapply(1:20,function(k) treasure(4,5,k))
852:卵の名無しさん
18/10/26 13:32:10.72 R30QBCXi.net
>789 や他の算数の問題の答え まだぁ
宅建国試浪人の中卒事務員さん
853:卵の名無しさん
18/10/26 15:19:53.27 +rLvZj8m.net
>>797
算数の問題ですらないじゃん、
おまえは問題作成すらできないドアホなのか？
最近のトピックをつかったこういう問題面白いだろ？
答える頭はお前にはないだろうけど。
ド底辺シリツ医大生が８人いてそのうち３人は裏口入学であることが判明している。
８人をA,B,C,D,E,F,G,Hとする。
A,B,C,D,E,F,G,Hの順に調査して３人みつかれば調査は終了とする。
A,B,C,D,E,F,G,Hの順の調査と比べて早くみつかる可能性の探索順を列挙せよ。
ある医大で合格率の男女比が1.2で男子優位という結果だったという。
定員100で男子800人女子200人が受験して合格率の男女比が
1.2であったときに統計的には有意差があると言えるか？
854:卵の名無しさん
18/10/26 16:34:55.99 +rLvZj8m.net
国立大学をでていれば、これくらいの計算ができて当然。
できない椰子はド底辺シリツ医大卒と同じ頭脳レベルといえる。
東京医大、本来合格者入学許可へ 今年の受験生５０人
2018年10月25日 02時06分
東京医科大＝８月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLﾘﾝｸ(www.nishinippon.co.jp)
URLﾘﾝｸ(www.tokyo-med.ac.jp)によると
学年 第1学年 第2学年
在学者数 133 113
昨年入学者の留年者や退学者が０として、
大学が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
855:卵の名無しさん
18/10/26 16:36:52.25 +rLvZj8m.net
library(rjags)
y=c(50,19)
N=c(133,113)
Ntotal=length(y)
a=1
b=1
dataList=list(y=y,N=N,Ntotal=Ntotal,a=a,b=b)
# JAGS model
modelString ="
model{
for(i in 1:Ntotal){
y[i] ~ dbin(theta,N[i])
}
theta ~ dbeta(a,b)
}
"
writeLines(modelString,'TEMPmodel.txt')
jagsModel=jags.model('TEMPmodel.txt',data=dataList)
codaSamples=coda.samples(jagsModel,var=c('theta'),n.iter=20000,na.rm=TRUE)
summary(codaSamples)
js=as.vector(as.matrix(codaSamples))
BEST::plotPost(js,,xlab="裏口確率")
BEST::plotPost(js,showMode = TRUE,xlab="裏口確率")
856:卵の名無しさん
18/10/26 16:43:26.58 +rLvZj8m.net
>>799
昨年と今年で裏口入学率に違いがあったといえるか？
裏口入学率の差の期待値とその９５％信頼区間を求めよ
857:卵の名無しさん
18/10/26 16:44:57.64 +rLvZj8m.net
>>801
解法はここに置いた
ｽﾚﾘﾝｸ(math板:208番)
858:卵の名無しさん
18/10/26 16:56:31.17 +rLvZj8m.net
>>800
裏口入学率の事前分布は一様分布よりJefferey分布の方がいいかなと思ってやってみたが、殆ど数値に影響なかったな。
859:卵の名無しさん
18/10/26 20:01:23.00 R30QBCXi.net
>789 や他の算数の問題の答え まだぁ
宅建国試浪人の中卒事務員さん
860:卵の名無しさん
18/10/26 20:05:38.40 46gTk/Hn.net
>>804
問題の数式はどれだよ？
861:卵の名無しさん
18/10/26 20:07:01.94 46gTk/Hn.net
こういう風にちゃんと数字を出して問題文を自作しろよ。
東京医大、本来合格者入学許可へ 今年の受験生５０人
2018年10月25日 02時06分
東京医科大＝８月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLﾘﾝｸ(www.nishinippon.co.jp)
URLﾘﾝｸ(www.tokyo-med.ac.jp)によると
学年 第1学年 第2学年
在学者数 133 113
昨年入学者の留年者や退学者が０として、
大学が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
862:卵の名無しさん
18/10/26 20:56:35.85 R30QBCXi.net
おらぁ 算数の問題だ
URLﾘﾝｸ(www.nicovideo.jp)
863:卵の名無しさん
18/10/26 21:44:41.27 +rLvZj8m.net
>>791
外れ値扱いせず、対数正規分布に従うデータとして検定してみた。
負の値になるシミュレーションデータの組み合わせは除外して対数化して検定
n1=17 ; n2=12 ;m1=36.2 ;m2=56.1 ; sd1=21.5 ; sd2=30.1
fx1 = function(){
tmp=scale(rnorm(n1))*sd1+m1
while(min(tmp)<0){
tmp=scale(rnorm(n1))*sd1+m1
}
return(tmp)
}
fx2 = function(){
tmp=scale(rnorm(n2))*sd2+m2
while(min(tmp)<0) tmp=scale(rnorm(n2))*sd2+m2
return(tmp)
}
k=1e4
T.p_value=replicate(k,t.test(log(fx1()),log(fx2()),var=T)$p.value)
W.p_value=replicate(k,t.test(log(fx1()),log(fx2()))$p.value)
mean(T.p_value) ; quantile(T.p_value,c(0.025,0.975))
mean(W.p_value) ; quantile(W.p_value,c(0.025,0.975))
> mean(T.p_value) ; quantile(T.p_value,c(0.025,0.975))
[1] 0.1453533
2.5% 97.5%
0.04098323 0.50619976
> mean(W.p_value) ; quantile(W.p_value,c(0.025,0.975))
[1] 0.1407244
2.5% 97.5%
0.0313094 0.5523005
全く、有意差なしの結論になった。
864:卵の名無しさん
18/10/26 21:45:35.06 +rLvZj8m.net
>>807
こういう風にちゃんと数字を出して問題文を自作しろよ。
東京医大、本来合格者入学許可へ 今年の受験生５０人
2018年10月25日 02時06分
東京医科大＝８月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLﾘﾝｸ(www.nishinippon.co.jp)
URLﾘﾝｸ(www.tokyo-med.ac.jp)によると
学年 第1学年 第2学年
在学者数 133 113
昨年入学者の留年者や退学者が０として、
大学が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
865:卵の名無しさん
18/10/26 21:58:38.54 R30QBCXi.net
まだ答えられんのか？
やはり事務員には無理だったか
算数の問題だ おらぁ
URLﾘﾝｸ(www.nicovideo.jp)
866:卵の名無しさん
18/10/26 22:25:41.39 QZUdd6Brt
>>810
ちゃんと数字を出して問題文を自作しろよ。問題自作する能力すらないの？
こいういうトピックでの自分で問題を考えて解決するのは楽しいぞ。
東京医大、本来合格者入学許可へ 今年の受験生５０人
2018年10月25日 02時06分
東京医科大＝８月、東京都新宿区
東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
URLﾘﾝｸ(www.nishinippon.co.jp)
URLﾘﾝｸ(www.tokyo-med.ac.jp)によると
学年 第1学年 第2学年
在学者数 133 113
昨年入学者の留年者や退学者が０として、
大学が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
>>796
数学板で同じようにプログラミングして解決しようとする人も現れて議論が深まって楽しいいな。
867:卵の名無しさん
18/10/26 22:29:20.99 QZUdd6Brt
シミュレーションで自己検算していたけれど、他のひとも独立にプログラミングしてその結果が一致すると( ・∀・)ｲｲ!!
Haskellの暗号のようなコードは俺には可読性が悪いんだが、Rの方が読みにくいという人もいる。
C：高速、Haskell:不定長整数が扱える R:結果がすぐにグラフ化できるで使い分けてる。
868:卵の名無しさん
18/10/26 23:20:04.96 +rLvZj8m.net
>>810
答えるもなにも問題を書けよ、リンクは踏まないから。
３ｘ４が完成
> (t34=sapply(1:12,function(k) treasure0(3,4,k)))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
短軸有利 5 26 73 133 167 148 91 37 9 1 0 0
長軸有利 5 27 76 140 176 153 92 37 9 1 0 0
同等 2 13 71 222 449 623 609 421 202 64 12 1
４ｘ５が完成
> (t45=sapply(1:20,function(k) treasure0(4,5,k)))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
短軸有利 9 84 463 1776 5076 11249 19797 28057 32243 30095 22749 13820
長軸有利 9 83 453 1753 5075 11353 20057 28400 32528 30250 22803 13831
同等 2 23 224 1316 5353 16158 37666 69513 103189 124411 122408 98319
[,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
短軸有利 6656 2486 695 137 17 1 0 0
長軸有利 6657 2486 695 137 17 1 0 0
同等 64207 33788 14114 4571 1106 188 20 1
869:卵の名無しさん
18/10/26 23:29:02.10 +rLvZj8m.net
これをCに移植するのは
870:俺には重荷だな。 treasure0 <- function(m=3,n=4,k=2){ y=1:(m*n) (z=matrix(y,ncol=n,byrow=T)) (P=as.vector(z)) (Q=as.vector(t(z))) PQ <- function(x){ p=q=numeric(k) for(i in 1:k){ p[i]=which(P==x[i]) q[i]=which(Q==x[i]) } min(p)-min(q) } tre=combn(m*n,k) re=apply(tre,2,PQ) return(c(短軸有利=sum(re<0),長軸有利=sum(re>0),同等=sum(re==0))) }
871:卵の名無しさん
18/10/27 01:16:19.77 RG7nMybg.net
算数の問題だぁ
国試浪人には無理だろうが、さっさと解かんか！
URLﾘﾝｸ(www.nicovideo.jp)
872:卵の名無しさん
18/10/27 02:47:24.95 0SWbBA8x.net
>>815
数値を出して自作問題を書けよ。
統計スレに相応しい問題にしろよ。
これ答えてみ！
昨年入学者の留年者や退学者が０として、
東京医大が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
873:卵の名無しさん
18/10/27 10:58:30.69 ZN2oZPvM.net
六面体のサイコロでP君のサイコロは３面が１、Q君のサイコロは２面が１とする。
サイコロを降ってどちらか一方が１であればそちらが勝者。
どちらも１であるときは引き分け
どちらも１でないならば少なくともどちらかが１がでるまでサイコロをふる。
P君、Q君の勝つ確率を求めよ。
p=1/2
q=1/3
q: win
(1-p)*q + (1-p)^2*(1-q)*q+(1-p)^3*(1-q)^2*q+(1-p)^4*(1-q)^3*q+....
=(1-p)*q *( 1 + (1-p)*(1-q) + ((1-p)*(1-q))^2 + ((1-p)*(1-q))^3+...
let r=(1-p)*(1-q)=1/2 * 2/3 =1/3
=(1-p)*q *(1 + r + r^2 + r^3 + ...) = (1-1/2)*1/3 * 3/2 = 1/4 = 0.25
p:win
(1-1/3) * 1/2 * 3/2 = 0.5
draw 1-1/4-1/2= 0.25
874:卵の名無しさん
18/10/27 11:01:51.29 ZN2oZPvM.net
引き分けなしの場合：
六面体のサイコロでP君のサイコロは３面が１、Q君のサイコロは２面が１とする。
サイコロを降ってどちらか一方が１であればそちらが勝者。
どちらも１であるときはもしくはどちらも１でないならば どちらか一方だけが１が出て勝者が決まるまで繰り返す。
P君、Q君の勝つ確率を求めよ。
875:卵の名無しさん
18/10/27 12:39:05.73 ZN2oZPvM.net
>>818
no_draw rule
Σ[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
で解ける（ド底辺シリツ医大卒の裏口バカには無理）が、シミュレーションプログラムの練習をしてみる。
P=c(1,1,1,3,4,5)
Q=c(1,1,2,3,4,5)
sub = function(){
p=sample(P,1)
q=sample(Q,1)
if(p==1 & q!=1) return(1)
if(p!=1 & q==1) return(-1)
return(0)
}
throw <- function(){
r=sub()
while(r==0) r=sub()
return(r)
}
re=replicate(1e3,mean(replicate(1e3,throw()==-1)))
summary(re)
876:卵の名無しさん
18/10/27 12:42:12.98 ZN2oZPvM.net
> re=replicate(1e3,mean(replicate(1e3,throw()==-1)))
> summary(re)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.2840 0.3230 0.3330 0.3331 0.3430 0.3830
理論値と合致。
シミュレーションのデバッグの方が面倒だった。
877:卵の名無しさん
18/10/27 16:35:18.47 kZ9du42U.net
おら ど底辺事務員
算数の問題 の答えは まだか？
やっぱり 解けないんだなＷ
878:卵の名無しさん
18/10/27 16:39:48.42 ZN2oZPvM.net
答えるもなにも問題を書けよ、リンクは踏まないから。
最近のトピックを題材にした問題
これ答えてみ！
昨年入学者の留年者や退学者が０として、
東京医大が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
ド底辺シリツ医大卒なら無理だろうが、国立大学でていて解けないなら恥ずかしいね。
879:卵の名無しさん
18/10/27 18:43:40.34 sR9utt0J.net
はよ 算数の問題 の答えだせや
し! _ -─ ‐- ､ , -─-､ -‐─_ﾉ
小 国 /／￣＞ ´ ￣ ￣ ｀ヽ Y , ´ ） 事 え
学 試 L_ ／ ／ ヽ 務 |
生 浪 / ' ' i員 !? マ
ま 人 / / く ジ
で が l ,ｨ/! / /l/!,l ／厶,
だ 許 i ,.lrH‐|'| /‐!-Lﾊ_ l /-!'|/l ／｀'ﾒ､_iヽ
よ さ l | |_|_|_|/| / /__!__ |/!ﾄi i/-- ､ ﾚ!/ / ,-- ﾚ､⌒Y⌒ヽ
ね れ _ゝ|/'/⌒ヽ ヽﾄ､|/ '/￣｀ヾ 、ヽﾄ､Ｎ'/⌒ヾ ,ｲ￣｀ヾ,ﾉ!
l る 「 l ′ 「1 /てヽ′| | | 「L! ' i'ひ} ﾘ
の ヽ | ヽ__Ｕ, ､ヽ シノ ﾉ! ! |ヽ_､ｿ, ヾシ _ノ _ﾉ
-┐ は ,√ ! ￣ ﾘ l ! ￣ ￣ ７/
ﾚ'⌒ヽ/ ! | 〈 _人__人ﾉ_ i く //!
人_,､ノL_,iノ! /! ヽ r─‐- ､ 「 Ｌ_ヽ r─‐- ､ u ノ/
/ / lﾄ､ ＼ ヽ, -‐┤ ﾉ キ 了＼ ヽ, -‐┤ /／
ハ キ { / ヽ,ﾄ､ヽ/!｀ｈﾉ ） モ |/! 「ヽ, `ｰ ／） _ ‐'
ハ ャ ヽ/ r-､‐' // / |-‐ く | ＞ / / ｀'／／-‐､ ／
ハ ハ ＞ ／＼＼// / /ヽ_ ! イ （ / / ／／ / `ｧ-‐ '
ハ ハ / /! ヽ レ'/ ﾉ ＞ ' ∠ -‐ ￣ノヽ /
{ i l ! / フ / -‐ /￣/〉 〈 ＼ /!
880:卵の名無しさん
18/10/27 18:55:20.24 NpJ0k/uw.net
答えるもなにも問題を書けよ、リンクは踏まないから。
最近のトピックを題材にした問題
これ答えてみ！
昨年入学者の留年者や退学者が０として、
東京医大が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
ド底辺シリツ医大卒なら無理だろうが、国立大学でていて解けないなら恥ずかしいね。
881:卵の名無しさん
18/10/27 19:14:52.46 kZ9du42U.net
このスレは事務員が日がな一日妄想を垂れ流し
見物人たちがそれを見てフルボッコにするスレである
事務員 とは？
同一者の別名として、薄汚いジジイ、国試浪人、統計野郎、自称医科歯科、がある
自称医科歯科卒、実際は九州の駅弁国立出身で、卒業後は実家の東京に帰り国試浪人となり23年間経過、家庭教師バイトなどを経て現在は事務員、とされている
本人は、勤務医でたまに開業医の手伝いと内視鏡バイト、専門医なし、と主張しているが、連日連夜の異常な書き込み数の多さからは、勤務医とは考え難く、彼の職業がたとえ何であったとしても、人としてド底辺なことは間違いがない
自ら事務員であると告白したこともある
ｽﾚﾘﾝｸ(hosp板:108番)-109
自作自演も頻繁に行なっている
事務仕事の際に手にした他人の給与明細や同窓会費振込票を盗撮し、自分のものとしてアップロードしたこともある、犯罪なんじゃねーの？
統計絶対主義でマウントを取ろうとするが、大学卒業後の人生の大部分においてそれが何の役に立つのかはなはだ疑問である
統計に基づいた重病の治療成績を患者に説明しても、患者の希望する治療は患者の事情によって異なることは一定年数以上の臨床経験があるものなら誰でも知っていることだ
それが理解できない彼は、彼自身が素性を明らかにしているようなものである
コンプ薬屋同様に、私立医への執拗な誹謗中傷を数年余にわたって続けている
裏口批判をしているが、実際にそれが可能な資産家の子弟がいるであろう右派スレには近寄らない、一度書き込んだ形跡があったが指摘したら別人だと言って逃げた
弱いものが多そうなスレしか狙わないのもコンプ薬屋と共通する異常者の特徴である
劣等コンプレックスの塊なんだろうな
自分の主張に対して都合の悪い突込みが入ると、相手を私立医と決めつけて、さてはシリツだな、裏口バカ、というセリフ
882:でワンパターンに返すことが多い 国試問題を挙げて簡単だとほざいてますが、本番で通らなければお話になりません いくらネットでわめいていても、医師国試の本番で通らなければお話になりません 大事なことなので２回言いました
883:卵の名無しさん
18/10/27 21:04:28.55 kZ9du42U.net
算数の問題の答え まだぁ？
それはそうと高学歴スレに宣伝して来てやったぞ
ありがたく思えド底辺め
ｽﾚﾘﾝｸ(newsplus板)
884:卵の名無しさん
18/10/28 06:31:16.86 Q5x4cT59.net
>>826
>>803
残念ながら俺、東大卒じゃないのよ。
1期校は滑り止めに受けた理一通ったけど行かなかった。
885:卵の名無しさん
18/10/28 07:07:44.53 Q5x4cT59.net
>>826
答えるもなにも問題を書けよ、リンクは踏まないから。
最近のトピックを題材にした問題
これ答えてみ！
昨年入学者の留年者や退学者が０として、
東京医大が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。
裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。
ド底辺シリツ医大卒なら無理だろうが、国立大学でていて解けないなら恥ずかしいね。
886:卵の名無しさん
18/10/28 09:59:54.67 Q5x4cT59.net
>>814
k=2
n=m+1
の特殊ケースを移植してくれた人がいた。
パラメータを変えるだけで宝の数を増やせる仕様でないのは残念。
int main(){
int i,j,P[5000],Q[5000],c,m,n,mn,Pwin,Qwin,Draw,pp,qq;
for(m=1;m<70;m++){
n=m+1;mn=m*n;
for(i=0,c=1;i<m;i++)for(j=0;j<n;j++)P[j*m+i+1]=Q[i*n+j+1]=c++;
for(i=1,Pwin=Qwin=Draw=0;i<mn;i++)for(j=i+1;j<=mn;j++){
pp=P[i]>P[j]?P[j]:P[i];
qq=Q[i]>Q[j]?Q[j]:Q[i];
if(pp>qq)Pwin++;else if(pp<qq) Qwin++;else Draw++;
}
printf("%3d*%3d: %d, %d, %d\n",m,n,Pwin,Qwin,Draw);
}
return 0;
}
887:卵の名無しさん
18/10/28 12:08:07.77 6treyEuT.net
算数の問題 と書いてあるから 間違いなく算数の問題なのだ
というわけで 回答まだぁ
それはそうとド底辺事務員の好きな高学歴スレに
宣伝して来てやったぞ
ありがたく思え
ｽﾚﾘﾝｸ(newsplus板:600番)-
入試の時だけ成績が良くても、その後の生活で人はダメになるんやね
国立私立問わず医学部に入学しても、その後の生活がダメダメで
国試に通らなかったり、専門医も取れなかったりすると
人としてダメダメになるんやね
888:卵の名無しさん
18/10/28 12:29:02.77 Q5x4cT59.net
treasure0 <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
return(c(短軸有利=sum(re<0),長軸有利=sum(re>0),同等=sum(re==0)))
}
結局、
combn(m*n,k)
ここがボトルネックなんだよなぁ。
12部屋から6部屋選ぶ組み合わせは924通りしかないのに
20部屋から10部屋だと184756通り、
30部屋から15部屋だと155117520通り、
という感じなのでどうしても時間やメモリを食う。
889:卵の名無しさん
18/10/28 13:22:21.24 Q5x4cT59.net
やはり、cは高速だな
#include<stdio.h>
int main(){
int i,j,P[10000],Q[10000],c,m,n,mn,Pwin,Qwin,Draw,pp,qq;
for(m=1;m<100;m++){
n=m+1;
mn=m*n;
for(i=0,c=1;i<m;i++)for(j=0;j<n;j++)P[j*m+i+1]=Q[i*n+j+1]=c++;
for(i=1,Pwin=Qwin=Draw=0;i<mn;i++)for(j=i+1;j<=mn;j++){
pp=P[i]>P[j]?P[j]:P[i];
qq=Q[i]>Q[j]?Q[j]:Q[i];
if(pp>qq)Pwin++;else if(pp<qq) Qwin++;else Draw++;
}
printf("%3d*%3d: %d, %d, %d\n",m,n,Pwin,Qwin,Draw);
}
return 0;
}
90* 91: 16820566, 16703175, 10214
91* 92: 17575980, 17454585, 10441
92* 93: 18356554, 18231065, 10671
93* 94: 19162841, 19033167, 10903
94* 95: 19995399, 19861448, 11138
95* 96: 20854793, 20716472, 11375
96* 97: 21741593, 21598808, 11615
97* 98: 22656376, 22509032, 11857
98* 99: 23599724, 23447725, 12102
99*100: 24572226, 24415475, 12349
890:卵の名無しさん
18/10/28 19:48:06.14 Q5x4cT59.net
>>829
このコードをもとに宝を増やすように改変したつもりで数学板投稿したらをバグを指摘された。勉強になって( ・∀・)ｲｲ!!
このスレのド底辺からは学ぶことがないなあ
891:。
892:卵の名無しさん
18/10/28 20:16:51.08 AAZUotga.net
## 第22回生命表（完全生命表）URLﾘﾝｸ(www.mhlw.go.jp)
F=c(178,32,20,12,8,8,8,8,7,7,7,7,7,7,8,10,12,13,15,16,17,19,20,22,23,24,25,27,28,30,31,32,34,36,39,41,42,45,50,56,
62,68,73,79,85,94,104,114,124,134,145,159,174,189,202,215,226,237,250,268,291,318,346,372,399,433,471,511,554,603,662,729,802,874,954,1053,1180,1332,1505,1699,1909,2143,2409,2701,3004,3310,3622,3938,4253,4531,4757,4918,5025,
5024,4876,4598,4132,3594,3025,2464,1941,1477,1085,769,525,345,217,132,76,42,22,11,5,2,1,0)
M=c(202,34,24,16,11,10,10,10,9,8,7,7,8,11,13,17,21,26,32,39,45,49,51,53,55,55,54,54,54,56,57,59,61,65,69,73,75,78,
84,93,103,113,122,131,144,159,176,195,215,236,257,283,310,340,373,411,450,488,525,568,620,688,764,839,910,994,
1081,1166,1256,1349,1450,1561,1675,1776,1885,2021,2185,2377,2594,2819,3046,3279,3504,3714,3900,4043,4116,4127,
4080,3973,3810,3580,3302,2967,2567,2123,1718,1352,1034,768,554,387,262,171,108,66,39,22,12,6,3,2,1)
LE <-function(ndx,Y,N0=10^5){
n=length(ndx)
lx=numeric(n)
lx[1]=N0
for(i in 1:(n-1))
lx[i+1] <- lx[i] - ndx[i]
nqx=ndx/lx
nLx=numeric(n)
for(i in 1:n)
nLx[i] <- mean(c(lx[i],lx[i+1]))
nLx[n]=0
Tx=rev(cumsum(rev(nLx)))
le=Tx/lx
return(round(le[Y+1],2))
}
yome = function(x){
le=LE(m,x)
round(uniroot(function(x,u0) LE(f,x)-le,c(1,100))$root)
}
y=sapply(18:80,yome)
plot(18:80,y-18:80)
893:卵の名無しさん
18/10/29 21:31:07.37 VzAgCvD5.net
さて 算数の問題です
URLﾘﾝｸ(m.youtube.com)
この動画は part１ が 2009年に投稿され
最終回まであと１話となっています
最終回の登校日はいつになるでしょう？
あらゆる手段を用いて求めましょう
894:卵の名無しさん
18/10/29 21:44:21.84 72tLNU2F.net
>>831
Haskellの遅延IOを使うととりあえずのメモリ不足は克服できる。
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
Prelude> take 5 $ combinations 15 [1..30]
[[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15],[1,2,3,4,5,6,7,8,9,10,11,12,13,14,16],[1,2,3,4,5,6,7,8,9,10,11,12,13,14,17],[1,2,3,4,5,6,7,8,9,10,11,12,13,14,18],[1,2,3,4,5,6,7,8,9,10,11,12,13,14,19]]
895:卵の名無しさん
18/10/29 21:46:20.38 72tLNU2F.net
Prelude> tre = combinations 15 [1..30]
Prelude> tre !! 100000
[1,2,3,4,5,6,7,8,11,12,17,20,21,25,30]
10万個目の組み合わせも出してはく
896:れるが、処理が終わる時間が予想がつかないなぁ。
897:卵の名無しさん
18/10/29 22:12:36.39 72tLNU2F.net
>>831
min(p) - min(q)までHaskellに移植できた。
import Data.List
import Data.List.Split
k = 2
q = [1..12]
matQ = chunksOf 3 q
matP = transpose matQ
p = concat matP
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
tre = combinations 2 q
ip x1 x2 = minimum [elemIndices x1 p, elemIndices x2 p] !! 0
iq x1 x2 = minimum (elemIndices x1 q, elemIndices x2 q) !! 0
pqi x1 x2 = ip x1 x2 - iq x1 x2
898:卵の名無しさん
18/10/29 23:42:39.13 72tLNU2F.net
Haskellへの移植プログラム完成！
import Data.List
import Data.List.Split
m = 5 -- 縦マス(短軸）
n = 6 -- 横マス(長軸)
k = 5 -- 宝の数
q = [0..m*n-1]
matQ = chunksOf m q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換
combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
treasure = combinations k q -- 宝の組み合わせ
ip y = minimum $ map(\x -> elemIndices x p!!0) y -- 宝の、配列ｐでの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 -- 短軸方向探索ｐが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q
main = do
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw
899:卵の名無しさん
18/10/30 00:23:28.38 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 = minimum $ map(\x -> elemIndices x p!!0) y -- 宝の、配列ｐでの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 -- 短軸方向探索ｐが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw
900:卵の名無しさん
18/10/30 00:28:02.10 Eer9XEMu.net
暴走せずに200万個を数えてくれる
>treasure 5 6 8
p1st = 1827737, q1st = 1825076, draw = 2200112
901:卵の名無しさん
18/10/30 01:00:07.00 Eer9XEMu.net
\hs>treasure 5 6 10
p1st = 7995426, q1st = 8023257, draw = 14026332
902:卵の名無しさん
18/10/30 07:08:46.76 Eer9XEMu.net
>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184
903:卵の名無しさん
18/10/30 07:49:18.11 BVJsV4IZ.net
やはり事務員には難しすぎたようだね
さて 算数の問題です
URLﾘﾝｸ(m.youtube.com)
この動画は part１ が 2009年に投稿され
最終回まであと１話となっています
最終回の登校日はいつになるでしょう？
あらゆる手段を用いて求めましょう
904:卵の名無しさん
18/10/30 08:53:02.81 kXLKssbE.net
>>840
５×６マスで宝の数を１０まで増やしていくと、
D:\bin>for %i in (1,2,3,4,5,6,7,8,9,10) do treasure 5 6 %i
D:\bin>treasure 5 6 1
p1st = 14, q1st = 14, draw = 2
D:\bin>treasure 5 6 2
p1st = 203, q1st = 197, draw = 35
D:\bin>treasure 5 6 3
p1st = 1801, q1st = 1727, draw = 532
D:\bin>treasure 5 6 4
p1st = 11418, q1st = 11008, draw = 4979
D:\bin>treasure 5 6 5
p1st = 55469, q1st = 54036, draw = 33001
D:\bin>treasure 5 6 6
p1st = 215265, q1st = 211894, draw = 166616
D:\bin>treasure 5 6 7
p1st = 685784, q1st = 680768, draw = 669248
D:\bin>treasure 5 6 8
p1st = 1827737, q1st = 1825076, draw = 2200112
D:\bin>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184
D:\bin>treasure 5 6 10
p1st = 7995426, q1st = 8023257, draw = 14026332
１：同等
１～８：短軸探索有利
９、１０：長軸探索有利
という結果になった。
久しぶりにバッチコマンドを書いた
905:卵の名無しさん
18/10/30 08:53:48.77 kXLKssbE.net
>>844
リンクは踏まないから問題文をかけよ。
906:卵の名無しさん
18/10/30 09:01:44.04 Eer9XEMu.net
>>844
＞あらゆる手段を用いて求めましょう
答： お前が答を書け
解答終了！！！
>>840
数学板の賢者がｃでプログラムするのを期待していたが、ｃだと最初にメモリを割り当てなくちゃいけないから難しいんじゃないかな？
動的メモリ確保にしてメモリ不足を回避する必要がでてくると思う。
Rで３０個から１５個取�
907:闖oす組み合わせをしようとしたら > nrow(combn(30,15)) Error: cannot allocate vector of size 8.7 Gb のメッセージがでて実行不能だった。
908:卵の名無しさん
18/10/30 10:25:52.99 Eer9XEMu.net
３０個から１５個取り出す組み合わせは155117520通りあるので
シミュレーションで頻度を出して確率を推測してみることにした。
m=5
n=6
mn=m*n
P=as.vector(matrix(1:mn,ncol=n,byrow=T))
Q=1:mn
sim <- function(k=15,all=FALSE){
tre=sample(mn,k)
p=q=numeric(k)
for(i in 1:k){
p[i]=which(tre[i]==P)
q[i]=which(tre[i]==Q)
}
p_order=ifelse(all,max(p),min(p))
q_order=ifelse(all,max(q),min(q))
p_order - q_order
}
choose(30,15)
re15=replicate(1e5,sim(15))
mean(re15<0) < mean(re15>0)
c(p1st_sim=mean(re15<0), q1st_sim=mean(re15>0) ,draw_sim=mean(re15==0))
> 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
909:卵の名無しさん
18/10/30 12:49:00.84 BVJsV4IZ.net
おらっ 新しい算数の問題だ
自分で回答が出せなければ０点だぞW
URLﾘﾝｸ(www.nicovideo.jp)
910:卵の名無しさん
18/10/30 14:03:59.31 Eer9XEMu.net
treasures <- function(m,n,k){
mn=m*n
Q=1:mn
(mat=matrix(Q,ncol=n,byrow=T))
(P=as.vector(mat))
tre=combn(mn,k)
PQ <- function(x){
p=q=p_q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
p_q[i]=p[i]-q[i]
}
return(p_q)
}
re=apply(tre,2,PQ)
result=matrix(rep(NA,k*3),ncol=3)
colnames(result)=c('p1st','q1st','even')
for(i in 1:k){
result[i,] = c(sum(re[i,]<0),sum(re[i,]>0),sum(re[i,]==0))
}
return(result)
}
treasures(4,5,5)
# URLﾘﾝｸ(tpcg.io)
911:卵の名無しさん
18/10/30 14:06:10.19 Eer9XEMu.net
D:\bin>treasure 5 6 11
p1st = 13346984, q1st = 13395944, draw = 27884372
D:\bin>treasure 5 6 12
p1st = 19312228, q1st = 19372871, draw = 47808126
912:卵の名無しさん
18/10/30 14:08:53.77 Eer9XEMu.net
>>850
この結果は面白いな
全体としては互角なんだが、勝負のルールによって勝者がかわるのがこれで分かった。
> treasures(4,5,5)
p1st q1st even
[1,] 1948 9680 3876
[2,] 5488 10016 0
[3,] 7752 7752 0
[4,] 10016 5488 0
[5,] 9680 1948 3876
913:卵の名無しさん
18/10/30 16:01:45.84 Eer9XEMu.net
>>851
メモリ不足でクラッシュするかと思ったが案外、頑強だな。
エラー処理皆無でコードを書いたのだが。
１億を超える演算をこなしている。
Haskellのコンパイラーは優秀で感心。
D:\bin>treasure 5 6 13
p1st = 24301031, q1st = 24358063, draw = 71100756
914:卵の名無しさん
18/10/30 16:05:37.28 BVJsV4IZ.net
おい、婚活に失敗したド底辺
算数の問題の答えはまだぁ
まぁオマエの頭じゃ無理だろうな
915:卵の名無しさん
18/10/30 16:06:30.58 BVJsV4IZ.net
このスレは事務員が日がな一日妄想を垂れ流し
見物人たちがそれを見てフルボッコにするスレである
事務員 とは？
同一者の別名として、薄汚いジジイ、国試浪人、統計野郎、自称医科歯科、がある
自称医科歯科卒、実際は九州の駅弁国立出身で、卒業後は実家の東京に帰り国試浪人となり23年間経過、家庭教師バイトなどを経て現在は事務員、とされている
本人は、勤務医でたまに開業医の手伝いと内視鏡バイト、専門医なし、と主張しているが、連日連夜の異常な書き込み数の多さからは、勤務医とは考え難く、彼の職業がたとえ何であったとしても、人としてド底辺なことは間違いがない
自ら事務員であると告白したこともある
ｽﾚﾘﾝｸ(hosp板:108番)-109
自作自演も頻繁に行なっている
事務仕事の際に手にした他人の給与明細や同窓会費振込票を盗撮し、自分のものとしてアップロードしたこともある、犯罪なんじゃねーの？
統計絶対主義でマウントを取ろうとするが、大学卒業後の人生の大部分においてそれが何の役に立つのかはなはだ疑問である
統計に基づいた重病の治療成績を患者に説明しても、患者の希望する治療は患者の事情によって異なることは一定年数以上の臨床経験があるものなら誰でも知っていることだ
それが理解できない彼は、彼自身が素性を明らかにしているようなものである
コンプ薬屋同様に、私立医への執拗な誹謗中傷を数年余にわたって続けている
裏口批判をしているが、実際にそれが可能な資産家の子弟がいるであろう右派スレには近寄らない、一度書き込んだ形跡があったが指摘したら別人だと言って逃げた
弱いものが多そうなスレしか狙わないのもコンプ薬屋と共通する異常者の特徴である
劣等コンプレックスの塊なんだろうな
自分の主張に対して都合の悪い突込みが入ると、相手を私立医と決めつけて、さてはシリツだな、裏口バカ、というセリフでワンパターンに返すことが多い
国試問題を挙げて簡単だとほざいてますが、本番で通らなければお話になりません
いくらネットでわめいていても、医師国試の本番で通らなければお話になりません
大事なことなので２回言いました
916:卵の名無しさん
18/10/30 16:12:20.65 Eer9XEMu.net
>>854
問題文がないのだから、無理だよ。
＞あらゆる手段を用いて求めましょう
答： お前が答を書け
解答終了！！！
917:卵の名無しさん
18/10/30 18:54:55.33 Eer9XEMu.net
treasure_order <- function(m,n,k){
mn=m*n
Q=1:mn
(mat=matrix(Q,ncol=n,byrow=T))
(P=as.vector(mat))
tre=combn(mn,k)
PQk <- function(x){
p=q=p_q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
pq=sort(unique(c(p,q)))
np=as.numeric(pq[1:k] %in% p)
nq=as.numeric(pq[1:k] %in% q)
10*np + nq
}
(re=apply(tre,2,PQk))
(result=matrix(NA,nrow=k,ncol=3))
for(i in 1:k){
result[i,]=c(sum(re[i,]==10),sum(re[i,]==1 ),sum(re[i,]==11))
}
colnames(result)=c('P_gets','Q_gets','both_get')
result
}
treasure_order(4,5,5)
918:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています