20/01/20 15:01:37.71 4hsZHhAP.net
>>132
論証も大変だが、プログラムも（）の対応で混乱しそう。
n=5
TE=gtools::permutations(2,n,v=0:1,rep=T) # TE : 第１行が00000 で始まり最終行が11111で終わる行列
colnames(TE)=LETTERS[1:n] # 各列の名前A~E
fn <- function(x){
if(sum(x)==n|sum(x)==0) return(FALSE)
cond <- function(P,Q) !(P & !Q) # P ⇒ Qの真偽を返す関数
all(c( # all : c()内の , で区切られた命題が全て正しいかTRUE/FALSEで返す
# c(Aが正直者で証言が正しい|Aが嘘つきで証言が嘘, Bが正直者で証言が正しい|Bが嘘つきで証言が嘘, ...)
(x[1]==1 & x[2]==1) | (x[1]==0 & x[2]==0),
(x[2]==1 & x[1]==1 ) | (x[2]==0 & x[1]==0),
(x[3]==1 & cond(x[2]==0,x[3]==0)) | (x[3]==0 & !cond(x[2]==0,x[3]==0)),
(x[4]==1 & cond(x[3]==1,x[4]==1)) | (x[4]==0 & !cond(x[3]==1,x[4]==1)),
(x[5]==1 & (cond(x[4]==0,x[5]==0)|cond(x[4]==1,x[5]==1)) )
| (x[5]==0 & !(cond(x[4]==0,x[5]==0) | cond(x[4]==1,x[5]==1)) )
))
}
TE[apply(TE,1,fn),] # 各行にfooを適用して返り値がTRUEのものを表示
151:卵の名無しさん
20/01/20 19:09:21 4hsZHhAP.net
冗長な記述を同値の命題に変更
n=5
TE=gtools::permutations(2,n,v=0:1,rep=T) # TE : 第１行が00000 で始まり最終行が11111で終わる行列
colnames(TE)=LETTERS[1:n] # 各列の名前A~E
fn <- function(x){
if(sum(x)==n|sum(x)==0) return(FALSE)
cond <- function(P,Q) !(P & !Q) # P ⇒ Qの真偽を返す関数
all(c( # all : c()内の , で区切られた命題が全て正しいかTRUE/FALSEで返す
# c(Aが正直者で証言が正しい|Aが嘘つきで証言が嘘, Bが正直者で証言が正しい|Bが嘘つきで証言が嘘, ...)
(x[1]==1 & x[2]==1) | (x[1]==0 & x[2]==0),
(x[2]==1 & x[1]==1 ) | (x[2]==0 & x[1]==0),
(x[3]==1 & cond(x[2]==0,x[3]==0)) | (x[3]==0 & !cond(x[2]==0,x[3]==0)),
(x[4]==1 & cond(x[3]==1,x[4]==1)) | (x[4]==0 & !cond(x[3]==1,x[4]==1)),
(x[5]==1 & (x[4]==x[5]) ) | (x[5]==0 & (x[4]!=x[5]) )
))
}
TE[apply(TE,1,fn),] # 各行にfnを適用して返り値がTRUEのものを表示
152:卵の名無しさん
20/01/20 20:53:54 4hsZHhAP.net
arg <- function(z) ifelse(Arg(z)<0,Arg(z)+2*pi,Arg(z))
arg(-1i)/pi
arg(1i)/pi
153:卵の名無しさん
20/01/20 21:06:36 4hsZHhAP.net
arg <- function(z) ifelse(Arg(z)<0,Arg(z)+2*pi,Arg(z))
arg(-1i)/pi
arg(1i)/pi
iArg <- function(B,A) arg(B)-arg(A) < pi # increase Arg?
iArg(-1+1i,2+1i)
iArg(-1-1i,2+1i)
154:卵の名無しさん
20/01/21 06:29:50 qPsKoWg0.net
signを使うと記述が簡略化できる。
iArg <- function(B,A){ # increase Arg?
d=arg(B)-arg(A)
ifelse(abs(d)<pi, sign(d),-sign(d))
}
iArg(1i,1)
iArg(1,1i)
iArg(1,-1i)
iArg(-1i,1)
155:卵の名無しさん
20/01/21 15:44:20 qPsKoWg0.net
"1～5の自然数が書かれた5枚のカードを、A～Eの生徒5人に先生が1枚ずつ配った。
5人はそれぞれ自分のカードの数は分かるが、他の人のカードの数はわからない。
また、先生は誰の数もわからない。
さて、先生とA～Eとの間で次のような会話があった。
なお、全員正直者であり、後から答える人は先の会話を聞いて参考にしている。
先生「Aさん、誰が1番大きい数ですか？」
A「わかりません」
先生「Bさん、あなたはCさんよりも大きい数ですか？」
B「わかりません」
先生「Cさん、あなたはDさんよりも大きい数ですか？」
C「わかりません」
先生「Dさん、あなたはBさんよりも大きい数ですか？」
D「○○○○○」
先生「Bさん、あなたはCさんよりも大きい数ですか？」
B「いいえ」
先生「たった今、皆さんの数がわかりました」
問1、○○○○○に入る言葉は「はい」「いいえ」「わかりません」のどれか？
問2、A～Eの数は何か？"
156:卵の名無しさん
20/01/21 15:45:28 qPsKoWg0.net
>>146
プログラムによる解法、場合分けがいるな。
pm=gtools::permutations(5,5)
Yes <- function(x){
all(c(
x[1]!=5,
x[2] %in% 2:4,
x[3] %in% 2:4,
x[4]==4,
x[2]<x[3]
))
}
pm[apply(pm,1,Yes)
157:,] No <- function(x){ all(c( x[1]!=5, x[2] %in% 2:4, x[3] %in% 2:4, x[4]==2, x[2]<x[3] )) } pm[apply(pm,1,No),]
158:卵の名無しさん
20/01/21 15:45:35 qPsKoWg0.net
DK <- function(x){
all(c(
x[1]!=5,
x[2] %in% 2:4,
x[3] %in% 2:4,
x[4]==3,
x[2]<x[3]
))
}
pm[apply(pm,1,DK),]
159:卵の名無しさん
20/01/21 20:38:57 qPsKoWg0.net
オリンピックサイズ・プール50m×25mの水の入ったプールの一つの角に監視員を置く。
この監視員は世界記録で移動するものとする。
水泳100m自由形 46秒91で水中を移動
陸上100m9秒58でプールサイドを」移動
この監視員がプールのどこへでも到達しうるには，最短で何秒必要か計算せよ。
160:卵の名無しさん
20/01/21 20:45:16 qPsKoWg0.net
>>149
答：10.77933776秒
161:卵の名無しさん
20/01/22 10:33:43.25 Y0XI94SC.net
hnr <- function(FUN,lo,up,n=100,maximum==FALSE,prec=1e-5){
xx=seq(lo,up,length=n)
by=(up-lo)/(n-1)
y=sapply(xx,FUN)
z=ifelse(maximum,max(y),min(y))
while(by>prec){
by=by/10
xx=seq(z-by,z+by,by)
y=sapply(xx,FUN)
z=ifelse(maximum,max(y),min(y))
}
}
162:卵の名無しさん
20/01/22 10:37:58.78 Y0XI94SC.net
数学板は雑な投稿すると速攻で突っ込みが入る。
ｽﾚﾘﾝｸ(math板:809番)
具体例を上げて詳述したら納得された。
>151は説明の手順を関数化してみた。
hnr は人間ニュートンラスソンから命名。
ド底辺頭脳医ってプログラム組めないの？
今どき小学生でもできるけど。
163:卵の名無しさん
20/01/22 10:43:14.97 Y0XI94SC.net
# under construction
hnr <- function(FUN,lo,up,n=100,maximum=FALSE,prec=1e-5){
xx=seq(lo,up,length=n)
by=(up-lo)/(n-1)
y=sapply(xx,FUN)
z=ifelse(maximum,max(y),min(y))
while(by>prec){
by=by/10
xx=seq(z-by,z+by,by)
y=sapply(xx,FUN)
z=ifelse(maximum,max(y),min(y))
}
return(z)
}
164:卵の名無しさん
20/01/22 14:18:42.69 L+skXSXt.net
>>153
工事終了
rm(list=ls())
graphics.off()
# 関数の極大極小値を返すlocal maximum and minimum
# FUN:関数 lo:範囲下限 up:範囲上限 n:範囲分割数 prec:精度
lmm <- function(FUN,lo,up,n=100,maximum=FALSE,prec=1e-7,print=TRUE){
xx=seq(lo,up,length=n) # 範囲をn個の数値に分割
by=(up-lo)/(n-1) # by:分割幅を計算
y=sapply(xx,FUN) # xxで実行して
if(print) plot(xx,y,bty='l',type='l') # 初期値範囲の関数を描画
z=ifelse(maximum,max(y),min(y)) # 極大値/極小値 zを調べて
x=xx[which(y==z)]# xに極値を与えた変数の値を返す
while(by>prec){ # 分割幅が精度より大きい限り
xx=seq(x-by,x+by,length=n) # xを中心に[x-by,x+by]で実行
y=sapply(xx,FUN)
z=ifelse(maximum,max(y),min(y))
x=xx[min(which(y==z))] # 仮極値を与えた変数値の小さい方
by=by/10 # 次回は分割幅byを1/10にして実行
}
re=c(value=z,at=xx[min(which(y==z))]) # 変数がatのときに極値valueをreにいれて
return(round(re,-log10(prec))) # 精度に応じた桁で返す
}
lmm(function(x) sin(x)+cos(x),-pi,pi,max=F)
lmm(function(x) x^3-6*x^2+10,-2,3,max=TRUE)
165:卵の名無しさん
20/01/22 14:24:54.99 L+skXSXt.net
検証
cosθ+sinθの最大値と最小値を求める
> lmm(function(x) sin(x)+cos(x),-pi,pi,max=T)
value at
1.4142136 0.7853981
> lmm(function(x) sin(x)+cos(x),-pi,pi,max=F)
value at
-1.4142136 -2.3561945
166:卵の名無しさん
20/01/22 14:47:44.55 L+skXSXt.net
実用に耐える精度が出た！
> options(digits=22)
> f=function(x) cos(x)+sin(x)
> lmm(f,-pi,pi,maximum=TRUE,prec=1e-16)
value at
1.41421356237309515 0.78539815960067616
> opt=optimise(f,c(-pi,pi),maximum=TRUE,tol=1e-7) ; c(opt[[2]],opt[[1]])
[1] 1.41421356237309492 0.78539816341094015
> c(f(pi/4), pi/4)
[1] 1.41421356237309492 0.78539816339744828
167:卵の名無しさん
20/01/22 14:54:28 L+skXSXt.net
Wolfram
1.414213562373095048801688724209698078569671875376948073176...
0.785398163397448309615660845819875721049292349843776455243...
168:卵の名無しさん
20/01/22 16:54:12 +i3wg50W.net
いい加減にしろ
169:卵の名無しさん
20/01/22 18:31:37 L+skXSXt.net
>>158
ド底辺頭脳医ってプログラム組めないの？
今どき小学生でもできるけど。
170:卵の名無しさん
20/01/22 18:38:17 L+skXSXt.net
こういうのって幼稚園児でも問題の意味はわかると思うよね。
AからEの5人の中に正直者と嘘つきの双方がいるとわかっている。
誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきを全員列挙せよ。
A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」
171:卵の名無しさん
20/01/22 23:02:13.43 dfSEVG8E.net
>>1認知症の原因、実は「処方薬」だった 認知症の疑いがあっても絶対に医者にはかかるな！ [268718286]
ｽﾚﾘﾝｸ(poverty板)
60以上のデパスやハルシオンでなるぞ
発ガン性もあり
と言うか親の金で塾やら私立行ってた奴が頭悪い
172:卵の名無しさん
20/01/23 07:39:15.87 8z/WT4oz.net
今日の水道の温度は7℃
65℃のお茶を98℃の湯を使って500ｍL作りたい。
yu = function(V=500,T=65,C=7,H=98) V*(T-C)/(H-C)
yu()
173:卵の名無しさん
20/01/23 12:20:31 +5jjGI8j.net
発がん性はかなりの薬であるだろ
デパスハルシオンは知らないけど
174:卵の名無しさん
20/01/23 15:06:38 7m5XnL0B.net
>>161
さてはシリツだな。
稀ではあるが、治療効果が期待できる認知症ってあるよね。
正常圧水頭症とか慢性硬膜下血腫とか。
175:卵の名無しさん
20/01/23 15:07:30 7m5XnL0B.net
いまやシリツ頭脳でも国立大学に入学できる時代なんだなぁ。
176:卵の名無しさん
20/01/23 15:59:39 +5jjGI8j.net
日医、自警なら新設は入れるだろ
177:卵の名無しさん
20/01/23 17:57:26 UWnMt3Wz.net
川崎医大
178:卵の名無しさん
20/01/23 18:18:57 tWC6ssYu.net
雨竜郡雨竜町出身介護福祉士やってます。
179:卵の名無しさん
20/01/23 18:35:17 tWC6ssYu.net
あの板おもしろいっすw
180:卵の名無しさん
20/01/23 18:50:45.37 tWC6ssYu.net
今は関東に在住してます。医者のなりたかったです
181:卵の名無しさん
20/01/23 20:18:01.31 8z/WT4oz.net
東京医大の入試不正、前文科省局長の息子だったのに、
医科大学が文科省の管轄であるという常識すら欠いているようなアホが裏口入学するのがド底辺シリツ医大。
182:卵の名無しさん
20/01/23 20:45:33.88 tWC6ssYu.net
雨竜は何にもなくダサかった やっぱ東京がいい
183:卵の名無しさん
20/01/23 21:46:32 tWC6ssYu.net
それにしてもあの人達はよく調べたよな。あーゆう人達が本当の頭の良い人なんだな
184:卵の名無しさん
20/01/23 22:05:24 tWC6ssYu.net
やはり医者を怒らせたら怖いよな 力持ってるし
185:卵の名無しさん
20/01/23 22:46:44 8z/WT4oz.net
決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
186:卵の名無しさん
20/01/23 22:48:27 8z/WT4oz.net
>>174
ド底辺シリツ医は特にそうみたいだね。
底辺私立医大卒が国立大学医学部卒を妬むスレ [転載禁止](c)2ch.net
ｽﾚﾘﾝｸ(hosp板:695番)
695 名前：卵の名無しさん[] 投稿日：2016/09/25(日) 22:04:25.53 ID:FzuOCpqV
匿名先生
精神科
▲
1
評価
私立出でも、素晴らしい医師がいるし、国立でも、どうしようもないのもいる。
そんなのは、あく
187:まで個人の資質差と思ってました。 が、 今いる病院（医師の大部分が、いわゆる底辺私立）に来て、その考えは変わりました。 とある、気に食わない医師の白衣やサンダルを隠す、机の上にゴミを置く、など、小学生じみた嫌がらせをしている医師というのを、はじめて見ました。 国立出の変人も多々みてきましたが、こんな、小学生のイタズラみたいなことをする人は、いませんでした。 やはり、底辺私立は、それなりなんだと思いました。 2016年09月02日 00時21分09秒
188:卵の名無しさん
20/01/23 22:52:15 tWC6ssYu.net
あんたは何も知らないんだねw
189:卵の名無しさん
20/01/23 23:15:03 8z/WT4oz.net
>>177
フェラしてもらうには何と言えばいい？
190:卵の名無しさん
20/01/23 23:20:36.13 tWC6ssYu.net
一般常識を知らない人が病院組織で上手く立ち回ることはできないぞよw人間関係のストレスがないバイトでもしなさいwまずは医者になってからだなwあ、医者になるのは無理かw
191:卵の名無しさん
20/01/23 23:21:44.58 8z/WT4oz.net
>>179
フェラしてもらうには何と言えばいい？
女子大生にフェラしてもらったことないの？
192:卵の名無しさん
20/01/23 23:43:32 tWC6ssYu.net
馴れ馴れしい奴だな
193:卵の名無しさん
20/01/24 00:37:03 0ZD3nIoE.net
聖マリ不正入試
URLﾘﾝｸ(news.yahoo.co.jp)
以前は偏差値が５０もなかったから入試の必要もなかった
今もほとんどが金とコネだけだな
194:卵の名無しさん
20/01/24 00:37:26 vNJL53lT.net
格安風俗性病罹患w
195:卵の名無しさん
20/01/24 00:52:42.36 vNJL53lT.net
雨竜は性病多かった
196:卵の名無しさん
20/01/24 04:50:22 p9O7hNNU.net
フェラには興味がないのか？
いや、単に頭が悪いだけなんだろうな。
このスレに合わせてみた。
入試前に寄付をしたド底辺シリツ受験生の親が
寄付金だけとられて不合格になるのでないかと心配して理事長と談判したところ、合格と寄付金は不可分であり、
お子様が面接で正しい主張をすれば
寄付金返却か
裏口合格になる
と言われた。
理事長は嘘をつかない。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？
197:卵の名無しさん
20/01/24 07:18:06 p9O7hNNU.net
これに答を出さずに逃げ回るってド底辺医はフェラ嫌いか、頭が悪いんだろうなぁ。
決して嘘をつかない女子大生から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
198:卵の名無しさん
20/01/24 07:53:19 p9O7hNNU.net
この方が実情にあってるかもな。
合格と寄付金は不可分を掲げるド底辺医大に入試前に寄付をした受験生の親が寄付金だけとられて不合格になるのでないかと心配して理事長と談判した。
合格と寄付金は不可分という大学の方針は常に正しいが、理事長は常に嘘をつく。
談判では、
「お子様が面接で正しい主張をすれば寄付金返却か裏口合格になる、間違った主張であれば寄付金没収で不合格になる」
と言われた。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？
199:卵の名無しさん
20/01/24 07:56:41 p9O7hNNU.net
￢（P⇒Q ∧ ￢P⇒￢Q）
￢（P⇒Q）∨ ￢（￢P⇒￢Q）where P⇒Q ≡ ￢（P∧￢Q）
￢￢（P∧￢Q）∨ ￢￢（￢P∧￢￢Q）
（P∧￢Q）∨ （￢P∧Q
200:卵の名無しさん
20/01/24 07:58:07.64 p9O7hNNU.net
￢（P⇒Q ∧ ￢P⇒￢Q）
￢（P⇒Q）∨ ￢（￢P⇒￢Q）where P⇒Q ≡ ￢（P∧￢Q）
￢￢（P∧￢Q）∨ ￢￢（￢P∧￢￢Q）
（P∧￢Q）∨ （￢P∧Q）
201:卵の名無しさん
20/01/24 09:01:22.19 p9O7hNNU.net
>>186
こういう風にすると少し難しくなるな。国立卒でもド底辺頭脳には無理と思う。
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
202:卵の名無しさん
20/01/24 10:14:06 9ccIONJj.net
ニュー速＋にも正しい意見が
49 名無しさん＠１周年 2020/01/23(木) 06:28:26.69 ID:fsvBnYrD0
私立医大はどうでもいい。
金で医者になるところだろ。
203:卵の名無しさん
20/01/24 13:44:27.94 p9O7hNNU.net
# URLﾘﾝｸ(detail.chiebukuro.yahoo.co.jp)
204:" ６人集まると、その中に「互いに知り合いの３人が存在する」か、または、「互いに知り合いでない３人が存在する」ことを示せ。 " rm(list=ls()) c2=gtools::combinations(6,2) # 6人の2人選ぶ組み合わせ 15通り p15=gtools::permutations(2,15,0:1,re=T) # 15通りが知り合い(1)か否（０）か 2^15通り c3=gtools::combinations(6,3) ; # 6人から3人を選ぶ組み合わせ 20通り ramsey <- function(x){ (acq=c2[x==1,]) # 知り合いの組み合わせ if(is.null(nrow(acq))) acq=t(as.matrix(acq)) is3acq=FALSE nr=nrow(c3) for(i in 1:nr){ j1=which(c3[i,1] == acq[,1]) # c3の組み合わせの1人目の1列目のindexを返す j2=which(c3[i,2] == acq[,1]) j3=which(c3[i,3] == acq[,1]) aq1=acq[c(j1,j2,j3),2] # 1列目の知り合い j1=which(c3[i,1] == acq[,2]) # c3の組み合わせの1人目の2列目のindexを返す j2=which(c3[i,2] == acq[,2]) j3=which(c3[i,3] == acq[,2]) aq2=acq[c(j1,j2,j3),1] # 2列目の知り合い su=sort(unique(c(aq1,aq2))) # 重複を除いてソート is3acq=all(c3[i,] %in% su) # 互いに知り合いかを返す if(is3acq) break } if(is3acq) return(1)
205:卵の名無しさん
20/01/24 13:44:54.01 p9O7hNNU.net
(acq=c2[x==0,]) # 他人の組み合わせ
if(is.null(nrow(acq))) acq=t(as.matrix(acq))
is3acq=FALSE
nr=nrow(c3)
for(i in 1:nr){
j1=which(c3[i,1] == acq[,1]) # c3の組み合わせの1人目の1列目のindexを返す
j2=which(c3[i,2] == acq[,1]) #
j3=which(c3[i,3] == acq[,1]) #
aq1=acq[c(j1,j2,j3),2] # 1列目の他人
j1=which(c3[i,1] == acq[,2]) # c3の組み合わせの1人目の2列目のindexを返す
j2=which(c3[i,2] == acq[,2]) #
j3=which(c3[i,3] == acq[,2]) #
aq2=acq[c(j1,j2,j3),1] # 2列目の他人
su=sort(unique(c(aq1,aq2))) # 重複を除いてソート
is3acq=all(c3[i,] %in% su) # 互いに他人かを返す
if(is3acq) break
}
if(is3acq) return(-1)
return(0)
}
sum(apply(p15,1,ramsey)==0)
206:卵の名無しさん
20/01/24 15:15:31.42 p9O7hNNU.net
>>187
数学板に題材を変えて投稿したら速攻で正解がきたなぁ。
207:卵の名無しさん
20/01/24 15:28:17 TuIVzM9K.net
問題は過去を克服することではありません。
さようなことができるわけはありません。
後になって過去を変えたり、
起こらなかったことにするわけにはまいりません。
しかし過去に目を閉ざす者は結局のところ現在にも盲目となります。(ヴァイツゼッカー)
URLﾘﾝｸ(imgur.com)
208:卵の名無しさん
20/01/24 15:30:15 TuIVzM9K.net
患者様からのありがたいお言葉です。
>>
膨大な情報を活用しなければいけない仕事だと認識しているからこそ、命のかかった病気においては特に、せめて自分よりは頭がいい先生に診てもらいたいと思ってしまいます。その場合に、出身校、入学した学校、は嘘をつかないとおもいます。
<<
>出身校、入学した学校、は嘘をつかないとおもいます。
>出身校、入学した学校、は嘘をつかないとおもいます。
>出身校、入学した学校、は嘘をつかないとおもいます。
大切な教えなので 繰り返しました。
209:卵の名無しさん
20/01/24 17:32:38 TuIVzM9K.net
##
"
N(=6)人集まると、その中に「互いに知り合いのM(=3)人が存在する」か、または、
「互いに知り合いでないM人が存在する」ことを示せ。
"
rm(list=ls())
N=6;M=3
(cN2=gtools::combinations(N,2)) # N人の2人選ぶ組み合わせ choose(N,2)(=15)通り
N2=nrow(cN2)
pN=gtools::permutations(2,N2,0:1,re=T) # N2通りが知り合い(1)か否（０）か 2^N2通り
(cNM=gtools::combinations(N,M)) ; # N人からM人を選ぶ組み合わせ choose(N,M)(=20)通り
(x=pN[10000,])
(x=pN[1,])
210:卵の名無しさん
20/01/24 17:32:54 TuIVzM9K.net
ramsey <- function(x){
(acq=cN2[x==1,]) # 知り合いの組み合わせ
if(is.null(nrow(acq))) acq=t(as.matrix(acq))
is3acq=FALSE
J=NULL
nr=nrow(cNM)
for(i in 1:nr){
for(j in 1:M){
J=c(J,which(cNM[i,j]==acq[,1])) # cN2の組み合わせのj人目の1列目のindexを返す
}
aq1=acq[J,2] # 1列目の知り合い
for(j in 1:M){
J=c(J,which(cNM[i,j]==acq[,2])) # cN2の組み合わせのj人目の2列目のindexを返す
}
aq2=acq[J,1] # 2列目の知り合い
su=sort(unique(c(aq1,aq2))) # 重複を除いてソート
is3acq=all(cNM[i,] %in% su) # 互いに知り合いかを返す
if(is3acq) break
}
if(is3acq) return(1)
211:卵の名無しさん
20/01/24 17:33:02 TuIVzM9K.net
(acq=cN2[x==0,]) # 他人の組み合わせ
if(is.null(nrow(acq))) acq=t(as.matrix(acq))
is3acq=FALSE
J=NULL
nr=nrow(cNM)
for(i in 1:nr){
for(j in 1:M){
J=c(J,which(cNM[i,j]==acq[,1])) # cN2の組み合わせのj人目の1列目のindexを返す
}
aq1=acq[J,2] # 1列目の知り合い
for(j in 1:M){
J=c(J,which(cNM[i,j]==acq[,2])) # cN2の組み合わせのj人目の2列目のindexを返す
}
aq2=acq[J,1] # 2列目の知り合い
su=sort(unique(c(aq1,aq2))) # 重複を除いてソート
is3acq=all(cNM[i,] %in% su) # 互いに知り合いかを返す
if(is3acq) break
}
if(is3acq) return(-1)
return(0)
}
sum(apply(pN,1,ramsey)==0)
212:卵の名無しさん
20/01/24 17:34:39 TuIVzM9K.net
一般化して書き直したので４人でも成立することが確認できた。
213:卵の名無しさん
20/01/24 20:28:09.43 TuIVzM9K.net
"A～Dの4人が受けた試験の結果について、4人はそれぞれ以下のように述べているが、
本当のことを述べているのは1人だけで、他の3人はすべてうそを述べている。
合格したのが2人であるとき、合格した2人と正直者は誰か。
A Bは合格した。
B AとCのうち1人だけが
214:合格した。 C Bが言っていることは正しい。 D Cが言っていることはうそである。 " TE=gtools::permutations(2,8,0:1,re=T) f <- function(x){ if(sum(x[1:4])!=1 | sum(x[5:8])!=2) return(FALSE) else{ all(c( x[1]==1 & x[6]==1 | x[1]==0 & x[6]==0, x[2]==1 & sum(x[c(5,7)])==1 | x[2]==0 & sum(x[c(5,7)])!=1, x[3]==1 & x[2]==1 | x[3]==0 & x[2]==0, x[4]==1 & x[3]==0 | x[4]==0 & x[3]==1 )) } } colnames(TE)=c(LETTERS[1:4],letters[1:4]) head(TE) # A~D:honest(1) or liar(0), a~d:success(1) or failure(0) TE[apply(TE,1,f),] > TE[apply(TE,1,f),] A B C D a b c d 0 0 0 1 1 0 1 0
215:卵の名無しさん
20/01/24 20:49:51.84 TuIVzM9K.net
"Ａ，Ｂ，Ｃ，Ｄの４人がマラソン競争をしました。そのときの順位について、
４人はそれぞれ次のように言っていますが、１人だけうそをついています。４人の正しい順位を答えなさい。ただし、順位は４人とも異なるものとします。
Ａ君：「私は、２位でも４位でもない」
Ｂ君：「私は、１位だった」
Ｃ君：「私は、Ａ君より速かった」
Ｄ君：「私は、４位だった」
"
TE=gtools::permutations(4,4)
f <- function(x){
sum(c(
x[1]!=2 & x[1]!=4,
x[2]==1,
x[3]<x[1],
x[4]==4
))
}
TE[apply(TE,1,f)==3,]
> TE[apply(TE,1,f)==3,]
[1] 3 2 1 4
216:卵の名無しさん
20/01/24 21:32:57.61 TuIVzM9K.net
A, B、Cは裏口医、正規医、認知症患者からなるが外見で区別できない。
裏口医は常に嘘をつき、正規医は常に真実を語り、認知症患者は嘘も真実も語る。
A 「Bは認知症患者です」
B 「Cは裏口ではありません」
C 「Aは認知症患者ではありません」
A、B、Cは誰かを答えよ。
TE=gtools::permutations(3,3,-1:1)
f <- function(x){
all(c(
x[1]==1 & x[2]==0 | x[1]==-1 & x[2]!=0 | x[1]==0,
x[2]==1 & x[3]!=-1 | x[2]==-1 & x[3]==-1 | x[2]==0,
x[3]==1 & x[1]!=0 | x[3]==-1 & x[1]==0 | x[3]==0
))
}
TE[apply(TE,1,f),]
217:卵の名無しさん
20/01/25 03:16:35 v8nYWWGU.net
ある幼稚園で、砂場で遊んでいたA,B,C,D 部屋で遊んでいたE,F,Gの７人の中に、逆上がりができる子が２人いることが分かっている。
そこで、A～Gに尋ねたところ、それぞれ以下の発言をした。
ただし、７人うち、本当のことを言っているのは２人だけで、あとの５人は間違ったことを言っていた。これらのことから確実にいえるのはどれか。
A：Bは逆上がりできるよ。
B：Aは間違ったことを言っているよ。
C：AもBも２人とも間違ったことを言っているよ。
D：砂場で遊んでいた子の中には逆上がりできる子はいないよ。
E：私は逆上がりできない。
F：逆上がりができるのは２人とも砂場で遊んでいた子だよ。
G：EとFの少なくともどっちかは本当のこと言っているよ。
f <- function(x){
x1=x[1] ; x2=x[2]
re=rep(0,7)
re[x1]=1; re[x2]=1
re
}
(cmb=t(combn(7,2,f)))
gr=as.matrix(expand.grid(1:21,1:21)) ; head(gr)
g <- function(ij){
i=ij[1] ; j=ij[2]
c(cmb[i,],cmb[j,])
}
TE=t(apply(gr,1,g)) ; tail(TE)
218:卵の名無しさん
20/01/25 03:20:58 v8nYWWGU.net
>>204
問題に応じてこういう行列を列挙させるプログラムを組むのがちょっと難しいな。
> tail(TE)
A B C D E F G a b c d e f g
[436,] 0 0 0 1 1 0 0 0 0 0 0 0 1 1
[437,] 0 0 0 1 0 1 0 0 0 0 0 0 1 1
[438,] 0 0 0 1 0 0 1 0 0 0 0 0 1 1
[439,] 0 0 0 0 1 1 0 0 0 0 0 0 1 1
[440,] 0 0 0 0 1 0 1 0 0 0 0 0 1 1
[441,] 0 0 0 0 0 1 1 0 0 0 0 0 1 1
219:卵の名無しさん
20/01/25 03:42:42 v8nYWWGU.net
問題は
1.本当のことを言っているは、１人は砂場で遊んでいた子であり、１人は部屋で遊んでいた子である。
2.Dは本当のことを言っており、Eは間違ったことを言っている。
3.Bは逆上がりができ、間違ったことは言っていない。
4.逆上がりできるのは、FとGの2人である。
5.逆上がりできるのは、１人は砂場で遊んでいた子であり、１人は部屋で遊んでいた子である。
プログラム解は
> TE[apply(TE,1,h),] # A-G:正直1嘘つき0 a-g:逆上がり可能1不能0
A B C D E F G a b c d e f g
[1,] 0 1 0 1 0 0 0 0 0 0 0 1 1 0
[2,] 0 1 0 1 0 0 0 0 0 0 0 1 0 1
よって正しいのは２
220:卵の名無しさん
20/01/25 03:44:35 v8nYWWGU.net
>>204
公務員試験の問題らしいな。医師国試の方が格段に易しいね。
こういう問題をいれて裏口入学の学生を排除すればいいのに。
221:卵の名無しさん
20/01/25 05:44:29 v8nYWWGU.net
"
A～Fの６人が100ｍ競争をした。
同順位の者はなく、それぞれ以下のように述べているが、
真実を述べているのは１位と６位の者のみである。
A：私は４位ではない。
B：私は４位ではない。
C：私はAより上位である。
D：私はBより上位である。
E：Bは１位である。
F：Bは４位である。
"
jun=gtools::permutations(6,6) ;
colnames(jun)=LETTERS[1:6]
tail(jun)
f <- function(x){
all(c(
(x[1]==1|x[1]==6) & x[1]!=4 | !(x[1]==1|x[1]==6) & x[1]==4,
(x[2]==1|x[2]==6) & x[2]!=4 | !(x[1]==1|x[2]==6) & x[2]==4,
(x[3]==1|x[3]==6) & x[3]<x[1] | !(x[3]==1|x[3]==6) & x[3]>=x[1],
(x[4]==1|x[4]==6) & x[4]<x[2] | !(x[4]==1|x[4]==6) & x[4]>=x[2],
(x[5]==1|x[5]==6) & x[2]==1 | !(x[5]==1|x[5]==6) & x[2]!=1,
(x[6]==1|x[6]==6) & x[2]==4 | !(x[6]==1|x[6]==6) & x[2]!=4)
)}
ans=jun[apply(jun,1,f),]
ans
222:卵の名無しさん
20/01/25 06:37:53.05 v8nYWWGU.net
>>208
> f <- function(x){
+ all(c(
+ (x[1]==1|x[1]==6) & x[1]!=4 | !(x[1]==1|x[1]==6) & x[1]==4,
+ (x[2]==1|x[2]==6) & x[2]!=4 | !(x[2]==1|x[2]==6) & x[2]==4,
+ (x[3]==1|x[3]==6) & x[3]<x[1] | !(x[3]==1|x[3]==6) & x[3]>=x[1],
+ (x[4]==1|x[4]==6) & x[4]<x[2] | !(x[4]==1|x[4]==6) & x[4]>=x[2],
+ (x[5]==1|x[5]==6) & x[2]==1 | !(x[5]==1|x[5]==6) & x[2]!=1,
+ (x[6]==1|x[6]==6) & x[2]==4 | !(x[6]==1|x[6]==6) & x[2]!=4)
+ )}
> ans=jun[apply(jun,1,f),]
> ans
A B C D E F
[1,] 1 4 2 5 3 6
[2,] 1 4 3 5 2 6
[3,] 4 1 5 2 6 3
[4,] 4 1 5 3 6 2
[5,] 4 6 5 1 2 3
[6,] 4 6 5 1 3 2
223:卵の名無しさん
20/01/25 06:59:03.02 v8nYWWGU.net
"
A,B,C,D,Eの５人の中におやつをつまみ食いした人が３人います。
事情を聞くと、５人は次のように答えました。
A 私は食べました。でもEも食べました。
B Cは食べました。一緒にDも食べました。
C 私は食べていません。Bも食べていません。
D 私はもちろん食べていません。同じようにCも食べていないですよ。
E 私は食べていません。でもCは食べました。
それぞれ半分は本当のことを言い、半分はうそを言っているとすると、つまみ食
いをした３人はだれですか。
"
f <- function(x){
x1=x[1] ; x2=x[2] ; x3=x[3]
re=rep(0,5)
re[x1]=1 ; re[x2]=1 ; re[x3]=1
re
}
(TE=t(combn(5,3,f)))
colnames(TE)=LETTERS[1:5]
g <- function(x){
all(c(
x[1]==1 & x[5]==0 | x[1]==0 & x[5]==1,
x[3]==1 & x[4]==0 | x[3]==0 & x[4]==1,
x[3]==0 & x[2]==1 | x[3]==1 & x[2]==0,
x[4]==0 & x[3]==1 | x[4]==1 & x[3]==0,
x[5]==0 & x[3]==0 | x[5]==1 & x[3]==1)
)
}
TE[apply(TE,1,g),]
224:卵の名無しさん
20/01/25 07:16:25.22 v8nYWWGU.net
頻繁に使うテンプレートを関数化
TEnr <- function(n,r,one=1){ # n(=5),r(=3)を指定して 0 0 1 1 1から1 1 1 0 0までの順列行列を作る
f=function(x){
re=numeric(n) # 容れ子
re[x]=one # 指定のindexにoneを代入
re
}
t(combn(n,r,f)) # oneを入れる組み合わせに上記関数fを実行して転置
}
225:卵の名無しさん
20/01/25 07:37:15.73 v8nYWWGU.net
A,B,Cの3人がいて、それぞれリンゴを1個か2個持っています。
1個持っている人は必ずウソをつきますが、2個持っている人はウソをつくかどうかわかりませ
ん。1個持っている人が少なくとも1人いるとき、3人がそれぞれ次のように言い
ました。
A BとCのリンゴの合計は２個ではありません。
B AとCのリンゴの合計は３個ではありません。
C AとBのリンゴの合計は４個ではありません。
このとき、A,B,Cが持っているリンゴの個数の合計と、ウソをついている人の組み合わせは？
"
TE=gtools::permutations(2,3,re=T)
f <- function(x){
if(sum(x==1)==0){ return(FALSE)
}else{
all(c(
x[1]==1 & sum(x[2:3])==2 | x[1]==2,
x[2]==1 & sum(x[c(1,3)])==3 | x[2]==2,
x[3]==1 & sum(x[1:2])==4 | x[3]==2)
)
}
}
TE[apply(TE,1,f),]
226:卵の名無しさん
20/01/25 07:38:09.98 v8nYWWGU.net
>>211
0 1 以外の場合もあるからこうした方がいいな
TEnr <- function(n,r,zero=0,one=1){ # n(=5),r(=3)を指定して 0 0 1 1 1から1 1 1 0 0までの順列行列を作る
f=function(x){
re=rep(zero,n) # 容れ子
re[x]=one # 指定のindexにoneを代入
re
}
t(combn(n,r,f)) # oneを入れる組み合わせに上記関数fを実行して転置
}
227:卵の名無しさん
20/01/25 07:42:33.30 v8nYWWGU.net
"
A、B、C、Dの４人の体重を調べました。BとDの平均体重が他の2人のどちらかの
体重と同じで、AはBより17㎏軽く、CとDとの体重の差は12㎏でした。
A 私はCとの体重の差が5㎏です。
B 私はCとの体重の差が22㎏です。
C 私は４人の体重の順番でAの直前または直後です。
D 私はCより重いです。
4人のうち、2人がうそをついているとすると、嘘を言っている２人は誰ですか。
"
228:卵の名無しさん
20/01/25 08:01:55.47 voAAA8hd.net
おまえとおまえのおかん
229:卵の名無しさん
20/01/25 09:14:00.29 v8nYWWGU.net
>>214
Aの体重を50kgとして可能な組み合わせで連立方程式を�
230:｢て体重の行列をつくって2人嘘つきがいる行列と組み合わせて関数でフィルターする。 A=50 # 設定値 B=67 C=D+12 | C=D-12 67+D=2*50 (D=33) | 67+D=2*C w8=rbind( c(50,67,45,33), c(50,67,55,43), c(50,67,11,33), c(50,67,79,91) ) (uso=TEnr(4,2)) (gr=as.matrix(expand.grid(1:4,1:6))) f <- function(i,j){ c(w8[i,],uso[j,]) } TE=t(mapply(f,gr[,1],gr[,2])) colnames(TE)=c(LETTERS[1:4],letters[1:4]) g <- function(x){ tf=c( x[5]==1 & abs(x[1]-x[3])==5 | x[5]==0 & abs(x[1]-x[3])!=5, x[6]==1 & abs(x[2]-x[3])==22 | x[6]==0 & abs(x[2]-x[3])!=22, x[7]==1 & abs(which(sort(x)==x[3])-which(sort(x)==x[1]))==1 | x[7]==0 & abs(which(sort(x)==x[3])-which(sort(x)==x[1]))!=1, x[8]==1 & x[4]>x[3] | x[8]==0 & x[4]<=x[3]) all(tf) } > TE[apply(TE,1,g),] # A-D:体重 a-d:正直1嘘つき0 A B C D a b c d 50 67 55 43 1 0 1 0
231:卵の名無しさん
20/01/25 09:15:27.09 v8nYWWGU.net
>>215
これ面白いだろ？
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
232:卵の名無しさん
20/01/25 11:37:09.68 v8nYWWGU.net
横に並んだいくつかのマスに○、△、×のいずれか１つを記入するとき、隣り合っ
たマスの少なくとも一方が○となるような記入の仕方を考えます。例えば２つの
マスであれば、下図の５通りあります。
○○ ○△ ○× △○ ×○
条件をみたすように７マスを記入するとき、その記入の仕方は何通りあるか。
> sapply(2:13,sim)
[1] 5 11 21 43 85 171 341 683 1365 2731 5461 10923
233:卵の名無しさん
20/01/25 12:22:34.98 bfvHv7+z.net
雨竜郡雨竜町出身介護福祉士でーす 雨竜は寒いでーす
234:卵の名無しさん
20/01/25 12:54:32.90 v8nYWWGU.net
ある人が第１ホールから順にゴルフをプレイして回ったところ、いずれのホールの成績もボギー（＋１）、ダブルボギー（＋２）又はトリプルボギー（＋３）で
あり、あるホールを終えてスコアがちょうど「＋７」となった時点で雨のため中断した。
このとき、各ホールの成績の並びの組合せは何通りあるか。
なお、ボギー、ダブルボギー、トリプルボギーとは、基準打数よりもそれぞれ１打、２打、３打多く打った成績を指しており、基準打数よりも多く打った打数の
総和がスコアである。例えば、第１ホールからの成績の並びがトリプルボギー、ダブルボギー、ボギーであった場合、スコアは（＋３）＋（＋２）＋（＋１）で「＋６」である。
（ I 種国家公務員試験の過去問 ）
数えるのが面倒なのでプログラムにカウントさせる。
for(i in 3:7){
pi=gtools::permutations(3,i,rep=T)
p7=pi[apply(pi,1,sum)==7,]
count[i]=ifelse(is.vector(p7),1,nrow(p7))
}
sum(count)
> sum(count)
[1] 44
235:卵の名無しさん
20/01/25 12:55:04.80 v8nYWWGU.net
>>219
女子大生にフェラをしてもらうには何と言えばいい？
236:卵の名無しさん
20/01/25 13:13:23.67 v8nYWWGU.net
>>220
階段登りの場合の数の計算に関数化しておこう。
# Σxi=n xi<=r 一度にｒ(=3)段まで上がれる時、総数ｎ(=7)段上がる場合の数
stairway <- function(n,3){
count=numeric(n)
m=ceiling(n/r)
for(i in m:n){
pi=gtools::permutations(r,i,rep=T)
pn=pi[apply(pi,1,sum)==n,]
count[i]=ifelse(is.vector(pn),1,nrow(pn))
}
sum(count)
}
237:卵の名無しさん
20/01/25 13:33:02.20 v8nYWWGU.net
１３段の階段を上るさい、１段上るか、２段上るか、３段上るかのいずれかを組合わせて上るこ
ととする。上り方は何通りあるか。
238:卵の名無しさん
20/01/25 13:35:18.76 v8nYWWGU.net
> stairway(13,3)
[1] 1705
おい、ド底辺。
書き出してカウントしてくれ。
オマエの底辺頭脳ではフェラ問題の答は無理。
数くらい数えられるだろ？
これ↓面白いだろ？
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
239:卵の名無しさん
20/01/25 14
240::21:14 ID:mWDUfEYx.net
241:卵の名無しさん
20/01/25 14:22:22 v8nYWWGU.net
考えるよりプログラムを入力して計算させた方が速いな。
"
Ａ～Ｆの６人が旅行に出かけた。宿泊所では、３人部屋と４人部屋に分かれて泊
まることになった。このとき、ＡとＢの２人が必ず同じ部屋に泊まる部屋割りは
何通りあるか。
"
pm=gtools::permutations(2,6,0:1,re=T)
f <- function(x){
if(sum(x)>4|sum(x==0)>3) return(FALSE)
else if(x[1]==x[2]) return(TRUE)
return(FALSE)
}
sum(apply(pm,1,f))
> sum(apply(pm,1,f))
[1] 15
>
242:卵の名無しさん
20/01/25 14:24:36 v8nYWWGU.net
>>225
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
女子大生にフェラしてもらった経験なければ、こっちに答えるのでもいいぞ。
合格と寄付金は不可分を掲げるド底辺医大に入試前に寄付をした受験生の親が寄付金だけとられて不合格になるのでないかと心配して理事長と談判した。
合格と寄付金は不可分という大学の方針は常に正しいが、理事長は常に嘘をつく。
談判では、
「お子様が面接で正しい主張をすれば寄付金返却か裏口合格になる、間違った主張であれば寄付金没収で不合格になる」
と言われた。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？
243:卵の名無しさん
20/01/25 14:33:18.07 mWDUfEYx.net
情報収集の圧倒的な差w
244:卵の名無しさん
20/01/25 14:43:02.34 mWDUfEYx.net
いいことを１つ教えてやろう Xビデオ で検索してみろ お前の欲しい情報がいっぱいあるぞ！w
245:卵の名無しさん
20/01/25 14:44:05.94 v8nYWWGU.net
１円玉３枚、５円玉１枚、１０円玉２枚、１００円玉２枚の合計２２８円を持っ
ている。いま、おつりをもらわないでちょうどぴったり支払いのできる金額は何
通りあるか。 ０円はカウントしない。
一行プログラムで答がでる。
> nrow((expand.grid(0:3,c(0,5),c(0,10,20),c(0,100,200))))-1
[1] 71
こういう解き方していると馬鹿になるな。
それすらできない馬鹿がド底辺シリツ医。
女子大生からフェラも引き出せないｗｗｗ
246:卵の名無しさん
20/01/25 14:46:11.67 v8nYWWGU.net
>>228
赤の他人であることもわからんドアホがオマエだよ。
フェラ問題はとけたのか？
247:卵の名無しさん
20/01/25 14:52:01.09 v8nYWWGU.net
>>228
圧倒的な頭脳の差！
シリツ医の頭脳は昆虫並み。
決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
248:卵の名無しさん
20/01/25 14:53:55.72 mWDUfEYx.net
なっ オナニーネタにちょうどいいだろ おまえにはw
249:卵の名無しさん
20/01/25 14:56:42.15 v8nYWWGU.net
AからFの６人が、横一列になって写真撮影をする。仲良しのAとBは必ず隣り同士
であり、不仲のAとCは必ず隣り合わないようにするような並べ方は何通りか。
３行で終わり
pm=gtools::permutations(6,6)
f <- function(x)abs(which(x==1)-which(x==2))==1 & abs(which(x==1)-which(x==3))!=101
nrow(pm[apply(pm,1,f),])
実行結果
> nrow(pm[apply(pm,1,f),])
[1] 192
250:卵の名無しさん
20/01/25 14:57:29.46 v8nYWWGU.net
>>233
いや、女子大生ファラの方がいいぞ。
251:卵の名無しさん
20/01/25 15:03:51 mWDUfEYx.net
よかったね！ずりネタ１つ増えてw
252:卵の名無しさん
20/01/25 15:15:54 mWDUfEYx.net
福井県越前市瓜生町 オナニー好きのあなたへ 気持ちよかったですか？
253:卵の名無しさん
20/01/25 15:27:26 v8nYWWGU.net
ある職場は、Ａ～Ｆの６人の職員で構成され、このうちＥおよびＦは新人職員で
ある。新人職員は、1人だけまたは２人だけでは外出または留守番をしないとす
るとき、外出する職員の組合せは、何通りあるか。
pm=gtools::permutations(2,6,0:1,re=T)
f <- function(x){
any(c(
(sum(x==1)==1 & (x[1]==1|x[2]==1)) | (sum(x==1)==2 & x[1]==1 & x[2]==1),
(sum(x==0)==1 & (x[1]==0|x[2]==0)) | (sum(x==0)==2 & x[1]==0 & x[2]==0)
))
}
pm[apply(pm,1,f),]
nrow(pm)- nrow(pm[apply(pm,1,f),]) - 1
実行結果
> nrow(pm)- nrow(pm[apply(pm,1,f),]) - 1
[1] 57
254:卵の名無しさん
20/01/25 15:29:50 v8nYWWGU.net
> stairway(13,3)
[1] 1705
おい、ド底辺。
書き出してカウントしてくれ。
オマエの底辺頭脳ではフェラ問題の答は無理。
数くらい数えられるだろ？
これ↓面白いだろ？
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
255:卵の名無しさん
20/01/25 15:30:56 v8nYWWGU.net
女子大生にフェラしてもらった経験なければ、こっちに答えるのでもいいぞ。
合格と寄付金は不可分を掲げるド底辺医大に入試前に寄付をした受験生の親が寄付金だけとられて不合格になるのでないかと心配して理事長と談判した。
合格と寄付金は不可分という大学の方針は常に正しいが、理事長は常に嘘をつく。
談判では、
「お子様が面接で正しい主張をすれば寄付金返却か裏口合格になる、間違った主張であれば寄付金没収で不合格になる」
と言われた。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？
256:卵の名無しさん
20/01/25 15:42:10 bfvHv7+z.net
気持ちよかったですか？
257:卵の名無しさん
20/01/25 15:43:39 bfvHv7+z.net
福井県越前市瓜生町
258:卵の名無しさん
20/01/25 15:56:00.38 v8nYWWGU.net
〇、△、×のいずれかのマークが１つずつかかれたカードがあります。〇３枚、
△２枚、×２枚の全７枚のカードを、同じマークのカードが隣り合わないように横
一列に並べます。
一番左に○を並べた時、残りの６枚の並べ方は何通りあるか。
割と面倒
v=rep(1:3,c(3,2,2))
pm=permutations(7,7,v,set=F)
pm=unique(pm)
f<-function(x){
if(x[1]!=1) return(FALSE)
else{
if(
1
259: %in% diff(which(x==1))| 1 %in% diff(which(x==2)) | 1 %in% diff(which(x==3)) ) return(FALSE) } return(TRUE) } nrow(pm[apply(pm,1,f),]) 実行結果 > nrow(pm[apply(pm,1,f),]) [1] 20
260:卵の名無しさん
20/01/25 15:58:02.20 v8nYWWGU.net
>>241
答がわからん馬鹿なの？
決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
261:卵の名無しさん
20/01/25 16:00:45.80 bfvHv7+z.net
そうかそうか 気持ちよかったかー 一日一善したw
262:卵の名無しさん
20/01/25 16:08:31 bfvHv7+z.net
それにしても彼らの情報収集能力は秀逸だ おまえはある偉い人々を怒らせたんだよ
しーらないっとw
263:卵の名無しさん
20/01/25 16:35:06 bfvHv7+z.net
福井県越前市瓜生町 おまえは通っていた病院で問題起こしたろw
264:卵の名無しさん
20/01/25 16:46:11 v8nYWWGU.net
>>246
未だに他人を同定しているどアホがオマエだよ。
まさに昆虫並みの頭脳だな。
ここまでアホだから、藤田か聖マだろ？
265:卵の名無しさん
20/01/25 16:47:07 v8nYWWGU.net
これにも答がだせないの？
馬鹿なんだろ。
決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
266:卵の名無しさん
20/01/25 18:10:05.47 v8nYWWGU.net
A～Gの７人を、３人、２人、２人の３グループに分ける分け方は何通りあるか。
267:卵の名無しさん
20/01/25 18:35:05 v8nYWWGU.net
赤、青、白の３色のカードが５枚ずつ、計１５枚あります。各色ともに、１～５
までの数が１つずつ書かれています。この１５枚の中から３枚を選びます。選ん
だ３枚に書かれた数の和が１２になるような選び方は何通りありますか。
v=rep(1:5,3)
cmb=gtools::combinations(15,3)
head(cmb)
mat=matrix(v[cmb],nrow=3,byrow=T)
dat=t(mat)
head(dat)
f=function(x) sum(x)==12
nrow(dat[apply(dat,1,f),])
> nrow(dat[apply(dat,1,f),])
[1] 37
>
268:卵の名無しさん
20/01/25 19:12:41 bfvHv7+z.net
福井県越前市瓜生町 おまえのことは全部バレバレなんだよw
269:卵の名無しさん
20/01/25 19:48:33 bfvHv7+z.net
国立私立どちらも含めて多くの医者を敵を作っちゃったね ケケケ
270:卵の名無しさん
20/01/25 19:53:06.19 v8nYWWGU.net
これにも答がだせないの？
馬鹿なんだろ。
決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？
271:卵の名無しさん
20/01/25 19:53:37.17 v8nYWWGU.net
>>246
未だに他人を同定しているどアホがオマエだよ。
まさに昆虫並みの頭脳だな。
ここまでアホだから、藤田か聖マだろ？
あたりみたいだな。
272:卵の名無しさん
20/01/25 20:02:55.47 bfvHv7+z.net
いや俺は長崎大学医学部看護学科出身の看護師だw
273:卵の名無しさん
20/01/25 20:16:08.55 bfvHv7+z.net
伝説の福井県越前市瓜生町 ケケケ
274:卵の名無しさん
20/01/25 20:23:52.34 bfvHv7+z.net
俺の知り合いに聞いたら あーアイツか あのクルクルパーだな だってさ 人気があるんだね うらやましいよ ケケケ
275:卵の名無しさん
20/01/25 20:40:43 bfvHv7+z.net
いろんな病院にかかって ずいぶんと医者や看護師を困らせているんだな もっと荒らしていいよw
276:卵の名無しさん
20/01/25 22:17:33 bfvHv7+z.net
精神障害者手帳持ってる？
277:卵の名無しさん
20/01/26 04:22:46 6QosvB6P.net
>>258
朱の他人を中傷しているオマエがアホ。
オマエの知り合いだから裏口入学容疑者だろ。
278:卵の名無しさん
20/01/26 05:37:56 6QosvB6P.net
f <- function(a) {
x=a-floor(a)
floor(a)+ifelse(x<0.5,floor(x),ceiling(x))
}
sim <- function(n=3,L=0,U=10){
x=runif(n,L,U)
y= numeric(n)
for(i in 1:n) y[i]=f(x[i])
f(sum(x)) != sum(y)
}
mean(replicate(1e3,sim()))
sim2 <- function(n){
pn=mean(replicate(1e3,sim(n)))
sqrt(n)*pn
}
sim3 <- function(n,L=0,U=100){
pn=mean(replicate(1e4,sim(n)))
pn
}
Sim3=Vectorize(sim3)
Sim3(3:10)
279:卵の名無しさん
20/01/26 05:43:46.84 6QosvB6P.net
f <- function(a) {
x=a-floor(a)
floor(a)+ifelse(x<0.5,floor(x),ceiling(x))
}
sim <- function(n=3,k=1e4){
sub <- function(n){
x=runif(n)
y= numeric(n)
for(i in 1:n) y[i]=f(x[i])
f(sum(x)) != sum(y)
}
mean(replicate(k,sub(n)))
}
sim()
280:卵の名無しさん
20/01/26 05:46:04.35 6QosvB6P.net
任意の3つの実数を四捨五入して整数にした上で和を取るのと、
和を取ってから四捨五入して整数にする場合で数値が異なる確率を求めよ。
ただし、3つの実数の小数部は[0,1)区間で独立な一様分布とする。
(アクチュアリー資格試験問題 - 一部改正)
281:卵の名無しさん
20/01/26 06:20:45 6QosvB6P.net
四捨五入で思い出した。
IEEE では任意の桁
282:位置における丸め処理法を次のように定めている（ RjpWiki の記事，『工学のためのデータサイエンス入門』 間瀬 茂 他 著 (数理工学社) から引用した）． 一番近い丸め結果候補が1つだけなら，その数に丸める． 一番近い丸め結果候補が2つある場合は，末尾が偶数のものに丸める（五入ばかりでなく五捨もあり得る！）． 丸め処理は1段階で行なわなければならない． 2.5の丸めは2 3.5の丸めは4 これって丸め誤差が小さくなると説明されていた。
283:卵の名無しさん
20/01/26 07:17:36 6QosvB6P.net
f45 <- function(a) {
x=a-floor(a)
floor(a)+ (x<=0.5)
}
f45(0.5)
round(0.5)
sim <- function(n=100,print=T){
x=runif(n,0,10)
M=mean(x)
tf <- abs(M-mean(round(x))) < abs(M-mean(f45(x)))
if(tf & print) cat(x,'\n')
return(tf)
}
tf=FALSE
k=1e3
while(tf==FALSE & k>0){
tf=sim(10)
k=k-1
}
tf
mean(replicate(k,sim()))
284:卵の名無しさん
20/01/26 07:33:47 T35MsftF.net
f45 <- function(a) {
x=a-floor(a)
floor(a)+ (x>=0)
}
f45(0.5) ; f45(1.5) ; f45(1.4999)
round(0.5) ; round(1.5) ; round(1.499)
sim <- function(n=100,print=T){
x=runif(n,0,10)
M=mean(x)
tf <- abs(M-mean(round(x))) < abs(M-mean(f45(x)))
if(tf & print) cat(x,'\n')
return(tf)
}
tf=FALSE
k=1e3
while(tf==FALSE & k>0){
tf=sim(10)
k=k-1
}
tf
mean(replicate(k,sim()))
285:卵の名無しさん
20/01/26 07:34:40 T35MsftF.net
f45 <- function(a) {
x=a-floor(a)
floor(a)+ (x>=0.5)
}
f45(0.5) ; f45(1.5) ; f45(1.4999)
round(0.5) ; round(1.5) ; round(1.499)
sim <- function(n=100,print=T){
x=runif(n,0,10)
M=mean(x)
tf <- abs(M-mean(round(x))) < abs(M-mean(f45(x)))
if(tf & print) cat(x,'\n')
return(tf)
}
tf=FALSE
k=1e3
while(tf==FALSE & k>0){
tf=sim(10)
k=k-1
}
tf
mean(replicate(k,sim()))
286:卵の名無しさん
20/01/26 07:35:57 T35MsftF.net
f45 <- function(a) {
x=a-floor(a)
floor(a)+ (x>=0.5)
}
f45(0.5) ; f45(1.5) ; f45(1.4999)
round(0.5) ; round(1.5) ; round(1.499)
sim <- function(n=100,print=T){
x=runif(n,0,10)
M=mean(x)
tf <- abs(M-mean(round(x))) < abs(M-mean(f45(x)))
if(tf & print) cat(x,'\n')
return(tf)
}
tf=FALSE
k=1e3
while(tf==FALSE & k>0){
tf=sim(10, print=F)
k=k-1
}
tf
mean(replicate(k,sim(print=F)))
287:卵の名無しさん
20/01/26 07:37:45 T35MsftF.net
f45 <- function(a) {
x=a-floor(a)
floor(a)+ (x>=0.5)
}
f45(0.5) ; f45(1.5) ; f45(1.4999)
round(0.5) ; round(1.5) ; round(1.499)
sim <- function(n=100,print=T){
x=runif(n,0,10)
M=mean(x)
tf <- abs(M-mean(round(x))) < abs(M-mean(f45(x)))
if(tf & print) cat(x,'\n')
return(tf)
}
tf=FALSE
k=1e3
while(tf==FALSE & k>0){
tf=sim(10, print=F)
k=k-1
}
tf
mean(replicate(k,sim(10, print=F)))
288:卵の名無しさん
20/01/26 09:40:02 LvdMCsAi.net
rm(list=ls())
f45 <- function(a) {
x=a-floor(a)
floor(a)+ (x>=0.5)
}
f45(0.5) ; f45(1.5)
round(0.5) ; round(1.5)
sim <- function(n=100,print=T){
x=runif(n,0,10)
M=mean(x)
tf <- abs(M-mean(round(x))) < abs(M-mean(f45(x)))
if(tf & print) cat(x,'\n')
return(tf)
}
tf=FALSE
k=1e5
while(tf==FALSE & k>0){
tf=sim(100, print=F)
k=k-1
}
k=1e5
mean(replicate(k,sim()))
289:卵の名無しさん
20/01/26 11:22:51 LvdMCsAi.net
>>264
rm(list=ls())
f45 <- function(a) { # 四捨五入
x=a-floor(a) # floor(a):aを超えない整数ガウス記号[x]と同じ,x:小数部分をxに入れる
floor(a)+ (x>=0.5) # xが0.5以上なら１をそうでないなら０を加える
}
f45(0.49) ; f45(1.5)
sim <- function(n=3,k=1e4){
sub <- function(n){
x=runif(n) # 一様分布乱数(実数)ｎ個の配列
y=numeric(n) # y：四捨五入での整数を入れる配列
for(i in 1:n) y[i]=f45(x[i]) # ｘの各実数を四捨五入してyに入れる
f45(sum(x))!=sum(y) # xの総和の四捨五入数とyの総和が異なればTRUEを返す
}
mean(replicate(k,sub(n))) # k個のシミュレーションでのTRUEの頻度を返す
}
#
n=3:100
pn=sapply(n,function(n) sim(n,1e5))
plot(n,pn,bty='l',pch=19)
290:卵の名無しさん
20/01/26 15:40:07.28 LvdMCsAi.net
四捨五入
C++ Java JavaScript Python2 PHP
IEEE
C# R Python3 Numpy
# 四捨五入 vs IEEE
comp <- function(n=10,print=T){
x=sample((1:9)/10,n,rep=T) # 0.1,0.2,..,0.9から重複を許してn個選ぶ
M=mean(x)
dif = abs(M-mean(round(x))) - abs(M-mean(f45(x))) # round後の平均と四捨五入後の平均の差
if(dif!=0 & print) cat(dif<0,' : ',x,'\n') # 差があれば表示 round後が小さければTRUE
return(dif)
}
291:卵の名無しさん
20/01/26 16:25:49.23 LvdMCsAi.net
>>265
実際に実行して確認してみる。
# 四捨五入 vs ＩＥＥＥ with FUN(=mean) for x(=c(-0.9,-0.8,...,0.8,0.9))
comp <- function(n=10,FUN=mean,x=sample((-9:9)/10,n,rep=T),print=T){
X=FUN(x)
dif = abs(X-FUN(round(x))) - abs(X-FUN(f45(x))) # round後の実行と四捨五入後の実行の差
if(dif!=0 & print) cat(dif<0,' : ',sort(x),'\n') # 差があれば表示 round後が小さければTRUE
return(dif) # 差を返す dif<0:round優位 dif>0:四捨五入優位
}
comp()
k=1e4
# mean
re=replicate(k,comp(print=F))
c(ＩＥＥＥ=mean(re<0),四捨五入=mean(re>0),引き分け=mean(re==0))
# prod
re=replicate(k,comp(FUN=prod,print=F))
c(ＩＥＥＥ=mean(re<0),四捨五入=mean(re>0),引き分け=mean(re==0))
# squared sum
comp(FUN=function(x) sum(x^2))
re=replicate(k,comp(FUN=function(x)sqrt(sum(x^2)),print=F))
c(ＩＥＥＥ=mean(re<0),四捨五入=mean(re>0),引き分け=mean(re==0))
292:卵の名無しさん
20/01/26 16:28:08.22 LvdMCsAi.net
> k=1e5
> # mean
> re=replicate(k,comp(print=F))
> c(ＩＥＥＥ=mean(re<0),四捨五入=mean(re>0),引き分け=mean(re==0))
ＩＥＥＥ 四捨五入 引き分け
0.25238 0.15829 0.58933
> # prod
> re=replicate(k,comp(
293:FUN=prod,print=F)) > c(ＩＥＥＥ=mean(re<0),四捨五入=mean(re>0),引き分け=mean(re==0)) ＩＥＥＥ 四捨五入 引き分け 0.00041 0.00000 0.99959 > # squared sum > comp(FUN=function(x) sum(x^2)) [1] 0 > re=replicate(k,comp(FUN=function(x)sqrt(sum(x^2)),print=F)) > c(ＩＥＥＥ=mean(re<0),四捨五入=mean(re>0),引き分け=mean(re==0)) ＩＥＥＥ 四捨五入 引き分け 0.40377 0.01346 0.58277 平均くらいでは大差はないが、統計の内部計算でつかう平方和に関しては IEEEの丸めの方が圧倒的に誤差が少ないなぁ。 国立卒とﾄﾞ底辺シリツ卒くらいの差がみてとれる。
294:卵の名無しさん
20/01/27 03:22:11 4ee1GOdM.net
>>264
rm(list=ls())
pn <- function(n=3,k=1e3){ # 1-pn = -1/2 < Σ(1→n)xi < 1/2, where xi~U(-1/2,1/2)
f <- function(n) abs(sum(replicate(n,runif(1,-1/2,1/2)))) < 1/2
1-mean(replicate(k,f(n)))}
pn(3,1e5)
y=sapply(1:100,pn)
plot(n,y)
295:卵の名無しさん
20/01/27 05:43:38.12 4ee1GOdM.net
>>260
おまえこそ、知恵遅れの障害者じゃないの？
296:卵の名無しさん
20/01/27 07:53:24 4ee1GOdM.net
"
任意の3つの実数を四捨五入して整数にした上で和を取るのと、
和を取ってから四捨五入して整数にする場合で数値が異なる確率1-pnを求めよ。
ただし、3つの実数の小数部は[0,1)区間で独立な一様分布とする。
"
fpn <- function(n=3,k=1e3){ # 1-pn = -1/2 < Σ(1→n)xi < 1/2, where xi~U[-1/2,1/2)
f <- function(n) abs(sum(replicate(n,runif(1,-1/2,1/2)))) < 1/2
mean(replicate(k,f(n)))
}
y=sapply(1:10,fpn)
plot(1:10,y,pch=19,bty='l')
fpn <- function(n=3,k=1e3){ # 1-pn = -1/2 < Σ(1→n)xi < 1/2, where xi~Unif[-1/2,1/2)
f <- function(n){ # 1 -pn = (n-1)/2 < Σ(1→n)Xi, where Xi~Unif[0,1)
s=sum(replicate(n,runif(1)))
(n-1)/2<s & s<(n+1)/2 # is round after sum(RAS) == sum after round(SAR)?
}
mean(replicate(k,f(n))) # proportion of RAS == SAR
}
n=3
pn=replicate(1e4,fpn(n,k=1e3))
hist(pn,freq=F,col='lightblue')
mean(pn) ; sd(pn)
s=replicate(1e5,sum(replicate(n,runif(1))))
mean(s)
sd(s)
m=n/2
sd=sqrt(n/12)
curve(dnorm(x,m,n),0,n,bty='l')
1- (pnorm((n+1)/2,m,sd)-pnorm((n-1)/2,m,sd))
297:卵の名無しさん
20/01/27 07:53:53 4ee1GOdM.net
# Qn=1-Pn Xi~Unif[-1/2,1/2) Pr[Σ(1:n)Xi∈[-1/2,1/2)
"
区間I=[-1/2,1/2)で一様分布する独立なn個の実数を足してもIに含まれる確率をP[n]とすると
P[n]=(1/π)∫[-∞,∞](sinx/x)^(n+1)dx
である。
"
Pn <- function(n){
if(n==0) return(1)
if(n==1) return(1)
f<-function(x) (sin(x)/x)^(n+1)/pi
integrate(f,-Inf,Inf)$value
}
sapply(0:10,Pn)
#
f <- function(x,n=3) (sin(x)/x)^(n+1)/pi
integrate(f,-Inf,Inf)
vf=Vectorize(f)
curve(vf,-pi,pi,bty='l')
298:卵の名無しさん
20/01/27 17:03:42 54kq03+D.net
"AからHの8人はそれぞれ正直者か嘘つきであり、誰が正直者か嘘つきかはお互いに知っている。
A,B,C,D,Eは嘘つきなら必ず嘘をつくが、F,G,Hは嘘つきでも正しいことを言う場合がある。
次の証言から誰を確実に正直者と断定できるか？
A「嘘つきの方が正直者より多い」
B「Hは嘘つきである」
C「Bは嘘つきである」
D「CもFも嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも正直者である」
"
TE=expand.grid(0:1,0:1,0:1,0:1,0:1,0:1,0:1,0:1)
colnames(TE)=LETTERS[1:8]
tail(TE)
TE=expand.grid(0:1,0:1,0:1,0:1,0:1,0:1,0:1,0:1)
colnames(TE)=LETTERS[1:8]
f <- function(x){
all(c(
x[1]==1 & sum(x==0)>sum(x==1) | x[1]!=1 & !(sum(x==0)>sum(x==1)),
x[2]==1 & x[8]==0 | x[2]!=1 & x[8]!=0,
x[3]==1 & x[2]==0 | x[3]!=1 & x[2]!=0 ,
x[4]==1 & (x[3]==0 & x[6]==0) | x[4]!=1 & !(x[3]==0 & x[6]==0),
x[5]==1 & sum(x==0)>=1 | x[5]!=1 & !(sum(x==0)>=1),
x[6]==1 & sum(x==0)>=2 | x[6]==0,
x[7]==1 & x[5]==0 | x[7]==0,
x[8]==1 & (x[1]==1 & x[6]==1) | x[8]==0
))
}
TE[apply(TE,1,f),]
299:卵の名無しさん
20/01/29 13:20:39.63 TTb+Qyhg.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
300:卵の名無しさん
20/01/29 13:20:46.52 TTb+Qyhg.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
301:卵の名無しさん
20/01/29 13:20:54.19 TTb+Qyhg.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
302:卵の名無しさん
20/01/29 17:45:15 qmYv2xxr.net
インバウンドとかやめるべき。リスクが高くてリターンが少い。中国様はSARSを起こして、またこれだよ。野性動物食うような連中が移動すれば常にこういうリスクを起こ�
303:ｷ 同じような習慣ある連中は、世界にはたくさんある。今までは貧乏国で旅行なんか出来なかっただけ。アフリカとかからも大量に旅行者が来るようになれば、奇病の感染とか有り得る 先進国通しが行き来すれば充分。リスクあるから無理にインバウンドなんかやる必要ない。大体、日本は封鎖したりで封じ込めとかできないしな。今まではラッキーだっただけ SARSの時は中国人は、まだ貧乏で旅行なんて殆んどいけなかったし、国内の移動も少なかったから広がらなかった。今は違うよ。あと、SARSの時は、飛行機消毒したりかなりやっていた 今回は、感染者がでずらかったSARSに対するよりも初動を失敗している。もう、かなり感染者が出てると考えた方がいいわ
304:卵の名無しさん
20/01/30 07:09:31 m/EdOA9B.net
まず入れない
増やさない
封じ込める
これができない政治家はパンデミックなったときの主犯です。
305:卵の名無しさん
20/01/30 07:39:24.72 m/EdOA9B.net
23日 中国政府 武漢市を封鎖
24日 安倍総理、春節中国人観光客歓迎祝辞
曰く
春節に際して、そしてまた、オリンピック・パラリンピック等の機会を通じて、更に多くの中国の皆様が訪日されることを楽しみにしています。その際、ぜひ東京以外の場所にも足を運び、その土地ならではの日本らしさを感じて頂ければ幸いです。
25日、26日安倍総理は私邸で待機
28日 中国政府、中国人の海外団体旅行を禁止
封じ込める気0だったんだからそりゃ拡散止まらんだろ
目先の春節インバウンド目当てでもっと長期の観光客利益失ったな
306:卵の名無しさん
20/01/31 01:13:23 24QiZJ0Y.net
f <- function(n,prec=10000){ # Σ 1/kを既約分数表示する
if(n==1){
cat(' ',1,'\n')
invisible(1)
}else{
GCD <- function(a,b){ # ユークリッドの互除法
r = a%%b # a=bq+r ⇒ a%%b=b%%rで最大公約数表示
while(r!=0){a = b ; b = r ; r = a%%b}
b }
library(Rmpfr)
one = mpfr(1, prec) # 1(one)を10000桁精度に設定
nn = 1:n # nn : 1 2 3 ... n
nume=numeric(n) # 分子の容器
for(i in nn) nume[i] = prod(nn[-i])*one # nnからi番目の要素を除いて乗算し精度アップ
nume = mpfr2array(nume, dim = length(nume)) # mpfr2arrayで加算を可能にする
Nume = sum(nume) # numeの総和を計算して分子に
Deno=factorialZ(n) # 分母 n! = factorial(n*one)
gcd = GCD(Nume,Deno) # Numerator/Denominator約分するため最大公約数を計算
res=list(nume = Nume/gcd,deno=Deno/gcd,ratio=as.numeric(Nume/Deno)) # 最大公約数で除算して
# 分数表示 give.head=FALSEでheader除去,digits.dで桁数を指定
# capture.outputで変数に取り込み
nm = capture.output(str(res$nume, give.head=FALSE,digits.d = prec))
dn = capture.output(str(res$deno, give.head=FALSE,digits.d = prec))
cat(pa
307:ste0(nm,'/',dn,'\n')) invisible(res) }} for(i in 1:50) f(i) f(77)
308:卵の名無しさん
20/01/31 04:18:10 Uud6Fjwu.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
309:卵の名無しさん
20/01/31 04:18:17 Uud6Fjwu.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
310:卵の名無しさん
20/01/31 04:18:26 Uud6Fjwu.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
311:卵の名無しさん
20/01/31 07:20:57.36 24QiZJ0Y.net
>>287
ｆ ＜－ ｆｕｎｃｔｉｏｎ（ｎ，ｐｒｅｃ＝１００００）｛ ＃ Σ １／ｋを既約分数表示する
ｉｆ（ｎ＝＝１）｛
ｃａｔ（＇ ＇，１，＇\ｎ＇）
ｉｎｖｉｓｉｂｌｅ（１）
｝ｅｌｓｅ｛
ＧＣＤ ＜－ ｆｕｎｃｔｉｏｎ（ａ，ｂ）｛ ＃ ユークリッドの互除法
ｒ ＝ ａ％％ｂ ＃ ａ＝ｂｑ＋ｒ ⇒ ａ％％ｂ＝ｂ％％ｒで最大公約数表示
ｗｈｉｌｅ（ｒ！＝０）｛ａ ＝ ｂ ； ｂ ＝ ｒ ； ｒ ＝ ａ％％ｂ｝
ｂ ｝
ｌｉｂｒａｒｙ（Ｒｍｐｆｒ）
ｏｎｅ ＝ ｍｐｆｒ（１， ｐｒｅｃ） ＃ １（ｏｎｅ）を１００００桁精度に設定
ｎｎ ＝ １：ｎ ＃ ｎｎ ： １ ２ ３ ．．． ｎ
ｎｕｍｅ＝ｎｕｍｅｒｉｃ（ｎ） ＃ 分子の容器
ｆｏｒ（ｉ ｉｎ ｎｎ） ｎｕｍｅ［ｉ］ ＝ ｐｒｏｄ（ｎｎ［－ｉ］）＊ｏｎｅ ＃ ｎｎからｉ番目の要素を除いて乗算
ｎｕｍｅ ＝ ｍｐｆｒ２ａｒｒａｙ（ｎｕｍｅ， ｄｉｍ ＝ ｌｅｎｇｔｈ（ｎｕｍｅ）） ＃ ｍｐｆｒ２ａｒｒａｙで加算を可能に
Ｎｕｍｅ ＝ ｓｕｍ（ｎｕｍｅ） ＃ ｎｕｍｅの総和を計算して分子に
Ｄｅｎｏ＝ｆａｃｔｏｒｉａｌＺ（ｎ） ＃ 分母 ｎ！ ＝ ｆａｃｔｏｒｉａｌ（ｎ＊ｏｎｅ）
ｇｃｄ ＝ ＧＣＤ（Ｎｕｍｅ，Ｄｅｎｏ） ＃ 最大公約数を計算
ｒｅｓ＝ｌｉｓｔ（ｎｕｍｅ ＝ Ｎｕｍｅ／ｇｃｄ，ｄｅｎｏ＝Ｄｅｎｏ／ｇｃｄ，ｒａｔｉｏ＝ａｓ．ｎｕｍｅｒｉｃ（Ｎｕｍｅ／Ｄｅｎｏ）） ＃ 最大公約数で除算して
＃ ｃａｐｔｕｒｅ．ｏｕｔｐｕｔで切り取りｓｕｂｓｔｒで［１］を除去
ｎｍ ＝ｓｕｂｓｔｒ（ｃａｐｔｕｒｅ．ｏｕｔｐｕｔ（ｒｅｓ＄ｎｕｍｅ）［２］，５，ｎｃｈａｒ（ｒｅｓ＄ｎｕｍｅ）） ＃ ＂［１］ １２３４．．＂ 文字列の５文字目から最後まで
ｄｎ ＝ｓｕｂｓｔｒ（ｃａｐｔｕｒｅ．ｏｕｔｐｕｔ（ｒｅｓ＄ｄｅｎｏ）［２］，５，ｎｃｈａｒ（ｒｅｓ＄ｄｅｎｏ））
ｃａｔ（ｐａｓｔｅ０（ｎｍ，＇／＇，ｄｎ，＇\ｎ＇））
ｉｎｖｉｓｉｂｌｅ（ｒｅｓ）
｝｝
312:卵の名無しさん
20/01/31 07:25:03.91 24QiZJ0Y.net
#Reed-Frost Model
ReedFrost=function(
p=0.04,
N=100,
T=40)
{
q=1-p
I=numeric(T)
S=numeric(T)
I[1]=1
S[1]=N-I[1]
for(t in 1:(T-1)){
I[t+1]=S[t]*(1-q^I[t])
S[t+1]=S[t]-I[t+1]
}
plot(1:T,I,type="l",lwd=2, ylim=c(0,N),xlab="time",ylab="persons",main=paste("Reed-Frost Model p= ",p))
lines(S,lty=2,col=2,lwd=2)
lines(N-S,lty=3,col=3,lwd=2)
legend("topright",bty="n",legend=c("Infected","Susceptible","Immunized"),lty=1:3,col=1:3,lwd=2)
}
313:卵の名無しさん
20/01/31 08:05:29 k8jciAtw.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
314:卵の名無しさん
20/01/31 08:05:37 k8jciAtw.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
315:卵の名無しさん
20/01/31 08:05:44 k8jciAtw.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
316:卵の名無しさん
20/01/31 12:55:23 24QiZJ0Y.net
＞安倍首相「これまで実施してきた水際作戦などのフェーズをもう一段引き上げる必要がある
安倍ぱよ三 春節祝辞抜粋
春節に際して、そしてまた、オリンピック・パラリンピック等の機会を通じて、更に多くの中国の皆様が訪日されることを楽しみにしています。
何とんでも無い事
317:抜かしてんだ？ゲリグソ野郎が https://www.cn.emb-japan.go.jp/itpr_ja/00_000351.html ここnot foundになって安倍がページを消して証拠隠滅したぞｗｗ
318:卵の名無しさん
20/01/31 16:02:18 24QiZJ0Y.net
まだ騙されている人がいそうだけど、
インバウンド、インバウンドと騒がれているけど、
インバウンドの消費額なんてＧＤＰの１％にも満たない、
家計の総消費額に対しても、たったの1.5%。
それを「観光立国」だの、
あたかも日本がインバウンドに支えられているかのごとく国民を騙して、
ごく一部の人達の金儲けのために
大量の中国人観光客を入国させて観光公害を撒き散らしてきたあげく、
いまや大勢のや日本人の生命が脅かされている。
319:卵の名無しさん
20/01/31 16:57:45 k8jciAtw.net
あへボケ三
320:卵の名無しさん
20/01/31 18:15:02 S7aQfrHj.net
アスペで人とのコミュニケーションが取れず、数学が得意であった。
某上位国立大学を卒業する。
医局入局する。上級医とのコミュニケーションも取れず、(旧研修医時代、徒弟制の)医局時代は厄介者扱いされ、地方の病院に飛ばされる。地方の病院でずっと独学で臨床に励む。尚、まともに教育されず独学のために系統だった臨床については…。
医局をいつの間にか辞めて、療養型病院に勤務する。
同期の医師が学会の要職、基幹病院の要職に就いている事をネットで見て屈折した感情を蓄積させる。
職場に私立卒の専門医を持っている上司が来る。
自分よりも学歴が下の癖に！と癇癪を起こす。
2ch(現5ch)で狂ったように鬱憤を晴らす。
321:卵の名無しさん
20/01/31 18:20:11 24QiZJ0Y.net
安倍は「対応が遅い」とかではない
●新型肺炎大流行になってるときに
どんどん観光客来て下さい！
●日本で人から人感染が起きてるときに
観光客減ると経済への影響が出るから！
と
患者をどんどん呼び込んだ
322:卵の名無しさん
20/02/01 14:02:48 hIisy8jC.net
180枚のくじを60人が一枚づつ引いて限りなく６人に近い人数に当たりが出るようにするには180枚のくじに何枚当たりを入れればよいでしょうか？
library(gmp)
n=180 ; r=60
chooseZ(180,60)
E <- function(x,n=180,r=60){
y=chooseZ(n,r)
z=0
for(i in 0:x) z = z + i*chooseZ(x,i)*chooseZ(n-x,r-i)
RE=z/y
cat(x,':',as.numeric(RE),'\n')
invisible(RE)
}
> for(i in 6:20) E(i)
6 : 2
7 : 2.333333
8 : 2.666667
9 : 3
10 : 3.333333
11 : 3.666667
12 : 4
13 : 4.333333
14 : 4.666667
15 : 5
16 : 5.333333
17 : 5.666667
18 : 6
19 : 6.333333
20 : 6.666667
323:卵の名無しさん
20/02/01 14:04:45.89 hIisy8jC.net
>>301
シミュレーションして確認。
sim <- function(x,n=180,r=60,k=1e4){
g <- function(x){
y=rep(1:0,c(x,n-x))
sum(sample(y,r))
}
mean(replicate(k,g(x)))
}
Sim=Vectorize(sim)
cbind(10:20,Sim(10:20))
> cbind(6:20,Sim(6:20))
[,1] [,2]
[1,] 6 2.0074
[2,] 7 2.3418
[3,] 8 2.6643
[4,]
324:9 3.0059 [5,] 10 3.3244 [6,] 11 3.6494 [7,] 12 3.9859 [8,] 13 4.3330 [9,] 14 4.6727 [10,] 15 5.0071 [11,] 16 5.3384 [12,] 17 5.6876 [13,] 18 5.9835 [14,] 19 6.3356 [15,] 20 6.7060 １８が近似している。
325:卵の名無しさん
20/02/01 14:41:44 eBpzYxtI.net
試験まであと1週間だな～受験生の健闘を祈る！
326:卵の名無しさん
20/02/01 14:44:30 eBpzYxtI.net
4月からはいよいよ医師としてのスタートだね！一緒に頑張ろう！
327:卵の名無しさん
20/02/01 15:17:51 hIisy8jC.net
H(n) = Σ[k=1,2,...,n] 1/k
とする。H(n)を既約分数で表したときの分子の整数をf(n)と表す。
f(1000)を求めよ。
5336291328229478504559104562404298040965247228038426009710134924845626888949710175750609790198503569140908873155046809837844
2172117885009464302344326566022502100278425632852081405544941210442510142672770294774712708917963967779610453224692426866468
8882815820719848971051107968732493191555293970175089315645199760857344730141832840117244122806490743077037366831700558002936
59235088589360235285852808160759574737836655413175508131522517
/
7128865274665093053166384155714272920668358861885893040452001991154324087581111499476444151913871586911717817019575256512980
2640676210092514658710043051310726862681432001966099748627459371883437050154344525237397452989631456749821282369562328237940
1106880926231770886197954079124775455804932647573782992335275179673524804246363805113703433121478174685087845348567802188807
5373249921995672056932029099390891687487672697950931603520000
おいド底辺、小学校でていたら既約分数くらい知っているだろう。計算があっているか確かめておいてくれ。
328:卵の名無しさん
20/02/01 15:19:38.51 hIisy8jC.net
バグ付きプログラムでの結果を数学板に投稿して撤回したから、デバッグして起きたかった。
329:卵の名無しさん
20/02/01 16:14:08.56 eBpzYxtI.net
試験は大学で学んだ事の集大成だよね。きちんと勉強してきた人は必ず合格する！
330:卵の名無しさん
20/02/01 16:18:33.36 ohxifYIE.net
外科系だとそうとも言えない、実例は山ほどある、
前の前の天皇の血管手術された先生は何処卒だっけ
331:卵の名無しさん
20/02/01 16:33:41 eBpzYxtI.net
この試験を乗り越えたらこれからの人生の自信になるだろう 医師の生活はたいへんだけれども充実したものになる！頑張ったら頑張った分戻ってくる
332:卵の名無しさん
20/02/01 19:38:15 hIisy8jC.net
【悲報】武漢からの帰国者192人、政府が用意したホテルが140室しかなく、2人1組の相部屋になることに怒号
library(gmp)
allocate.rooms <- function(m,n){ # m:rooms n:people
if(m==n
333:) return(factorialZ(m)) else if(m==1) return(1) else m*Recall(m,n-1) + m*Recall(m-1,n-1) } allocate.rooms(140,192)
334:卵の名無しさん
20/02/01 19:40:08 hIisy8jC.net
>>308
シリツ医大にも手先の器用な人はいるが、頭が器用ならシリツ医大なんぞには行かない。
まあ、多浪と女医は裏口入学の可能性が少ないらしい。
335:卵の名無しさん
20/02/01 19:41:44 hIisy8jC.net
東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。
336:卵の名無しさん
20/02/01 19:50:19 hIisy8jC.net
頭の不器用なシリツ医がそのうち事故を起こすだろうなと思っていたら
URLﾘﾝｸ(i.imgur.com)
これも表にでた分だからもっとひどいのは隠蔽されていると俺は思っている。
337:卵の名無しさん
20/02/02 07:06:52 X/HgwjDC.net
こういう投稿をみるとできると分かっていても解かずにいると気持ちわるい。
"
P,Q,Rを倫理(ﾏﾏ)命題とするとき、次の2つの倫理(ﾏﾏ)式は同値であるか、同値でないか。
真偽表を作成し、判定せよ
￢R⇒(Q ∨ P) と R ∨￢(Q ∨ P)
"
f <- function(p,q,r){
naraba<-function(x,y) !(x&!y)
naraba(!r,q|p)
}
g <- function(p,q,r) r | !(q|p)
sg = expand.grid(c(TRUE,FALSE),c(TRUE,FALSE),c(TRUE,FALSE))
colnames(sg)=c('P','Q','R')
d1=mapply(f,sg[,1],sg[,2],sg[,3])
d2=mapply(g,sg[,1],sg[,2],sg[,3])
cbind(sg,d1,d2)
P Q R d1 d2
1 TRUE TRUE TRUE TRUE TRUE
2 FALSE TRUE TRUE TRUE TRUE
3 TRUE FALSE TRUE TRUE TRUE
4 FALSE FALSE TRUE TRUE TRUE
5 TRUE TRUE FALSE TRUE FALSE
6 FALSE TRUE FALSE TRUE FALSE
7 TRUE FALSE FALSE TRUE FALSE
8 FALSE FALSE FALSE FALSE TRUE
338:卵の名無しさん
20/02/02 09:48:27.72 m4KFXXFo.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
339:卵の名無しさん
20/02/02 09:48:34.77 m4KFXXFo.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
340:卵の名無しさん
20/02/02 09:48:40.84 m4KFXXFo.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
341:卵の名無しさん
20/02/02 13:34:48.37 +5dNqMpE.net
315 名前：あぼ～ん[NGWord:あほシリツあほシリツ] 投稿日：あぼ～ん
342:卵の名無しさん
20/02/02 15:25:48 +5dNqMpE.net
３D世界の問題
一辺の長さが10mの立方体のプールの一つの角に監視員を置く.この監視員は水中は秒速1mで,プー ルの縁上は秒速 2m で移動するものとする.この監視員がプールのどこへでも到達しうるには，最短で何秒 必要か計算せよ.
343:卵の名無しさん
20/02/02 15:26:25 +5dNqMpE.net
di<-function(X,Y) sqrt(sum((X-Y)^2)) # X,Yの距離
Lx=Ly=Lz=10 # プールの大きさ
vs=2 # プールサイドの移動速度velocity along side
vw=1 # 水中速度 velocity in water
O=c(0,0,0) # 監視員の位置
S=c(x,y,z) # Swimmerの位置
sim <- function(S){
# O-S (Oから水没)
O=di(O,S)/vw
# O-X(X軸上から水没)
f=function(x) x/vs+di(c(x,0,0),S)/vw
OX=optimise(f,c(0,Lx))$obj
# O-Y(Y軸上から水没)
f=function(y) y/vs+di(c(0,y,0),S)/vw
OY=optimise(f,c(0,Ly))$obj
# O-Z(Z軸上から水没)
f=function(z) z/vs+di(c(0,0,z),S)/vw
OZ
344:=optimise(f,c(0,Lz))$obj
345:卵の名無しさん
20/02/02 15:26:41 +5dNqMpE.net
# O-X-Y(Xを全長走行してY軸上から水没)
f=function(y) Lx/vs+y/vs+di(c(Lx,y,0),S)/vw
OXY=optimise(f,c(0,Ly))$obj
# O-Y-X(Yを全長走行してX軸上から水没)
f=function(x) Ly/vs+x/vs+di(c(x,Ly,0),S)/vw
OYX=optimise(f,c(0,Lx))$obj
# O-Y-Z(Yを全長走行してZ軸上から水没)
f=function(z) Ly/vs+z/vs+di(c(0,Ly,z),S)/vw
OYZ=optimise(f,c(0,Lz))$obj
# O-Z-Y(Zを全長走行してY軸上から水没)
f=function(y) Lz/vs+y/vs+di(c(0,y,Lz),S)/vw
OZY=optimise(f,c(0,Ly))$obj
# O-Z-X(Zを全長走行してX軸上から水没)
f=function(x) Lz/vs+x/vs+di(c(x,0,Lz),S)/vw
OZX=optimise(f,c(0,Lx))$obj
# O-X-Z(Xを全長走行してZ軸上から水没)
f=function(z) Lx/vs+z/vs+di(c(Lx,0,z),S)/vw
OXZ=optimise(f,c(0,Lz))$obj
346:卵の名無しさん
20/02/02 15:26:47 +5dNqMpE.net
# O-X-Y-Z(X,Yを全長走行してZ軸上から水没)
f=function(z) Lx/vs+Ly/vs+di(c(Lx,Ly,z),S)/vw
OXYZ=optimize(f,c(0,Lz))$obj
# O-X-Z-Y(X,Zを全長走行してY軸上から水没)
f=function(y) Lx/vs+Lz/vs+di(c(Lx,y,Lz),S)/vw
OXZY=optimize(f,c(0,Ly))$obj
# O-Y-Z-X(Y,Zを全長走行してx軸上から水没)
f=function(x) Ly/vs+Lz/vs+di(c(x,Ly,Lz),S)/vw
OYZX=optimize(f,c(0,Lx))$obj
route=c('O','OX','OY','OZ',
'OXY','OYX','OYZ','OZY','OZX','OXZ',
'OXYZ','OXZY','OYZX')
t=c(O,OX,OY,OZ,OXY,OYX,OYZ,OZY,OZX,OXZ,OXYZ,OXZY,OYZX)
m=min(t)
cat(route[which(t==m)],':',m,'\n')
invisible(m)
}
347:卵の名無しさん
20/02/02 15:38:17 +5dNqMpE.net
options(digits=7)
(opt=optim(c(5,5,5),sim,method='L',control=list(fnscale=-1)))
opt[1]
opt[2]
> opt[1]
$par
[1] 7.691099 7.691099 7.691099
> opt[2]
$value
[1] 13.26518
348:卵の名無しさん
20/02/02 15:38:36 o+g6Fp1r.net
最後のラストスパートだね！悔いの残らないように頑張れよ！
349:卵の名無しさん
20/02/02 15:47:50 o+g6Fp1r.net
今まで勉強してきたことは無駄ではないことを味わって欲しい！まとまった勉強をすることはこの先あまり無いよ！
350:卵の名無しさん
20/02/02 15:50:22 m4KFXXFo.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
351:卵の名無しさん
20/02/02 15:50:32.69 m4KFXXFo.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
352:卵の名無しさん
20/02/02 15:50:41.06 m4KFXXFo.net
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
あほシリツあほシリツあほシリツあほシリツあほ
353:卵の名無しさん
20/02/02 16:42:29 OISskyPV.net
私立に落ちた人の恨みは怖いねｗ
354:卵の名無しさん
20/02/02 16:43:26 +5dNqMpE.net
人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。
これな！
不朽の名投稿
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８－９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<
355:卵の名無しさん
20/02/02 16:43:33 +5dNqMpE.net
東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。
356:卵の名無しさん
20/02/02 16:55:29 fRKGw+nR.net
>>319
これ東工大の入試問題を３次元に拡張した問題。
数理で解くなら偏微分が必要になる。
初等幾何学だけで解けたらネ申だな。
357:卵の名無しさん
20/02/02 16:56:04 fRKGw+nR.net
>>329
私立を受けただけでも汚らわしいよ。
358:卵の名無しさん
20/02/02 16:59:12 o+g6Fp1r.net
卒後20数年の無職だよーん ある掲示板では有名だよーん
359:卵の名無しさん
20/02/02 17:04:11 o+g6Fp1r.net
福井県越前市瓜生町
360:卵の名無しさん
20/02/02 17:12:41 fRKGw+nR.net
学習しないアホさからは藤田か聖魔と思われるな。
361:卵の名無しさん
20/02/02 17:16:14 o+g6Fp1r.net
よく国立医学部入れたね。偉いね。多浪だけどね。卒業までが限界だったんだね w w w
362:卵の名無しさん
20/02/02 17:19:56 fRKGw+nR.net
現役で国立大学進学が親との暗黙の了解。
一期は滑り止めに理?を選択。
363:卵の名無しさん
20/02/02 17:20:51 fRKGw+nR.net
>>323
コンピュータ使うなら、こういうのよりこっちの方が興味深いな。
数学の公式なら数値をいれれば一意に面積とか出てくるけど、疫学を含めた臨床は確率事象を扱っているから
当然ばらつきがでる。 SEIRモデルに従うシミュレーションプログラムを組んでみた。
サイコロの１の目の出る確率が1/6だからといって6回投げたら必ず１が１回でるわけではない。
１人の感染者を含む帰国者２０人を相部屋にいれたときのシミュレーションがこれ。
広がりが収束する場合と遷延する場合があるのがみてとれる。
URLﾘﾝｸ(i.imgur.com)
多重加算平均をとればSEIRモデルのグラフに収束していくのだろうが、現実的には何度も繰り返せるわけではない。