20/01/14 12:40:34 HETN1jKh.net
回の配置だと8でもダメ
612:132人目の素数さん
20/01/14 12:42:50 HETN1jKh.net
対角線が一致しないようにチョットずらして
613:132人目の素数さん
20/01/14 12:44:08 Kk48Gk/U.net
>>579
n>8は確認できた。
URLリンク(i.imgur.com)
小数だと見にくいので整数になる8点を探索させた。
複素平面での値は
> seek.concave.int(8)
[1] 0+ 0i 345+358i 165+362i 58+ 82i 460+ 16i 264+130i 334+342i 238+345i
614:132人目の素数さん
20/01/14 12:44:17 HETN1jKh.net
あーずらすとダメか
615:132人目の素数さん
20/01/14 12:46:44 HETN1jKh.net
回の外の口を横に縮めて目に近くしたら反例になるから
やっぱりn>8
616:132人目の素数さん
20/01/14 12:49:28 HETN1jKh.net
5角形ということで外郭は三角形か四角形かしかないわけだし
四角形2つの入れ子が反例になるのが最後じゃないかな
たぶんn=9
617:132人目の素数さん
20/01/14 12:54:41 HETN1jKh.net
さて証明はどうしよう
おそらく引き出し論法使うんだろうって気がするけど
見込み角の選択をどう表現するかな
618:132人目の素数さん
20/01/14 13:59:22.56 Kk48Gk/U.net
>>585
俺もその予感。
n=9だとシミュレーションプログラムがwhileループから抜け出せない。
プログラムのバグで抜け出せない可能性もあるけど、n=8までは上手くいった。
>見込み角の選択をどう表現するかな
シミュレーションだと外積つかって点が三角形の内部にあるかどうか判断した。
619:132人目の素数さん
20/01/14 15:11:49 HETN1jKh.net
9点をすべて含む凸領域の共通部分を考えると
点のうちいくつかを頂点として持つ凸多角形となることはすぐ証明できるはず
それらの点を最外殻と呼ぶことにしよう
とするとその凸多角形が5角形であればだめだから
最外殻は四角形か三角形
それを取り除くと内部に残るもので同様に考えて最外殻は四角形か三角形
それも取り除くと内部には・・・
四角形が四角形を含み内部に1点
四角形が三角形を含み内部に2点
三角形が四角形を含み内部に2点
三角形が三角形を含み内部に三角形
の4パターンしかない
620:イナ ◆/7jUdUKiSM
20/01/14 16:54:16 +pxsRseD.net
>>578俺の辞書には肉薄と出てた。肉薄の説明の中に肉迫と書いてあったからどっちでもいいとわかった。肉がうすいより肉に迫るほうが勢いあっていい。
/∥__∥∩∩]∥ |゚○。
|∩∩((-。-)。∥ ∩∩ ゚
( (`)(っ[ ̄]∥(`) )゚
( ̄ ̄)「 ̄ ̄]∥(_υ_)~
(__)_∩∩UU□∥∩∩~ ~
~ ~~(____)~~~ (`)__)~ ~ ~~~ ~~~ ~~~ ~~~ ~~~~ ~~~ ~~~ ~ ~ ~ ~ ~~前>>567
621:132人目の素数さん
20/01/14 17:51:16 Kk48Gk/U.net
>>589
俺が高校生の頃は肉迫は不正解だったなぁ。
試験に出る英単語 という本には constructの対義語はdestroy
destructとなんて単語は存在しないと書かれていたけど
今じゃ、起爆装置detonatorでロケットとかを爆破するときはdestructが使われる。
622:132人目の素数さん
20/01/14 18:01:31 ITApl1nr.net
>>568
M=a_2020、S={a_1,‥}とおく。
(i)M:oddのとき。
和がMになる組み合わせ(1,M-1),‥,((M-1)/2,(M+1)/2)の中の高々ひと組みしか両方Sに属するものがないとSの元数は高々(M-1)/2+1+1=(M+3)/2≦2019となり#S=2020に矛盾する。
∴ 上記ペア中2組以上が共にSに属する。
(i)M:evenのとき。
和がMになる組み合わせ(1,M-1),‥,(M/2-1,M/2+1)の中の高々ひと組みしか両方Sに属するものがないとSの元数は高々M/2-1+3=(M+4)/2≦2019となり#S=2020に矛盾する。
∴ 上記ペア中2組以上が共にSに属する。
623:イナ ◆/7jUdUKiSM
20/01/14 19:11:26 +pxsRseD.net
前>>589
>>590デストロイといえば仮面ライダーV3の敵は地獄のデストロン。原案によるとデストロイヤーになってたかもしれなんだって今知った。
624:イナ ◆/7jUdUKiSM
20/01/14 21:23:57 +pxsRseD.net
前>>592
>>565をなんとか解読しないな。なんとか解読できないものか。
625:132人目の素数さん
20/01/14 22:23:15.37 Kbe89GpS.net
ヒント
答えが
(2 - 7 t + 6 t^2 - 3 t^3)/((1 - t + t^2)^2 sqrt(4 - 4 t + 3 t^2))
なのでt=0のとき2/√2、t=1のとき-2/√3など計算しやすい値を自分の計算結果に入れてみて検算していけばどこで間違ったか見つけられる。
626:132人目の素数さん
20/01/14 22:44:25 kjiwvSoy.net
>>593
(√f)/g の微分 ただし、f=3x^2-4x+4、g=x^2-x+1
((√f)/g
627:)' ={(√f)’*g-(√f)*g’}/g^2 ={(1/2)(1/√f)*f’*g-(√f)*g’}/g^2 ={(1/2)*f’*g-f*g’}/{g^2*√f} ={(1/2)*(6x-4)*(x^2-x+1)-(3x^2-4x+4)*(2x-1)}/{(x^2-x+1)^2*√(3x^2-4x+4)} ={(3x-2)*(x^2-x+1)-(3x^2-4x+4)*(2x-1)}/{(x^2-x+1)^2*√(3x^2-4x+4)} ={(3x^3-5x^2+5x-2)-(6x^3-11x^2+12x-4)}/{(x^2-x+1)^2*√(3x^2-4x+4)} ={-3x^3+6x^2-7x+2}/{(x^2-x+1)^2*√(3x^2-4x+4)}
628:132人目の素数さん
20/01/14 23:48:24.90 Kk48Gk/U.net
>>593
ここにはプレーンテキストでしか貼れないけど
Wolframのサイトだと√とか読みやすいから、Wolframのサイトで解読した方が良いのではと思う。
629:132人目の素数さん
20/01/15 00:19:46.03 jJaC6yxD.net
>>559
へぇー
すごいなあ
630:イナ ◆/7jUdUKiSM
20/01/15 02:17:16 eGPXp69J.net
前>>593
>>595
√f'gのほうはあってる。
√fg'のほうに疑問が2つある。
1つは√fとgに分け、gを1/√の形ととらえたことで、前微分×後+前×後微分にしたはずなのに真ん中の符号が-であること。
もう1つはgが(-1)乗で、微分すると-2倍の-2乗になるはずなのに(-2)が掛かってないこと。
計算結果が違うのに、いくら簡単な数字を代入しようとも、値が異なるのは当然です。
631:イナ ◆/7jUdUKiSM
20/01/15 02:29:33 eGPXp69J.net
前>>598あ、わかったかも。
gを1/√の形ととらえたことで、微分したg'は(-1)が掛かって真ん中の符号が-になるんだ。
あとはなんで(-1)乗を微分して(-2)が掛かってないかだ。
前半は(3x-2)が因数だけども、後半は符号はともかくとして、因数の(2x-1)が(4x-2)になるんじゃないかと思った。
632:イナ ◆/7jUdUKiSM
20/01/15 02:35:25 eGPXp69J.net
前>>599
>>598後半もわかった。
-1乗を微分すると-1が掛かって-2乗になるわ。
x^3を微分すると3x^2だわ。
解決を見ました。
633:132人目の素数さん
20/01/15 02:44:36 tC7vGnci.net
三角形付け足して作っていったらいいんじゃねと思ったけど
違うわな
634:132人目の素数さん
20/01/15 04:11:02.51 V75Cib9r.net
円に内接する凸n角形(n≧4)を、1つの頂点からn-3本の対角線を引くことによって、n-2個の三角形に分割する。
このとき、各三角形の内接円の半径の和は分割の仕方に関係なく一定であることを示せ。
635:イナ
20/01/15 06:00:30.78 eGPXp69J.net
ポッと入れてパッと出るもんじゃねえんだな、
~∩∩前>>600 ∩∩
(-.-))カルダノ (`) )
[ ̄]_) って。 U⌒U、
 ̄ ̄]/\___∩∩ノ(γ)
__/\/,,(`.`))⌒゙,|
 ̄ ̄\/彡`-`ミυ`υυ|
 ̄ ̄|\_U⌒U、___/| |
□ | ∥~U~U~ ̄∥ | /
__| ∥ □ □ ∥ |/
_____`∥_______∥/
 ̄ ̄ ̄ ̄ ̄ ̄ ̄ ̄
636:132人目の素数さん
20/01/15 08:21:07.91 0Yai4QdM.net
(√f)/g の最大値を求めたいなら、 f/gg の最大値を求めて平方根とればいいんぢゃね?
0 = 2{(√f)/g}(d/dx){(√f)/g} = (d/dx)(f/gg),
637:132人目の素数さん
20/01/15 12:39:31.69 pa0Zl04A.net
>>602
四角形の場合に示せば十分である。
半径を1の円に四角形ABCDが内接するとしてAB,BC,CD,DAの円周角をx,y,z,wとする。
△ABCの内接円の半径は
4sin(x/2)sin(y/2)sin((π-x-y)/2)=cosx+cosy-1-cos(x+y)。
よって△ABCと△CDAの内接円の半径の和はx+y+z+w=πによりcosx+cosy+cosz+cosw-2。
内接円の半径の公式
URLリンク(mathtrain.jp)
638:132人目の素数さん
20/01/15 13:07:30.27 Y/R0zC7W.net
そうか、書いた後気づいたけど各辺の円周角をxiとすると
Σ内接円の半径 / 外接円の半径 = Σcos xi-n+2
になるのか。
美しいな。
639:132人目の素数さん
20/01/15 15:05:45 l229ykjv.net
>>602
いつものようにシミュレーションプログラムで確認してみた
"
円に内接する凸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系の剰余0のとき 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)
640:132人目の素数さん
20/01/15 15:06:52 l229ykjv.net
百角形での結果(まぁ、当然といえば当然。
> sim(100)
[1] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[12] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[23] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[34] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[45] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[56] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[67] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[78] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[89] 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977 1.914977
[100] 1.914977
>
641:132人目の素数さん
20/01/15 15:27:20 l229ykjv.net
何の役にたつかはわからんが、
せっかくシミュレーションプログラムを作ったので正方形から500角形までの分割三角内接円半径総和をグラフにしてみた。
URLリンク(i.imgur.com)
ある値に向かって収束しているみたいだな。
これは答のある問題ではありませんので、実務家の俺に解説や理論を求めてはなりませんwww
642:132人目の素数さん
20/01/15 15:39:41 l229ykjv.net
>>609
10万角形だと
> Soric(100000)
[1] 1.999902
2に収束するようにみえるなぁ。
643:132人目の素数さん
20/01/15 15:41:17 GJSc+40l.net
何コレ?
内接円の半径の総和が、nのみによって決まるわけないじゃん?
644:132人目の素数さん
20/01/15 16:40:13.62 KbwiaZJW.net
>>602
問題がおかしくないか?
>> 円に内接する凸n角形(n≧4)を、1つの頂点からn-3本の対角線を引くことによって、n-2個の三角形に分割する。
>> このとき、各三角形の内接円の半径の和は分割の仕方に関係なく一定であることを示せ。
「分割の仕方に関係なく一定」と書かれているが、
「1つの頂点からn-3本の対角線を引くことによって、n-2個の三角形に分割する」という操作は、ユニークに定まる。
「分割の仕方」は関与し得ない。
従って、本来の問題の意図は、「円に内接する凸n角形(n≧4)を、n-3本の交わらない対角線を引いてn-2個の三角形に分割する」
ではないのか? それならば、「分割の仕方」が関与し得る。
カタラン数C[n-2]=(2n-4)!/{(n-2)!(n-1)!}通りの分割方法があるが、そのどの分割であっても、
一定値になる、ということの証明が求められているのでは?
いずれにしろ、>>605 でほとんど解決(もう一言二言、補足があった方がよいとは思う)していることには、間違いない。
645:132人目の素数さん
20/01/15 16:52:44 l229ykjv.net
>>611
n→∞で2に収束する予感
646:132人目の素数さん
20/01/15 16:54:32 l229ykjv.net
>>612
どの頂点をつかって分割しても内接円半径の総和は一定という意味と理解したんだが。
647:132人目の素数さん
20/01/15 16:57:42 k/Eco663.net
>>614
そんなわけないやん?
n=3で考えたらわかる。
半径1の円に内接する三角形の内接円の半径が一定なわけないでしょ?
648:132人目の素数さん
20/01/15 17:08:10 l229ykjv.net
>>615
ある一つの多角形について、どの頂点を使って分割してもその多角形については内接半径の総和は一定ということじゃないの?
649:132人目の素数さん
20/01/15 17:15:30 k/Eco663.net
>>616
既に指摘されてる通り問題文には誤りがあるけど、エスパーすると
半径の総和は多角形のみによって決まり、分割の仕方にはよらない。
だと思う。
もちろん多角形の取り
650:方には依存し、nだけで決まるなどという事があるハズはない。
651:132人目の素数さん
20/01/15 18:35:33 l229ykjv.net
>>617
一様分布で乱数発生させたから、nを無限大にしていくと収束するように見えたんだろうな。
652:132人目の素数さん
20/01/15 20:06:14 l229ykjv.net
>>605
すいません。
>四角形の場合に示せば十分である。
というのは何故なんでしょうか?
653:132人目の素数さん
20/01/15 21:14:20 VFF1h4ZA.net
>>619
まぁもうすでに>>606で一歩先ゆく方針が出てるから今更になるけど。
四角形の場合に示せたとする。
以下を示せばよい。
頂点Aを任意に選ぶとき、任意の分割の時の値は全ての対角線の一端がAである分割の時の値に等しい。
n=4では示せている。
n<Nで示せたとしてn=Nのときを考える。
分割に用いた対角線の中にABが入っているときはABで分けた各々について帰納法の仮定を用いてよい。
そうでないときを考える。
分割の中で用いられている対角線CDをとる。
ただしCDで分けた二つの多角形のうちAを含む側は四角形以上とする。
(n≧5のとき分割で3つ以上の三角形が出てくるからそのようなCDは必ずとれる)
CDで分けられた2つの多角形に対して帰納法の仮定を用いることにより、CDで分けられた内、Aを含む側の分割の対角線は全て一端がAとしてよい。
その中にはAを一端とする対角線が少なくとも一つは出てくるのですでに示された場合に帰着できた。□
654:132人目の素数さん
20/01/15 22:24:13 VFF1h4ZA.net
>>606の方針による別解。
多角形の各辺の円周角をxi、外接円の半径を1とするとき
Σ内接円の半径=Σcos xi - n +2
を帰納法で示す。
n=3のとき各辺の円周角とはすなわち対角であり、この場合
URLリンク(mathtrain.jp)
などにより
r=4sin(x1/2)sin(x2/2)sin(x3/2)
であるから、積和公式を2回用いて主張を得る。
n<Nで示せたとしてn=Nとする。
与えられた分割に対して、その中で用いられている対角線ABをとる。
ABによって多角形がk角形とl角形に分かれたとする。
k角形の方に出てくる円周角をABに対するyとそれ以外のyp、l角形の方に出てくる円周角をABに対するzとそれ以外のzqとする。
帰納法の仮定により
Σk角形の側の半径の和=cos y + Σcos yp - k +2,
Σl角形の側の半径の和=cos z + Σcos zq - l + 2。
辺々足してk+l=n+2とy+z=πを用いて主張を得る。□
こちらの方が良い気がする。
655:イナ
20/01/15 23:44:53.30 eGPXp69J.net
前>>603
>>602
四角形ABCDに対角線ACを引き、
△ABCの内接円の半径をr,△CDAの内接円の半径をRとすると、
四角形ABCDの面積=(AB+BC)r/2+(CD+DA)R/2
四角形ABCDに対角線BDを引き、
△BCDの内接円の半径をя,△DABの内接円の半径を㌃とすると、
四角形ABCDの面積=(BC+CD)я/2+(DA+AB)㌃/2
BCとこれに垂直な半径rとの接点をE,CDとこれに垂直な半径Rとの接点をFとすると、
四角形ABCD=AC(r+R)+BEr+FDR
CDとこれに垂直な半径Rとの接点をG,DAとこれに垂直な半径㌃との接点Hとすると、
四角形ABCD=BD(я+㌃)+CDя+DA㌃
お膳立ては完璧。
あとはうまく消去して、
r+R=я+㌃を示す。
656:132人目の素数さん
20/01/16 05:35:20.59 aOBtX3Rl.net
3本の平行線が与えられているとき、各平行線上に各頂点があるような正三角形を定規とコンパスを使って作図する方法を説明せよ。
657:132人目の素数さん
20/01/16 06:34:28.03 07zfQGLo.net
>>620
解説ありがとうございました。
658:132人目の素数さん
20/01/16 08:04:16.33 552XyIPx.net
密度が一定の球形の惑星がある。
この惑星の表面点Pでの重力を最大にするために
体積と密度一定の条件で惑星の形を変形する。
点Pでの重力は球形のときの最大何倍になるか。
またこのときの惑星の形はどうなるか。
659:132人目の素数さん
20/01/16 09:03:43.37 ALo92Jf3.net
>>623
平行線を順にl,m,nとし、巾の比をa:bとする。
正三角形ABCを作図し、辺ABをa:bに内分する点Dを作図する。
(直線AB外の点XとAX:XY=a:bなるYをXが直線ABと重ならないようにと�
660:閨ABX//DYとなる点を作図すれば良い。) m上にC'D'をCD=C'D'となるようにとりB'をnの側に△B'C'D'が△BCDと合同になるようにとる。 直線B'D'とl,nの交点をA",B"とすればA"とB"を一辺とする正三角形の一方が求めるもの。
661:132人目の素数さん
20/01/16 09:21:53 YUeZdYQq.net
>>626
>正三角形ABCを作図し、辺ABをa:bに内分する点Dを作図する。
>(直線AB外の点XとAX:XY=a:bなるYをXが直線ABと重ならないようにとり、BX//DYとなる点を作図すれば良い。)
平行線の垂線の垂直二等分線で正三角形作れば良いじゃん
662:132人目の素数さん
20/01/16 09:39:28 PPdqHAdt.net
>>627
だな。
l,n上にABを好きにとればよかった。
663:132人目の素数さん
20/01/16 09:40:02 YUeZdYQq.net
定規とコンパスである2点間の長さを半径とする円をその2点とは別の点を中心に描くのは面倒くさい(幾何学原論ではコンパスは中心からある点までの距離を半径とする円を描けるだけで長さを保存して中心を移動させてはいけない)
664:132人目の素数さん
20/01/16 09:43:18 552XyIPx.net
>>623
中心の平行線をl、他の平行線をm,nとする。
l上に点Oを取り、O点を通りlとのなす角が60°,120°の直線p,qを引く。
直線pとmの交点をA、直線qとnの交点をBとすれば、AB=AC=BCとなるl上の点Cが取れる。
665:132人目の素数さん
20/01/16 09:45:23 2aq5oe6O.net
まぁもっと楽な方法があったら教えて下さい。
自分は解けたのでもう満足。
666:132人目の素数さん
20/01/16 12:17:23 /nA7kUIy.net
>>625
半円かな?
667:哀れな素人
20/01/16 12:50:07 7jomMr1V.net
>>602
とりあえず四角形の場合について初等幾何的に証明すると-
四辺をa、b、c、d、対角線をe、f、内心半径をr1、r2、r3、r4とする。
eと、Aと二つの内心とCを結ぶ線で分けられる四つの図形の面積を
S1、S2、S3、S4とすると、S2+S3:S1+S4=2e:a+b+c+d
同様にfと、BとDを結ぶ線で分けられる四つの図形の面積を
S5、S6、S7、S8とすると、S6+S7:S5+S8=2f:a+b+c+d
ゆえにS2+S3:S6+S7=e:f
ところでS2+S3=e(r1+r2)、S6+S7=f(r3+r4)
ゆえにe(r1+r2):f(r3+r4)=e:f
ゆえにr1+r2=r3+r4
668:132人目の素数さん
20/01/16 13:37:53 /nA7kUIy.net
びっくりするくらいわからん
669:イナ ◆/7jUdUKiSM
20/01/16 14:02:50 izvR1SI9.net
前>>622で言いたかったことをn≧4に拡張すると、
>>633そういうことです。
670:イナ ◆/7jUdUKiSM
20/01/16 14:09:07 izvR1SI9.net
前>>635
拡張はしてないか。
四角形ABCDそのものか。
n=4のときどちらの対角線で分割しても内接円の半径の和が等しいことを示したかった。
671:132人目の素数さん
20/01/16 16:08:24 2LVdbF4s.net
>>633
S1、S2、S3、S4とすると、S2+S3:S1+S4=2e:a+b+c+d
これ内接円の半径等しいという仮定ないとダメなのでは?
672:哀れな素人
20/01/16 17:28:33.93 7jomMr1V.net
>>637
なるほど、うっかりしていた(笑
比の計算で、単純に足したり引いたりしてはいけない、
という基本的なことを忘れていた。
673:イナ ◆/7jUdUKiSM
20/01/16 17:46:42 izvR1SI9.net
前>>636 なんだ6>>33は違うのか。
>>622まではいいと思う。
674:132人目の素数さん
20/01/16 22:56:53 aOBtX3Rl.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のこの行動は最善なのか?
675:132人目の素数さん
20/01/16 23:35:32.99 l4Y8hLMv.net
>>625
等周問題の変形版だね
x^2+y^2+z^2-h^(4/3)*z^(2/3)=0 (hは高さ)
URLリンク(imgur.com)
676:132人目の素数さん
20/01/16 23:37:04.96 l4Y8hLMv.net
>>625
重力は3*5^(1/3)/5=1.026倍
677:イナ
20/01/16 23:47:02.74 izvR1SI9.net
前>>639
>>640
A「(Cに)Bの野郎、最善尽くして右へ右へと逸れていきよる思うんよ」
B「(Aに)π倍の速さで走ればいいじゃねえか」
A「いや、それだと逃げられる。おめぇさんがまっすぐ走ってるあいだに、円周という拘束がある俺は、まっすぐ走ってるおめぇさんに対してr走ってから直角に曲がってπr走るようなもんだ」
Aの心の声「ピタゴラスの定理より、
√{r^2+(πr)^2}=r√(1+π^2)
つまりr(1+π)/r√(1+π^2)倍走らな追いつかん。
約分し、
(1+π)/√(1+π^2)
分母を有理化し」
∴(1+π)√(1+π^2)/(1+π^2)倍以上の速さで走ればいい。
678:132人目の素数さん
20/01/16 23:47:37.27 552XyIPx.net
>>641 >>642
正解です。
679:イナ ◆/7jUdUKiSM
20/01/16 23:55:12 izvR1SI9.net
前>>643
(1+π)√(1+π^2)/(1+π^2)
=1.25620499……
1.256205倍なら追いつく。
680:132人目の素数さん
20/01/17 00:11:50 m+4frRir.net
>>641
kwsk
681:132人目の素数さん
20/01/17 00:23:11 gpQhzbqr.net
>>640
某サイトに答えあるの見つけた。
コピペするのもあれなのでROMします。
682:132人目の素数さん
20/01/17 00:26:36 pPPuMatU.net
>>646
点Pでの重力を体積一定の下で最大化するような断面の形状を変分法で求める
惑星は軸対称(軸対称でないとすると点Pから少し動いたところにより重心に近い点がある)だから断面の形状だけ考えれば良いのがポイント
↓のpart II の 4 に等長問題のオイラーラグランジュの方程式を求める方法が書いてある
URLリンク(citeseerx.ist.psu.edu)
683:132人目の素数さん
20/01/17 00:30:13 wItShquN.net
ちなみに速度比が2倍では捕まえられない。
Bの速度は1とする。
速度比が2倍ならBは半径1/2の地点まではAの角速度を上回っているのでここまでBはAとの偏角の差が180°の状態をキープできる。
この地点から円周の最短地点までBが要する時間は1。
Aが要する時間はπなのでBの逃げ切り成功になる。
684:132人目の素数さん
20/01/17 00:45:02 XDQcvfLQ.net
>>648
この問題に最大解がある事の証明はどうすればいいのでしょう?
変分法使う以上は別に解の存在示しとかないとダメだとおもうんですが。
685:132人目の素数さん
20/01/17 00:54:20.56 pPPuMatU.net
>>650
リンク先にも書いてあるけど,厳密には第二変分を考えなきゃいけない
686:132人目の素数さん
20/01/17 00:56:48.70 Vq8XRGW9.net
>>659
そんなのがあるんですか?
なんか言葉のイメージだとHessianみたいなものですか?
でもそれだと結局極大である事の保証しか得られないのでは?
687:132人目の素数さん
20/01/17 01:30:08 pPPuMatU.net
>>652
具体的に求めてみると分かるけど汎関数を与える関数(被積分関数)が凸関数だから結局極大値が最大値になる
688:132人目の素数さん
20/01/17 01:46:52.25
689: ID:0cuf6Wtp.net
690:132人目の素数さん
20/01/17 01:53:19 pufwLXqa.net
>>653
そうなんですか。
つまりf、gがともにV(f)≦V0、V(g)≦V0を満たすとき任意の0<t<1に対してV(tf+(1-t)g)、S(tf+(1-t)g)について凸不等式が使えるという事でしょうか?
体積の方は行けそうですね。
Sの方はまだ手動かしてないのでやってみます。
ありがとうございました。
691:132人目の素数さん
20/01/17 02:11:57 B81I4URB.net
今チェックしてみたら行けそうですね。
素晴らしい。
692:132人目の素数さん
20/01/17 08:18:19 F8/vjJrP.net
別スレより
1歩で1段または2段のいずれかで階段を昇るとき、1歩で2段昇ることは連続しないものとする。15段の階段を昇る昇り方は何通りあるかを求めよ。
693:132人目の素数さん
20/01/17 09:06:38 Ax5k6DPo.net
2段昇りの回数をnとすると、連続できないので最低限
間に1段昇りがはいることから、2n+(n-1)≦15 よって、
0≦n≦5
n=0の場合、一通り
n≠0の場合、15-2n個の1段昇りの隙間と両端に2段昇
りが入りうるので、その組み合わせは C(15-2n +1, n)
よって、全部で1+Σ(n=1~5)C(16-2n,n)
694:132人目の素数さん
20/01/17 10:15:09 CPtT5WKR.net
n段登る方法の数をAnとする。
A1=1, A2=2,A(n+2)=A(n+1)+An
によりFをフィボナッチ数列とすればAn=Fn+F(n-1)。
695:132人目の素数さん
20/01/17 10:16:22 CPtT5WKR.net
あ、一歩で二段連続なしかorz
696:132人目の素数さん
20/01/17 10:20:39 F8/vjJrP.net
プログラム組んで数えさせたら49になったんだけど合ってる?
[[1]]
[1] 1 2 1 2 1 2 1 2 1 2
から
[[49]]
[1] 2 2 2 2 2 2 2 1
697:132人目の素数さん
20/01/17 10:25:37 CPtT5WKR.net
となると
A1=1,A2=1,A3=2,A(n+3)=A(n+1)+Anですな。
steps=map head $ iterate (\[x,y,z]->[y,z,x+y]) [1,1,2]
main = do
print $ take 16 steps
[1,1,2,2,3,4,5,7,9,12,16,21,28,37,49,65]
698:132人目の素数さん
20/01/17 11:22:45 F8/vjJrP.net
同じ結果
> sapply(1:16,sim)
[1] 1 1 2 2 3 4 5 7 9 12 16 21 28 37 49 65
699:132人目の素数さん
20/01/17 14:08:51 F8/vjJrP.net
>>657
2歩を連続させてはいけないというルールにしてみた。
1歩で1段または2段のいずれかで階段を昇るとき、2歩で続けて昇ることは連続しないものとする。
15段の階段を昇る昇り方は何通りあるかを求めよ。
パソコンに数えてもらったら、次のようになったけど、合ってる?
> sapply(1:16,function(x) sim(x,2, print=FALSE))
[1] 1 2 3 4 6 9 13 19 28 41 60 88 129 189 277 406
700:132人目の素数さん
20/01/17 14:11:58 2HphDdoJ.net
最後のステップが一段で、合計n段になる登り方をA[n]、
最後のステップが二段で、合計n段になる登り方をB[n]とすると、
求めるものは、C[n]=A[n]+B[n]。
それぞれの漸化式は
A[n]=A[n-1]+B[n-1]
B[n]=A[n-2]
つまり、A[n]=A[n-1]+A[n-3]
初期条件として、A[1]=1,A[2]=1,A[3]=2 を用いて、n=1から、列挙すると、
1,1,2,3,4,6,9,13,19,28,41,60,88,129,189
C[15]=A[15]+B[15]=189+88=277
これは、>>658さん、>>664さんの結果とも一致
701:132人目の素数さん
20/01/17 15:34:25 k0YCxJCN.net
受験数学でクソ頻出の3項間の項の数が増えてるだけですがな。
これが即答できないのは受験レベルの3項間の勉強の時に解き方だけ覚えて理屈がわかってなかった証拠。
702:132人目の素数さん
20/01/17 16:34:49 ElKm6FvG.net
2007年の京大の問題か
703:132人目の素数さん
20/01/17 17:04:55 F8/vjJrP.net
>>667
検索したらそうみたいですね。
704:俺は、ニュー速+に書いてあったを転載したんだけど。 【勉強】数をかぞえられない学生たち ありもしない「公式」に頼り…算数教育の珍現象 https://asahi.5ch.net/test/read.cgi/newsplus/1579164227/
705:132人目の素数さん
20/01/17 17:50:57 F8/vjJrP.net
>>640
Bは常にAと正反対の方向に向かって逃げてればいいような気がするな。
こんな感じ URLリンク(imgur.com)
706:132人目の素数さん
20/01/17 18:54:31 F8/vjJrP.net
>>669
Bが必勝のような気がしてきたな。
707:イナ ◆/7jUdUKiSM
20/01/17 19:56:54 3S1g3K0T.net
前>>645再考。
>>669が見えないからなんとも言えないんだが、
1.256205倍は直感的すぎるし、πrで追いつかないものがπr+rぐらいで追いつくとも思えない。
もっと基本的な話、泥警といっしょで、泥棒Bは警察Aに対して90°の方向に逃げる。仮にBが90°以上の方向に逃げようとした場合、Aは逆方向に先回りする。ていうか最善を尽くすという題意だからその時点でそうする。
BはAのせいで自然と円軌道をとることになり、
スタートは12時の方向にいるAを背後に見て6時の方向に進みだすが、Aが1時2時と時計回りに進むに従って、Bも時計回りに最短距離で外周に達しようとする。警官Aが泥棒Bを捕まえる外周の位置は6時から9時のあいだにある。
6時のとき、π倍。
7時のとき、Aが(7/6)πr走るあいだにBは9時の外周を中心とした6分円、
2πr(60°/360°)=πr/3を走る。
Aの速さはBの速さの、 (7/6)πr÷πr/3=3.55
9時の外き、3倍だから、
ABA外周の7./2走っったとき追Bにいつくたた5倍
文字化けで書き直し。
∴AがBの3.5倍の速さで外周の7/12走ったとき、AがBに追いつく。
708:132人目の素数さん
20/01/17 20:11:26 F8/vjJrP.net
>>671
画像のURL
URLリンク(i.imgur.com)
709:132人目の素数さん
20/01/17 20:37:41.55 F8/vjJrP.net
3.5倍の速度でも常にAと反対向きに進路変更すれば逃げ切れる気がするんだが。
URLリンク(i.imgur.com)
710:132人目の素数さん
20/01/17 20:49:50.57 F8/vjJrP.net
>>673
これ撤回。 Bが加速している。
711:イナ ◆/7jUdUKiSM
20/01/17 21:21:42 3S1g3K0T.net
前>>671
計算上、泥棒B(ルパン)が警官A(銭形)の2/7以上の速さで走れば、わずかにAに捕まる前に外周に逃げのびることができる。
逆走についてはAはBの逆走にただちに反応して逆走するルールだからむだな抵抗。袋のねずみ状態と考えられる。ただし穴あき袋、ある程度の速さで逃げれば外周に達する。
不二子「銭形警部の2/7ぐらいのスピード、ルパンならわけないわ」
ただ3.5倍という速さはそうとう速いので、感覚的に逃げきれそうなのはAのBに対する身体的優位性が実感できないからだと思う。
あと、3.5倍前後が少し怪しいといえば怪しい。Bが逃げる円軌道の中心をAの円軌道の中心と9時方向の端点を結ぶ直線上にとって半径rの6分円の外周(中心角60°の扇形の孤)を7時方向の端点まで逃げると考えた。
微分して極値を求めたわけじゃないから確証が持てないけど、確実にπ倍は超えてる。
712:132人目の素数さん
20/01/17 21:43:57 k0YCxJCN.net
>>673
BがAの2/7倍なら逃げられるよ。
Aの速度が1,Bが2/7として半径2/7の地点まではBの角速度が勝るのでこの地点までBはAと偏角の差をπに保てる。
そこからBは円周の最短地点まで距離は5/7,所用時間は5/2。
一方Aはその地点まで道のりπ、所用時間πなのでBの勝ち。
713:イナ
20/01/17 23:29:10.25 3S1g3K0T.net
前>>675
AがBを捕まえる外周の位置は6時方向から9時方向のあいだにあると考えられる。Bが逆走したらAも逆走するから、一周することはない。
6時のとき、
Aが外周の半分を走るあいだにBはまっすぐr逃げる。
2πr/2÷r=π(倍)
7時のとき、
Aが(7/6)πr追うあいだにBは9時の外周を中心とした6分円、
2πr(60°/360°)=πr/3を走る。
Aの速さはBの速さの、
(7/6)πr÷πr/3=3.5(倍)
9時のとき、Aの速さはBの速さの、
3πr/2÷2π(r/2)(1/2)=3(倍)
8時のとき、Aの速さはBの速さの、
8πr/6÷2π(r/√3)(1/3)
=(4πr/3)(3√3/2πr)
=2√3
=3.4641016……
7時と8時のあいだ方向が怪しいといえば怪しい。
そのとき、AはBの3.5倍を超える速さ以上の速さで走らなければならない可能性があるといえばある。
714:132人目の素数さん
20/01/17 23:31:45.00 8XEo1F0J.net
A[n] = A[n-1] + A[n-3] より
A[n] = c・α^n + d・|β|^n・cos(nω+θ),
ここに
α、β、β~ は特性値。 (特性方程式 t^3 - t^2 -1 = 0 の根)
α = {1 + [(29-3√93)/2]^(1/3) + [(29+3√93)/2]^(1/3)}/3
= 1.465571231876768
β = |β| e^(iω),
β~ = |β| e^(-iω),
|β| = 1/(√α) = 0.826031357654187
Re{β} = -(α-1)/2 = -0.232785615938384
cos(ω) = Re{β}/|β| = -0.281812081080629
ω = arg(β) = 1.856478541471303
また、
c^3 - c^2 +(9/31)c -(1/31) = 0,
c = 1/3 + {[4(√31-√27)]^(1/3) + [4(√31 +√27)]^(1/3)}/(3√31)
= 0.6114919919508125
d = (1-c)/cosθ,
715:132人目の素数さん
20/01/18 05:44:37.53 2fzpXJcR.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人嘘つきがいる」
716:132人目の素数さん
20/01/18 08:24:40.93 OC/W80tJ.net
H/L A B C D E F G H
0/8 × × × × ◯ ◯ ◯ ◯
1/7 ◯ × × × ◯ ◯ ◯ ◯
2/6 ◯ ◯ × × ◯ ◯ ◯ ◯
3/5 ◯ ◯ ◯ × ◯ ◯ ◯ ◯
4/4 ◯ ◯ ◯ ◯ ◯ ◯ ◯ ◯
5/3 ◯ ◯ ◯ ◯ ◯ ◯ ◯ ×
6/2 ◯ ◯ ◯ ◯ ◯ ◯ × ×
7/1 ◯ ◯ ◯ ◯ ◯ × × ×
8/0 ◯ ◯ ◯ ◯ × × × ×
GとHが嘘つき。
717:132人目の素数さん
20/01/18 08:49:58.52 aX+4W30D.net
いつもの総当りプログラムでも最後の2人が嘘つき
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
718:132人目の素数さん
20/01/18 08:57:10.68 2fzpXJcR.net
正解です。
719:132人目の素数さん
20/01/18 09:39:19 aX+4W30D.net
>>679
こういう問題にすると
A「8人の中に、1人正直者がいる」
B「8人の中に、2人正直者がいる」
C「8人の中に、3人正直者がいる」
D「8人の中に、4人正直者がいる」
E「8人の中に、1人嘘つきがいる」
F「8人の中に、2人嘘つきがいる」
G「8人の中に、3人嘘つきがいる」
H「8人の中に、4人嘘つきがいる」
答が変わるね。
720:132人目の素数さん
20/01/18 09:54:54 aX+4W30D.net
こういうバリエーションも可能。
A「8人の中に、1人嘘つきがいる」
B「8人の中に、2人嘘つきがいる」
C「8人の中に、3人嘘つきがいる」
D「8人の中に、4人嘘つきがいる」
E「8人の中に、5人嘘つきがいる」
F「8人の中に、6人嘘つきがいる」
G「8人の中に、7人嘘つきがいる」
H「8人の中に、8人嘘つきがいる」
721:132人目の素数さん
20/01/18 09:58:04 PXyx+OwH.net
>>678
Re{β} = -(1/2)|β|^4 = -α^(-2)/2,
cos(ω) = -(1/2)|β|^3,
722:132人目の素数さん
20/01/18 10:57:32.36 aX+4W30D.net
>>679
これって嘘つきの証言も正しいという前提?
723:132人目の素数さん
20/01/18 11:27:02.31 aX+4W30D.net
>>679
Aが嘘つきの場合は
「8人の中に、少なくとも1人正直者がいる」は嘘の証言で
正直者は誰もいないとして計算?
724:132人目の素数さん
20/01/18 11:43:01.17 EsuXaTYh.net
この問題で問題文の意味の取りようなんて議論の余地ないだろ?
725:132人目の素数さん
20/01/18 11:50:29.99 PXyx+OwH.net
>>669
>>670
円の半径を1, Bの速さを1, Aの速さをaとする。
Aがθ(t)で追ったとする。
|dθ/dt| ≦ a,
最初、Bは Aと逆方向 π+θ(t) を保ちながら速さ1で逃げる。
中心~Bの距離を r(t) とする。
dr/dt = √{1 - rr(dθ/dt)^2} ≧ √{1-(ar)^2},
r ≧ sin(at)/a,
時刻 π/2a までに r=1/a に到達する。
次に、Bは円周Cまで直進する。所要時間: (a-1)/a,
AがCに到着するまでの所要時間: π/a
a<π+1 ならば逃げ切れる。
r=1/a で直角に曲がらずに丸く曲がった方が短くなる。
aはもっと大きい可能性・・・
726:132人目の素数さん
20/01/18 12:09:00.92 aX+4W30D.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も嘘つきである」
727:658
20/01/18 12:56:09.45 E9RqLHbc.net
>>665
漸化式かっこいい!でも、一般項は求まらないか、、、。
>658まで書いたところで、急にウンチしたくなって、出勤時間も
迫ってたんで、値まで計算せずに書き込みして出かけちゃいまし
たが、検算ありがとう。
728:132人目の素数さん
20/01/18 12:56:14.78 EsuXaTYh.net
全員正直はCと矛盾するのであり得ない。
よってE=H。G=L。
CによりB又はCが嘘つきなのでG以外の嘘つきがいる。
よってF=H、H=L。
E,Fが正直なのでB=H、C=L、D=H。
ここまでで4人以上正直がいるのでA=L。
729:132人目の素数さん
20/01/18 15:05:40 zvEtuxN+.net
>>691
>一般項は求まらないか、、、。
3項漸化式で特性方程式が3次
その解を使えば?
730:132人目の素数さん
20/01/18 15:08:21 kUGqA2Ne.net
嘘つき問題を論理式で解いたらどうなるんだろ
731:132人目の素数さん
20/01/18 15:34:40 EsuXaTYh.net
普通に命題論理の問題になるけど例えば>>679の問題の
B「8人の中に、少なくとも2人正直者がいる」
を命題論理の式にすると
B=(A∩B)∪(A∩C)‥∪(G∩H)
になって手計算では無理な話になってしまう。
計算機使えれば一撃だけど。
732:132人目の素数さん
20/01/18 15:46:56 zvEtuxN+.net
>>694
自己言及があるが?
733:132人目の素数さん
20/01/18 16:54:38.67 dWbEz1QY.net
>>690
プログラム組んで解いたら
> TE[apply(TE,1,g),]
[1] 0 1 0 1 1 1 0 0
嘘つきは A C G H
734:132人目の素数さん
20/01/18 17:24:05 dWbEz1QY.net
プログラム解を前提にすれば、ファジーな嘘つきを入れてこんな問題もできる。
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も嘘つきである」
735:132人目の素数さん
20/01/18 17:27:51 dWbEz1QY.net
どこかの中学入試の問題らしい
A、B、C、D、Eの5人のうち2人は常に本当のことを言う正直者です。あとの3人は嘘つきですが、その発言内容は本当のときもあります。
彼らに誰が嘘つきか尋ねたところ次のように答えました。
A 「BとEは嘘つきではない」
B 「Cは嘘つきだ」
C 「Dは嘘つきだ」
D 「Eは嘘つきだ」
E 「BとCは嘘つきだ」
さて、正直者は誰と誰でしょうか?
736:132人目の素数さん
20/01/18 18:21:23 EsuXaTYh.net
>>699
A=HとするとE=Hとなって矛盾するのでA=L。
E=Hとすると残り1人の正直者はDしかあり得ないが、するとするとE=Lとなって矛盾するのでE=L。
B=C=HはBの証言に矛盾するのであり得ないからE=H。
するとD=L、C=Hとわかる。
737:132人目の素数さん
20/01/18 18:22:46 EsuXaTYh.net
>>698
この手の問題計算機使って解かないと解けないのでは単なる計算機の演習にしかならない。
738:132人目の素数さん
20/01/18 21:12:52 zvEtuxN+.net
>>701
試行錯誤で解けるが
739:132人目の素数さん
20/01/18 21:17:38 EsuXaTYh.net
>>702
解答プレーズ
740:イナ ◆/7jUdUKiSM
20/01/18 21:20:38 hIqEloY9.net
前>>679挑戦者暫定1位。
>>640問題�
741:B 逃走ルートをなるべく最短にしたいBは、7時半の方向に逃げるはずだ。 AがBを追った距離は、 2πr{(7+30/60)/12} =5πr/4 Bが逃げた距離は、半径が直角になりこれはかなり最短だと思うんだけど半径r/√2の四分円孤、 2π(r/√2)(1/4)=πr/2√2 Aの速さはBの速さの、 5πr/4÷πr/2√2=5√2/2 =3.53553391……(倍) まさかの3.5倍超え。さすがに4倍は超えないと思うけど、まだわからんな。
742:イナ ◆/7jUdUKiSM
20/01/18 21:27:02 hIqEloY9.net
前>>704
前々>>677
アンカー訂正。
m(__)m
743:132人目の素数さん
20/01/18 21:31:04 dWbEz1QY.net
>>701
2^8=256だから、手作業でやれなくもないとは思うけど
プログラムを組む方が早い。
256通りを以下のような関数でTRUEを返すものを選ぶだけ。まさに無思考www
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), # Cが正直者で証言が正しい | Cが嘘つきで証言が嘘
(x[4]==1 & x[3]==0) | (x[4]==0 & x[3]==1), # Dが正直者で証言が正しい | Dが嘘つきで証言が嘘
(x[5]==1 & L>=1) | (x[5]==0 & L<1 ), # Eが正直者で証言が正しい | Eが嘘つきで証言が嘘
(x[6]==1 & L>=2) | (x[6]==0 ), # Fが正直者で証言が正しい | Fが嘘つき
(x[7]==1 & x[5]==0) | (x[7]==0 ), # Gが正直者で証言が正しい | Gが嘘つき
(x[8]==1 & x[1]==0 & x[6]==0) | (x[8]==0 ) ) ) # Hが正直者で証言が正しい | Hが嘘つき
可能な組み合わせ(1:正直者 0:嘘つきと
> TE[apply(TE,1,k),]
A B C D E F G H
[1,] 0 1 0 1 1 0 0 1
[2,] 0 1 0 1 1 1 0 0
どちらでも正直者なのは、B D E
744:132人目の素数さん
20/01/18 21:40:14 EsuXaTYh.net
あれ?
オレの計算結果と違う?
[True,False,False,True,True,False,False,False]
[True,False,False,True,False,True,False,False]
[True,False,False,True,False,False,True,False]
[True,False,False,False,True,True,False,False]
[True,False,False,False,False,True,True,False]
[False,False,False,True,True,True,False,False]
[False,False,False,True,True,False,False,True]
[False,False,False,True,False,True,True,False]
[False,False,False,True,False,False,True,True]
になった。
まぁこんな計算できてもどうでもいいけど。
745:イナ ◆/7jUdUKiSM
20/01/19 01:44:23 6fYcVEZg.net
前>>705挑戦者1位。
>>640問題。
Bが7時20分の方向に逃げたとき、
AがBを追った距離は、
2πr{(7+1/3)/12}
=11πr/9
Bが逃げた距離は、半径R,中心角80°の扇形円孤、
2πR(80°/360°)=4πR/9
Aの速さはBの速さの、
11πr/9÷4πR/9=11r/4R
=3.……(倍)
正十八角形の1辺の長さを外接円の半径rで表すことができれば、三角形の相似比からRをrで表せて、倍率が出る。と考えたがとりあえず値だけ出して、
Rsin40°=r/2より、
R=r/2sin40°
11r/4R=(11/2)sin40°
=3.53533185……(倍)
7時30分の短針の方向にわずかに及ばない。
円軌道で7時半の短針の方向に出るBにAが追いつくとき、AはBの3.53553391倍以上の速さが必要が最大で、これ以上の値はみつからない。
Bの逃走ルートとして螺旋のルートを考えると、先に求めた円軌道の最大値、
πr√2/4を保ったままコイルをひっぱるように螺旋状にのばすと、
Bがπr√2/4(>r)逃げるあいだにAが一周2πr走って追いつくことも可能で、
Aの速さはBの速さの、
2πr/(πr√2/4)=4√2
=5.65685……(倍)
というべらぼうな速さを出せるが、
かならずAが最善の走りに気づいて逆回りするからアウトだと思う。
746:132人目の素数さん
20/01/19 01:50:22 0RKATGn1.net
>>698
> 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も嘘つきである」
Dの発言からCとDの属性(正直者or嘘つき)は一致しないので、
Eの発言は真。ゆえにEは正直者。
C,Dのうち片方とEを合わせれば合計二人の正直者がいることになるから、
Bの発言は真。ゆえにBは正直者。
CとDの発言から、C嘘つきとD正直者が確定。
この時点で正直者が三人発覚しているが、
Aが正直者なら正直者が四人以上いることになり矛盾。ゆえにAは嘘つき。
Gの発言は偽であるからGは嘘つき。
Aが嘘つきであることから、Hの発言の真偽はFの属性(正or嘘)と一致。
Fの発言は真であるが、Fの属性については何もわからない。以上から、HとFについて
(H,F)=(正,正), (正,嘘), (嘘,嘘)
の可能性がある以外は、これまでの議論の通り全て確定している。
以上から、答えはB,D,E.
747:イナ
20/01/19 04:45:53.99 6fYcVEZg.net
前>>708念のため7時40分の方向にBが逃げたときを考える。
>>640問題。
Bが逃げる軌跡は、
底角40°,頂角100°の二等辺三角形を含む扇形の円孤になり、
Aの速さはBの速さの、
2πr(7+2/3)/12÷2π(r/2cos40°)(100°/360°)
=(23/3)/12÷5/18・2cos40°
=23・2cos40°/10
=(23/5)cos40°
=3.52380444……(倍)
<3.53553391……
あとは7時25,27,29分あたりがちょっと気になる。
7時台の速さ比(A/B)のグラフが見たい。
7時半の方向が最大だと思うけど。
748:132人目の素数さん
20/01/19 06:34:16.10 EEjXWdAC.net
>>709
レスありがとうございます。
(H,F)=(正,正)は成立しません。 Hの発言: H「AもFも嘘つきである」 からFが嘘つきになり矛盾します。
(H,F)=(嘘,嘘)も成立しません。 Aは必ず嘘をつく嘘つきと確定しているので正直者は3人ではありえません。
(H,F)=(嘘,嘘)ならば正直者がB,D,Eの3人になるのでAが真実を語っていることになり、矛盾します。
749:132人目の素数さん
20/01/19 07:04:20 EEjXWdAC.net
>>700
Aが正直者とすると正直者が3人になるのでAは嘘つき。
B,C,D,Eに二人の正直者がいる。
その組み合せは6通り。
相互に矛盾しないのはB,Dの組合せだけ。
750:イナ ◆/7jUdUKiSM
20/01/19 08:36:20 6fYcVEZg.net
前>>710
>>699
正直者はAとC。
Aが正直者とすると矛盾が生じてAは嘘つきとなったが、うまくいかない。
もう一度、嘘つきが本当のことを言うこともあるという広い心で題意を読みかえしてみた。
おのずとAとCが正直者だったとわかった。
751:132人目の素数さん
20/01/19 08:52:07 MNtmCELY.net
1cm刻みの数直線上の隣り合った3点に●がある(3つの●は区別がない)。
1秒ごとに●1つが他の●1つを飛び越す。
その場合、飛んだあとの●の位置は、飛び越された●からの距離がもとの位置からと同じになるようにする。
飛んだあとの位置にすでに他の●があるときは飛べない。
2021秒後に、●●●全てが開始位置に戻ることはできるか?
例、
●●●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ●●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ● ●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ● ●
╋╋╋╋╋╋╋╋╋╋╋
↓
:
↓
●●●
╋╋╋╋╋╋╋╋╋╋╋
752:132人目の素数さん
20/01/19 10:17:19 p/jcyrav.net
>>713
正解が投稿されてから誤答を繰り返すいつもの芸風 乙
753:132人目の素数さん
20/01/19 10:18:00 0RKATGn1.net
>>711
うわ、二つもアウトな場合取りこぼすとかぼーっとしてた
なるほど一度使った条件もまた後から効いてくる場合があったか
754:132人目の素数さん
20/01/19 10:22:05 mAiNEiBE.net
いくらなんでも2021がムリクリすぎん?
755:132人目の素数さん
20/01/19 10:31:11 flshczPE.net
>>703
2^8=256通り試すだけ
756:132人目の素数さん
20/01/19 10:36:09 Ld5ZhyuY.net
>>714
A B C
╋╋╋╋╋╋╋╋╋╋╋
↓
(2021回後)
↓
C A B
╋╋╋╋╋╋╋╋╋╋╋
は元に戻った事になりますか?
757:132人目の素数さん
20/01/19 10:37:14 gvXkKI2o.net
ズレたorz
でも疑問の意味はわかってもらえると信じて
758:132人目の素数さん
20/01/19 10:56:30 U/f9+Pty.net
とりあえず3個の石がそれぞれ全部元位置に戻るのは不可。
3個の位置を(a,b,c)としてAがBを飛び越すと(-a+2b,b,c)になる。
この変換を行列で表した時のdetは-1。
他も同様。
もし最初(0,1,2)から始めて(0,1,2)に戻ったとすると(n,n+1,n+2)から始めると(n,n+1,n+2)に戻り、よって変換は単位行列にならなければならない。
しかし一個の行列の行列式が-1なのでどのように2021個組み合わせてかけてもその行列式は-1にしかならない。
759:132人目の素数さん
20/01/19 11:28:59.40 Yx0eiL2o.net
あ、相互位置交換も不可かな?
対称性から(a,b,c)→(b,a,c)が不可である事を示せば十分。
(0,1,2)→(1,0,2)が不可である事を示せば十分だけど各変換はmod2で恒等写像なのでコレは明らか。
760:132人目の素数さん
20/01/19 11:49:47 MNtmCELY.net
>>719
●は区別がないので、なります。
761:132人目の素数さん
20/01/19 12:04:17 0RKATGn1.net
(0,1,2)→(0,3,2)→(0,3,4)→(6,3,4)→(2,3,4)→(2,3,0)→(2,1,0)
偶数回ではあるけど位置交換自体は可能みたい
762:イナ ◆/7jUdUKiSM
20/01/19 12:55:24 6fYcVEZg.net
前>>713
>>699
BとDが正直者。
>>715読みなおした。やっぱり矛盾が生じてAは嘘つきとなった。
Bが正直者とすると、
C嘘つき,D正直者,E嘘つきで矛盾しないと思う。
763:132人目の素数さん
20/01/19 12:58:23 EEjXWdAC.net
シミュレーションしていたら、こういうジャンプで元に戻る
> sim()
0 : 1 2 3
1 : 1 3 4
2 : -1 1 4
3 : -1 4 7
4 : -6 -1 7
5 : -1 4 7
6 : -1 1 4
7 : 1 3 4
8 : -1 1 4
9 : 1 3 4
10 : 1 2 3
最短だと
> sim()
0 : 1 2 3
1 : 0 1 3
2 : 1 2 3
2021回は奇数回だから 元に戻れないと思う。
根拠は直感のみw
764:132人目の素数さん
20/01/19 13:57:20 0RKATGn1.net
>>714
こんな感じかな
~~~~~~~~
互いに素な正整数の組(a,b)全体からなる集合CPに対して関数 f:CP→{1,-1} を
f(1,1)=1,
a>b ならば f(a,b)=f(b,a)
a<b ならば f(a,b)=-f(a,b-a)
と定める。well-defined性は、互除法の成立に関する命題と同様にして示せる。
盤面で一番左にある●と真ん中の●の距離a(cm)と、真ん中の●と一番右の●の距離b(cm)を用いて、
盤面の状態を正整数の組(a,b)で表すことにする。
(複数の盤面で同じ正整数の組になることはあるが、そのような複数の盤面は
平行移動で互いに移り合えるものに限られる)
盤面の状態(a,b)から一秒後に移ることができる盤面は、a>bの時
(a+b,b), (b,a+b), (a-b,b), (b,a-b), (a,a+b), (a+b,b)
のみであり、いずれの場合も新たな整数の組はCPに属する。
更にfの定義を用いると、いずれの場合もfの値が元のf(a,b)と異なるものになることが導かれる。
a<bやa=b(つまりa=b=1)の場合も同様。
以上より、2021秒後の盤面(a,b)のfによる値は
(-1)^2021=-1 であり、(a,b)=(1,1) となることはあり得ない。
~~~~~~~~
つまり互いに素な正整数の組に対して定まる、ある種の『符号』みたいなものが
存在するということなのね、知らなかった
765:132人目の素数さん
20/01/19 13:59:39 BJcaYRJm.net
偶数回でいいなら逆回ししたらいいだけでは?
766:イナ ◆/7jUdUKiSM
20/01/19 15:25:14 6fYcVEZg.net
前>>725
>>719がなりますだったらなんだってなりますだよ。>>714だってなりますだ。
左に一個ずつズレていいならいくつズレようが三個並べばいいってことじゃないか。
正解は2021回だと元の位置にはならない、だと思う。
767:132人目の素数さん
20/01/19 15:32:23 EEjXWdAC.net
>>674
Aと円の中心と対称となる点を、目指して常に方向転換するという方針で、定速で作図してみた。
URLリンク(i.imgur.com)
Bの曲線の接線と円周の交点がその時点でのAと対称になっている(はず)。
図はAの速度がBの1.2倍のとき。
768:132人目の素数さん
20/01/19 15:45:57 EEjXWdAC.net
Bは円周ギリギリで円運動しておいて、タイミングを見計らって円周に到達すればいいんじゃないかな?
Aがいくら速くても円周との距離を限りなく0に近づければBは捕まらないと思う。
769:132人目の素数さん
20/01/19 17:26:30 umPR//wq.net
>>721は撤回します。
吊ってくる
770:イナ ◆/7jUdUKiSM
20/01/19 18:14:18 6fYcVEZg.net
前>>729
>>730袋のねずみじゃないか。もう捕まってるじゃないか。Bが棄権負けだよ。�
771:s要な時間稼ぎ。外周から遠ざかってる。Bが反則負け。 だいたいAがそんな遅くて捕まえられるわけがない。スタート時点でAの進行方向に対して90°の方向にBが逃げると、π倍以上の速さが必要なのは問題にあるとおり。 それじゃ捕まるってんでBが逃げるんだから、Aはもっと速くないと捕まえられないじゃないか。 簡単な話、同じ時間走ってるわけだから、AとBの速さは、AとBが走った距離に比例する。 それより7時25分から7時29分までの短針の方向にBが逃げてAが捕まえたときのグラフが見たいよ。
772:132人目の素数さん
20/01/19 18:28:48 EEjXWdAC.net
>>728
逆回しでなくてもいい経路があるみたい。
> sim()
0 : 1 2 3
1 : 1 3 4
2 : 1 4 5
3 : 1 5 6
4 : 1 4 5
5 : 1 3 4
6 : -1 1 4
7 : 1 3 4
8 : 1 2 3
773:132人目の素数さん
20/01/19 18:38:03 EEjXWdAC.net
>>734
> sim()
[,1] [,2] [,3]
[1,] 0 1 3
[2,] -1 0 3
[3,] 0 1 3
[4,] 0 3 5
[5,] -3 0 5
[6,] 0 3 5
[7,] 0 1 3
[8,] -1 0 3
[9,] -1 3 6
[10,] -1 0 3
[11,] -1 3 6
[12,] -1 0 3
[13,] -2 -1 3
[14,] -2 3 7
[15,] -2 -1 3
[16,] -1 0 3
[17,] 0 1 3
[18,] 1 2 3
774:132人目の素数さん
20/01/19 18:49:27 EEjXWdAC.net
●が長旅してみました。
[,1] [,2] [,3]
[1,] 1 3 4
[2,] -1 1 4
[3,] -1 4 7
[4,] -1 1 4
[5,] -1 4 7
[6,] -6 -1 7
[7,] -6 7 15
[8,] -19 -6 15
[9,] -32 -19 15
[10,] -19 -6 15
[11,] -6 7 15
[12,] -19 -6 15
[13,] -32 -19 15
[14,] -45 -32 15
[15,] -32 -19 15
[16,] -19 -6 15
[17,] -6 7 15
[18,] -6 -1 7
[19,] -6 7 15
[20,] -6 -1 7
[21,] -6 7 15
[22,] -6 15 23
[23,] -6 7 15
[24,] -6 -1 7
[25,] -6 7 15
[26,] -6 -1 7
[27,] -6 7 15
[28,] -19 -6 15
[29,] -6 7 15
[30,] -6 -1 7
[31,] -1 4 7
[32,] -6 -1 7
[33,] -6 7 15
[34,] -6 15 23
[35,] -6 7 15
[36,] -6 -1 7
[37,] -1 4 7
[38,] -6 -1 7
[39,] -1 4 7
[40,] -1 1 4
[41,] -3 -1 4
[42,] -1 1 4
[43,] -1 4 7
[44,] -1 1 4
[45,] 1 3 4
[46,] 1 4 5
[47,] 1 5 6
[48,] 1 4 5
[49,] 1 3 4
[50,] 1 2 3
775:132人目の素数さん
20/01/19 19:23:37 EEjXWdAC.net
>>733
反則負け??
目標に達するのと相手から離れて捕獲を逃れるというとの兼ね合いじゃないの?
776:イナ
20/01/19 20:22:20.24 6fYcVEZg.net
前>>733記録更新した。最善を尽くす、という題意に則って。
>>640問題。
7時25分の短針の方向の外周でAがBを捕まえたとすると、
AがBを追った距離は、
2πr(7+5/12)/12=89πr/72
Bが逃げた距離は、
7時から9時の60°のうちの19/24を底角に持つ二等辺三角形を含む扇形の弧で、その中心角は、
180°-2(19/24)60°=85°だから、
2πR(85°/360°)
二等辺三角形を二等分した直角三角形についてピタゴラスの定理より、
Rcos47.5°=r/2
Rを消してBが逃げた距離は、
2πR(85°/360°)
=17πr/72cos47.5°
速さの倍率は距離に比例し、
A/B=89cos47.5°/17
=3.53691344……(倍)
超えた。
7時半の方向の少し手前で捕まえればAのBに対する速さの倍率は少し大きくなる。
7時25分のあたりだけど、まだ確定じゃない。
777:132人目の素数さん
20/01/19 21:06:37 umPR//wq.net
>>640は何故か出題者が全然顔出さないな。
オレはこの問題出題した元サイトも答えも知ってるので答え書けないけど、Aの最小値を求めるための方程式はまぁまぁシンプル。
色々考えさせて結局これかい?みたいな。
ただし解は明示的には出せないようなので元のパズルサイトでは有効数字6桁まで求めよになってる。
ちなみに>>689はいい線行ってる。
r=1/aより外をどう逃げるか?
778: 【大凶】
20/01/20 00:08:01 jCCIPOX7.net
前>>738最大値はこれだ!!
>>640問題。
7時x分の短針の方向の外周でAがBを捕まえたとすると、
AがBを追った距離は、
2πr(7+x/60)/12=(420+x)πr/360
Bが逃げた距離は、
7時から9時の60°のうちの(1-x/120)を底角に持つ二等辺三角形を含む扇形の弧で、
その中心角は、
180°-2(1-x/120)60°
=(60+x)°だから、
2πR(60+x)/360
=πr(60+x)/180
二等辺三角形を二等分した直角三角形についてピタゴラスの定理より、
Rcos(60-x/2)°=r/2
Rを消してBが逃げた距離は、
2πR{(60+x)°/360°)
=2πr(60+x)/360・2cos(60-x/2)°
=πr(60+x)/360cos(60-x/2)°
速さの倍率は距離に比例し、
A/B=(420+x)πr/360÷πr(60+x)/360cos(60-x/2)°
={(420+x)/(60+x)}cos(60-x/2)°
={7-6x/(60+x)}cos(60-x/2)°
x=26のとき、
A/B=3.53687522……(倍)
x=27のとき、
A/B=3.53671834……(倍)
x=28のとき、
A/B=3.53644261……(倍)
x=29のとき、
A/B=3.53604786……(倍)
逆だ。減ってる。
x=25.5のとき、
A/B=3.53690915……(倍)
x=25.4のとき、
A/B=3.53691238……(倍)
x=25.3のとき、
A/B=3.53691442……(倍)
x=25.2のとき、
A/B=3.53691528……(倍)
x=25.19のとき、
A/B=3.5369153……(倍)
x=25.18のとき、
A/B=3.53691531……(倍)
x=25.17のとき、
A/B=3.5369153……(倍)
x=25.16のとき、
A/B=3.53691529……(倍)
出ましたな、最大値。
7時25.18分すなわち、
7時25分10秒8の短針の方向の外周で捕まえればAのBに対する速さの倍率は最大になる。
そのときAの速度はBの速度の、
3.53691531……倍。
779:132人目の素数さん
20/01/20 00:13:53 vpDptCaR.net
すでにπ+1で脱出可能ってレスが出てるんだけどねぇ。
780:イナ
20/01/20 02:04:33.14 jCCIPOX7.net
前>>740なぜか冒頭が欠けてたみたい。加筆。
>>640問題。
7時x分の短針の方向の外周でAがBを捕まえたとすると、AがBを追った距離は、
2πr(7+x/60)/12=(420+x)πr/360
Bが逃げた距離は、7時から9時の60°のうちの(1-x/120)を底角に持つ二等辺三角形を含む扇形の弧で、その中心角は、
180°-2(1-x/120)60°
=(60+x)°だから、
2πR(60+x)/360
=πr(60+x)/180
二等辺三角形を二等分した直角三角形についてピタゴラスの定理より、
Rcos(60-x/2)°=r/2
Rを消してBが逃げた距離は、
2πR{(60+x)°/360°}
=2πr(60+x)/360・2cos(60-x/2)°
=πr(60+x)/360cos(60-x/2)°
速さの倍率は距離に比例し、
A/B=(420+x)πr/360÷πr(60+x)/360cos(60-x/2)°
={(420+x)/(60+x)}cos(60-x/2)°
={7-6x/(60+x)}cos(60-x/2)°
x=26のとき、
A/B=3.53687522……(倍)
x=27のとき、
A/B=3.53671834……(倍)
x=28のとき、
A/B=3.53644261……(倍)
x=29のとき、
A/B=3.53604786……(倍)
逆だ。減ってる。
x=25.5のとき、
A/B=3.53690915……(倍)
x=25.4のとき、
A/B=3.53691238……(倍)
x=25.3のとき、
A/B=3.53691442……(倍)
x=25.2のとき、
A/B=3.53691528……(倍)
x=25.19のとき、
A/B=3.5369153……(倍)
x=25.18のとき、
A/B=3.53691531……(倍)
x=25.17のとき、
A/B=3.5369153……(倍)
x=25.16のとき、
A/B=3.53691529……(倍)
最大値は、
7時25.18分すなわち、
7時25分10秒8の短針の方向の外周で捕まえればAのBに対する速さの倍率は最大になる。そのときAの速度はBの速度の、
3.53691531……倍。
中心角でいうと、
360°{(7+25.18/60)}/12
=222.59°の方向の外周。
速度4倍もは必要ない。
781:132人目の素数さん
20/01/20 02:20:46.64 vpDptCaR.net
4倍ならAは捕まえられないという正しい答えがすでに出てんのに、なんで自信満々におかしな答え書き込めるん?
782:イナ ◆/7jUdUKiSM
20/01/20 05:48:52 jCCIPOX7.net
前>>742
>>640問題。
速度比A/Bをさらに大きくできないかと考えた。
中心角でいうと、
360°{(7+25.18/60)}/12
=222.59°
の方向の外周が、AがBを捕まえる地点。だという読みだった。
が、もう逃げるBが大きく螺旋状に内側にカーブして222.59°よりも先に逃げようとしたところで、Aは逆回りしてスタート地点に戻るより、そのまま外周付近で接近する。
中心方向に逃げるのは題意に反するが、Bの逃走距離はカーブしてるぶんrを確実に超えていて、Aがスタートした地点まではコイルをのばすようにして行けるんじゃないか。
Bの逃走距離を7時25分10秒8の短針の方向の外周に達するときと変えることなく、Aがスタートした地点までコイルをのばすようにBの逃走ルートをのばすと、
速度比A/B=2πr/{πr(60+25.18)/360cos(60-25.18/2)°}
=720/85.18cos47.91°
=5.66581262……(倍)
どうだ。この速度比。最大だろう。
783:132人目の素数さん
20/01/20 06:45:36 eVrMliO5.net
嘘つき問題に条件文を加味してみた。
AからEの5人はそれぞれ正直者か嘘つきのどちらかであり、誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきの可能性があるのは誰か?
A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」
784:132人目の素数さん
20/01/20 07:10:21 eVrMliO5.net
>>745
無思考解のRでのプログラム(未検証)
n=5
TE=gtools::permutations(2,n,v=0:1,rep=T) # TE : 第1行が00000 で始まり最終行が11111で終わる行列
colnames(TE)=LETTERS[1:n] # 各列の名前A~E
foo <- function(x){ # TEの各行を判定する関数
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,foo),] # 各行にfooを適用して返り値がTRUEのものを表示
785:イナ ◆/7jUdUKiSM
20/01/20 07:55:45 jCCIPOX7.net
前>>744
>>640問題。
Aが外周を一周してスタート地点で最短距離を逃げたBを捕まえたとしたら究極、
速度比A/B=2πr/r
=2π
=6.2831853……(倍)
これが最大値か。
786:132人目の素数さん
20/01/20 09:11:31 jCm5AHw1.net
>>747
逃げられるわけないwww
787:イナ
20/01/20 14:01:00.13 jCCIPOX7.net
>>640前>>747
─
788:「Aの速度がBの速度の何倍以上でなければならないか」 この題意、最大値じゃなくて最小値か? 速度比A/Bの最小値か。 Aが一周するあいだに、 Bは3/4周か? いやなるべく早く外周に到達したいはず。270°より手前で外周に達することができる。 3.5369153倍と3.5369153倍のあいだにある3.53691531倍。これが最速だ。少数第7位までだと決着がつかない。少数第8位を比べる必要がある。元ネタと大きく違う点。これが面白い問題たる由縁。 ∴Aの速度がBの速度の3.53691531倍以上でなければならない。
789:132人目の素数さん
20/01/20 15:10:05 M+lgfEpa.net
何ですでに上がってるπ+1を無視するの?
790:イナ ◆/7jUdUKiSM
20/01/20 16:10:56 jCCIPOX7.net
前>>749
>>750どういう式と計算でπ+1が出たか、AがBを捕まえた地点は外周のどこなのか、式と計算と言葉で示さないと。
Bは円軌道で外周に逃げるのが最善と思って半径Rをrで表して円弧の中心角から逃走距離を求めた。
Bが螺旋か楕円かハートか、スタート方向をAに対して90°にしたまま最短で外周に逃げる形がほかにあるのかどうか。
カブトガニの例もある。途中から直線なんてのもあなどれない。
791:132人目の素数さん
20/01/20 16:25:22 jCm5AHw1.net
>>751
示されてる。
理解できてないのは君だけ。
自分の読解力がないのを他人のせいに平気でするから嫌われるんだよ。
いくら便所の落書きでもそういう最低限の礼節は守らないといけない。
ましてやコテハンで自分の素性も明らかにしてるんでしょ?
自分の知り合いに自分が掲示板でそういう言動してるのバレたりしたときの事とか考えないの?
792:イナ
20/01/20 16:51:17.85 jCCIPOX7.net
前>>751
AがBの外周到達地点を予測して逆回りしてBを捕まえられる限界折り返し地点は、
1時18分10秒8の短針の方向。
Aは残り半周でBに追いつくからどっちを回ってもいい。
もうAが逆回りしないとなった瞬間、Bは最善の方法をとる。
すなわちBは円軌道の必要がなくなり、直線軌道に変える。円弧より内側をえぐったほうが速い。
円弧と直線の交点をつきとめれば速度比は決まる。
793:132人目の素数さん
20/01/20 17:40:59 62Dbolk5.net
オレは出題者じゃないからあんまり書き込むのも何だとは思うけどちょい書いてみる。
この問題はいわゆるゲームの理論なんだけどちゃんと数学的に記述するのはかなり難しい。
一例としてはa,bを正の定数としてA,Bがとりうる選択肢は
F={f : [0,∞) → ∂D | f(0)=(1,0), d(f(t),f(u)) ≦ a|t-u|}
G={f : [0,∞) → D | g(0)=(0,0), d(g(t),g(u)) ≦ b|t-u|}。
Aの戦略とはS:G→Fで
( 0≦∀t≦T g1(t)=g2((t) )⇒ ( 0≦∀t≦T S(g1)(t)=S(g2)(t) )
(Bの選択した関数に対し時刻Tまでの出方で時刻Tまでの対抗行動は決まる。)
Bの戦略T:F→Gも同様に定める。
見つけるべき定数Cとは
(1)∀a/b>C ∃S ∀g∈G ∃t0 S(g)(t0)∈∂D, S(g)(t0)=g(t0), S(g(t)) ∈int(D) (∀t<t0) or ∀t g(t) ∈ int(D)
(2)∀a/b<C ∃T ∀f∈F ∃t0 T(f)(t0)∈∂D, T(f)(t0)≠f(t0), T(f(t)) ∈int(D) (∀t<t0)
の両方を満たす定数。
この戦略関数SとTをCと抱き合わせで見つけないといけないのが難しい。
Aの戦略はまぁそりゃそうだというもの、簡単。
基本Bの動きに応じて右に回るか左に回るかしかないんだから。
Bの脱出戦略Tが難しい。
AとBのその時点での相対位置からどっちに向かうのが得か?
半径b/a-εの地点まで�
794:フ戦略は簡単なのだけれどその先がムズイ。 答え聞くとまぁそりゃそうなのかもなと思えるけど。 ちなみに出題してる元サイトではCの値出せたら正解みたいだった。 それだけなら方程式を勘で当てれなくもない。
795:132人目の素数さん
20/01/20 18:20:54 62Dbolk5.net
Aの戦略関数Sは簡単なので一例として書いてみる。
Dの点pに対し∂D上のベクトル場X(p)をX(0,0)=0、p=(0,0)以外に対してはφをpに最も近い円周上の点として
X(p)(θ)
=正の方向に向かう大きさaのベクトル(θからφへは正の方向に向かう方が近いときか、θとφが原点対称のとき)
=負の方向に向かう大きさaのベクトル(θからφへは負の方向に向かう方が近いとき)
=0(θ=φの時)
て定めてS(g)(t)=exp(X(g(t)))(1,0)で定める。
つまりは常にg(t)との偏角差をなくす方向に速度aで向かう。
偏角差0なら動かない、偏角差πなら正の方向。
まぁこれが最適戦略なのはそりゃそうだと思える。
この戦略で任意のgを捕まえられるa/bの下限がC。
796:イナ ◆/7jUdUKiSM
20/01/20 18:37:08 jCCIPOX7.net
前>>753
>>754数字は3~9も使ったほうがいい値が出せると思う。
計算がまだ追いついてないだけで、Bの逃走経路の、
(1+25.18/60)/(7+25.18/60)を円弧のまま、
残りを222.59°地点まで直線としたら、
速度比A/Bは3.53691531をわずかに超えるはず。
797:132人目の素数さん
20/01/20 18:45:29.00 MTgckgfx.net
>>746
論証を考えていたらバグを発見した。関数を訂正。
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)) )
))
}
798:132人目の素数さん
20/01/20 19:19:04.08 MTgckgfx.net
>>746
Eの記述を簡略化
n=5
TE=gtools::permutations(2,n,v=0:1,rep=T) # TE : 第1行が00000 で始まり最終行が11111で終わる行列
colnames(TE)=LETTERS[1:n] # 各列の名前A~E
fn <- function(x){
if(sum(x)==n|sum(x)==0) return(FALSE) # 全員が正直か嘘つきなら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のものを表示
実行結果
A B C D E
1 1 1 1 0
Eが嘘つき
799:イナ ◆/7jUdUKiSM
20/01/20 22:37:06 jCCIPOX7.net
前>>756
>>752ここで面白い問題を解いたり教えてもらったりしてることは俺たちだけの秘密だよ。
800:132人目の素数さん
20/01/20 23:07:06.30 YsJCrV7U.net
一辺の長さが 10m の正方形のプールの一つの角に監視員を置く。
この監視員は水中は秒速 1m で,プールの縁上は秒速 2m で移動するものとする。
この監視員がプールのどこへでも到達しうるには,最短で何秒必要か計算せよ。
(某AO入試問題)
801:132人目の素数さん
20/01/21 00:19:12 sGzwnwIV.net
>>760
監視員の位置を原点、プールを0≦x≦10, 0≦y≦10とする。
プール内の0≦y≦xにある地点に到達する所用時間の最大値を求めればよい。
この場合監視員の陸路はy=0とx=10を移動する場合のみを考えればよい。
この時時刻tまでに監視員が到達できる領域は
x+√3y≦2t‥?、-√3x+y≦2t-10√3-10‥?
である。
最後まで残る点はy=x上でありy=x上の?に含まれる点は
x≦2t/(√3+1)の部分であり、?のそれはx≧(-2t+10√3+10)/(√3-1)を満たす部分である。よって?、?で全て覆われる時間は
2t/(√3+1) = (-2t+10√3+10)/(√3-1)
の時でありt =5/√3の時である。
かな。自信なし。
802: 【吉】
20/01/21 00:31:04 WFhY3+vZ.net
前>>759ほんとはプールサイドから斜めに飛びこむわなぁ。
>>760
プールの対角線上のじゅうぶん遠いところに、救出に最長時間が必要な地点と考えられる。
x秒かかるとすると最初に監視員がいるコーナーから、
10x(m)離れている。
縁を端まで行くと10/12=5
803:/6(秒)かかる。 コーナーで折れて最長時間が必要な救出地点に縁上で最接近するため10x/√2(m)縁を行く。 水に入って10-10x/√2(m)泳ぐ。 救出時間で等式を作ると、 x=10/12+(10x/√2)/12+(10-10x/√2)/10 分母を払って、 12x√2=10√2+10x+12√2-12x 12x+x√2=22 x=22/(12+√2) =22(12-√2)/(144-2) =11(12-√2)/71 =1.64005142……(秒) ただ優秀な監視員なら縁から斜めに飛びこんで1.6秒ぐらいで救出する可能性がある。
804:132人目の素数さん
20/01/21 01:11:21.37 i1LLIpbZ.net
>>761
> 2t/(√3+1) = (-2t+10√3+10)/(√3-1)
> の時でありt =5/√3の時である。
式はあってると思うけど、計算間違いかな
明らかに10秒ちょっとかかるし
805:132人目の素数さん
20/01/21 01:12:26.37 i1LLIpbZ.net
t=5+10/√3
806:132人目の素数さん
20/01/21 01:32:59.13 sGzwnwIV.net
>>761
あ、対角の位置から飛び込む方が早い可能性抜けてた。
対角の位置から飛び込んだ場合にカバーされる領域は
(x-10)^2+(y-10)^2≦(t-10)^2
x=y上ではx=10-(t-10)/√2。
2t/(√3+1)=10-(t-10)/√2
wolfram大先生によると
t=10 (2 + sqrt(2)))/(4 + sqrt(2) + 4 sqrt(3)
=2. 76624393725438801
だそうな。立式まちがってんのかな?
807:132人目の素数さん
20/01/21 01:35:50.13 sGzwnwIV.net
>>763
計算が間違ってるのはwolfram先生にも教えてもらった。
立式も②の領域は対角のコーナーから飛び込んだ時に負けるみたい。
808:イナ
20/01/21 01:47:56.36 WFhY3+vZ.net
前>>762
>>760ごめん、目を疑うほど間違えた。一行飛ばして速さ二桁にしてた。
x秒とする。
(x/√2)/2=x/2√2秒
10-x/√2秒
5+x/2√2+10-x/√2=x
分母を払って、
10√2+x+20√2-2x=2x√2
30√2=(2√2+1)x
x=30√2/(2√2+1)
=30√2(2√2-1)/7
=(120-30√2)/7
=11.0819419……(秒)
斜めに飛びこむときがはずだけど、とりあえず。
809:132人目の素数さん
20/01/21 02:06:14 i1LLIpbZ.net
>>766
立式は結局
>>761
> x+√3y≦2t‥?、-√3x+y≦2t-10√3-10‥?
> 2t/(√3+1) = (-2t+10√3+10)/(√3-1)
で合っていると思うけど?
wolfram大先生も
t=5+10/√3≒10.774
と言ってくれている
> 立式も?の領域は対角のコーナーから飛び込んだ時に負けるみたい。
対角に着いてから飛び込むよりも、対角に着く前に飛び込んだ方が速いよ
> 2t/(√3+1)=10-(t-10)/√2
を解くと
t=10*(2+√2)(1+√3)/(4+√2+√6)≒11.862
810:132人目の素数さん
20/01/21 02:07:54 sGzwnwIV.net
またまた訂正。
?は?に負けない。
ので?と?をy=x上で解いた>>764さんが正解。
811:132人目の素数さん
20/01/21 02:11:53 sGzwnwIV.net
>>768
うん、立式合ってた。
wolfram先生に教えてもらう時/が一個抜けてた。
そりゃそうだよな。
対角から飛び込んで勝つハズない。
最初はありえないと思って無視したんだけど一応と思ってwolfram先生に聞く時打ち間違えた。
812:イナ
20/01/21 03:32:58.98 WFhY3+vZ.net
前>>767
>>760
縁から50°ぐらいもかなり速いと思うけど、斜め45°に飛びこむときで解く。
プールの対角線上のじゅうぶん遠いところに、救出に最長時間が必要な地点と考える。
x秒かかるとすると最初に監視員がいるコーナーから、
x(m)離れている。
縁を端まで行くと10/2=5(秒)かかる。
コーナーで折れて最長時間が必要な救出地点を斜め45°に見る縁まで、
(1/2)(10-x√2/2)
=5-x√2/4(秒)
水に入って(10-x√2/2)√2(m)泳ぐ。
救出時間で等式を作ると、
5+5-x√2/4+(10-x√2/2)√2=x
10-x√2/4+10√2-x=x
10+10√2=(2+√2/4)x
分母を払って、
40(1+√2)=(8+√2)x
x=40(1+√2)/(8+√2)
=40(1+√2)(8-√2)/(64-2)=20(8-√2+8√2-2)/31
=20(6+7√2)/31
=10.2577387(秒)
813:イナ
20/01/21 03:47:08.74 WFhY3+vZ.net
前>>771
斜め45°に飛びこむとき、じゅうぶん速くてびっくりした。
斜め40°から斜め50°のとき、意外な極値があるかも。
>>760
814:132人目の素数さん
20/01/21 05:24:17.87 4VohdIcv.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の数は何か?
815:イナ
20/01/21 06:41:12.31 WFhY3+vZ.net
>>773前>>772
1 いいえ はい
2 A 3 4 3 2
B 2 2 4 4
C 4 3 2 3
D 1 1
E 5 5 出番なし
816:132人目の素数さん
20/01/21 06:42:13.56 Y0gh5JcA.net
x>0で
0^x=0
x^0=1
0^0=1とするのはその方が0でも辻褄が合う法則が多いから?
0の偏角は不定?それとも0?
偏角が0なら実数という法則を成立させるには0の偏角=0と定義でいいと思うけど。
817:132人目の素数さん
20/01/21 06:46:16.71 Y0gh5JcA.net
>>759
イナさんの芸風は楽しみにしています。
お気になさらず続けてください。
読みたくない人はコテハンをNGに設定すればいいだけですから。
818:132人目の素数さん
20/01/21 08:29:21.13 udoX+djG.net
>>775
0個の物から重複を許して0個取り出して並べる順列は1通りだけど,これは0^0通りとも計算できるから,0^0=1
実数に対しては色々定義がありうるけど基数としては明確に定まる
819:132人目の素数さん
20/01/21 09:47:53 Y0gh5JcA.net
>>760
プログラムを組んで
Oの位置にいる監視員がZで溺れている人に到達する時間を 経路OZ, OXZ, OPYZ, OPQRZ で計算してみると。
URLリンク(i.imgur.com)
場所によって最短到達経路に違いがでる。
(5,5)だとOXZで6.8秒,(8,9)だとOPYZで11.2秒が最短になった。
820:哀れな素人
20/01/21 09:56:57 dWPrQnYr.net
>>773
問1 「いいえ」
問2 A=3か4 B=2 C=3か4 D=1 E=5
Dは1か5。なぜなら明確に答えられるのは1か5のカードを持っている生徒だけだから。
しかしBが「いいえ」と答えたということはD=1、B=2。
A、B、Cがいずれも「分りません」と答えたということはE=5。
今のところ、AとCのカードは不明。
821:132人目の素数さん
20/01/21 10:08:46 udoX+djG.net
>>775
原点を通って偏角一定の曲線(=直線)を考えると0の偏角を0にしちゃうと原点で偏角が不連続になるからまずい
822:哀れな素人
20/01/21 10:12:29.11 dWPrQnYr.net
>>773
やや訂正。次のような場合もある。
問1「分りません」
問2 A=1 B=2 C=? D=? E=5
823:132人目の素数さん
20/01/21 10:24:39 Y0gh5JcA.net
>>778
バグ発見したので図以外は>778は撤回します。
824:哀れな素人
20/01/21 10:44:36 dWPrQnYr.net
>>773
分った。
問1 「はい」
問2 A=1 B=2 C=3 D=4 E=5
Bは1でも5でもないと分るから、Bは2か3か4。
Dが4を持っていれば確実に「はい」と答えることができる。
825:132人目の素数さん
20/01/21 12:57:25 NlSt5Qji.net
>>739
r=1/a まで来たらあとは直進ですか。
>>751 >>753 にありますね。
Aが逆転すれば、Bはその時のOBに垂直な向きに進む。
逆転がない場合は
Aの進む距離 (弧長) はπ+θ、Bの進む距離 (半弦) は sinθ
ここに、中心角θ = arccos(1/a)
逃げ切り条件:
tanθ - θ = a・sinθ - θ < π (0<θ<π/2)
から
θ < 1.35181680431927
a = 1/cosθ < 4.6033388487517
π+1 = 4.1416 より大きい 。
826:132人目の素数さん
20/01/21 13:11:35 9Sn3mJld.net
おぉ、出ましたね。
a/bの臨界値のための方程式。
元サイトではその数値出せば正解です。
のでそれで終わりでもいいし、興味ある人は
a/b>4.6033‥のときのAの補足戦略と
a/b<4.6034..のときのBの逃走戦略
に挑戦してみてはどうでしょうか?
元サイト
URLリンク(www.research.ibm.com)
827:132人目の素数さん
20/01/21 13:23:03 Y0gh5JcA.net
作図と計算をやり直してみた。
原点Oの監視員がZにまで達する時間と経路別に計算
URLリンク(i.imgur.com)
Zの座標から各経路での最短時間を計算させて表示。
> sim(7.9+7.9i,print=F)
OZ PZ QZ RZ OXZ OUZ ORWZ OPYZ OPQWZ
11.17229 13.17435 12.96985 13.17435 10.79160 10.79160 10.76865 10.76865 12.86865
座標を0.1区切りで組み合わせたら最短時間がもっとも大きいのが上記であった。
経路は座標が(7.9,7.9)のときOPYZ(またはORWZ)の経路で最短でも10.76秒かかるという結果になった。
828:132人目の素数さん
20/01/21 13:53:12.95 NlSt5Qji.net
>>773
A<5, B≠1,5, C≠1,5
・D=5 なら D「はい」
・D=1 なら D「いいえ」
・2≦D≦4 のとき
A=1, {B,C,D} = {2,3,4} E=5
D=4 ならD「はい」
D=2 ならD「いいえ」
∴D「分かりません」はD=3のみ。
B「いいえ」 より B=2, C=4
>>779
D「いいえ」の場合
A=1, {B,C}={3,4}, D=2, E=5
もある。
>>781
A=1, B=2, C=3, D=4, E=5
ならD「はい」
>>783
D「はい」の場合
D=5 もある。
829:イナ
20/01/21 14:02:00.34 WFhY3+vZ.net
前>>774先生は「わかりません」と言えなかったんだね。
>>771斜め45°に飛びこんだほうが速いと思う。
20(6+7√2)/31
=10.2577387(秒)
830:132人目の素数さん
20/01/21 14:10:50.45 9Sn3mJld.net
プールのやつは>>764さんの
5+10/√3=10.773502691896...
だろ?
831:イナ ◆/7jUdUKiSM
20/01/21 14:30:17 WFhY3+vZ.net
20(6+7√2)/31
=10.2577387(秒)
>>789見てないのか? コンマ5以上速いぞ? まだもっと速い角度で飛ぶ奴いる気がして探してるけど。前>>788
832:132人目の素数さん
20/01/21 14:34:03 9Sn3mJld.net
ふっと考えたんだけど>>764さんの数値と>>786さんの数値がまぁまぁ離れてるのはなるほどですな。
本物はその地点までの最短到達時間を三次元的なグラフにした場合を考えると各格しだピラミッドみたいな形になる。
いわゆる微分可能な関数の極直ではないからモンテカルロやメッシュがあまりいい数値を出せないんだな。
833:132人目の素数さん
20/01/21 14:37:45 9Sn3mJld.net
>>790
何度で飛び込むのが最適かすら間違ってるのに読む気になどならない。
834:132人目の素数さん
20/01/21 15:42:18 Y0gh5JcA.net
> pm[apply(pm,1,Yes),]
[1] 1 2 3 4 5
> pm[apply(pm,1,No),]
[1] 1 3 4 2 5
> pm[apply(pm,1,DK),]
[1] 1 2 4 3 5
はい で 1 2 3 4 5
いいえ で 1 3 4 2 5
分からんで 1 2 4 3 5
835:イナ ◆/7jUdUKiSM
20/01/21 15:43:35 WFhY3+vZ.net
前>>790
45°は勘だけど、じゅうぶん速かった。
ほかに10.25秒台は出てない。
初め90°出して次に30°出していっしょだな、と。
コンピューターの図があるレスもそこの数値は同じだと出てる。
あいだだ。10.2577387秒が今のところ最速。
3:4:5は11.0819419
836:秒かかる。 4:3:5は10.868秒かかる。
837:イナ ◆/7jUdUKiSM
20/01/21 16:27:37 WFhY3+vZ.net
前>>794
√3:1:2のとき、60°で飛ぶ奴は10.53849秒。一様分布じゃないみたい。
あいだにある4:3:5が遅くて45°が逆に速い。なぜかはわからん。
46°~50°があるいは。
たぶん距離の影響と速さの影響の兼ね合いではないかと思う。
838:132人目の素数さん
20/01/21 16:39:01 Y0gh5JcA.net
>>761
お手数ですが、この不等式の導入法を解説していただけませんか?
839:132人目の素数さん
20/01/21 16:43:00 Y0gh5JcA.net
>>760
グリッドつくって等高線表示させてみた
URLリンク(i.imgur.com)
840:132人目の素数さん
20/01/21 16:45:11 uUzv/iS9.net
>>785
あ、元サイトってパズルの国のアリスじゃなかったのか。
URLリンク(www.nikkei-science.com)
841:132人目の素数さん
20/01/21 17:04:17 Y0gh5JcA.net
>>797
3D グラフにしてみた。
URLリンク(i.imgur.com)