20/01/10 09:43:17 1nBZfk4O.net
タイトルまんまです
すなわち、タイトルの正しさを検証するスレです。
※前スレ
底辺私立医大を卒業した頭悪いよね？
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part2
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part3
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part4
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part5
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part6
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part7
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part8
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part9
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part10
URLﾘﾝｸ(egg.2ch.net)
底辺私立医大を卒業した医者って頭悪いよね？ Part11
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part12
ｽﾚﾘﾝｸ(hosp板)
底辺私立医大を卒業した医者って頭悪いよね？ Part13
ｽﾚﾘﾝｸ(hosp板)
2:卵の名無しさん
20/01/10 12:29:20 JemsuqLY.net
【参考】
五浪で三流医大に入った古根高の末路↓
病的な虚言癖と妄想癖の精神科医 古根高の病名を診断するスレ
ｽﾚﾘﾝｸ(hosp板)
古根高（サッポロファクトリーメンタルクリニック院長）のご尊顔
URLﾘﾝｸ(i.imgur.com)
3:卵の名無しさん
20/01/10 16:28:31 1nBZfk4O.net
ド底辺シリツ医大卒ってこれだよな
底辺私立医大卒が国立大学医学部卒を妬むスレ [転載禁止](c)2ch.net
ｽﾚﾘﾝｸ(hosp板:695番)
695 名前：卵の名無しさん[] 投稿日：2016/09/25(日) 22:04:25.53 ID:FzuOCpqV
匿名先生
精神科
▲
1
評価
私立出でも、素晴らしい医師がいるし、国立でも、どうしようもないのもいる。
そんなのは、あくまで個人の資質差と思ってました。
が、
今いる病院（医師の大部分が、いわゆる底辺私立）に来て、その考えは変わりました。
とある、気に食わない医師の白衣やサンダルを隠す、机の上にゴミを置く、など、小学生じみた嫌がらせをしている医師というのを、はじめて見ました。
国立出の変人も多々みてきましたが、こんな、小学生のイタズラみたいなことをする人は、いませんでした。
やはり、底辺私立は、それなりなんだと思いました。
2016年09月02日 00時21分09秒
4:卵の名無しさん
20/01/10 16:28:58 1nBZfk4O.net
人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。
これな！
不朽の名投稿
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８－９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<
東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。
5:卵の名無しさん
20/01/10 16:29:40 1nBZfk4O.net
Last but not least, three laws of Do-Teihen Medical School, currently called Gachi'Ura by its graduates.
最後にド底辺医大の三法則を掲げましょう。
1: It is not the bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。
2: The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater which they have gone through.
ド底辺シリツ医大卒は恥ずかしくて、学校名を皆さま言いません。
3: The Do-Teihen graduates are so ashamed of having bought their way into the exclusively bottom-leveled medical school
that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
ド底辺特殊シリツ医大卒は裏口入学の負い目から裏口馬鹿を暴く人間を偽
6:医者扱いしたがる。
7:卵の名無しさん
20/01/10 20:03:40 1nBZfk4O.net
これをコンピュータなしで解ける頭脳には感服する。
赤色のカメレオン３匹と、青色のカメレオン５匹と、黄色のカメレオンが８匹いる。
もし、異なる色の２匹のカメレオンが出会えば、２匹とも３番目の色に変わる。
例えば、青色と黄色が出会えば2匹とも赤色になる。
同じ色の２匹のカメレオンが出会っても色は変わらない。
出会いを重ねて、全てのカメレオンが同じ色になることは可能か？
可能ならば何色になるか？
8:卵の名無しさん
20/01/10 23:45:44 5zvjNhcE.net
糞河野と糞犬今すぐ死ねや
毎日糞犬吠えさせんなカス
9:卵の名無しさん
20/01/11 07:06:03 KGOpuKV/.net
>>6
こんな数になるとコンピュータも処理しきれないが理論値は正解を出してくれるんだよなぁ。
赤色のカメレオン９９匹と、青色のカメレオン１００匹と、黄色のカメレオンが１０１匹の場合はどうか？
赤色のカメレオン７７７匹と、青色のカメレオン８８８匹と、黄色のカメレオンが９９９匹の場合はどうか？
こういうのをサラッと答えられる頭脳が理学部や工学部にいくんだな。
考えようともしない奴がいくのが裏口でド底辺シリツ医大であるのは申すまでもない。
10:卵の名無しさん
20/01/11 07:49:28.69 KGOpuKV/.net
改竄したり隠蔽する安倍が原因だろ。
官僚はアクロバット擁護で安倍に貸しを作って役人が望む政策を実現させてる。
官僚のことだから安倍を脅すネタとして内密に名簿は保管しているだろう。
共産党の議員にまでミンシュガーをやってたのが安倍。
立憲の議員にも「共産党！」とやったのも安倍。
俺は国賓招聘した習近平を指差しして「共産党！！」とやってほしいんだが。
11:卵の名無しさん
20/01/11 08:28:59.74 KGOpuKV/.net
これをプログラムしたいなぁ。
赤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
12:卵の名無しさん
20/01/11 08:29:08.55 KGOpuKV/.net
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
にすることが可能
13:卵の名無しさん
20/01/11 10:20:12 taaoxAl3.net
一体、いつから日本はこうなってしまったのでしょうか？
若者は正社員でも月15万円ほどの手取りでボーナス無し。
子どもの7人に1人が貧困状態。(先進国の中でも最悪な水準)
少子高齢化で人手不足になり、やっと給料が上がろうとしている所に移民を入れ、給与の上昇を抑えようとする政府。
災害で国民の命が奪われ続け、巨万の富を抱えているにも関わらず十分な防災対策をしない国。
デメリットしかないとわかっているにも関わらず実行された消費税の増税。
世紀の大発明であるiPS細胞にですら、支援(年10億円)を打ち切ろうとする政府・財務省。
バカがこの国を運営している。
14:卵の名無しさん
20/01/11 10:37:11 taaoxAl3.net
少子化だし、シリツ大学廃止でいいんじゃね？
1 すらいむ ★ 2020/01/10(金) 16:47:33.70 ID:yEOurkSc9
名城大学で学生が教員刺す 学生の身柄確保
愛知県警察本部によりますと10日午後４時すぎ、名古屋市天白区の名城大学で学生が教員を刺したと通報があったということです。
教員は病院に運ばれましたが、詳しい容体はわかっていません。
警察は学生の身柄を確保し詳しい状況を調べています。
NHK NEWS WEB 2020年1月10日 16時45分
URLﾘﾝｸ(www3.nhk.or.jp)
15:卵の名無しさん
20/01/11 22:58:47 2ClnjBAU.net
"赤色のカメレオン13匹と、青色のカメレオン15匹と、黄色のカメレオンが17匹いる。
もし、異なる色の2匹のカメレオンが出会えば、2匹とも3番目の色に変わる。
例えば、青色と黄色が出会えば2匹とも赤色になる。
同じ色の2匹のカメレオンが出会っても色は変わらない。
このとき、全てのカメレオンが同じ色になることは可能か？"
ct <- function(x,y){ # camelleon transformation
if(x==y)return(0) # 同色なら0を返す
else return((1:3)[-c(x,y)]) # 異色なら第３色を返す
}
ce <- function(x){ # camelleon encounter x=c(1,2,2,3,3,3)
N=length(x)
ab=sample(N,2) # indexから２個選ぶ
a=x[ab[1]] # indexに相当するxの要素
b=x[ab[2]]
z=ct(a,b) # camelleon transformation
if(z!=0){ # 同色でなければ
x[ab[1]]=z # 第３の色に入れ替える
x[ab[2]]=z
}
return(x)
}
is.uc <- function(x){ # is uniform color?
y=unique(x)
if(length(y)==1) return(y) # 一色ならその色番号を返す
else return(0) # それ以外は0を返す
}
ca <- function(r,b,y) rep(1:3,c(r,b,y)) # camelleon array
16:卵の名無しさん
20/01/11 22:58:52 2ClnjBAU.net
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)) # 統一不能 無限ループ
17:卵の名無しさん
20/01/12 07:27:39 aRcuq7iu.net
直感を理詰めでどう証明するかの問題だな。
【問1】宇宙空間に2020個の星が、どの2つの星の距離も異なるように配置されている。
ただし、星は宇宙空間に固定されているとする（公転などは考えない）。
各星には天文学者がいて、自分と最も近い星だけを観測している。
このとき、「どの天文学者にも観測されない星が必ず存在する」は真か偽か。
【問2】2021個のときではどうか。
18:卵の名無しさん
20/01/12 09:09:25.20 aRcuq7iu.net
>>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 がすべての星を含んでいるか
}
}
19:卵の名無しさん
20/01/12 14:11:10 aRcuq7iu.net
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
}
20:卵の名無しさん
20/01/12 14:15:43.57 aRcuq7iu.net
> 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
21:卵の名無しさん
20/01/12 16:50:24 aRcuq7iu.net
星１６個の場合
URLﾘﾝｸ(i.imgur.com)
> 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
22:卵の名無しさん
20/01/12 16:50:48 aRcuq7iu.net
# 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を返す
}
23:卵の名無しさん
20/01/12 16:51:37 aRcuq7iu.net
# すべて観察される配列があればその配列を返す、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
24:卵の名無しさん
20/01/12 21:31:15 aRcuq7iu.net
# 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を返す
}
25:卵の名無しさん
20/01/12 21:31:24 aRcuq7iu.net
# すべて観察される配列があればその配列を返す、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)
# URLﾘﾝｸ(i.imgur.com)
sim(15) # Endless
26:卵の名無しさん
20/01/12 21:33:09 aRcuq7iu.net
すべての星が最短距離の星から観察されている配置が存在するとして、その星の集合をAとする。
Aからa->b b<-aのように相互監視されている星を除く集合をBとする。除かれた星の集合をCと呼ぶ。
Bの要素が1つの星であればその星はどこからも観察されていないからAの前提に反する。
Bの要素の数をｂとする。b人のうち一人でもCに属する星を観察しているとBの中に観察されない星が出現する。
それはAの前提に反するから、b人はすべてBに属する星を観察していることになる。
Bに属する星で星間の距離が最短の2個の星は相互監視していることになりBの前提に反する。
ゆえにBは空集合である。
27:卵の名無しさん
20/01/12 21:39:09 aRcuq7iu.net
The Disappearance of Bambi Woods
URLﾘﾝｸ(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以上の情報はなくてつまらなかった。
28:卵の名無しさん
20/01/13 08:05:25.34 mlDScjUX.net
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
>
29:卵の名無しさん
20/01/13 10:16:11 Z39Ioy5G.net
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
30:卵の名無しさん
20/01/13 11:31:39.04 FPasfEXg.net
治療ハンドブックの最新版が到着。
今年も通読予定。毎年通読していると読むのが楽しみ。
31:卵の名無しさん
20/01/13 19:39:10 FPasfEXg.net
三次方程式の解の公式は面白くないんで三次方程式の解の公式に行く前に、
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
>
32:卵の名無しさん
20/01/13 19:39:45 FPasfEXg.net
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))
33:卵の名無しさん
20/01/14 06:48:44.23 77soHecG.net
平面上に、どの3点も同一直線上にないようにn個の点を配置するとき、それらの中の5点を頂点とする凸五角形が少なくとも1つ存在するためのnの最小値を求めよ。
例、n=6のとき
・ ・
・・
・ ・
のように配置すると、どこにも凸五角形ができないので、n=6は不適。
34:卵の名無しさん
20/01/14 08:28:35.63 77soHecG.net
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)
}
35:卵の名無しさん
20/01/14 08:30:43.28 77soHecG.net
# 一直線上にあるか？外積=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軸で全て正か全て負かで三角形の内部と確認
}
36:卵の名無しさん
20/01/14 08:31:20.10 77soHecG.net
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
37:卵の名無しさん
20/01/14 10:41:45 77soHecG.net
>>32
n=7が不適な例をパソコンで出してみた。
URLﾘﾝｸ(i.imgur.com)
複素平面での７つの座標は
[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
38:卵の名無しさん
20/01/14 12:37:12 77soHecG.net
>>32
こう問題に変えてみた。
１００万を持っているチンパンジーと賭けをする
壁に向かって目をつむって無作為にダーツ矢を６本投げる。
・ ・
・・
・ ・
の配置のように、どの５点を選んでも凸五角形ができない場合はチンパンジーの勝ちであなたは掛け金を全て失う。
凸五角形ができる５点を選べる配置ならチンパンジーから１００万円がもらえる。
掛け金がいくらまでなら有利な賭けといえるか？
39:卵の名無しさん
20/01/14 13:50:37 77soHecG.net
sim.n <- function(n=7,draw=TRUE,concave.draw=TRUE){
ve=c(0,runif(n-1)+runif(n-1)*1i) # vertex 一様分布乱数で頂点をn個発生させて
v3=matrix(ve[combn(n,3)],nrow=3) # 3点の組み合わせ
m=choose(n,3) # m通りに
for(i in 1:n){
aln=align(v3[,i]) # いずれかの3点が直線上にあれば
if(aln) return(NA) # NAを返す
}
m5=choose(n,5) # n個から５個選ぶ組み合わせ
v5=matrix(ve[combn(n,5)],nrow=5) # の5行m5列を作成
convx=FALSE # 凸か?のflag
for(i in 1:m5){
ver5=v5[,i] # i列目が凸か調べる
convx=convex(ver5)
if(convx){ # 凸ならば
if(draw){ # 描画して
plot(ver5,xlim=c(0,1),ylim=c(0,1),bty='l',pch=19,
axes=TRUE, ann=FALSE,asp=1)
combn(5,2,function(x) coep(x,ver5))
points(ve) # n点描画
}
break # for loopから抜ける
}
}
if(!convx & concave.draw){ # 凸が無い場合
print(ve) # ｎ点を表示して赤で描画
plot(ve,xlim=c(0,1),ylim=c(0,1),bty='l',pch=19,col='red',
axes=TRUE, ann=FALSE,asp=1)
}
return(convx)
}
40:卵の名無しさん
20/01/14 13:51:02 77soHecG.net
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)
41:卵の名無しさん
20/01/14 16:58:13 77soHecG.net
これって名投稿だと感心する。
>>(quote)
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８－９割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<(unquote)
＞同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く
同業者の発言：
【ウハも】 開業医達の集い 8診 【粒も】 [無断転載禁止]©2ch.net
>>
670 名前：卵の名無しさん[] 投稿日：2017/05/20(土) 11:15:40.12 ID:46exOAAP
学会で川崎の医者が発表してたら、「馬鹿が何偉そうにしゃべってる。」と思う自分が嫌になるが、
これだけは学生時代から続く反射なので止められない
<<
患者の発言：
【医療】医者は患者にコレを言われると、内心ものすごくムッとする★4 [無断転載禁止]©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..
ド底辺シリツ医大が悪いんじゃない、本人の頭が悪いんだ。
42:卵の名無しさん
20/01/14 17:32:12.29 77soHecG.net
"次の性質(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))
43:卵の名無しさん
20/01/14 17:32:34.24 77soHecG.net
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)
44:卵の名無しさん
20/01/14 18:41:59.76 77soHecG.net
療養病床勤務医・老健施設勤務医・回復期病棟勤務医が集うスレ☆ 復活版 Part.4
ｽﾚﾘﾝｸ(hosp板:992番)
992 名前：卵の名無しさん[] 投稿日：2019/07/20(土) 19:03:01.20 ID:WRt9tsMe
当直と外来の穴埋めは看護師も事務も断れないやつとか変な正義感持っているやつを狙って頼んでる
表面的には助かりますー、先生のおかげですーとかおだててるが内心はまたコイツにやらせとけばいいだろって下に見てる
感謝なんかしてない、そいうバカが変にフォローするせいで外来の医師数増やさなかったり、当直医を本気で確保しない（当直料安い）
キレイなキラキラした病院でしか働いたことないやつはこのことをわかってない。ブラック病院では標的にされて良いように使われるだけ。
大抵は俺だけ働いている、俺だけちゃんとした医療してるって頭おかしくなってやめていく
なんにもしない爺医や患者を適当にあしらう医師は人生経験でそのこと知ってんだよ
わかってて無能、バカのフリしてんだよ、気づけよ
45:卵の名無しさん
20/01/15 07:22:43 kW1YHlTi.net
>>41
存在証明だけなのだから
整数に拘らず有理数解があることを証明するだけでいいと数学板で教わった。無駄なプログラムを書いていたもんだ。
46:卵の名無しさん
20/01/15 07:27:35 kW1YHlTi.net
１００万円を持っている女子大生と賭けをする
壁に向かって目をつむって無作為にダーツ矢を６本投げる。
・ ・
・・
・ ・
の配置のように、どの５点を選んでも凸五角形ができない場合は女子大生の勝ちであなたは掛け金を全て失う。
凸五角形ができる５点を選べる配置なら女子大生から１００万円がもらえる。
掛け金がいくらまでなら有利な賭けといえるか？
47:卵の名無しさん
20/01/15 10:51:52.08 OeJJ0/zs.net
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))
48:卵の名無しさん
20/01/15 15:36:37.08 kW1YHlTi.net
簡単にシミュレーションできると思ったがわりとデバッグに苦労した。配列にA[0]を許さないRに対応してようやく完成。
プログラミングはロジカル思考の訓練になるなぁ。
知恵遅れのド底辺シリツには無理だけど。
##
rm(list=ls())
"
円に内接する凸n角形(n≧4)を、1つの頂点からn-3本の対角線を引く�
49:ｱとによって、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)
50:卵の名無しさん
20/01/15 18:10:41 kW1YHlTi.net
巨大数を扱えるというふれこみの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:
URLﾘﾝｸ(www.wolframalpha.com)
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
51:卵の名無しさん
20/01/15 18:25:32.39 iAhmqtla.net
藤林丈司
52:卵の名無しさん
20/01/15 22:07:31 kW1YHlTi.net
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)
53:卵の名無しさん
20/01/16 06:57:58.74 NA2F+h8M.net
災害が発生していたるところに重症被災者がいる。消防署から出動して救急センターに患者を搬送する
消防署から救急センターへの距離は100km 救急車のガソリンは20L、患者を乗せない状態では燃費は10km/L、患者を乗せての燃費は5km/Lである
患者を救える地域の面積はいくらになるか？
砂漠の基地Aからもうひとつの基地Bに向かって出かけた戦車がGPSの故障で進路を見失ってさまよった挙げ句にガス欠で止まってしまった。別の戦車で基地Aからこの戦車へ救助に向かい、燃料を分け与えて一緒に基地Bに行くことになった。しかし、戦車にはAB間をちょうど往復できるだけしか燃料は積めない。AB間の距離をRとして、救出可能な領域の面積をRを使って表しなさい。
ただし、戦車の燃費はいずれも同じものとする。
54:卵の名無しさん
20/01/16 14:26:58 NA2F+h8M.net
このシミュレーションプログラミングを考えているんだが、手がつかん。
密度が一定の球形の惑星がある。
この惑星の表面点Pでの重力を最大にするために体積と密度一定の条件で惑星の形を
55:変形する。 点Pでの重力は球形のときの最大何倍になるか。 またこのときの惑星の形はどうなるか。
56:卵の名無しさん
20/01/16 19:03:25 BUr06wt4.net
私立医から医師免許を剥奪せよ。
私立医学部を全て解体せよ。
本来なら到底医師になれず、
カネで医者になったニセ医者を迫害仲間外れにして糾弾しょう！
奴らの行動は患者殺しと国立医への逆差別、嫌がらせである。
私立医の患者殺しを許すな！
57:卵の名無しさん
20/01/16 23:28:50 NA2F+h8M.net
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)
58:卵の名無しさん
20/01/16 23:47:17 NA2F+h8M.net
>>52
シミュレーションの準備に まず、質点を一様分布配置しなくちゃならん。
URLﾘﾝｸ(i.imgur.com)
じゃだめで
URLﾘﾝｸ(i.imgur.com)
だな
59:卵の名無しさん
20/01/16 23:52:33 NA2F+h8M.net
改憲以前に安倍の議員資格に疑問符がついているんですが…
まず桜を見る会での公費による有権者買収は
公職選挙法違反なんじゃないんですか？
60:卵の名無しさん
20/01/17 00:30:38 HL0McwQD.net
# 単位円内に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)
61:卵の名無しさん
20/01/17 10:22:18 HL0McwQD.net
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)
62:卵の名無しさん
20/01/17 10:23:54 HL0McwQD.net
1歩で1段または2段のいずれかで階段を昇るとき、1歩で2段昇ることは連続しないものとする。15段の階段を昇る昇り方は何通りあるかを求めよ。
63:卵の名無しさん
20/01/17 10:26:40.96 HL0McwQD.net
>>59
13段の方が興味をもたれるかな。
64:卵の名無しさん
20/01/17 12:07:43 HL0McwQD.net
>>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
65:卵の名無しさん
20/01/17 12:21:43 HL0McwQD.net
>>61
不定長さの扱えるHaskell だと
72428926648244409475025064434360344869036057631812487293633154973026123860267808630768185105127650714516006185052479679199
になった。
66:卵の名無しさん
20/01/17 12:46:59 HL0McwQD.net
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の結果と一致
67:卵の名無しさん
20/01/17 13:32:20.32 HL0McwQD.net
ｍ３から配信されてきたのだが、最近の国試って馬鹿を医者にする問題だらけだな。
第112回医師国試問題◆
胸やけの誘因となりにくいのはどれか。
国試112-第112回医師国家試験問題解説書（医学評論社）より出題 学生正答率80.2％
A 過食
B 運動
C 肥満
D 高脂肪食
E 前屈姿勢
こういう問題を出して、裏口入学のシリツ医学生を排除すべきだと思うよなぁ。
1歩で1段または2段のいずれかで階段を昇るとき、1歩で2段昇ることは連続しないものとする。15段の階段を昇る昇り方は何通りあるかを求めよ。
68:卵の名無しさん
20/01/17 15:00:52 HL0McwQD.net
こういう非定型な問題で場合分けして系統的にカウントできるかを問う問題は他人の結果と合致すると自分の計算に確信がもてる。
1歩で1段または2段のいずれかで階段を昇るとき、2歩で続けて昇ることは3回以上は連続しないものとする。（2回続けて2歩で昇るのは可）
15段の階段を昇る昇り方は何通りあるかを求めよ。
69:卵の名無しさん
20/01/17 18:51:10.19 HL0McwQD.net
床に描かれた、ある程度大きな円の円周上に人物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のこの行動は最善なのか？
70:卵の名無しさん
20/01/17 18:51:24.86 HL0McwQD.net
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:2
71:0){ t=0.05 pt(p2d(pi+v*t*n),n) pt(p2d(v*t*n,t*n),n,col=2) } } sim(pi)
72:卵の名無しさん
20/01/17 19:16:56.48 HL0McwQD.net
2018年2月5日月曜日
自殺者と変死体
最近では，遺書や目撃証言といった具体的な証拠がない限り，自殺とは認めず，変死体として処理するようになったと（上記記事）。
統計上の自殺者数というのは，現場のさじ加減で動くものです。統計には表れない「暗数」に思いを巡らせないといけません。
URLﾘﾝｸ(tmaita77.blogspot.com)
URLﾘﾝｸ(2.bp.blogspot.com)
73:卵の名無しさん
20/01/17 20:33:03 HL0McwQD.net
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])
74:卵の名無しさん
20/01/17 20:38:48 HL0McwQD.net
>>66
こういうふうに逃げれば
URLﾘﾝｸ(i.imgur.com)
75:卵の名無しさん
20/01/18 06:33:53.78 WIxYzvkW.net
トランプだって在日米軍を引き上げたがり、日本の核武装を容認発言をしていたわけで、
日本が自立する最高のチャンスだった。
しかしそれすらも活かせず何の戦略もなくアメリカの犬であり続けたのが安倍自民政権。
それだけでなく安倍がネオリベグローバリストで保守的政策は何一つやらずに売国政策だけは猛スピードで次々に成立させている現実を見れば、
安倍の本性がなにかなんて馬鹿でもわかる問題だろうが。
明らかに安倍はDSのスパイなわけだが、100歩譲って仮にそうではなく仕方なく売国をやっていると仮定しても
日本自立の最大のチャンスすら何も活かせずにアメリカの犬をやり続け、
保守的政策は何一つやらずに売国政策だけは猛スピードで次々に成立させているわけで、
保守にとって何のメリットもない総理で、今すぐ倒閣し日本の破壊を止めさせる以外の選択肢はない。
水島は安倍売国に加担し続けている売国奴として我々保守は絶対に許すわけにはいかん。
76:卵の名無しさん
20/01/18 06:45:16.37 WIxYzvkW.net
>>415
自民党憲法改正草案(本文抜粋)
・第12条
憲法が国民に保障する自由及び権利は、国民の不断の努力により、保持されなければならない。
国民は、これを濫用してはならず、自由及び権利には責任及び義務が伴うことを自覚し、常に公益及び公の秩序に反してはならない。
・第13条
生命、自由及び幸福追求に対する国民の権利については、公益及び公の秩序に反しない限り、立法その他の国政の上で、最大限に尊重されなければならない。
・第14条-3
(現行)栄誉、勲章その他栄典の授与は、いかなる特典も伴わない。
↓
(改正)
栄誉、勲章その他栄典の授与は、現にこれを有し、又は将来これを受ける者の一代に限り、その効力を有する。
・第18条
(現行)何人も、いかなる奴隷的拘束も受けない。
↓
(改正)社会的又は経済的関係において身体を拘束�
77:ｳれない。 ・第29条 財産権の内容は、公益及び公にけ公の秩序に適合するように、法律で定める。 ・第97条 (現行)この憲法が日本国民に保障する基本的人権は、人類の多年にわたる自由獲得の努力の成果であつて、 これらの権利は、過去幾多の試錬に堪へ、現在及び将来の国民に対し、侵すことのできない永久の権として信託されたものである。 ↓ (改正)全文削除
78:卵の名無しさん
20/01/18 06:49:25 WIxYzvkW.net
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人嘘つきがいる」
79:卵の名無しさん
20/01/18 07:20:50 WIxYzvkW.net
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
80:卵の名無しさん
20/01/18 07:28:00.76 WIxYzvkW.net
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
81:卵の名無しさん
20/01/18 07:54:16 WIxYzvkW.net
"
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
82:卵の名無しさん
20/01/18 07:58:05 WIxYzvkW.net
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
83:卵の名無しさん
20/01/18 08:43:49.79 yI6kNSEQ.net
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)
84:卵の名無しさん
20/01/18 08:51:00.86 yI6kNSEQ.net
通勤中に解決してしまった。
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
85:卵の名無しさん
20/01/18 09:46:55 yI6kNSEQ.net
こういう問題にすると
A「8人の中に、1人正直者がいる」
B「8人の中に、2人正直者がいる」
C「8人の中に、3人正直者がいる」
D「8人の中に、4人正直者がいる」
E「8人の中に、1人嘘つきがいる」
F「8人の中に、2人嘘つきがいる」
G「8人の中に、3人嘘つきがいる」
H「8人の中に、4人嘘つきがいる」
面白い。
86:卵の名無しさん
20/01/18 09:57:30 yI6kNSEQ.net
こっちの方が面白い。
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人嘘つきがいる」
87:卵の名無しさん
20/01/18 10:42:43.57 rf9LEbFI.net
頭の良さはこういう事では無いと思うのですが
88:卵の名無しさん
20/01/18 10:50:28.63 yI6kNSEQ.net
>>82
んで>81の答は？
89:卵の名無しさん
20/01/18 11:23:13 yI6kNSEQ.net
# 嘘つきは必ず嘘の証言をする設定に変更
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),]
90:卵の名無しさん
20/01/18 12:07:11.98 yI6kNSEQ.net
嘘つきを人数でなくて特定する証言を含めてみる。
AからHの8人はそれぞれ正直者か嘘つきであり、次のように証言している。証言でも嘘つきは必ず嘘をつく。
嘘つきは誰か、全ての組合せを答えよ。
A「8人の中に、正直者は3人いる」
B「8人の中に、少なくとも2人正直者がいる」
C「Bは嘘つきである」
D「Cは嘘つきである」
E「8人の中に、少なくとも1人嘘つきがいる」
F「8人の中に、少なくとも2人嘘つきがいる」
G「Eは嘘つきである」
H「AもFも嘘つきである」
URLﾘﾝｸ(o.5ch.net)
91:卵の名無しさん
20/01/18 13:41:30 tceOUuEc.net
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),]
92:卵の名無しさん
20/01/18 14:41:10.06 tceOUuEc.net
>>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),]
93:卵の名無しさん
20/01/18 15:17:30.87 tceOUuEc.net
どこかに、問題がないかと探っていたら、中学入試の問題があった
"Ａ、Ｂ、Ｃ、Ｄ、Ｅの５人のうち２人は常に本当のことを言う正直者です。\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"
うまく動作している！
94:卵の名無しさん
20/01/18 15:17:50.49 tceOUuEc.net
Ａ、Ｂ、Ｃ、Ｄ、Ｅの５人のうち２人は常に本当のことを言う正直者です。あとの３人は嘘つきですが、その発言内容は本当のときもあります。
彼らに誰が嘘つきか尋ねたところ次のように答えました。
Ａ 「ＢとＥは嘘つきではない」
Ｂ 「Ｃは嘘つきだ」
Ｃ 「Ｄは嘘つきだ」
Ｄ 「Ｅは嘘つきだ」
Ｅ 「ＢとＣは嘘つきだ」
さて、正直者は誰と誰でしょうか？
95:卵の名無しさん
20/01/18 15:57:50 tceOUuEc.net
>>82
おい、ド底辺。
これに答えてみ！小学生向けの問題だけど。
あるド底辺シリツ医大の新入生が101人とする。
101人は学力考査で入学した正規学生と金で入った裏口学生のどちらかである。
正規学生は常に本当のことを語り、裏口学生は虚しか言わない。
全ての学生は誰が裏口学生かを知っている。
学生全員に裏口学生は何人いるかを尋ねたところ
学生A1「1人以上の裏口学生がいる」
学生A2「2人以上の裏口学生がいる」
学生A3「3人以上の裏口学生がいる」
...
学生A99「99人以上の裏口学生がいる」
学生A100「100人以上の裏口学生がいる」
学生A101「全員が裏口学生である」
との証言を得た。
以上の情報からこのド底辺シリツ医大の裏口学生は誰かを述べよ。
96:卵の名無しさん
20/01/18 17:26:15.10 tceOUuEc.net
>>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),]
97:卵の名無しさん
20/01/18 20:59:59 hxUSgTgG.net
↑荒らし工作員
底辺私立医の実態バラされると死人がでるからな
98:卵の名無しさん
20/01/19 06:38:09 xnoCWj3b.net
>>92
どの問題にも解答の出せないアホ乙！
99:卵の名無しさん
20/01/19 06:40:13 xnoCWj3b.net
>>91
数学板に書いたら、プログラムで解いたと俺のとは違う解答が出てきたので困惑していたが
論証で解答してくれた人がいて、俺の解答と一致していたので俺のプログラムの方が正しいと確信できてほっとした。
このスレにはド底辺頭脳しか、いないようだな。
100:卵の名無しさん
20/01/19 06:59:05.78 THI8CzOO.net
>>93
ズボシで必死だな
死ねキチガイ
101:卵の名無しさん
20/01/19 07:02:39.34 xnoCWj3b.net
>>95
どれか答えてみ！
ド底辺頭脳クン
102:卵の名無しさん
20/01/19 07:06:19.38 xnoCWj3b.net
お～い、ド底辺。
これに答えてみ！小学生向けの問題だけど。
あるド底辺シリツ医大の新入生が101人とする。
101人は学力考査で入学した正規学生と金で入った裏口学生のどちらかである。
正規学生は常に本当のことを語り、裏口学生は虚しか言わない。
全ての学生は誰が裏口学生かを知っている。
学生全員に裏口学生は何人いるかを尋ねたところ
学生A1「1人以上の裏口学生がいる」
学生A2「2人以上の裏口学生がいる」
学生A3「3人以上の裏口学生がいる」
...
学生A99「99人以上の裏口学生がいる」
学生A100「100人以上の裏口学生がいる」
学生A101「全員が裏口学生である」
との証言を得た。
以上の情報からこのド底辺シリツ医大の裏口学生は誰かを述べよ。
103:卵の名無しさん
20/01/19 07:23:42 xnoCWj3b.net
国立卒でも>97の裏口学生が誰か答えられないなら裏口学生レベルの頭脳だね。
104:卵の名無しさん
20/01/19 07:42:52 xnoCWj3b.net
フェルマーの最終定理の反例が見つかる
これをみつけられるのは、それはそれで凄いな。
URLﾘﾝｸ(livedoor.blogimg.jp)
105:卵の名無しさん
20/01/19 08:12:15.97 aWxCzLG/.net
↑あほ
106:卵の名無しさん
20/01/19 10:15:45 Vy3HtBW+.net
外来やっていても単語しか会話できない知恵遅れっているよな。
おい、ド底辺。
これに答えてみ！小学生向けの問題だけど。
あるド底辺シリツ医大の新入生が101人とする。
101人は学力考査で入学した正規学生と金で入った裏口学生のどちらかである。
正規学生は常に本当のことを語り、裏口学生は虚しか言わない。
全ての学生は誰が裏口学生かを知っている。
学生全員に裏口学生は何人いるかを尋ねたところ
学生A1「1人以上の裏口学生がいる」
学生A2「2人以上の裏口学生がいる」
学生A3「3人以上の裏口学生がいる」
...
学生A99「99人以上の裏口学生がいる」
学生A100「100人以上の裏口学生がいる」
学生A101「全員が裏口学生である」
との証言を得た。
以上の情報からこのド底辺シリツ医大の裏口学生は誰かを述べよ。
107:卵の名無しさん
20/01/19 10:38:32.87 eyf/xZES.net
死ね
108:卵の名無しさん
20/01/19 10:52:19.70 Vy3HtBW+.net
どれか答えてみ！
ド底辺頭脳クン
109:卵の名無しさん
20/01/19 11:28:23 FHzWRk73.net
へー
110:卵の名無しさん
20/01/19 11:28:42 FHzWRk73.net
ノータリン
111:卵の名無しさん
20/01/19 11:29:07 xnoCWj3b.net
>>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,]
112:卵の名無しさん
20/01/19 11:29:42 xnoCWj3b.net
>>106
３桁までの数字では、画像にある１組だけだな。
> cmb[idx,]
[1] 70 212
113:卵の名無しさん
20/01/19 12:59:49 xnoCWj3b.net
1cm刻みの数直線上の隣り合った3点に●がある(3つの●は区別がない)。
1秒ごとに●1つが他の●1つを飛び越す。
その場合、飛んだあとの●の位置は、飛び越された●からの距離がもとの位置からと同じになるようにする。
飛んだあとの位置にすでに他の●があるときは飛べない。
2021秒後に、●●●全てが開始位置に戻ることはできるか？
例、
●●●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ●●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ● ●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ● ●
╋╋╋╋╋╋╋╋╋╋╋
↓
：
↓
●●●
114:卵の名無しさん
20/01/19 13:00:13 xnoCWj3b.net
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)
115:卵の名無しさん
20/01/19 13:00:58 xnoCWj3b.net
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()
116:卵の名無しさん
20/01/19 13:01:16 xnoCWj3b.net
このジャンプで復帰するね。
> sim()
0 : 1 2 3
1 : 1 3 4
2 : 1 4 5
3 : -2 1 5
4 : 1 4 5
5
117: : -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
118:卵の名無しさん
20/01/19 16:37:41 xnoCWj3b.net
URLﾘﾝｸ(www.nippon-foundation.or.jp)
日本財団「18歳意識調査」第20回 テーマ：「国や社会に対する意識」（9カ国調査）国や社会に対する意識 各国1,000人に聞く
日本の若者の数字の低さ際立つ
URLﾘﾝｸ(www.nippon-foundation.or.jp)
URLﾘﾝｸ(www.nippon-foundation.or.jp)
119:卵の名無しさん
20/01/19 16:53:55 xnoCWj3b.net
英語に変更します。
> 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"
120:卵の名無しさん
20/01/19 17:11:43.50 xnoCWj3b.net
正解が一つの方が問題としては美しいな。
"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),]
121:卵の名無しさん
20/01/19 18:55:45.90 xnoCWj3b.net
>>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
122:卵の名無しさん
20/01/19 18:56:02.25 xnoCWj3b.net
[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
123:卵の名無しさん
20/01/19 18:56:21.96 xnoCWj3b.net
[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
124:卵の名無しさん
20/01/19 18:56:34.24 xnoCWj3b.net
[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
ようやく、復帰！
125:卵の名無しさん
20/01/19 20:21:15 3Op6J520.net
シリツ医はやくざと同じ
大半がチョン
税金還付で実質脱税
126:卵の名無しさん
20/01/19 20:29:50 xnoCWj3b.net
論理的思考が試される問題
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も正直者である」
127:卵の名無しさん
20/01/19 21:16:09 xnoCWj3b.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も正直者である」
128:卵の名無しさん
20/01/19 21:23:30.33 xnoCWj3b.net
2013年の発言
129:http://imgur.com/MGbp1xx.jpg ↓ http://imgur.com/LU42EXK.jpg
130:卵の名無しさん
20/01/19 23:05:49 xnoCWj3b.net
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)
131:卵の名無しさん
20/01/19 23:08:31 xnoCWj3b.net
>>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
だが、処理速度が落ちるな。
132:卵の名無しさん
20/01/19 23:27:51 xnoCWj3b.net
イランイラク戦争のあとイラクに傀儡政権としてフセイン政権を作ったのはCIA
同じように敗戦後日本に自民党という傀儡政権をCIAが作った
CIAはA級戦犯の児玉誉士夫と岸信介に目をつけた
自民党を作ったのは緒方竹虎自由党と鳩山一郎民主党
その緒方の後継者となり自民党総裁となり60年日米安保を改正した
緒方→CIA工作員、岸、鳩山を総裁につけた児玉誉士夫→CIA工作員、
岸信介→CIA工作員、佐藤栄作→CIA工作員、読売オーナー正力松太郎→CIA工作員
かくして自民党は1958年からCIAからの7000万～1億円(現在の10億円)の資金を十数回受け取り政界工作をして55年体制を維持した
50年経過しアメリカ情報公開法により2007年CIAからの資金提供やコードネームなどが明らかとなった
現代でもアメリカの意向に逆らえば政界から失脚する
田中角栄や小沢一郎がいい例
133:卵の名無しさん
20/01/20 05:00:41 oNRiQofW.net
>>124
４桁までを探索させた結果は
[,1] [,2]
[1,] 70 212
[2,] 700 2120
134:卵の名無しさん
20/01/20 05:20:38 oNRiQofW.net
>>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
135:卵の名無しさん
20/01/20 06:42:53 oNRiQofW.net
AからEの5人はそれぞれ正直者か嘘つきのどちらかであり、誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきの可能性があるのは誰か？
A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」
136:卵の名無しさん
20/01/20 07:24:27 oNRiQofW.net
馬鹿は、死ななきゃ治らない の 待遇を 馬鹿�
137:ｪ治れば死んでいる と答えるのがド底辺シリツ医だから >128に正答できないだろうな。 ド底辺頭脳の国立卒も似たようなものだが。
138:卵の名無しさん
20/01/20 07:25:10 oNRiQofW.net
馬鹿は、死ななきゃ治らない の 対偶を 馬鹿が治れば死んでいる と答えるのがド底辺シリツ医だから
>128に正答できないだろうな。 ド底辺頭脳の国立卒も似たようなものだが。
139:卵の名無しさん
20/01/20 07:47:14 ICG8EQVz.net
高校数学おもろいか？w
次元の違う高等数学ゴマンとあるぞ
あほやろおまえ
140:卵の名無しさん
20/01/20 07:52:44 oNRiQofW.net
可能性より特定できた方が問題らしいな。
"
AからEの5人の中に正直者と嘘つきの双方がいるとわかっている。
誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきを全員列挙せよ。
A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」
"
141:卵の名無しさん
20/01/20 07:54:12 oNRiQofW.net
>>131
じゃあ、大学教養レベルの問題に答えてみ！
これはシリツ医や底辺国立頭脳には無理だろうな。
ド底辺シリツ裏口調査団が100人を順次調査した。
裏口判明人数をそのまま公表はヤバすぎる結果であったため、
連続して裏口がみつかった最大数は4人であったとだけ公表した。
公表結果が正しいとして裏口入学人数の期待値、最頻値、及び95%信頼区間を述べよ。
142:卵の名無しさん
20/01/20 07:57:43 oNRiQofW.net
臨床医学は確率事象を扱う。
少数データから断定するのはアホ。
信頼区間を設定して議論するのが臨床医の正しい姿。
これ答えてみ。高校数学の知識では信頼区間算出は無理。
ド底辺シリツ医にも当然、無理。
ド底辺シリツ医大にある年の裏口受験票が保存されていることが判明した。
枚数は１０００枚以下であることは判明しているが、その数を推定したい。
裏口入学調査員が無作為に１０枚選んで番号を記録して元に戻した。
別の裏口入学調査員が無作為に２０枚選んで番号を記録した。
二人の調査員の記録した番号を照合すると３枚の番号が一致していた。
この情報からこの年の裏口受験証の数の最頻値、期待値とその９５％信頼区間を求めよ。
143:卵の名無しさん
20/01/20 08:01:56 oNRiQofW.net
>>66
Aと円の中心と対称となる点を、目指して常に方向転換するという方針で、定速で作図してみた。
偏角を使った方が楽だった。数値積分を組み込むは煩瑣。
URLﾘﾝｸ(i.imgur.com)
Bの曲線の接線と円周の交点がその時点でのAと対称になっている（はず）。
図はAの速度がBの１．２倍のとき。
144:卵の名無しさん
20/01/20 08:17:09 oNRiQofW.net
これ生データがほしいなぁ
URLﾘﾝｸ(link.springer.com)
【医学】胸の大きな女性は小さな女性より風邪の治りが2倍以上も遅い！ポーランドの最新研究で明らかに
ｽﾚﾘﾝｸ(newsplus板)
標準回帰係数が
respiratory infections (standardized β?=?0.20, p?=?.011), their duration (standardized β?=?0.18, p?=?.027),
145:卵の名無しさん
20/01/20 08:21:36 oNRiQofW.net
>>131
別に実用上必要な数値がだせればいい。不定積分より数値積分。
これ、計算できる？
災害が発生していたるところに重症被災者がいる
消防署から出動して救急センターに患者を搬送する
消防署から救急センターへの距離は100km
救急車のガソリンは50L、患者を乗せない状態では燃費は10km/L
患者を乗せての燃費は５km/Lである
患者を救える地域の面積はいくらになるか？
146:卵の名無しさん
2020/01/
147:20(月) 08:27:16 ID:oNRiQofW.net
148:卵の名無しさん
20/01/20 12:24:58 NsqRuhKW.net
今日は13件と内視鏡が多め。
サテライト、行きつけのステーキハウスで食事して休もうっと。
149:卵の名無しさん
20/01/20 12:37:50 Ss+GOu8N.net
嘘つき問題はブール代数使うよりプログラム組んだ方が楽だと分かったのは寝当直での収獲だったな。
150:卵の名無しさん
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