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:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています