底辺私立医大を卒業した医者って頭悪いよね？ Part14

タイトルまんまです

すなわち、タイトルの正しさを検証するスレです。

※前スレ
底辺私立医大を卒業した頭悪いよね？
http://egg.2ch.net/test/read.cgi/hosp/1472273747/
底辺私立医大を卒業した医者って頭悪いよね？ Part2
http://egg.2ch.net/test/read.cgi/hosp/1520807375/
底辺私立医大を卒業した医者って頭悪いよね？ Part3
http://egg.2ch.net/test/read.cgi/hosp/1522967067/
底辺私立医大を卒業した医者って頭悪いよね？ Part4
http://egg.2ch.net/test/read.cgi/hosp/1526677794/
底辺私立医大を卒業した医者って頭悪いよね？ Part5
http://egg.2ch.net/test/read.cgi/hosp/1529666055/
底辺私立医大を卒業した医者って頭悪いよね？ Part6
http://egg.5ch.net/test/read.cgi/hosp/1531305684/
底辺私立医大を卒業した医者って頭悪いよね？ Part7
http://egg.2ch.net/test/read.cgi/hosp/1532472024/
底辺私立医大を卒業した医者って頭悪いよね？ Part8
http://egg.2ch.net/test/read.cgi/hosp/1533510399/
底辺私立医大を卒業した医者って頭悪いよね？ Part9
http://egg.2ch.net/test/read.cgi/hosp/1535547275/
底辺私立医大を卒業した医者って頭悪いよね？ Part10
https://egg.2ch.net/test/read.cgi/hosp/1537867846/
底辺私立医大を卒業した医者って頭悪いよね？ Part11
https://egg.5ch.net/test/read.cgi/hosp/1544320163/
底辺私立医大を卒業した医者って頭悪いよね？ Part12
https://egg.5ch.net/test/read.cgi/hosp/1550969602/
底辺私立医大を卒業した医者って頭悪いよね？ Part13
https://egg.5ch.net/test/read.cgi/hosp/1575242106/

【参考】
五浪で三流医大に入った古根高の末路↓

病的な虚言癖と妄想癖の精神科医　古根高の病名を診断するスレ
https://egg.5ch.net/test/read.cgi/hosp/1529634250/

古根高（サッポロファクトリーメンタルクリニック院長）のご尊顔
https://i.imgur.com/WMWcuRx.jpg

ド底辺シリツ医大卒ってこれだよな

底辺私立医大卒が国立大学医学部卒を妬むスレ [転載禁止](c)2ch.net
http://potato.5ch.net/test/read.cgi/hosp/1424846816/695

695 名前：卵の名無しさん[] 投稿日：2016/09/25(日) 22:04:25.53 ID:FzuOCpqV
匿名先生
精神科
▲
1
評価

私立出でも、素晴らしい医師がいるし、国立でも、どうしようもないのもいる。
そんなのは、あくまで個人の資質差と思ってました。

が、

今いる病院（医師の大部分が、いわゆる底辺私立）に来て、その考えは変わりました。

とある、気に食わない医師の白衣やサンダルを隠す、机の上にゴミを置く、など、小学生じみた嫌がらせをしている医師というのを、はじめて見ました。

国立出の変人も多々みてきましたが、こんな、小学生のイタズラみたいなことをする人は、いませんでした。

やはり、底辺私立は、それなりなんだと思いました。

2016年09月02日 00時21分09秒

人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。

これな！

不朽の名投稿
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８−９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<

東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。

これをコンピュータなしで解ける頭脳には感服する。


赤色のカメレオン３匹と、青色のカメレオン５匹と、黄色のカメレオンが８匹いる。

もし、異なる色の２匹のカメレオンが出会えば、２匹とも３番目の色に変わる。

例えば、青色と黄色が出会えば2匹とも赤色になる。

同じ色の２匹のカメレオンが出会っても色は変わらない。

出会いを重ねて、全てのカメレオンが同じ色になることは可能か？

可能ならば何色になるか？

糞河野と糞犬今すぐ死ねや
毎日糞犬吠えさせんなカス

>>6
こんな数になるとコンピュータも処理しきれないが理論値は正解を出してくれるんだよなぁ。

赤色のカメレオン９９匹と、青色のカメレオン１００匹と、黄色のカメレオンが１０１匹の場合はどうか？

赤色のカメレオン７７７匹と、青色のカメレオン８８８匹と、黄色のカメレオンが９９９匹の場合はどうか？

こういうのをサラッと答えられる頭脳が理学部や工学部にいくんだな。

考えようともしない奴がいくのが裏口でド底辺シリツ医大であるのは申すまでもない。

改竄したり隠蔽する安倍が原因だろ。
官僚はアクロバット擁護で安倍に貸しを作って役人が望む政策を実現させてる。
官僚のことだから安倍を脅すネタとして内密に名簿は保管しているだろう。

共産党の議員にまでミンシュガーをやってたのが安倍。
立憲の議員にも「共産党！」とやったのも安倍。
俺は国賓招聘した習近平を指差しして「共産党！！」とやってほしいんだが。

これをプログラムしたいなぁ。

赤0 青1 黄2の数字が書いてあるとして
Z/3で考えたら
0+1=2+2 赤と青が黄2個に変わる
1+2=0+0 青と黄が赤2個に変わる
2+0=1+1 黄と赤が青2個に変わる
で変化の前後で合計は変わらない

0*13+1*15+2*17=1
13+15+17=0
0*x=1は決して成立しない
逆に
0x+1y+2z=d(x+y+z)
となるdが存在するとき
1x+2y+0z=(d+1)(x+y+z)
2x+0y+1z=(d+2)(x+y+z)
であるので
一般性を失わずd=0としてよい
すなわち
0x+1y+2z=0
z=y+3n

n=0なら
0(x+y+z)+1*0+2*0=0
とできる
n>0なら
0(x+y)+1*0+2(3n)=0
0(x+y)+1(3n)+2*0=0
n=-m<0なら
0(x+z)+1(3m)+2*0=0
とできるので
そもそも
0n+1(3m)+2*0=0
m>0
としてよい
n=0なら終了なのでn>0とすると
0(n-1)+1(3m-1)+2*2=0
0(n+3)+1(3(m-1))+2*0=0
とできるので次第に減らして
0(n+3m)+1*0+2*0=0
にすることが可能

一体、いつから日本はこうなってしまったのでしょうか？

若者は正社員でも月15万円ほどの手取りでボーナス無し。
子どもの7人に1人が貧困状態。(先進国の中でも最悪な水準)

少子高齢化で人手不足になり、やっと給料が上がろうとしている所に移民を入れ、給与の上昇を抑えようとする政府。

災害で国民の命が奪われ続け、巨万の富を抱えているにも関わらず十分な防災対策をしない国。
デメリットしかないとわかっているにも関わらず実行された消費税の増税。
世紀の大発明であるiPS細胞にですら、支援(年10億円)を打ち切ろうとする政府・財務省。

バカがこの国を運営している。

少子化だし、シリツ大学廃止でいいんじゃね？


1 すらいむ ★ 2020/01/10(金) 16:47:33.70 ID:yEOurkSc9
名城大学で学生が教員刺す 学生の身柄確保

　愛知県警察本部によりますと10日午後４時すぎ、名古屋市天白区の名城大学で学生が教員を刺したと通報があったということです。
　教員は病院に運ばれましたが、詳しい容体はわかっていません。

　警察は学生の身柄を確保し詳しい状況を調べています。

NHK NEWS WEB　2020年1月10日 16時45分
https://www3.nhk.or.jp/news/html/20200110/k10012241491000.html

sim <- function(r,b,y,k=1e3,verbose=FALSE){
x=ca(r,b,y) # 色別の数からカメレオン配列形成
U=FALSE # 一色か？
i=0
while(!U){ # camelleon encounterを一色になるまで
# while(!U & k>0){ # camelleon encounterを最大k回繰り返す
if(verbose) cat(i,':',x,'\n') ; i=i+1 #　配列表示
x=ce(x)
U=is.uc(x)
# k=k-1
if(U) {
print(x)
break # 途中で一色ならloopを抜ける
}
}
return(U) # 一色になったかを返す
}

sim(1,3,4) # 小数剰余が統一色 1 0 1　ゆえ2番め
replicate(10,sim(1,4,7)) # 剰余が等しいから１ 1 1ゆえ3色統一色になれる
try(sim(1,2,3)) # 統一不能 無限ループ

直感を理詰めでどう証明するかの問題だな。

【問1】宇宙空間に2020個の星が、どの2つの星の距離も異なるように配置されている。
ただし、星は宇宙空間に固定されているとする（公転などは考えない）。
各星には天文学者がいて、自分と最も近い星だけを観測している。
このとき、「どの天文学者にも観測されない星が必ず存在する」は真か偽か。

【問2】2021個のときではどうか。

>>16
シミュレーション

rm(list=ls())

gsc <- function(N,M,cmplx=TRUE){ # M×Mの大きさの複素平面にN個の星を配置
f=function(x){
ab=sample(x,2,replace=TRUE)
ifelse(cmplx,ab[1]+1i*ab[2],ab[1])}
replicate(N,f(M))
}

dista <- function(x,y) abs(x-y) # 星の距離

snest <- function(i){ # i番目の星に最も近い星を選びそのindexを返す
dis=dista(cons[i],cons) # 距離0の自身を含めてすべての星との距離を出して
which(dis==min(dis[-i])) # ０を除いて最小の値のindexを返す
}
Snest <- Vectorize(snest) # ベクトル対応化

sim <- function(N,M,cmp){
cons=gsc(N,M,cmp) # constellation 星を配置
d=combn(N,2,function(z) dista(cons[z[1]],cons[z[2]])) # すべての2星間の距離
if(length(unique(d))==length(d)){ # 星の間の距離が異なることを確認
os=Snest(1:N) # os : observed star 各星から観察している星のindex配列を返す
all(1:N %in% os) # os がすべての星を含んでいるか
}
}

snest <- function(i,cons){ # 星の配列consのi番目の星に最も近い星を選びそのindexを返す
dis=abs(cons-cons[i]) # 距離0の自身を含めてすべての星との距離を出して
which(dis==min(dis[-i])) # ０を除いて最小の値のindexを返す
}

als <- function(N=4,M=16){ # allocate stars
cons=sample(M,N) # constellation 星を数直線1,2,..,M上に配置
d=combn(N,2,function(z) abs(cons[z[1]] - cons[z[2]])) # すべての2星間の距離
if(length(unique(d))==length(d)){ # 星の間の距離が異なることを確認
os=sapply(1:N,function(x) snest(x,cons)) # os:observed star 観察している星のindex配列を返す
as=all(1:N %in% os) # all stars : os がすべての星を含んでいるか
if(as) return(sort(cons)) # すべて含んでいれば星の配列をソートして返す
}
else return(0) # 同一距離がある、非観察星があれば0を返す
}

# すべて観察される配列があればその配列を返す、k回までになければ0を返す
sim <- function(N,M,k=NULL){ # k:上限回数、NULLなら制限なし
AS=0 # all star observed
if(is.null(k)){ # みつかるまで繰り返す
while(sum(AS)==0){
AS=als(N,M)
}
}else{
while(sum(AS)==0 & k > 0){ # 上限k回
AS=als(N,M)
k=k-1
}
}
AS
}

> sim(4,16)
[1] 8 9 13 16
> sim(4,16)
[1] 3 5 11 12
> sim(6,24)
[1] 1 2 12 16 21 24
> sim(8,64)
[1] 5 6 33 36 40 42 50 55
> sim(10,128)
[1] 2 6 16 17 41 44 53 60 112 125

星１６個の場合

https://i.imgur.com/PPfJ1pT.jpg

> sim(16,256)
[1] 7+198i 18+ 73i 36+ 16i 36+ 59i 43+212i 61+ 12i 80+102i 103+ 86i 105+127i
[10] 115+122i 140+230i 143+235i 194+ 39i 210+228i 219+ 85i 241+235i

# search nearest star from constellation
snest <- function(i,cons){ # 星の配列consのi番目の星に最も近い星を選びそのindexを返す
dis=abs(cons-cons[i]) # 距離0の自身を含めてすべての星との距離を出して
which(dis==min(dis[-i])) # ０を除いて最小の値のindexを返す
}

# allocate stars # 星間距離がすべて異なり、全て観察されている星の配列を探す
als <- function(N=4,M=16){ # 乱数でN個の星を1,2,...Mの数直線に配置し、適合しているか検証する
consRe=sample(M,N,replace=TRUE) # constellation 星を数直線1,2,..,M上に配置
consIm=sample(M,N,replace=TRUE)
cons=complex(N)
for(i in 1:N) cons[i]=consRe[i]+1i*consIm[i]
d=combn(N,2,function(z) abs(cons[z[1]] - cons[z[2]])) # すべての2星間の距離
if(length(unique(d))==length(d)){ # 星の間の距離が異なれば
os=sapply(1:N,function(x) snest(x,cons)) # os:observed star 観察している星のindex配列を返す
as=all(1:N %in% os) # all stars : os がすべての星を含んでいるか
if(as){
plot(cons,bty='l',pch=19,asp=1)
return(sort(cons)) # すべて含んでいれば星の配列をソートして返す
}
}
else return(0) # 同一距離がある、非観察星があれば0を返す
}

# すべて観察される配列があればその配列を返す、k回までになければ0を返す
sim <- function(N,M,k=NULL){ # k:上限回数、NULLなら制限なし
AS=0 # all star observed
if(is.null(k)){ # みつかるまで繰り返す
while(sum(AS)==0){
AS=als(N,M)
}
}else{
i=0
while(sum(AS)==0 & i < k){ # 上限k回
AS=als(N,M)
i=i+1
}
cat('i =',i,'\n')
}
AS
}
> sim(4,16)
[1] 7+ 1i 10+ 6i 15+12i 16+ 9i
> sim(6,24)
[1] 1+ 3i 1+24i 3+18i 4+ 8i 21+ 9i 24+19i
> sim(8,64)
[1] 15+41i 21+37i 28+41i 28+49i 34+26i 44+24i 51+16i 58+18i
> sim(10,128)
[1] 1+ 41i 9+ 4i 10+ 98i 26+ 2i 30+ 51i 45+ 88i 86+ 31i 112+122i 118+124i
[10] 123+ 9i
> sim(12,256)
[1] 10+ 16i 45+ 47i 56+249i 65+227i 93+178i 128+231i 143+161i 169+ 45i 174+ 43i
[10] 176+252i 201+ 4i 232+ 25i
> sim(16,256)
[1] 2+184i 46+ 53i 76+194i 99+ 27i 125+108i 167+134i 170+ 63i 184+ 26i 184+ 66i
[10] 184+223i 189+ 21i 216+ 44i 216+236i 222+166i 240+ 40i 249+199i

# 1 × 1の大きさの複素平面にN個の星を一様分布で発生させて (1)星間距離不同 (2)最短星観察 (3)全星被観察　を満たすものが、みつかるまで繰り返す

# search nearest star from constellation
snest <- function(i,cons){ # 星の配列consのi番目の星に最も近い星を選びそのindexを返す
dis=abs(cons-cons[i]) # 距離0の自身を含めてすべての星との距離を出して
which(dis==min(dis[-i])) # ０を除いて最小の値のindexを返す
}

# draw arrow from complex a to complex b 矢印描出
arrow <- function(a,b,...){
arrows(Re(a),Im(a),Re(b),Im(b),...)
}

# allocate stars # (1)星間距離不同 (2)最短星観察 (3)全星被観察の配列を探す
als <- function(N=16){ # 乱数でN個の星を複素平面に配置し、適合しているか検証する
consRe=runif(N) # 実部:一様分布[0,1]で乱数発生
consIm=runif(N) # 虚部:一様分布[0,1]で乱数発生
cons=complex(N)　# constellationの容れ物(複素数配列)
for(i in 1:N) cons[i]=consRe[i]+1i*consIm[i] # 複素数化
d=combn(N,2,function(z) abs(cons[z[1]] - cons[z[2]])) # すべての2星間の距離
if(length(unique(d))==length(d)){ # (1)星間距離不同ならば
os=sapply(1:N,function(x) snest(x,cons)) # 最短被観察星のindex配列をos(observed star)に代入
as=all(1:N %in% os) # (3)全星被観察 os がすべての星を含んでいるかをas(all stars)に代入
if(as){ # 全被観察の条件を満たせばグラフ化
plot(cons,bty='l',pch='★',asp=1) # 星描画
for(i in 1:N) arrow(cons[i],cons[os[i]],length=0.1,col='blue') # 観察矢印描画
return(sort(cons)) # 星の配列を実部でソートして返す
}
}
else return(0) # 同一距離があるか、非観察星があれば0を返す
}

# すべて観察される配列があればその配列を返す、k回までになければ0を返す
sim <- function(N,k=NULL){ # k:上限回数、NULLなら制限なし
AS=0 # all star observed 初期値
if(is.null(k)){ # k=NULLならみつかるまで繰り返す
while(sum(AS)==0){ # ASは0か配列なのでsumで数値にして比較(max(AS)でも可)
AS=als(N)
}
}else{
i=0
while(sum(AS)==0 & i < k){ # 上限k回
AS=als(N)
i=i+1
}
cat('i =',i,'\n') # 上限k回の何回目の試行で条件を満たしたか
}
AS
}

sim(16)
# https://i.imgur.com/WOwnFwI.jpg

sim(15) # Endless

すべての星が最短距離の星から観察されている配置が存在するとして、その星の集合をAとする。
Aからa->b b<-aのように相互監視されている星を除く集合をBとする。除かれた星の集合をCと呼ぶ。
Bの要素が1つの星であればその星はどこからも観察されていないからAの前提に反する。
Bの要素の数をｂとする。b人のうち一人でもCに属する星を観察しているとBの中に観察されない星が出現する。
それはAの前提に反するから、b人はすべてBに属する星を観察していることになる。
Bに属する星で星間の距離が最短の2個の星は相互監視していることになりBの前提に反する。
ゆえにBは空集合である。

The Disappearance of Bambi Woods

https://www.ama●zon.com/Disappearance-Bambi-Woods-anthology-Crime-ebook/dp/B07L159FNN/ref=sr_1_1?keywords=The+Disappearance+of+Bambi+Woods&qid=1578832618&s=books&sr=1-1

Wikipedia以上の情報はなくてつまらなかった。

1〜100000までの自然数の中から、「どの3個を選んでも等差数列を成さない2020個の数」が選べることを示せ

プログラムして100個作ってみた。

100個でも朝食の時間に計算終わってた。

> a
[1] 1 2 4 5 10 11 13 14 28 29 31 32 37 38 40 41 82 83 85 86 91 92 94 95 109
[26] 110 112 113 118 119 121 122 244 245 247 248 253 254 256 257 271 272 274 275 280 281 283 284 325 326
[51] 328 329 334 335 337 338 352 353 355 356 361 362 364 365 730 731 733 734 739 740 742 743 757 758 760
[76] 761 766 767 769 770 811 812 814 815 820 821 823 824 838 839 841 842 847 848 850 851 973 974 976 977
>

rm(list=ls())

p3c <- function(x){ # pick 3 numbers and check if arithmatic sequence
is.as3 <- function(x) diff(x)[1]==diff(x)[2]
is.as=function(y) is.as3(c(x[y[1]],x[y[2]],x[y[3]]))
n=length(x)
any(combn(n,3,is.as)) #
}

M=20
a=c(1,1,2)
i=1
AS=FALSE

while(length(a) < M){
a=append(a,rev(a)[1])
AS=p3c(a)
if(AS){
a=a[-length(a)]
}else{
a=append(a,i)
if(AS){
a=a[-length(a)]
}
}
i=i+1
}

a

治療ハンドブックの最新版が到着。
今年も通読予定。毎年通読していると読むのが楽しみ。

三次方程式の解の公式は面白くないんで三次方程式の解の公式に行く前に、
a=√(4-4t+3t^2)/(1-t+t^2) をゆっくり微分してもらえませんか？

Rの答
> f <- expression(sqrt(4-4*t+3*t^2)/(1-t+t^2))
> D(f , 't')
0.5 * ((3 * (2 * t) - 4) * (4 - 4 * t + 3 * t^2)^-0.5)/(1 - t +
t^2) - sqrt(4 - 4 * t + 3 * t^2) * (2 * t - 1)/(1 - t + t^2)^2
>

Possible derivation:
d/dt(sqrt(4 - 4 t + 3 t^2)/(1 - t + t^2))
Use the quotient rule, d/dt(u/v) = (v ( du)/( dt) - u ( dv)/( dt))/v^2, where u = sqrt(3 t^2 - 4 t + 4) and v = t^2 - t + 1:
= (-sqrt(4 - 4 t + 3 t^2) (d/dt(1 - t + t^2)) + (1 - t + t^2) (d/dt(sqrt(4 - 4 t + 3 t^2))))/(1 - t + t^2)^2
Differentiate the sum term by term and factor out constants:
= ((1 - t + t^2) (d/dt(sqrt(4 - 4 t + 3 t^2))) - d/dt(1) - d/dt(t) + d/dt(t^2) sqrt(4 - 4 t + 3 t^2))/(1 - t + t^2)^2
The derivative of 1 is zero:
= ((1 - t + t^2) (d/dt(sqrt(4 - 4 t + 3 t^2))) - sqrt(4 - 4 t + 3 t^2) (-(d/dt(t)) + d/dt(t^2) + 0))/(1 - t + t^2)^2
Simplify the expression:
= (-sqrt(4 - 4 t + 3 t^2) (-(d/dt(t)) + d/dt(t^2)) + (1 - t + t^2) (d/dt(sqrt(4 - 4 t + 3 t^2))))/(1 - t + t^2)^2
The derivative of t is 1:
= ((1 - t + t^2) (d/dt(sqrt(4 - 4 t + 3 t^2))) - sqrt(4 - 4 t + 3 t^2) (d/dt(t^2) - 1))/(1 - t + t^2)^2
Use the power rule, d/dt(t^n) = n t^(n - 1), where n = 2.
d/dt(t^2) = 2 t:
= ((1 - t + t^2) (d/dt(sqrt(4 - 4 t + 3 t^2))) - sqrt(4 - 4 t + 3 t^2) (-1 + 2 t))/(1 - t + t^2)^2
Using the chain rule, d/dt(sqrt(3 t^2 - 4 t + 4)) = ( dsqrt(u))/( du) ( du)/( dt), where u = 3 t^2 - 4 t + 4 and d/( du)(sqrt(u)) = 1/(2 sqrt(u)):
= (-(-1 + 2 t) sqrt(4 - 4 t + 3 t^2) + (1 - t + t^2) (d/dt(4 - 4 t + 3 t^2))/(2 sqrt(3 t^2 - 4 t + 4)))/(1 - t + t^2)^2
Differentiate the sum term by term and factor out constants:
= (-(-1 + 2 t) sqrt(4 - 4 t + 3 t^2) + d/dt(4) - 4 d/dt(t) + 3 d/dt(t^2) (1 - t + t^2)/(2 sqrt(4 - 4 t + 3 t^2)))/(1 - t + t^2)^2
The derivative of 4 is zero:
= (-(-1 + 2 t) sqrt(4 - 4 t + 3 t^2) + ((1 - t + t^2) (-4 (d/dt(t)) + 3 (d/dt(t^2)) + 0))/(2 sqrt(4 - 4 t + 3 t^2)))/(1 - t + t^2)^2
Simplify the expression:
= (-(-1 + 2 t) sqrt(4 - 4 t + 3 t^2) + ((1 - t + t^2) (-4 (d/dt(t)) + 3 (d/dt(t^2))))/(2 sqrt(4 - 4 t + 3 t^2)))/(1 - t + t^2)^2
The derivative of t is 1:
= (-(-1 + 2 t) sqrt(4 - 4 t + 3 t^2) + ((1 - t + t^2) (3 (d/dt(t^2)) - 1 4))/(2 sqrt(4 - 4 t + 3 t^2)))/(1 - t + t^2)^2
Use the power rule, d/dt(t^n) = n t^(n - 1), where n = 2.
d/dt(t^2) = 2 t:
= (-(-1 + 2 t) sqrt(4 - 4 t + 3 t^2) + ((1 - t + t^2) (-4 + 3 2 t))/(2 sqrt(4 - 4 t + 3 t^2)))/(1 - t + t^2)^2
Simplify the expression:
= (((-4 + 6 t) (1 - t + t^2))/(2 sqrt(4 - 4 t + 3 t^2)) - (-1 + 2 t) sqrt(4 - 4 t + 3 t^2))/(1 - t + t^2)^2
Simplify the expression:
Answer: |
| = (2 - 7 t + 6 t^2 - 3 t^3)/((1 - t + t^2)^2 sqrt(4 - 4 t + 3 t^2))

rm(list=ls())
graphics.off()

# draw segment of complex a to complex b
seg <- function(a,b,...){
segments(Re(a),Im(a),Re(b),Im(b),...)
}
# draw text y at complex x
pt <- function(x,y=NULL,...){
text(Re(x),Im(x), ifelse(is.null(y),'+',y), ...)
}
# connect every point(vertex)
coep <- function(ij,ve){
seg(ve[ij[1]],ve[ij[2]],col='gray',lty=3)
invisible(0)
}

# 一直線上にあるか？外積=0と同等
align <- function(xyz){ # y-y1 = m(x-x1), m=(y2-y1)/(x2-x1) をz1､z2が満たすか?
x=xyz[1];y=xyz[2];z=xyz[3]
(Im(z)-Im(x))*(Re(x)-Re(y)) == (Im(x)-Im(y))*(Re(z)-Re(x))
}

# ある点が三角形の内部にあるか？

oup <- function(A,B){ # 外積 outer prodct
Ax=A[1];Ay=A[2];Az=A[3]
Bx=B[1];By=B[2];Bz=B[3]
c(Ay*Bz-Az*By, Az*Bx-Ax*Bz, Ax*By-Ay*Bx)
}
# 二次元なので Z軸の値の正負のみが必要

opc <- function(a,b){ # outer product on complex plane
Re(a)*Im(b)-Im(a)*Re(b)
}

align <- function(ABC){ # ベクトルABとベクトルACの外積が0なら3点は直線上にある
opc(ABC[2]-ABC[1],ABC[3]-ABC[1])==0
}


in3 <- function(P,A,B,C){ # is P inside triangle ABC?
sum(opc(B-A,P-A)>0,opc(C-B,P-B)>0,opc(P-A,C-A)>0)%%3==0
#　右ねじ外積の方向がZ軸で全て正か全て負かで三角形の内部と確認
}

concave <- function(ve){
for(i in 2:5){
tri=ve[-i] # i番目をPとする
tm=matrix(tri[combn(4,3)],nrow=3) # triangl matrix for ABC (4候補)
for(j in 1:4){
inside=in3(ve[i],tm[1,j],tm[2,j],tm[3,j]) # 頂点がABCの内部にあるか？
if(inside) return(TRUE)
}
}
return(FALSE)
}
sim <- function(draw=TRUE){
ve=c(0,runif(4)+runif(4)*1i) # vertex 一様分布乱数で頂点を発生
v3=matrix(ve[combn(5,3)],nrow=3) # 3点の組み合わせ10通り
for(i in 1:10){
aln=align(v3[,i]) # いずれかの3点が直線上にあれば
if(aln) return(NA) # NAを返す
}
if(draw){
plot(ve,xlim=c(0,1),ylim=c(0,1),bty='l',pch=as.character(1:5), cex=0.75,
axes=TRUE, ann=FALSE,asp=1)
combn(5,2,function(x) coep(x,ve))}
if(concave(ve)) return(FALSE) # 凹ならばFALSE
else　return(TRUE) # 凸ならTRUEを返す
}
sim()

k=1e6
cvx5=replicate(k,sim(draw=F))
sum(is.na(cvx5))
> mean(cvx5)
[1] 0.270576

>>32
n=7が不適な例をパソコンで出してみた。

https://i.imgur.com/nKvKEpM.jpg

複素平面での７つの座標は

[1] 0.0000000+0.0000000i 0.8328849+0.3868275i 0.9219985+0.4125254i 0.8392956+0.9398254i
[5] 0.8249211+0.8054688i 0.7128166+0.5062512i 0.3850924+0.3886969i

graphics.off()
k=1e3
sim.n(6)
mean(replicate(k,sim.n(6,F,F)))
sim.n(7)


# 凹が見つかるまで探す(上限k回)
seek.concave <- function(n=8,k=1e4){
convex5=TRUE
i=0
while(convex5 & k >0){
convex5=sim.n(n,draw=F)
k=k-1
i=i+1
}
return(i)
}
seek.concave(n=8,k=1e5)
seek.concave(n=9,k=1e6)

これって名投稿だと感心する。
>>(quote)
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８−９割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<(unquote)

＞同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く

同業者の発言：
【ウハも】　開業医達の集い　8診　【粒も】 [無断転載禁止]&#169;2ch.net
>>
670 名前：卵の名無しさん[] 投稿日：2017/05/20(土) 11:15:40.12 ID:46exOAAP
学会で川崎の医者が発表してたら、「馬鹿が何偉そうにしゃべってる。」と思う自分が嫌になるが、
これだけは学生時代から続く反射なので止められない
<<

患者の発言：
【医療】医者は患者にコレを言われると、内心ものすごくムッとする★4 [無断転載禁止]&#169;2ch.net
>>
810 名前：名無しさん＠１周年[] 投稿日：2017/05/21(日) 00:11:22.04 ID:+h+2h2fq0
旧帝医卒の医者が（患者としては嫌だが）
多少偉そうにしているのはわからんでもないが

底辺私立に偉そうにされたら
そりゃ患者としてはむかつくだろｗｗ
<<

Last but not least,
it is not the kickass bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron..
ド底辺シリツ医大が悪いんじゃない、本人の頭が悪いんだ。

"次の性質(1)(2)(3)をすべて持つ四面体が存在することを証明せよ。(1)どの辺の長さも整数 (2)どの面の面積も整数 (3)体積が整数"
rm(list=ls())
library(numbers)

# while loop version for dec2n
dec2nw <- function(num, N, digit = 4){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

V=function(a,Print=TRUE){
a1=a[1];a2=a[2];a3=a[3];a4=a[4];a5=a[5];a6=a[6]
V2=
a1^2*a5^2*(a2^2+a3^2+a4^2+a6^2-a1^2-a5^2)+
a2^2*a6^2*(a1^2+a3^2+a4^2+a5^2-a2^2-a6^2)+
a3^2*a4^2*(a1^2+a2^2+a5^2+a6^2-a3^2-a4^2)-
a1^2*a2^2*a4^2-a2^2*a3^2*a5^2-a1^2*a3^2*a6^2-a4^2*a5^2*a6^2
if(V2<0) return(0)
else{
v=sqrt(V2)/12
if(Print) cat(v,'\n')
isNatural(v)
}
}
V(c(6,7,8,9,10,11))

S=function(a,Print=TRUE){
a1=a[1];a2=a[2];a3=a[3];a4=a[4];a5=a[5];a6=a[6]
s1=(a1+a2+a4)/2
s2=(a2+a3+a5)/2
s3=(a3+a6+a1)/2
s4=(a4+a5+a6)/2

S1=s1*(s1-a1)*(s1-a2)*(s1-a4)
S2=s2*(s2-a2)*(s2-a3)*(s2-a5)
S3=s3*(s3-a3)*(s3-a6)*(s3-a1)
S4=s4*(s4-a4)*(s4-a5)*(s4-a6)
S2=c(S1,S2,S3,S4)
if(any(S2<0)) return(0)
else{
if(Print) cat(sum(sqrt(S2)),'\n')
all(isNatural(sqrt(S2)))
}
}
S(c(6,7,8,9,10,11))
N=50 # 各辺の最大長
Nmin=sum(N^(5:0)) # 1,1,1,1,1,1から開始
dec2nw(Nmin,N=N)
Nmax=N^6-1
i=Nmin
SV=FALSE # 面積体積がすべて整数か？
while(!SV & i<Nmax){
i=Nmin
side=dec2nw(i,N)
SV=S(side,F) & V(side,F)
i=i+1
}
dec2nw(i-1,N)

療養病床勤務医・老健施設勤務医・回復期病棟勤務医が集うスレ☆　復活版 Part.4
https://egg.5ch.net/test/read.cgi/hosp/1550459754/992

992 名前：卵の名無しさん[] 投稿日：2019/07/20(土) 19:03:01.20 ID:WRt9tsMe
当直と外来の穴埋めは看護師も事務も断れないやつとか変な正義感持っているやつを狙って頼んでる
表面的には助かりますー、先生のおかげですーとかおだててるが内心はまたコイツにやらせとけばいいだろって下に見てる
感謝なんかしてない、そいうバカが変にフォローするせいで外来の医師数増やさなかったり、当直医を本気で確保しない（当直料安い）

キレイなキラキラした病院でしか働いたことないやつはこのことをわかってない。ブラック病院では標的にされて良いように使われるだけ。
大抵は俺だけ働いている、俺だけちゃんとした医療してるって頭おかしくなってやめていく

なんにもしない爺医や患者を適当にあしらう医師は人生経験でそのこと知ってんだよ
わかってて無能、バカのフリしてんだよ、気づけよ

>>41
存在証明だけなのだから
整数に拘らず有理数解があることを証明するだけでいいと数学板で教わった。無駄なプログラムを書いていたもんだ。

sim <- function(n=7){
th=runif(n,-pi ,pi)
p2d=function(th,r=1) r*(cos(th)+1i*sin(th))
f <- function(x){
x1=x[1];x2=x[2];x3=x[3]
v1=p2d(th[x1]);v2=p2d(th[x2]);v3=p2d(th[x3])
a=abs(v1-v2);b=abs(v2-v3);c=abs(v3-v1)
s=(a+b+c)/2
S=sqrt(s*(s-a)*(s-b)*(s-c))
2*S/(a+b+c) # S/s
}
sum(combn(n,3,f))
}

replicate(10,sim(7))

簡単にシミュレーションできると思ったがわりとデバッグに苦労した。配列にA[0]を許さないRに対応してようやく完成。
プログラミングはロジカル思考の訓練になるなぁ。
知恵遅れのド底辺シリツには無理だけど。


##
rm(list=ls())
"
円に内接する凸n角形(n≧4)を、1つの頂点からn-3本の対角線を引くことによって、n-2個の三角形に分割する。
このとき、各三角形の内接円の半径の和は分割の仕方に関係なく一定であることを示せ。
"
sim <- function(n,r=1){
th=c(0,sort(runif(n-1,0,2*pi))) # n角形の偏角θ(th)を[0,2π]で乱数発生させてソートして並べる
p2d=function(x) r*(cos(x)+1i*sin(x)) # 極形式を複素数に
c2r <- function(v1,v2,v3){ # complex number to radius of inscribed circle
a=abs(v1-v2);b=abs(v2-v3);c=abs(v3-v1) # 複素数間の距離を公式にいれる
s=(a+b+c)/2
S=sqrt(s*(s-a)*(s-b)*(s-c)) # 面積:ヘロンの公式
return(S/s) #　公式 r = 2*S/(a+b+c): 内接円の半径
}
ver=p2d(th) # vertex(複素数表示)
sor=numeric(n) # sum of raius of inscribed circle 頂点ver[i]からの対角線で分割した時
z2n <- function(x,m=n) ifelse(x%%m,x%%m,m) # ｎ系の剰余０のとき n を返す
for(i in 1:n){
for(j in 1:(n-2)){
sor[i]=sor[i] + c2r(ver[z2n(i)],ver[z2n(i+j)],ver[z2n(i+j+1)])
}
}
　return(sor)
}
sim(100)

巨大数を扱えるというふれこみのRmpfrって正確じゃないな。５０の階乗を計算させてみた

R with Rmpfr
> mpfr(factorial(50),1e5)
1 'mpfr' number of precision 100000 bits
30414093201713018969967457666435945132957882063457991132016803840

Haskell:
Prelude> product[1..50]
30414093201713378043612608166064768844377641568960512000000000000

Python:
import math
print(math.factorial(50))
30414093201713378043612608166064768844377641568960512000000000000

Wolfram:
https://www.wolframalpha.com/input/?i=50!
30414093201713378043612608166064768844377641568960512000000000000

と思っていたが、

library(Rmpfr)

one = mpfr(1, 5000)
factorial(one*50)

> one = mpfr(1, 5000)
> factorial(one*50)
1 'mpfr' number of precision 5000 bits
[1] 30414093201713378043612608166064768844377641568960512000000000000

藤林丈司

library(Rmpfr)
one = mpfr(1, 1000)

a2n <- function(x) as.numeric(paste(as.character(x),collapse=''))
nsplit <- function(x) as.numeric(unlist(strsplit(as.character(x),'')))

n=15*one
( x=sample(n) )
( y=a2n(x) )
( z=na.omit(nsplit(y)) )
sum(z)

災害が発生していたるところに重症被災者がいる。消防署から出動して救急センターに患者を搬送する
消防署から救急センターへの距離は100km　救急車のガソリンは20L、患者を乗せない状態では燃費は10km/L、患者を乗せての燃費は5km/Lである
患者を救える地域の面積はいくらになるか？


砂漠の基地Aからもうひとつの基地Bに向かって出かけた戦車がGPSの故障で進路を見失ってさまよった挙げ句にガス欠で止まってしまった。別の戦車で基地Aからこの戦車へ救助に向かい、燃料を分け与えて一緒に基地Bに行くことになった。しかし、戦車にはAB間をちょうど往復できるだけしか燃料は積めない。AB間の距離をRとして、救出可能な領域の面積をRを使って表しなさい。
ただし、戦車の燃費はいずれも同じものとする。

このシミュレーションプログラミングを考えているんだが、手がつかん。


密度が一定の球形の惑星がある。

この惑星の表面点Pでの重力を最大にするために体積と密度一定の条件で惑星の形を変形する。

点Pでの重力は球形のときの最大何倍になるか。

またこのときの惑星の形はどうなるか。

私立医から医師免許を剥奪せよ。
私立医学部を全て解体せよ。
本来なら到底医師になれず、
カネで医者になったニセ医者を迫害仲間外れにして糾弾しょう！
奴らの行動は患者殺しと国立医への逆差別、嫌がらせである。
私立医の患者殺しを許すな！

rm(list=ls())
graphics.off()

curve( sqrt(1-x^2), xlim=c(-1,1), ylim=c(-1,1), asp=1)
curve(-sqrt(1-x^2), xlim=c(-1,1), ylim=c(-1,1), add=T)


cir <- function(ce,r,...){
a=ce[1];b=ce[2] # (x-a)^2+(y-b)^2=r^2 ,y = b ± √(r^2-(x-a)^2)
curve(b+sqrt(r^2-(x-a)^2),xlim=c(a-r,a+r),ylim=c(b-r,b+r),asp=1,bty='l',...)
curve(b-sqrt(r^2-(x-a)^2),add=TRUE)
}

circle <- function(z,r,...){
if(is.complex(z)){ a=Re(z) ; b=Im(z)
}else{ a=z[1] ; b=z[2] }
x=seq(a-r,a+r,length=100)
y=b+sqrt(r^2-(x-a)^2)
plot(x,y,ylim=c(b-r,b+r),type='l',bty='l',...)
lines(x,2*b-y,...)
}
circle(c(1,3),5,lwd=2,col='maroon')


plot.circle <- function(x, y, r, ...){
theta <- seq(-pi, pi, length=100)
plot(x + r*cos(theta), y + r*sin(theta), type="l", asp=1, bty='l', ...)
}
plot.circle(1,2,3)

>>52
シミュレーションの準備に　まず、質点を一様分布配置しなくちゃならん。
https://i.imgur.com/O4o6aWt.jpg
じゃだめで
https://i.imgur.com/zN0GVeb.jpg
だな

改憲以前に安倍の議員資格に疑問符がついているんですが…
まず桜を見る会での公費による有権者買収は
公職選挙法違反なんじゃないんですか？

# 単位円内にN個の点を一様分布で発生させる
N=1e3
zt=0i
z=NULL
while(length(z) < N){
zt=sum(runif(2,-1,1)*c(1,1i))
if(abs(zt)<=1) z=append(z,zt)
}
plot(z,asp=1,bty='l',ann=F) ; points(p,pch='+',cex=2)

rm(list=ls())
library(gtools)

a=permutations(2,15,rep=T)

no21 <- function(x) { # 1の連続は無いときTRUEを返す No Two one
( y=rle(x) ) # run length endcoding で連続頻度リスト
( z=matrix(unlist(y),nrow=2,byrow=T) ) # 行列にして
( i1=which(z[2,]==1) ) # 1の頻度を配列化
all(z[1,i1]<2) # すべて２未満であればTRUE
}

no21i=which(apply(a,1,no21)) # no21i : 1が連続しないaのindex（行番号)
b=a[no21i,] # b:それに該当する行列

step15 <- function(x){ # 累積和が途中で15になればTRUEを返す step15
15 %in% cumsum(x)
}

s15i=which(apply(b,1,step15)) # s15i累積和が途中で15になるbのindex（行番号)
c=b[s15i,] # b:それに該当する行列

stp15 <- function(x){ # １５段になるまでのstepを返す
i=which(cumsum(x)==15)
if(length(i)!=0) return(x[1:i])
else return(NULL)
}

d=apply(c,1,stp15)
unique(d)

1歩で1段または2段のいずれかで階段を昇るとき、1歩で2段昇ることは連続しないものとする。15段の階段を昇る昇り方は何通りあるかを求めよ。

>>59

13段の方が興味をもたれるかな。

>>60
999段のとき　72428926648244374994035719791113496129988312942838688570843879389313216865724883592391343074367517004225503378684709437440　通りになったんだが、
数値計算精度の制約から誤答かもしれん。

ド底辺シリツ医でも数くらい数えられるだろうから、紙に書いて数えてくれ。

13段の場合だとこんな感じ。
[[1]]
[1] 1 2 1 2 1 2 1 2 1

[[2]]
[1] 1 2 1 2 1 2 2 2

[[3]]
[1] 1 2 1 2 2 1 2 2

....

[[27]]
[1] 2 2 2 2 2 1 2

[[28]]
[1] 2 2 2 2 2 2 1

>>61
不定長さの扱えるHaskell　だと

72428926648244409475025064434360344869036057631812487293633154973026123860267808630768185105127650714516006185052479679199

になった。

Pythonでも

N=1000
f = [0]*N
f[0]=1
f[1]=1
f[2]=2
for n in range(0,N-3):
f[n+3]=f[n+1]+f[n]
print (f[999-1])


72428926648244409475025064434360344869036057631812487293633154973026123860267808630768185105127650714516006185052479679199

とHaskellの結果と一致

ｍ３から配信されてきたのだが、最近の国試って馬鹿を医者にする問題だらけだな。


第112回医師国試問題◆
胸やけの誘因となりにくいのはどれか。

国試112-第112回医師国家試験問題解説書（医学評論社）より出題　学生正答率80.2％

A 過食
B 運動
C 肥満
D 高脂肪食
E 前屈姿勢

こういう問題を出して、裏口入学のシリツ医学生を排除すべきだと思うよなぁ。

1歩で1段または2段のいずれかで階段を昇るとき、1歩で2段昇ることは連続しないものとする。15段の階段を昇る昇り方は何通りあるかを求めよ。

こういう非定型な問題で場合分けして系統的にカウントできるかを問う問題は他人の結果と合致すると自分の計算に確信がもてる。

1歩で1段または2段のいずれかで階段を昇るとき、2歩で続けて昇ることは3回以上は連続しないものとする。（2回続けて2歩で昇るのは可）
15段の階段を昇る昇り方は何通りあるかを求めよ。

床に描かれた、ある程度大きな円の円周上に人物Aがいて、円の中心に人物Bがいる。
さて、次のようなゲームをする。
●AはBを捕まえるのが目的。しかし、円周上しか動けない。
●BはAに捕まらずに円周に到着するのが目的。その際、円周に向かってまっすぐ走ってもいいし、途中で向きを変えてもいい。また、無意味な時間稼ぎはしない(スタート地点に何時間も動かずにいるetc.)。
【問題】AがBを捕まえるためには、Aの速度がBの速度の何倍以上でなければならないか。ただし、AとBの速度は一定とし、お互い最善を尽くすものとする。
【注】BがAのいる場所とは正反対の方向に向かってまっすぐ行くと、Bの移動距離は半径r、Aの移動距離は半円周πrなので、AはBのπ倍の速さで行くと、Bが円周に到着した瞬間にBを捕まえられてAの勝ち、π倍未満ならばBの勝ちとなるが、果たしてBのこの行動は最善なのか？

circle <- function(z,r,...){
if(is.complex(z)){ a=Re(z) ; b=Im(z)
}else{ a=z[1] ; b=z[2] }
x=seq(a-r,a+r,length=100)
y=b+sqrt(r^2-(x-a)^2)
plot(x,y,ylim=c(b-r,b+r),asp=1,type='l',bty='l',...)
lines(x,2*b-y,...)}
# draw segment of complex a to complex b
seg <- function(a,b,...){
segments(Re(a),Im(a),Re(b),Im(b),col=2,...)}
# draw text y at complex x
pt <- function(x,y=NULL,...){
text(Re(x),Im(x), ifelse(is.null(y),'+',y), ...)}
p2d <- function(th,r=1) r*(cos(th)+1i*sin(th))

sim <- function(v=pi){
circle(0i,1) ; points(0+0i,pch='B',col=2) ; points(-1+0i,pch='A')
for(n in 1:20){
t=0.05
pt(p2d(pi+v*t*n),n)
pt(p2d(v*t*n,t*n),n,col=2)
}
}
sim(pi)

2018年2月5日月曜日
自殺者と変死体

最近では，遺書や目撃証言といった具体的な証拠がない限り，自殺とは認めず，変死体として処理するようになったと（上記記事）。
統計上の自殺者数というのは，現場のさじ加減で動くものです。統計には表れない「暗数」に思いを巡らせないといけません。


http://tmaita77.blogspot.com/2018/02/blog-post_5.html

https://2.bp.blogspot.com/-v86LnfvXaNE/Wne5R--T3hI/AAAAAAAAMSQ/a1mMae2gYEsPTUry7hr2_w2A-ipJlqQVgCLcBGAs/s1600/%25E5%25A4%2589%25E6%25AD%25BB%25E4%25BD%2593%25E2%2591%25A0.png

N=100
A = [0]*N
A[0] = 1
A[1] = 1
A[2] = 2
for n in range(3,N):
A[n] = (A[n-1] + A[n-3])
print (A[N-1]+A[N-3])

>>66
こういうふうに逃げれば
https://i.imgur.com/hZk4Buv.jpg

トランプだって在日米軍を引き上げたがり、日本の核武装を容認発言をしていたわけで、
日本が自立する最高のチャンスだった。
しかしそれすらも活かせず何の戦略もなくアメリカの犬であり続けたのが安倍自民政権。

それだけでなく安倍がネオリベグローバリストで保守的政策は何一つやらずに売国政策だけは猛スピードで次々に成立させている現実を見れば、
安倍の本性がなにかなんて馬鹿でもわかる問題だろうが。

明らかに安倍はDSのスパイなわけだが、100歩譲って仮にそうではなく仕方なく売国をやっていると仮定しても
日本自立の最大のチャンスすら何も活かせずにアメリカの犬をやり続け、
保守的政策は何一つやらずに売国政策だけは猛スピードで次々に成立させているわけで、
保守にとって何のメリットもない総理で、今すぐ倒閣し日本の破壊を止めさせる以外の選択肢はない。

水島は安倍売国に加担し続けている売国奴として我々保守は絶対に許すわけにはいかん。

>>415
自民党憲法改正草案(本文抜粋)
・第12条
憲法が国民に保障する自由及び権利は、国民の不断の努力により、保持されなければならない。
国民は、これを濫用してはならず、自由及び権利には責任及び義務が伴うことを自覚し、常に公益及び公の秩序に反してはならない。

・第13条
生命、自由及び幸福追求に対する国民の権利については、公益及び公の秩序に反しない限り、立法その他の国政の上で、最大限に尊重されなければならない。

・第14条-3
(現行)栄誉、勲章その他栄典の授与は、いかなる特典も伴わない。
↓
(改正)
栄誉、勲章その他栄典の授与は、現にこれを有し、又は将来これを受ける者の一代に限り、その効力を有する。

・第18条
(現行)何人も、いかなる奴隷的拘束も受けない。
↓
(改正)社会的又は経済的関係において身体を拘束されない。

・第29条
財産権の内容は、公益及び公にけ公の秩序に適合するように、法律で定める。

・第97条
(現行)この憲法が日本国民に保障する基本的人権は、人類の多年にわたる自由獲得の努力の成果であつて、
これらの権利は、過去幾多の試錬に堪へ、現在及び将来の国民に対し、侵すことのできない永久の権として信託されたものである。
↓
(改正)全文削除

AからHの8人はそれぞれ正直者か嘘つきであり、次のように証言している。このとき、嘘つきは誰か、全て答えよ。
A「8人の中に、少なくとも1人正直者がいる」
B「8人の中に、少なくとも2人正直者がいる」
C「8人の中に、少なくとも3人正直者がいる」
D「8人の中に、少なくとも4人正直者がいる」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「8人の中に、少なくとも3人嘘つきがいる」
H「8人の中に、少なくとも4人嘘つきがいる」

rm(list=ls())

dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

ah=sapply(0:255,function(x) dec2nw(x,2,8))
AH=t(ah)
AH=rbind(c(8,0),t(apply(AH[c(-1,-256),],1,table)),c(0,8))
AH

rm(list=ls())

dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

n=8

ah=sapply(0:(2^n-1),function(x) dec2nw(x,2,n))
AH=t(ah)
tf=rbind(c(n,0),t(apply(AH[c(-1,-(2^n-1),],1,table)),c(0,n))
tf

"
AからHの8人はそれぞれ正直者か嘘つきであり、次のように証言している。このとき、嘘つきは誰か、全て答えよ。
A「8人の中に、少なくとも1人正直者がいる」 tf[,2]
B「8人の中に、少なくとも2人正直者がいる」
C「8人の中に、少なくとも3人正直者がいる」
D「8人の中に、少なくとも4人正直者がいる」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「8人の中に、少なくとも3人嘘つきがいる」
H「8人の中に、少なくとも4人嘘つきがいる」
"
rm(list=ls())

dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

n=8

ah=sapply(0:(2^n-1),function(x) dec2nw(x,2,n))
AH=t(ah)
AH
tf=rbind(c(n,0),t(apply(AH[c(-1,-(2^n)),],1,table)),c(0,n))
tf

A「8人の中に、少なくとも1人正直者がいる」 tf[,2]>=1
B「8人の中に、少なくとも2人正直者がいる」 tf[,2]>=2
C「8人の中に、少なくとも3人正直者がいる」 tf[,2]>=3
D「8人の中に、少なくとも4人正直者がいる」 tf[,2]>=4
E「8人の中に、少なくとも1人嘘つきがいる」 tf[,1]>=1
F「8人の中に、少なくとも2人嘘つきがいる」 tf[,1]>=2
G「8人の中に、少なくとも3人嘘つきがいる」 tf[,1]>=3
H「8人の中に、少なくとも4人嘘つきがいる」 tf[,1]>=4

dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

n=8
te=sapply(0:(2^n-1),function(x) dec2nw(x,2,n)) # testimony
TE=t(te)

f <- function(x){
H=sum(x)
L=n-H
all(c(H>0,H>1,H>2,H>3,L>0,L>1,L>2,L>3)==x)
}

apply(TE,1,f)

通勤中に解決してしまった。


dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

n=8
te=sapply(0:(2^n-1),function(x) dec2nw(x,2,n)) # testimony
TE=t(te)

f <- function(x){
H=sum(x)
L=n-H
all(c(H>0,H>1,H>2,H>3,L>0,L>1,L>2,L>3)==x)
}

TE[apply(TE,1,f),]

> TE[apply(TE,1,f),]
[1] 1 1 1 1 1 1 0 0

こういう問題にすると
A「8人の中に、1人正直者がいる」
B「8人の中に、2人正直者がいる」
C「8人の中に、3人正直者がいる」
D「8人の中に、4人正直者がいる」
E「8人の中に、1人嘘つきがいる」
F「8人の中に、2人嘘つきがいる」
G「8人の中に、3人嘘つきがいる」
H「8人の中に、4人嘘つきがいる」
面白い。

こっちの方が面白い。

AからHの8人はそれぞれ正直者か嘘つきであり、次のように証言している。このとき、嘘つきは誰か、全ての組合せを答えよ。

A「8人の中に、1人嘘つきがいる」
B「8人の中に、2人嘘つきがいる」
C「8人の中に、3人嘘つきがいる」
D「8人の中に、4人嘘つきがいる」
E「8人の中に、5人嘘つきがいる」
F「8人の中に、6人嘘つきがいる」
G「8人の中に、7人嘘つきがいる」
H「8人の中に、8人嘘つきがいる」

頭の良さはこういう事では無いと思うのですが

>>82
んで>81の答は？

# 嘘つきは必ず嘘の証言をする設定に変更

dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

n=8
te=sapply(0:(2^n-1),function(x) dec2nw(x,2,n)) # testimony
TE=t(te)

f <- function(x){
H=sum(x)
L=n-H
all(c(
ifelse(x[1],H>=1,H<1),
ifelse(x[2],H>=2,H<2),
ifelse(x[3],H>=3,H<3),
ifelse(x[4],H>=4,H<4),

ifelse(x[5],L>=1,L<1),
ifelse(x[6],L>=2,L<2),
ifelse(x[7],L>=3,L<3),
ifelse(x[8],L>=4,L<4)

)==x)
}

TE[apply(TE,1,f),]

嘘つきを人数でなくて特定する証言を含めてみる。

AからHの8人はそれぞれ正直者か嘘つきであり、次のように証言している。証言でも嘘つきは必ず嘘をつく。
嘘つきは誰か、全ての組合せを答えよ。

A「8人の中に、正直者は3人いる」
B「8人の中に、少なくとも2人正直者がいる」
C「Bは嘘つきである」
D「Cは嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも嘘つきである」

dec2nw <- function(num, N=2, digit){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

n=8
te=sapply(0:(2^n-1),function(x) dec2nw(x,2,n)) # testimony
TE=t(te)

f <- function(x){
H=sum(x) ; cat(H,'\n')
L=n-H
all(c(
(x[1]==1 & H>=1) | x[1]==0 & H<1,
(x[2]==1 & H>=2) | x[2]==0 & H<2,
(x[3]==1 & H>=3) | x[3]==0 & H<3,
(x[4]==1 & H>=4) | x[4]==0 & H<4,

(x[5]==1 & L>=1) | x[5]==0 & L<1,
(x[6]==1 & L>=2) | x[6]==0 & L<2,
(x[7]==1 & L>=3) | x[7]==0 & L<3,
(x[8]==1 & L>=4) | x[8]==0 & L<4

))
}
TE[apply(TE,1,f),]

>>85
デバッグ終了（のつもり）
##
"AからHの8人はそれぞれ正直者か嘘つきであり、誰が正直者か嘘つきかはお互いに知っている。
次のように証言している。証言でも嘘つきは必ず嘘をつく。
嘘つきは誰かすべて答えよ。
A「8人の中に、正直者は3人いる」
B「8人の中に、少なくとも2人正直者がいる」
C「Bは嘘つきである」
D「Cは嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも嘘つきである」
"
TE=gtools::permutations(2,8,v=0:1,rep=T)
g <- function(x){
H=sum(x) # H:正直者の数
L=n-H # H:嘘つきの数
all(c( # 全て正しければTRUEを返す
(x[1]==1 & H==3) | (x[1]==0 & H!=3 ), # Aが正直者で証言が正しい | Aが嘘つきで証言が嘘
(x[2]==1 & H>=2) | (x[2]==0 & H <2 ), # Bが正直者で証言が正しい | Bが嘘つきで証言が嘘
(x[3]==1 & x[2]==0) | (x[3]==0 & x[2]==1),
(x[4]==1 & x[3]==0) | (x[4]==0 & x[3]==1),
(x[5]==1 & L>=1) | (x[5]==0 & L<1 ),
(x[6]==1 & L>=2) | (x[6]==0 & L<2 ),
(x[7]==1 & x[5]==0) | (x[7]==0 & x[5]==1),
(x[8]==1 & x[1]==0 & x[6]==0) | (x[8]==0 & (x[1]==1 | x[6]==1))
))
}
TE[apply(TE,1,g),]

どこかに、問題がないかと探っていたら、中学入試の問題があった

"Ａ、Ｂ、Ｃ、Ｄ、Ｅの５人のうち２人は常に本当のことを言う正直者です。\nあとの３人は嘘つきですが、その発言内容は本当のときもあります。\n彼らに誰が嘘つきか尋ねたところ次のように答えました。\n\nＡ　「ＢとＥは嘘つきではない」\nＢ　「Ｃは嘘つきだ」\nＣ　「Ｄは嘘つきだ」\nＤ　「Ｅは嘘つきだ」\nＥ　「ＢとＣは嘘つきだ」\nさて、正直者は誰と誰でしょうか？"
> TE=gtools::permutations(2,5,0:1,rep=T)
> f <- function(x){
+ if(sum(x)!=2){ return(FALSE)
+ }else{
+ all(c(
+ (x[1]==1&x[2]==1&x[5]==1) | x[1]==0,
+ (x[2]==1&x[3]==0) | x[2]==0,
+ (x[3]==1&x[4]==0) | x[3]==0,
+ (x[4]==1&x[5]==0) | x[4]==0,
+ (x[5]==1&x[2]==0&x[3]==0) | x[5]==0))
+ }
+ }
> LETTERS[1:5][TE[apply(TE,1,f),]==1]
[1] "B" "D"

うまく動作している！

Ａ、Ｂ、Ｃ、Ｄ、Ｅの５人のうち２人は常に本当のことを言う正直者です。あとの３人は嘘つきですが、その発言内容は本当のときもあります。

彼らに誰が嘘つきか尋ねたところ次のように答えました。

Ａ　「ＢとＥは嘘つきではない」

Ｂ　「Ｃは嘘つきだ」

Ｃ　「Ｄは嘘つきだ」

Ｄ　「Ｅは嘘つきだ」

Ｅ　「ＢとＣは嘘つきだ」

さて、正直者は誰と誰でしょうか？

>>82
おい、ド底辺。

これに答えてみ！小学生向けの問題だけど。

あるド底辺シリツ医大の新入生が101人とする。
101人は学力考査で入学した正規学生と金で入った裏口学生のどちらかである。
正規学生は常に本当のことを語り、裏口学生は虚しか言わない。
全ての学生は誰が裏口学生かを知っている。
学生全員に裏口学生は何人いるかを尋ねたところ

学生A1「1人以上の裏口学生がいる」
学生A2「2人以上の裏口学生がいる」
学生A3「3人以上の裏口学生がいる」
...
学生A99「99人以上の裏口学生がいる」
学生A100「100人以上の裏口学生がいる」
学生A101「全員が裏口学生である」
との証言を得た。

以上の情報からこのド底辺シリツ医大の裏口学生は誰かを述べよ。

>>87
更に、嘘をつくかどうかわからないファジーな人物も参加させて問題にしてみた。

AからHの8人はそれぞれ正直者か嘘つきであり、誰が正直者か嘘つきかはお互いに知っている。
A,B,C,D,Eは嘘つきなら必ず嘘をつくが、F,G,Hは嘘つきでも正しいことを言う場合がある。
次の証言から確実に正直者と断定できるものを全て挙げよ。

A「8人の中に、正直者は3人いる」
B「8人の中に、少なくとも2人正直者がいる」
C「Bは嘘つきである」
D「Cは嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも嘘つきである」

TE=gtools::permutations(2,n,v=0:1,rep=T)
colnames(TE)=LETTERS[1:n]
k <- function(x){
H=sum(x) # H:正直者の数
L=n-H # H:嘘つきの数
all(c( # 全て正しければTRUEを返す
(x[1]==1 & H==3) | (x[1]==0 & H!=3 ), # Aが正直者で証言が正しい | Aが嘘つきで証言が嘘
(x[2]==1 & H>=2) | (x[2]==0 & H <2 ), # Bが正直者で証言が正しい | Bが嘘つきで証言が嘘
(x[3]==1 & x[2]==0) | (x[3]==0 & x[2]==1),
(x[4]==1 & x[3]==0) | (x[4]==0 & x[3]==1),
(x[5]==1 & L>=1) | (x[5]==0 & L<1 ),
(x[6]==1 & L>=2) | (x[6]==0 ),
(x[7]==1 & x[5]==0) | (x[7]==0 ),
(x[8]==1 & x[1]==0 & x[6]==0) | (x[8]==0 ) ) )
}
TE[apply(TE,1,k),]

↑荒らし工作員
底辺私立医の実態バラされると死人がでるからな

>>92
どの問題にも解答の出せないアホ乙！

>>91
数学板に書いたら、プログラムで解いたと俺のとは違う解答が出てきたので困惑していたが

論証で解答してくれた人がいて、俺の解答と一致していたので俺のプログラムの方が正しいと確信できてほっとした。

このスレにはド底辺頭脳しか、いないようだな。

>>93
ズボシで必死だな
死ねキチガイ

>>95
どれか答えてみ！
ド底辺頭脳クン

お〜い、ド底辺。

これに答えてみ！小学生向けの問題だけど。

あるド底辺シリツ医大の新入生が101人とする。
101人は学力考査で入学した正規学生と金で入った裏口学生のどちらかである。
正規学生は常に本当のことを語り、裏口学生は虚しか言わない。
全ての学生は誰が裏口学生かを知っている。
学生全員に裏口学生は何人いるかを尋ねたところ

学生A1「1人以上の裏口学生がいる」
学生A2「2人以上の裏口学生がいる」
学生A3「3人以上の裏口学生がいる」
...
学生A99「99人以上の裏口学生がいる」
学生A100「100人以上の裏口学生がいる」
学生A101「全員が裏口学生である」
との証言を得た。

以上の情報からこのド底辺シリツ医大の裏口学生は誰かを述べよ。

国立卒でも>97の裏口学生が誰か答えられないなら裏口学生レベルの頭脳だね。

フェルマーの最終定理の反例が見つかる

これをみつけられるのは、それはそれで凄いな。

https://livedoor.blogimg.jp/math_stady/imgs/b/2/b28b082c-s.jpg

↑あほ

外来やっていても単語しか会話できない知恵遅れっているよな。

おい、ド底辺。

これに答えてみ！小学生向けの問題だけど。

あるド底辺シリツ医大の新入生が101人とする。
101人は学力考査で入学した正規学生と金で入った裏口学生のどちらかである。
正規学生は常に本当のことを語り、裏口学生は虚しか言わない。
全ての学生は誰が裏口学生かを知っている。
学生全員に裏口学生は何人いるかを尋ねたところ

学生A1「1人以上の裏口学生がいる」
学生A2「2人以上の裏口学生がいる」
学生A3「3人以上の裏口学生がいる」
...
学生A99「99人以上の裏口学生がいる」
学生A100「100人以上の裏口学生がいる」
学生A101「全員が裏口学生である」
との証言を得た。

以上の情報からこのド底辺シリツ医大の裏口学生は誰かを述べよ。

死ね

どれか答えてみ！
ド底辺頭脳クン

へー

ノータリン

>>99
これを探すスクリプトを書いてみた。
options(scipen = 10)
n=1000
cmb=gtools::combinations(n,2)
cmb=cmb[cmb[,1]<cmb[,2],]

dsplit <- function(x){ # digit split for natural number
y=unlist(strsplit(as.character(x),''))
if('.' %in% y) return(NULL)
else as.numeric(unlist(y))
}
a2d <- function(x) as.numeric(paste(as.character(x),collapse=''))

f <- function(xy){
x=xy[1] ; y=xy[2]
x3=x^3 ; y3=y^3
a=dsplit(x3) ; b=dsplit(y3)
y3 ; b
dl=length(b)-length(a)
a0=append(rep(0,dl),a)
c=a0+b
c2=a2d(c)
c3=c2^(1/3)
d=round(c3,10)
if(numbers::isNatural(d)) return(d)
else return(0)
}

idx=which(apply(cmb,1,f)>0)
cmb[idx,]

>>106
３桁までの数字では、画像にある１組だけだな。

> cmb[idx,]
[1] 70 212

AB <- function(abc){
a=abc[1];b=abc[2];c=abc[3]
y=b+b-a
if(y>=c) return(abc)
else return(c(b,y,c))
}

BC <- function(abc){
a=abc[1];b=abc[2];c=abc[3]
c(a,c,2*c-b)
}

BA <- function(abc){
a=abc[1];b=abc[2];c=abc[3]
c(a-(b-a),a,c)
}

CB <- function(abc){
a=abc[1];b=abc[2];c=abc[3]
y=b-(c-b)
if(y<=a) return(abc)
else return(c(a,y,b))
}

J=c(AB,BC,BA,CB)

jump <- function(x){
j=J[[sample(4,1)]](x)
while(all(j==x)){
j=J[[sample(4,1)]](x)
}
return(j)
}

lim=1e5
sim <- function(x=c(1,2,3)){
cat(0, ': ',x,'\n')
j=jump(x)
i=1
cat(i, ': ',j,'\n')
while(any(j!=x) & -lim<min(j) & max(j)<lim){
j=jump(j)
i=i+1
cat(i, ': ',j,'\n')
}
}

sim()

このジャンプで復帰するね。
> sim()
0 : 1 2 3
1 : 1 3 4
2 : 1 4 5
3 : -2 1 5
4 : 1 4 5
5 : -2 1 5
6 : 1 4 5
7 : 1 5 6
8 : 1 4 5
9 : 1 3 4
10 : -1 1 4
11 : 1 3 4
12 : 1 2 3

https://www.nippon-foundation.or.jp/who/news/pr/2019/20191130-38555.html

日本財団「18歳意識調査」第20回 テーマ：「国や社会に対する意識」（9カ国調査）国や社会に対する意識 各国1,000人に聞く
日本の若者の数字の低さ際立つ

https://www.nippon-foundation.or.jp/app/uploads/2019/11/new_pr_20191130_01.png

https://www.nippon-foundation.or.jp/app/uploads/2019/11/new_pr_20191130_02.png

英語に変更します。

> Sys.setenv(LANGUAGE="en")

> fisher.test(c(1,2,3,4))

Error in fisher.test(c(1, 2, 3, 4)) :

if 'x' is not a matrix, 'y' must be given

> Sys.getenv("LANGUAGE")

[1] "en"

>
Sys.setenvでLANGUAGEをenに設定すると、メッセージが英語になります。
エラーメッセージの表示も英語になっています。
先ほどは空だったSys.getenv("LANGUAGE")の結果で"en"が返っています。

日本語に戻します。

> Sys.setenv(LANGUAGE="ja")

> fisher.test(c(1,2,3,4))

以下にエラー fisher.test(c(1, 2, 3, 4)) :

'x' が行列でなければ、'y' を与えなければなりません

> Sys.getenv("LANGUAGE")

[1] "ja"

正解が一つの方が問題としては美しいな。

"AからHの8人はそれぞれ正直者か嘘つきであり、誰が正直者か嘘つきかはお互いに知っている。
A,B,C,D,Eは嘘つきなら必ず嘘をつくが、F,G,Hは嘘つきでも正しいことを言う場合がある。
次の証言から確実に正直者と断定できるものを全て挙げよ。

A「8人の中に、正直者は3人いる」
B「8人の中に、少なくとも2人正直者がいる」
C「Bは嘘つきである」
D「CもFも嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも嘘つきである」
"
n=8
TE=gtools::permutations(2,n,v=0:1,rep=T)
colnames(TE)=LETTERS[1:n]
k <- function(x){
H=sum(x) # H:正直者の数
L=n-H # H:嘘つきの数
all(c( # 全て正しければTRUEを返す
(x[1]==1 & H==3) | (x[1]==0 & H!=3 ) , # Aが正直者で証言が正しい | Aが嘘つきで証言が嘘
(x[2]==1 & H>=2) | (x[2]==0 & H <2 ) , # Bが正直者で証言が正しい | Bが嘘つきで証言が嘘
(x[3]==1 & x[2]==0) | (x[3]==0 & x[2]==1) ,
(x[4]==1 & x[3]==0 & x[6]==0) | (x[4]==0 & !(x[3]==0 & x[6]==0)) ,
(x[5]==1 & L>=1) | (x[5]==0 & L<1 ) ,
(x[6]==1 & L>=2) | (x[6]==0 ),
(x[7]==1 & x[5]==0) | (x[7]==0 ),
(x[8]==1 & x[1]==0 & x[6]==0) | (x[8]==0 ) ) )
}
TE[apply(TE,1,k),]

>>108

> ans
[,1] [,2] [,3]
[1,] 1 3 4
[2,] -1 1 4
[3,] -1 4 7
[4,] -1 1 4
[5,] -3 -1 4
[6,] -3 4 9
[7,] -10 -3 9
[8,] -3 4 9
[9,] -3 9 14
[10,] -3 4 9
[11,] -3 -1 4
[12,] -5 -3 4
[13,] -5 4 11
[14,] -5 11 18
[15,] -5 18 25
[16,] -28 -5 25
[17,] -5 18 25
[18,] -5 11 18
[19,] -5 4 11
[20,] -5 11 18
[21,] -5 18 25
[22,] -5 11 18
[23,] -5 4 11
[24,] -14 -5 11
[25,] -14 11 27
[26,] -14 -5 11
[27,] -5 4 11

[28,] -5 -3 4
[29,] -3 -1 4
[30,] -3 4 9
[31,] -3 9 14
[32,] -3 14 19
[33,] -20 -3 19
[34,] -3 14 19
[35,] -3 9 14
[36,] -3 4 9
[37,] -3 -1 4
[38,] -5 -3 4
[39,] -5 4 11
[40,] -5 -3 4
[41,] -3 -1 4
[42,] -3 4 9
[43,] -3 -1 4
[44,] -5 -3 4
[45,] -7 -5 4
[46,] -9 -7 4
[47,] -7 -5 4
[48,] -7 4 13
[49,] -18 -7 13
[50,] -18 13 33
[51,] -18 -7 13
[52,] -7 4 13
[53,] -18 -7 13
[54,] -7 4 13
[55,] -7 13 22

[56,] -27 -7 22
[57,] -7 13 22
[58,] -7 22 31
[59,] -7 31 40
[60,] -7 22 31
[61,] -7 13 22
[62,] -7 22 31
[63,] -36 -7 31
[64,] -7 22 31
[65,] -36 -7 31
[66,] -7 22 31
[67,] -36 -7 31
[68,] -7 22 31
[69,] -7 13 22
[70,] -27 -7 22
[71,] -7 13 22
[72,] -7 4 13
[73,] -18 -7 13
[74,] -18 13 33
[75,] -18 33 53
[76,] -18 53 73
[77,] -89 -18 73
[78,] -89 73 164
[79,] -251 -89 164
[80,] -251 164 417
[81,] -251 -89 164
[82,] -89 73 164
[83,] -89 -18 73
[84,] -89 73 164
[85,] -89 164 255

[86,] -89 73 164
[87,] -89 -18 73
[88,] -160 -89 73
[89,] -160 73 235
[90,] -160 -89 73
[91,] -89 -18 73
[92,] -18 53 73
[93,] -18 33 53
[94,] -18 53 73
[95,] -18 33 53
[96,] -18 53 73
[97,] -89 -18 73
[98,] -18 53 73
[99,] -18 73 93
[100,] -18 53 73
[101,] -18 33 53
[102,] -18 13 33
[103,] -18 -7 13
[104,] -7 4 13
[105,] -7 -5 4
[106,] -5 -3 4
[107,] -3 -1 4
[108,] -1 1 4
[109,] 1 3 4
[110,] 1 2 3

ようやく、復帰！

シリツ医はやくざと同じ
大半がチョン
税金還付で実質脱税

論理的思考が試される問題

AからHの8人はそれぞれ正直者か嘘つきであり、誰が正直者か嘘つきかはお互いに知っている。
A,B,C,D,Eは嘘つきなら必ず嘘をつくが、F,G,Hは嘘つきでも正しいことを言う場合がある。
この中から確実に正直者と言えるのは誰か？

A「8人の中に、正直者は3人いる」
B「Hは嘘つきである」
C「Bは嘘つきである」
D「CもFも嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも正直者である」

類題

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も正直者である」

2013年の発言

http://imgur.com/MGbp1xx.jpg
↓
http://imgur.com/LU42EXK.jpg

require(stats); require(graphics)
summary(anscombe)

##-- now some "magic" to do the 4 regressions in a loop:
ff <- y ~ x
mods <- setNames(as.list(1:4), paste0("lm", 1:4))
for(i in 1:4) {
ff[2:3] <- lapply(paste0(c("y","x"), i), as.name)
## or ff[[2]] <- as.name(paste0("y", i))
## ff[[3]] <- as.name(paste0("x", i))
mods[[i]] <- lmi <- lm(ff, data = anscombe)
print(anova(lmi))
}

## See how close they are (numerically!)
sapply(mods, coef)
lapply(mods, function(fm) coef(summary(fm)))

## Now, do what you should have done in the first place: PLOTS
op <- par(mfrow = c(2, 2), mar = 0.1+c(4,4,1,1), oma = c(0, 0, 2, 0))
for(i in 1:4) {
ff[2:3] <- lapply(paste0(c("y","x"), i), as.name)
plot(ff, data = anscombe, col = "red", pch = 21, bg = "orange", cex = 1.2,
xlim = c(3, 19), ylim = c(3, 13))
abline(mods[[i]], col = "blue")
}
mtext("Anscombe's 4 Regression data sets", outer = TRUE, cex = 1.5)
par(op)

>>107
メモリー消費を抑えるには

N=1e3
Fer=NULL
for(i in 1:N){
for(j in i:N){
if(f(c(i,j))>0) Fer=rbind(Fer,c(i,j))
}
}
Fer

だが、処理速度が落ちるな。

イランイラク戦争のあとイラクに傀儡政権としてフセイン政権を作ったのはCIA
同じように敗戦後日本に自民党という傀儡政権をCIAが作った

CIAはA級戦犯の児玉誉士夫と岸信介に目をつけた
自民党を作ったのは緒方竹虎自由党と鳩山一郎民主党
その緒方の後継者となり自民党総裁となり60年日米安保を改正した

緒方→CIA工作員、岸、鳩山を総裁につけた児玉誉士夫→CIA工作員、
岸信介→CIA工作員、佐藤栄作→CIA工作員、読売オーナー正力松太郎→CIA工作員

かくして自民党は1958年からCIAからの7000万〜1億円(現在の10億円)の資金を十数回受け取り政界工作をして55年体制を維持した
50年経過しアメリカ情報公開法により2007年CIAからの資金提供やコードネームなどが明らかとなった

現代でもアメリカの意向に逆らえば政界から失脚する
田中角栄や小沢一郎がいい例

>>124
４桁までを探索させた結果は

[,1] [,2]
[1,] 70 212
[2,] 700 2120

>>124
処理に時間がかかるから、途中経過を表示するようにした方がいいな。

N=1e3
Fer=NULL
for(i in 1:N){
for(j in i:N){
if(f(c(i,j))>0){
Fer=rbind(Fer,c(i,j))
print(Fer)
}
}
}
Fer

AからEの5人はそれぞれ正直者か嘘つきのどちらかであり、誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきの可能性があるのは誰か？

A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」

馬鹿は、死ななきゃ治らない　の　待遇を　馬鹿が治れば死んでいる　と答えるのがド底辺シリツ医だから

>128に正答できないだろうな。　ド底辺頭脳の国立卒も似たようなものだが。

馬鹿は、死ななきゃ治らない　の　対偶を　馬鹿が治れば死んでいる　と答えるのがド底辺シリツ医だから

>128に正答できないだろうな。　ド底辺頭脳の国立卒も似たようなものだが。

高校数学おもろいか？w
次元の違う高等数学ゴマンとあるぞ
あほやろおまえ

可能性より特定できた方が問題らしいな。

"
AからEの5人の中に正直者と嘘つきの双方がいるとわかっている。
誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきを全員列挙せよ。

A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」
"

>>131
じゃあ、大学教養レベルの問題に答えてみ！

これはシリツ医や底辺国立頭脳には無理だろうな。

ド底辺シリツ裏口調査団が100人を順次調査した。
裏口判明人数をそのまま公表はヤバすぎる結果であったため、
連続して裏口がみつかった最大数は4人であったとだけ公表した。
公表結果が正しいとして裏口入学人数の期待値、最頻値、及び95%信頼区間を述べよ。

臨床医学は確率事象を扱う。
少数データから断定するのはアホ。
信頼区間を設定して議論するのが臨床医の正しい姿。

これ答えてみ。高校数学の知識では信頼区間算出は無理。
ド底辺シリツ医にも当然、無理。

ド底辺シリツ医大にある年の裏口受験票が保存されていることが判明した。
枚数は１０００枚以下であることは判明しているが、その数を推定したい。
裏口入学調査員が無作為に１０枚選んで番号を記録して元に戻した。
別の裏口入学調査員が無作為に２０枚選んで番号を記録した。
二人の調査員の記録した番号を照合すると３枚の番号が一致していた。
この情報からこの年の裏口受験証の数の最頻値、期待値とその９５％信頼区間を求めよ。

>>66

Aと円の中心と対称となる点を、目指して常に方向転換するという方針で、定速で作図してみた。

偏角を使った方が楽だった。数値積分を組み込むは煩瑣。

https://i.imgur.com/DedHINQ.jpg

Bの曲線の接線と円周の交点がその時点でのAと対称になっている（はず）。

図はAの速度がBの１．２倍のとき。

これ生データがほしいなぁ

https://link.springer.com/article/10.1007/s10508-019-1464-z?shared-article-renderer

【医学】胸の大きな女性は小さな女性より風邪の治りが2倍以上も遅い！ポーランドの最新研究で明らかに
https://asahi.5ch.net/test/read.cgi/newsplus/1579432638/

標準回帰係数が
respiratory infections (standardized β?=?0.20, p?=?.011), their duration (standardized β?=?0.18, p?=?.027),

>>131
別に実用上必要な数値がだせればいい。不定積分より数値積分。

これ、計算できる？

災害が発生していたるところに重症被災者がいる
消防署から出動して救急センターに患者を搬送する
消防署から救急センターへの距離は100km
救急車のガソリンは50L、患者を乗せない状態では燃費は10km/L
患者を乗せての燃費は５km/Lである
患者を救える地域の面積はいくらになるか？

こういう問題の意味が小学生にでもわかるような方が難しい。

俺はプログラムして解いたけど最短解なのか確信がもてん。

問題 : 11 リットルと 7 リットルの容器を使って 9 リットルの水を測るにはどうすればいい？

今日は13件と内視鏡が多め。
サテライト、行きつけのステーキハウスで食事して休もうっと。

嘘つき問題はブール代数使うよりプログラム組んだ方が楽だと分かったのは寝当直での収獲だったな。

>>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のものを表示

冗長な記述を同値の命題に変更

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のものを表示

arg <- function(z) ifelse(Arg(z)<0,Arg(z)+2*pi,Arg(z))
arg(-1i)/pi
arg(1i)/pi

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)

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)

"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の数は何か？"

>>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),]

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),]

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),]

オリンピックサイズ・プール50m×25mの水の入ったプールの一つの角に監視員を置く。
この監視員は世界記録で移動するものとする。
水泳100m自由形 46秒91で水中を移動
陸上100m9秒58でプールサイドを」移動
この監視員がプールのどこへでも到達しうるには，最短で何秒必要か計算せよ。

>>149

答：10.77933776秒

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))
}
}

数学板は雑な投稿すると速攻で突っ込みが入る。

https://rio2016.5ch.net/test/read.cgi/math/1572866819/809

具体例を上げて詳述したら納得された。

>151は説明の手順を関数化してみた。
hnr は人間ニュートンラスソンから命名。

ド底辺頭脳医ってプログラム組めないの？
今どき小学生でもできるけど。

# 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)
}

>>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)

検証

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

実用に耐える精度が出た！

> 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

Wolfram

1.414213562373095048801688724209698078569671875376948073176...

0.785398163397448309615660845819875721049292349843776455243...

いい加減にしろ

>>158
ド底辺頭脳医ってプログラム組めないの？
今どき小学生でもできるけど。

こういうのって幼稚園児でも問題の意味はわかると思うよね。

AからEの5人の中に正直者と嘘つきの双方がいるとわかっている。
誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきを全員列挙せよ。

A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」

>>1認知症の原因、実は「処方薬」だった　認知症の疑いがあっても絶対に医者にはかかるな！ [268718286]
https://leia.5ch.net/test/read.cgi/poverty/1579697596/

60以上のデパスやハルシオンでなるぞ
発ガン性もあり

と言うか親の金で塾やら私立行ってた奴が頭悪い

今日の水道の温度は7℃
65℃のお茶を98℃の湯を使って500ｍL作りたい。

yu = function(V=500,T=65,C=7,H=98) V*(T-C)/(H-C)
yu()

発がん性はかなりの薬であるだろ
デパスハルシオンは知らないけど

>>161
さてはシリツだな。

稀ではあるが、治療効果が期待できる認知症ってあるよね。
正常圧水頭症とか慢性硬膜下血腫とか。

いまやシリツ頭脳でも国立大学に入学できる時代なんだなぁ。

日医、自警なら新設は入れるだろ

川崎医大

雨竜郡雨竜町出身介護福祉士やってます。

あの板おもしろいっすw

今は関東に在住してます。医者のなりたかったです

東京医大の入試不正、前文科省局長の息子だったのに、
医科大学が文科省の管轄であるという常識すら欠いているようなアホが裏口入学するのがド底辺シリツ医大。

雨竜は何にもなくダサかった　やっぱ東京がいい

それにしてもあの人達はよく調べたよな。あーゆう人達が本当の頭の良い人なんだな

やはり医者を怒らせたら怖いよな　力持ってるし

決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。

女子大生にフェラをしてもらうには何と言えばいいか？

>>174
ド底辺シリツ医は特にそうみたいだね。

底辺私立医大卒が国立大学医学部卒を妬むスレ [転載禁止](c)2ch.net
http://potato.5ch.net/test/read.cgi/hosp/1424846816/695

695 名前：卵の名無しさん[] 投稿日：2016/09/25(日) 22:04:25.53 ID:FzuOCpqV
匿名先生
精神科
▲
1
評価

私立出でも、素晴らしい医師がいるし、国立でも、どうしようもないのもいる。
そんなのは、あくまで個人の資質差と思ってました。

が、

今いる病院（医師の大部分が、いわゆる底辺私立）に来て、その考えは変わりました。

とある、気に食わない医師の白衣やサンダルを隠す、机の上にゴミを置く、など、小学生じみた嫌がらせをしている医師というのを、はじめて見ました。

国立出の変人も多々みてきましたが、こんな、小学生のイタズラみたいなことをする人は、いませんでした。

やはり、底辺私立は、それなりなんだと思いました。

2016年09月02日 00時21分09秒

あんたは何も知らないんだねw

>>177
フェラしてもらうには何と言えばいい？

一般常識を知らない人が病院組織で上手く立ち回ることはできないぞよw人間関係のストレスがないバイトでもしなさいwまずは医者になってからだなwあ、医者になるのは無理かw

>>179
フェラしてもらうには何と言えばいい？
女子大生にフェラしてもらったことないの？

馴れ馴れしい奴だな　

聖マリ不正入試

https://news.yahoo.co.jp/pickup/6348838

以前は偏差値が５０もなかったから入試の必要もなかった

今もほとんどが金とコネだけだな

格安風俗性病罹患w

雨竜は性病多かった

フェラには興味がないのか？
いや、単に頭が悪いだけなんだろうな。

このスレに合わせてみた。

入試前に寄付をしたド底辺シリツ受験生の親が
寄付金だけとられて不合格になるのでないかと心配して理事長と談判したところ、合格と寄付金は不可分であり、
お子様が面接で正しい主張をすれば
寄付金返却か
裏口合格になる
と言われた。
理事長は嘘をつかない。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？

これに答を出さずに逃げ回るってド底辺医はフェラ嫌いか、頭が悪いんだろうなぁ。

決して嘘をつかない女子大生から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

この方が実情にあってるかもな。

合格と寄付金は不可分を掲げるド底辺医大に入試前に寄付をした受験生の親が寄付金だけとられて不合格になるのでないかと心配して理事長と談判した。
合格と寄付金は不可分という大学の方針は常に正しいが、理事長は常に嘘をつく。
談判では、
「お子様が面接で正しい主張をすれば寄付金返却か裏口合格になる、間違った主張であれば寄付金没収で不合格になる」
と言われた。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？

¬（P⇒Q　∧　¬P⇒¬Q）
¬（P⇒Q）∨ ¬（¬P⇒¬Q）where P⇒Q　≡　¬（P∧¬Q）
¬¬（P∧¬Q）∨ ¬¬（¬P∧¬¬Q）
（P∧¬Q）∨ （¬P∧Q

¬（P⇒Q　∧　¬P⇒¬Q）
¬（P⇒Q）∨ ¬（¬P⇒¬Q）where P⇒Q　≡　¬（P∧¬Q）
¬¬（P∧¬Q）∨ ¬¬（¬P∧¬¬Q）
（P∧¬Q）∨ （¬P∧Q）

>>186
こういう風にすると少し難しくなるな。国立卒でもド底辺頭脳には無理と思う。

某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。

この女子大の学生（嘘つきかどうかは不明）から

「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。

女子大生にフェラをしてもらうには何と言えばいいか？

ニュー速＋にも正しい意見が

49 名無しさん＠１周年 2020/01/23(木) 06:28:26.69 ID:fsvBnYrD0
私立医大はどうでもいい。
金で医者になるところだろ。

# https://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q14116594724
"
６人集まると、その中に「互いに知り合いの３人が存在する」か、または、「互いに知り合いでない３人が存在する」ことを示せ。
"
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)

(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)

>>187
数学板に題材を変えて投稿したら速攻で正解がきたなぁ。

問題は過去を克服することではありません。
さようなことができるわけはありません。
後になって過去を変えたり、
起こらなかったことにするわけにはまいりません。
しかし過去に目を閉ざす者は結局のところ現在にも盲目となります。(ヴァイツゼッカー)

http://imgur.com/Ku3BD42.jpg

患者様からのありがたいお言葉です。
>>
膨大な情報を活用しなければいけない仕事だと認識しているからこそ、命のかかった病気においては特に、せめて自分よりは頭がいい先生に診てもらいたいと思ってしまいます。その場合に、出身校、入学した学校、は嘘をつかないとおもいます。
<<


>出身校、入学した学校、は嘘をつかないとおもいます。

>出身校、入学した学校、は嘘をつかないとおもいます。

>出身校、入学した学校、は嘘をつかないとおもいます。


大切な教えなので　繰り返しました。

##
"
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,])

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)

(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)

一般化して書き直したので４人でも成立することが確認できた。

"A〜Dの4人が受けた試験の結果について、4人はそれぞれ以下のように述べているが、
本当のことを述べているのは1人だけで、他の3人はすべてうそを述べている。
合格したのが2人であるとき、合格した2人と正直者は誰か。

A　Bは合格した。
B　AとCのうち1人だけが合格した。
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

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),]

ある幼稚園で、砂場で遊んでいた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)

>>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

問題は

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

よって正しいのは２

>>204
公務員試験の問題らしいな。医師国試の方が格段に易しいね。
こういう問題をいれて裏口入学の学生を排除すればいいのに。

"
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

>>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

"
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),]

頻繁に使うテンプレートを関数化

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を実行して転置
}

>>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を実行して転置
}

"
A、B、C、Dの４人の体重を調べました。BとDの平均体重が他の2人のどちらかの
体重と同じで、AはBより17�s軽く、CとDとの体重の差は12�sでした。
A　私はCとの体重の差が5�sです。
B　私はCとの体重の差が22�sです。
C　私は４人の体重の順番でAの直前または直後です。
D　私はCより重いです。
4人のうち、2人がうそをついているとすると、嘘を言っている２人は誰ですか。
"

おまえとおまえのおかん

>>214
Aの体重を50kgとして可能な組み合わせで連立方程式を解いて体重の行列をつくって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

>>215
これ面白いだろ？

某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

雨竜郡雨竜町出身介護福祉士でーす　雨竜は寒いでーす

ある人が第１ホールから順にゴルフをプレイして回ったところ、いずれのホールの成績もボギー（＋１）、ダブルボギー（＋２）又はトリプルボギー（＋３）で
あり、あるホールを終えてスコアがちょうど「＋７」となった時点で雨のため中断した。
このとき、各ホールの成績の並びの組合せは何通りあるか。

なお、ボギー、ダブルボギー、トリプルボギーとは、基準打数よりもそれぞれ１打、２打、３打多く打った成績を指しており、基準打数よりも多く打った打数の
総和がスコアである。例えば、第１ホールからの成績の並びがトリプルボギー、ダブルボギー、ボギーであった場合、スコアは（＋３）＋（＋２）＋（＋１）で「＋６」である。

（　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

>>219
女子大生にフェラをしてもらうには何と言えばいい？

>>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)
}

１３段の階段を上るさい、１段上るか、２段上るか、３段上るかのいずれかを組合わせて上るこ
ととする。上り方は何通りあるか。

> stairway(13,3)
[1] 1705

おい、ド底辺。
書き出してカウントしてくれ。
オマエの底辺頭脳ではフェラ問題の答は無理。
数くらい数えられるだろ？

これ↓面白いだろ？
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

福井県越前市瓜生町

考えるよりプログラムを入力して計算させた方が速いな。

"
Ａ〜Ｆの６人が旅行に出かけた。宿泊所では、３人部屋と４人部屋に分かれて泊
まることになった。このとき、ＡとＢの２人が必ず同じ部屋に泊まる部屋割りは
何通りあるか。

"
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
>

>>225

某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

女子大生にフェラしてもらった経験なければ、こっちに答えるのでもいいぞ。

合格と寄付金は不可分を掲げるド底辺医大に入試前に寄付をした受験生の親が寄付金だけとられて不合格になるのでないかと心配して理事長と談判した。
合格と寄付金は不可分という大学の方針は常に正しいが、理事長は常に嘘をつく。
談判では、
「お子様が面接で正しい主張をすれば寄付金返却か裏口合格になる、間違った主張であれば寄付金没収で不合格になる」
と言われた。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？

情報収集の圧倒的な差w

いいことを１つ教えてやろう　Xビデオ　で検索してみろ　お前の欲しい情報がいっぱいあるぞ！w

１円玉３枚、５円玉１枚、１０円玉２枚、１００円玉２枚の合計２２８円を持っ
ている。いま、おつりをもらわないでちょうどぴったり支払いのできる金額は何
通りあるか。　０円はカウントしない。

一行プログラムで答がでる。
> nrow((expand.grid(0:3,c(0,5),c(0,10,20),c(0,100,200))))-1
[1] 71

こういう解き方していると馬鹿になるな。

それすらできない馬鹿がド底辺シリツ医。

女子大生からフェラも引き出せないｗｗｗ

>>228
赤の他人であることもわからんドアホがオマエだよ。
フェラ問題はとけたのか？

>>228
圧倒的な頭脳の差！

シリツ医の頭脳は昆虫並み。

決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

なっ　オナニーネタにちょうどいいだろ　おまえにはw

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

>>233
いや、女子大生ファラの方がいいぞ。

よかったね！ずりネタ１つ増えてw

福井県越前市瓜生町　オナニー好きのあなたへ　気持ちよかったですか？

ある職場は、Ａ〜Ｆの６人の職員で構成され、このうちＥおよびＦは新人職員で
ある。新人職員は、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

> stairway(13,3)
[1] 1705

おい、ド底辺。
書き出してカウントしてくれ。
オマエの底辺頭脳ではフェラ問題の答は無理。
数くらい数えられるだろ？

これ↓面白いだろ？
某女子大には決して嘘をつかない女子大生と必ず嘘をつく女子大生がいることがわかっている。
この女子大の学生（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、間違っていれば何もしてあげない」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

女子大生にフェラしてもらった経験なければ、こっちに答えるのでもいいぞ。

合格と寄付金は不可分を掲げるド底辺医大に入試前に寄付をした受験生の親が寄付金だけとられて不合格になるのでないかと心配して理事長と談判した。
合格と寄付金は不可分という大学の方針は常に正しいが、理事長は常に嘘をつく。
談判では、
「お子様が面接で正しい主張をすれば寄付金返却か裏口合格になる、間違った主張であれば寄付金没収で不合格になる」
と言われた。
裏口入学したいド底辺シリツ受験生は面接でなんと主張すればよいか？

気持ちよかったですか？

福井県越前市瓜生町

〇、△、×のいずれかのマークが１つずつかかれたカードがあります。〇３枚、
△２枚、×２枚の全７枚のカードを、同じマークのカードが隣り合わないように横
一列に並べます。
一番左に○を並べた時、残りの６枚の並べ方は何通りあるか。

割と面倒

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 %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

>>241
答がわからん馬鹿なの？

決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

そうかそうか　気持ちよかったかー　一日一善したw

それにしても彼らの情報収集能力は秀逸だ　おまえはある偉い人々を怒らせたんだよ
しーらないっとw

福井県越前市瓜生町　おまえは通っていた病院で問題起こしたろw

>>246
未だに他人を同定しているどアホがオマエだよ。
まさに昆虫並みの頭脳だな。

ここまでアホだから、藤田か聖マだろ？

これにも答がだせないの？
馬鹿なんだろ。

決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

A〜Gの７人を、３人、２人、２人の３グループに分ける分け方は何通りあるか。

赤、青、白の３色のカードが５枚ずつ、計１５枚あります。各色ともに、１〜５
までの数が１つずつ書かれています。この１５枚の中から３枚を選びます。選ん
だ３枚に書かれた数の和が１２になるような選び方は何通りありますか。

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
>

福井県越前市瓜生町　おまえのことは全部バレバレなんだよw

国立私立どちらも含めて多くの医者を敵を作っちゃったね　ケケケ

これにも答がだせないの？

馬鹿なんだろ。

決して嘘をつかない女子大生から「あなたのいうことが間違っていなければ手コキか、フェラをしてあげる」と言われた。
女子大生にフェラをしてもらうには何と言えばいいか？

>>246
未だに他人を同定しているどアホがオマエだよ。
まさに昆虫並みの頭脳だな。

ここまでアホだから、藤田か聖マだろ？

あたりみたいだな。

いや俺は長崎大学医学部看護学科出身の看護師だw

伝説の福井県越前市瓜生町　ケケケ

俺の知り合いに聞いたら　あーアイツか　あのクルクルパーだな　だってさ　人気があるんだね　うらやましいよ　ケケケ

いろんな病院にかかって　ずいぶんと医者や看護師を困らせているんだな　もっと荒らしていいよw

精神障害者手帳持ってる？

>>258
朱の他人を中傷しているオマエがアホ。
オマエの知り合いだから裏口入学容疑者だろ。

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)

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()

任意の3つの実数を四捨五入して整数にした上で和を取るのと、
和を取ってから四捨五入して整数にする場合で数値が異なる確率を求めよ。
ただし、3つの実数の小数部は[0,1)区間で独立な一様分布とする。

(アクチュアリー資格試験問題 - 一部改正)

四捨五入で思い出した。


IEEE では任意の桁位置における丸め処理法を次のように定めている（ RjpWiki の記事，『工学のためのデータサイエンス入門』 間瀬 茂 他 著 (数理工学社) から引用した）．

一番近い丸め結果候補が1つだけなら，その数に丸める．
一番近い丸め結果候補が2つある場合は，末尾が偶数のものに丸める（五入ばかりでなく五捨もあり得る！）．
丸め処理は1段階で行なわなければならない．

2.5の丸めは2
3.5の丸めは4

これって丸め誤差が小さくなると説明されていた。

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()))

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()))

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()))

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)))

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)))

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()))

>>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)

四捨五入
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)
}

>>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))

> 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(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の丸めの方が圧倒的に誤差が少ないなぁ。
国立卒とﾄﾞ底辺シリツ卒くらいの差がみてとれる。