20/01/16 08:04:16.33 552XyIPx.net
密度が一定の球形の惑星がある。
この惑星の表面点Pでの重力を最大にするために
体積と密度一定の条件で惑星の形を変形する。
点Pでの重力は球形のときの最大何倍になるか。
またこのときの惑星の形はどうなるか。
659:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/16 09:39:28 PPdqHAdt.net
>>627
だな。
l,n上にABを好きにとればよかった。
663:１３２人目の素数さん
20/01/16 09:40:02 YUeZdYQq.net
定規とコンパスである2点間の長さを半径とする円をその2点とは別の点を中心に描くのは面倒くさい（幾何学原論ではコンパスは中心からある点までの距離を半径とする円を描けるだけで長さを保存して中心を移動させてはいけない）
664:１３２人目の素数さん
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:１３２人目の素数さん
20/01/16 09:45:23 2aq5oe6O.net
まぁもっと楽な方法があったら教えて下さい。
自分は解けたのでもう満足。
666:１３２人目の素数さん
20/01/16 12:17:23 /nA7kUIy.net
>>625
半円かな？
667:哀れな素人
20/01/16 12:50:07 7jomMr1V.net
>>602
とりあえず四角形の場合について初等幾何的に証明すると－
四辺をａ、ｂ、ｃ、ｄ、対角線をｅ、ｆ、内心半径をr1、r2、r3、r4とする。
ｅと、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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/17 00:11:50 m+4frRir.net
>>641
kwsk
681:１３２人目の素数さん
20/01/17 00:23:11 gpQhzbqr.net
>>640
某サイトに答えあるの見つけた。
コピペするのもあれなのでROMします。
682:１３２人目の素数さん
20/01/17 00:26:36 pPPuMatU.net
>>646
点Pでの重力を体積一定の下で最大化するような断面の形状を変分法で求める
惑星は軸対称（軸対称でないとすると点Pから少し動いたところにより重心に近い点がある）だから断面の形状だけ考えれば良いのがポイント
↓のpart II の 4 に等長問題のオイラーラグランジュの方程式を求める方法が書いてある
URLﾘﾝｸ(citeseerx.ist.psu.edu)
683:１３２人目の素数さん
20/01/17 00:30:13 wItShquN.net
ちなみに速度比が2倍では捕まえられない。
Bの速度は1とする。
速度比が2倍ならBは半径1/2の地点まではAの角速度を上回っているのでここまでBはAとの偏角の差が180°の状態をキープできる。
この地点から円周の最短地点までBが要する時間は1。
Aが要する時間はπなのでBの逃げ切り成功になる。
684:１３２人目の素数さん
20/01/17 00:45:02 XDQcvfLQ.net
>>648
この問題に最大解がある事の証明はどうすればいいのでしょう？
変分法使う以上は別に解の存在示しとかないとダメだとおもうんですが。
685:１３２人目の素数さん
20/01/17 00:54:20.56 pPPuMatU.net
>>650
リンク先にも書いてあるけど,厳密には第二変分を考えなきゃいけない
686:１３２人目の素数さん
20/01/17 00:56:48.70 Vq8XRGW9.net
>>659
そんなのがあるんですか？
なんか言葉のイメージだとHessianみたいなものですか？
でもそれだと結局極大である事の保証しか得られないのでは？
687:１３２人目の素数さん
20/01/17 01:30:08 pPPuMatU.net
>>652
具体的に求めてみると分かるけど汎関数を与える関数（被積分関数）が凸関数だから結局極大値が最大値になる
688:１３２人目の素数さん
20/01/17 01:46:52.25
689: ID:0cuf6Wtp.net
690:１３２人目の素数さん
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:１３２人目の素数さん
20/01/17 02:11:57 B81I4URB.net
今チェックしてみたら行けそうですね。
素晴らしい。
692:１３２人目の素数さん
20/01/17 08:18:19 F8/vjJrP.net
別スレより
1歩で1段または2段のいずれかで階段を昇るとき、1歩で2段昇ることは連続しないものとする。15段の階段を昇る昇り方は何通りあるかを求めよ。
693:１３２人目の素数さん
20/01/17 09:06:38 Ax5k6DPo.net
2段昇りの回数をnとすると、連続できないので最低限
間に1段昇りがはいることから、2n+(n-1)≦15 よって、
0≦n≦5
n=0の場合、一通り
n≠0の場合、15-2n個の１段昇りの隙間と両端に２段昇
りが入りうるので、その組み合わせは C(15-2n +1, n)
よって、全部で1+Σ(n=1～5）C(16-2n,n)
694:１３２人目の素数さん
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:１３２人目の素数さん
20/01/17 10:16:22 CPtT5WKR.net
あ、一歩で二段連続なしかorz
696:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/17 15:34:25 k0YCxJCN.net
受験数学でクソ頻出の3項間の項の数が増えてるだけですがな。
これが即答できないのは受験レベルの3項間の勉強の時に解き方だけ覚えて理屈がわかってなかった証拠。
702:１３２人目の素数さん
20/01/17 16:34:49 ElKm6FvG.net
2007年の京大の問題か
703:１３２人目の素数さん
20/01/17 17:04:55 F8/vjJrP.net
>>667
検索したらそうみたいですね。
704:俺は、ニュー速＋に書いてあったを転載したんだけど。 【勉強】数をかぞえられない学生たち ありもしない「公式」に頼り…算数教育の珍現象 https://asahi.5ch.net/test/read.cgi/newsplus/1579164227/
705:１３２人目の素数さん
20/01/17 17:50:57 F8/vjJrP.net
>>640
Bは常にAと正反対の方向に向かって逃げてればいいような気がするな。
こんな感じ URLﾘﾝｸ(imgur.com)
706:１３２人目の素数さん
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:１３２人目の素数さん
20/01/17 20:11:26 F8/vjJrP.net
>>671
画像のURL
URLﾘﾝｸ(i.imgur.com)
709:１３２人目の素数さん
20/01/17 20:37:41.55 F8/vjJrP.net
3.5倍の速度でも常にAと反対向きに進路変更すれば逃げ切れる気がするんだが。
URLﾘﾝｸ(i.imgur.com)
710:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/18 08:49:58.52 aX+4W30D.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
718:１３２人目の素数さん
20/01/18 08:57:10.68 2fzpXJcR.net
正解です。
719:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/18 09:58:04 PXyx+OwH.net
>>678
Re{β} = -(1/2)|β|^4 = -α^(-2)/2,
cos(ω) = -(1/2)|β|^3,
722:１３２人目の素数さん
20/01/18 10:57:32.36 aX+4W30D.net
>>679
これって嘘つきの証言も正しいという前提？
723:１３２人目の素数さん
20/01/18 11:27:02.31 aX+4W30D.net
>>679
Aが嘘つきの場合は
「8人の中に、少なくとも1人正直者がいる」は嘘の証言で
正直者は誰もいないとして計算？
724:１３２人目の素数さん
20/01/18 11:43:01.17 EsuXaTYh.net
この問題で問題文の意味の取りようなんて議論の余地ないだろ？
725:１３２人目の素数さん
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) を保ちながら速さ１で逃げる。
中心～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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/18 15:05:40 zvEtuxN+.net
>>691
>一般項は求まらないか、、、。
3項漸化式で特性方程式が3次
その解を使えば？
730:１３２人目の素数さん
20/01/18 15:08:21 kUGqA2Ne.net
嘘つき問題を論理式で解いたらどうなるんだろ
731:１３２人目の素数さん
20/01/18 15:34:40 EsuXaTYh.net
普通に命題論理の問題になるけど例えば>>679の問題の
B「8人の中に、少なくとも2人正直者がいる」
を命題論理の式にすると
B=(A∩B)∪(A∩C)‥∪(G∩H)
になって手計算では無理な話になってしまう。
計算機使えれば一撃だけど。
732:１３２人目の素数さん
20/01/18 15:46:56 zvEtuxN+.net
>>694
自己言及があるが？
733:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/18 17:27:51 dWbEz1QY.net
どこかの中学入試の問題らしい
Ａ、Ｂ、Ｃ、Ｄ、Ｅの５人のうち２人は常に本当のことを言う正直者です。あとの３人は嘘つきですが、その発言内容は本当のときもあります。
彼らに誰が嘘つきか尋ねたところ次のように答えました。
Ａ 「ＢとＥは嘘つきではない」
Ｂ 「Ｃは嘘つきだ」
Ｃ 「Ｄは嘘つきだ」
Ｄ 「Ｅは嘘つきだ」
Ｅ 「ＢとＣは嘘つきだ」
さて、正直者は誰と誰でしょうか？
736:１３２人目の素数さん
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:１３２人目の素数さん
20/01/18 18:22:46 EsuXaTYh.net
>>698
この手の問題計算機使って解かないと解けないのでは単なる計算機の演習にしかならない。
738:１３２人目の素数さん
20/01/18 21:12:52 zvEtuxN+.net
>>701
試行錯誤で解けるが
739:１３２人目の素数さん
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:１３２人目の素数さん
20/01/18 21:31:04 dWbEz1QY.net
>>701
2＾8＝256だから、手作業でやれなくもないとは思うけど
プログラムを組む方が早い。
256通りを以下のような関数でTRUEを返すものを選ぶだけ。まさに無思考ｗｗｗ
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が嘘つき
可能な組み合わせ（１：正直者 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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/19 06:34:16.10 EEjXWdAC.net
>>709
レスありがとうございます。
(H,F)=(正,正)は成立しません。 Hの発言： H「AもFも嘘つきである」 からFが嘘つきになり矛盾します。
(H,F)=(嘘,嘘)も成立しません。 Aは必ず嘘をつく嘘つきと確定しているので正直者は３人ではありえません。
(H,F)=(嘘,嘘)ならば正直者がB,D,Eの３人になるのでAが真実を語っていることになり、矛盾します。
749:１３２人目の素数さん
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:１３２人目の素数さん
20/01/19 08:52:07 MNtmCELY.net
1cm刻みの数直線上の隣り合った3点に●がある(3つの●は区別がない)。
1秒ごとに●1つが他の●1つを飛び越す。
その場合、飛んだあとの●の位置は、飛び越された●からの距離がもとの位置からと同じになるようにする。
飛んだあとの位置にすでに他の●があるときは飛べない。
2021秒後に、●●●全てが開始位置に戻ることはできるか？
例、
●●●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ●●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ● ●
╋╋╋╋╋╋╋╋╋╋╋
↓
● ● ●
╋╋╋╋╋╋╋╋╋╋╋
↓
：
↓
●●●
╋╋╋╋╋╋╋╋╋╋╋
752:１３２人目の素数さん
20/01/19 10:17:19 p/jcyrav.net
>>713
正解が投稿されてから誤答を繰り返すいつもの芸風 乙
753:１３２人目の素数さん
20/01/19 10:18:00 0RKATGn1.net
>>711
うわ、二つもアウトな場合取りこぼすとかぼーっとしてた
なるほど一度使った条件もまた後から効いてくる場合があったか
754:１３２人目の素数さん
20/01/19 10:22:05 mAiNEiBE.net
いくらなんでも2021がムリクリすぎん？
755:１３２人目の素数さん
20/01/19 10:31:11 flshczPE.net
>>703
2＾8=256通り試すだけ
756:１３２人目の素数さん
20/01/19 10:36:09 Ld5ZhyuY.net
>>714
A B C
╋╋╋╋╋╋╋╋╋╋╋
↓
(2021回後)
↓
C A B
╋╋╋╋╋╋╋╋╋╋╋
は元に戻った事になりますか？
757:１３２人目の素数さん
20/01/19 10:37:14 gvXkKI2o.net
ズレたorz
でも疑問の意味はわかってもらえると信じて
758:１３２人目の素数さん
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:１３２人目の素数さん
20/01/19 11:28:59.40 Yx0eiL2o.net
あ、相互位置交換も不可かな？
対称性から(a,b,c)→(b,a,c)が不可である事を示せば十分。
(0,1,2)→(1,0,2)が不可である事を示せば十分だけど各変換はmod2で恒等写像なのでコレは明らか。
760:１３２人目の素数さん
20/01/19 11:49:47 MNtmCELY.net
>>719
●は区別がないので、なります。
761:１３２人目の素数さん
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:１３２人目の素数さん
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回は奇数回だから 元に戻れないと思う。
根拠は直感のみｗ
764:１３２人目の素数さん
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:１３２人目の素数さん
20/01/19 13:59:39 BJcaYRJm.net
偶数回でいいなら逆回ししたらいいだけでは？
766:イナ ◆/7jUdUKiSM
20/01/19 15:25:14 6fYcVEZg.net
前>>725
>>719がなりますだったらなんだってなりますだよ。>>714だってなりますだ。
左に一個ずつズレていいならいくつズレようが三個並べばいいってことじゃないか。
正解は2021回だと元の位置にはならない、だと思う。
767:１３２人目の素数さん
20/01/19 15:32:23 EEjXWdAC.net
>>674
Aと円の中心と対称となる点を、目指して常に方向転換するという方針で、定速で作図してみた。
URLﾘﾝｸ(i.imgur.com)
Bの曲線の接線と円周の交点がその時点でのAと対称になっている（はず）。
図はAの速度がBの１．２倍のとき。
768:１３２人目の素数さん
20/01/19 15:45:57 EEjXWdAC.net
Bは円周ギリギリで円運動しておいて、タイミングを見計らって円周に到達すればいいんじゃないかな？
Aがいくら速くても円周との距離を限りなく０に近づければBは捕まらないと思う。
769:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/20 06:45:36 eVrMliO5.net
嘘つき問題に条件文を加味してみた。
AからEの5人はそれぞれ正直者か嘘つきのどちらかであり、誰が正直者か嘘つきかはお互いに知っている。
嘘つきなら必ず嘘をつく。嘘つきの可能性があるのは誰か？
A「Bは正直者である」
B「Aは正直者である」
C「Bが嘘つきなら私も嘘つきである」
D「Cが正直なら私も正直である」
E「Dが嘘つきなら私も嘘つきであるし、Dが正直ものなら私も正直者である」
784:１３２人目の素数さん
20/01/20 07:10:21 eVrMliO5.net
>>745
無思考解のRでのプログラム（未検証)
n=5
TE=gtools::permutations(2,n,v=0:1,rep=T) # TE : 第１行が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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/20 19:19:04.08 MTgckgfx.net
>>746
Eの記述を簡略化
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) # 全員が正直か嘘つきなら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:１３２人目の素数さん
20/01/20 23:07:06.30 YsJCrV7U.net
一辺の長さが 10m の正方形のプールの一つの角に監視員を置く。
この監視員は水中は秒速 1m で，プールの縁上は秒速 2m で移動するものとする。
この監視員がプールのどこへでも到達しうるには，最短で何秒必要か計算せよ。
(某AO入試問題)
801:１３２人目の素数さん
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:１３２人目の素数さん
20/01/21 01:11:21.37 i1LLIpbZ.net
>>761
> 2t/(√3+1) = (-2t+10√3+10)/(√3-1)
> の時でありt =5/√3の時である。
式はあってると思うけど、計算間違いかな
明らかに10秒ちょっとかかるし
805:１３２人目の素数さん
20/01/21 01:12:26.37 i1LLIpbZ.net
t=5+10/√3
806:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/21 02:07:54 sGzwnwIV.net
またまた訂正。
?は?に負けない。
ので?と?をy=x上で解いた>>764さんが正解。
811:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/21 06:46:16.71 Y0gh5JcA.net
>>759
イナさんの芸風は楽しみにしています。
お気になさらず続けてください。
読みたくない人はコテハンをNGに設定すればいいだけですから。
818:１３２人目の素数さん
20/01/21 08:29:21.13 udoX+djG.net
>>775
0個の物から重複を許して0個取り出して並べる順列は1通りだけど,これは0^0通りとも計算できるから,0^0=1
実数に対しては色々定義がありうるけど基数としては明確に定まる
819:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
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:１３２人目の素数さん
20/01/21 14:34:03 9Sn3mJld.net
ふっと考えたんだけど>>764さんの数値と>>786さんの数値がまぁまぁ離れてるのはなるほどですな。
本物はその地点までの最短到達時間を三次元的なグラフにした場合を考えると各格しだピラミッドみたいな形になる。
いわゆる微分可能な関数の極直ではないからモンテカルロやメッシュがあまりいい数値を出せないんだな。
833:１３２人目の素数さん
20/01/21 14:37:45 9Sn3mJld.net
>>790
何度で飛び込むのが最適かすら間違ってるのに読む気になどならない。
834:１３２人目の素数さん
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:１３２人目の素数さん
20/01/21 16:39:01 Y0gh5JcA.net
>>761
お手数ですが、この不等式の導入法を解説していただけませんか？
839:１３２人目の素数さん
20/01/21 16:43:00 Y0gh5JcA.net
>>760
グリッドつくって等高線表示させてみた
URLﾘﾝｸ(i.imgur.com)
840:１３２人目の素数さん
20/01/21 16:45:11 uUzv/iS9.net
>>785
あ、元サイトってパズルの国のアリスじゃなかったのか。
URLﾘﾝｸ(www.nikkei-science.com)
841:１３２人目の素数さん
20/01/21 17:04:17 Y0gh5JcA.net
>>797
3D グラフにしてみた。
URLﾘﾝｸ(i.imgur.com)
842:１３２人目の素数さん
20/01/21 17:10:23 B4OFa3kR.net
>>640
>>785
π＋1を正解にしてるサイトもある
URLﾘﾝｸ(www.arp-nt.co.jp)
843:哀れな素人
20/01/21 17:11:22 dWPrQnYr.net
>>787
なるほど。
844:１３２人目の素数さん
20/01/21 17:22:59 Y0gh5JcA.net
グリッド幅を狭くしていったら
(7.88675135,7.88675135)のときに10.77350269秒が最大という計算になった。
5+10/√3＝10.77350269189625764509148780501957455647601751270126876018...
に一致していて、プログラムは正確みたいでほっとした。
俺には理論はわからないけどｗｗｗｗ
845:１３２人目の素数さん
20/01/21 17:28:50 Y0gh5JcA.net
>>795
>なぜかはわからん。
多分、風が吹いているんじゃないの？
馬耳東風という風がｗｗｗ
846:１３２人目の素数さん
20/01/21 18:20:00 Y0gh5JcA.net
オリンピックのプールと世界最高記録を使って計算してみた。
オリンピックサイズ・プール50m*25m
水泳100m自由形 46秒91
陸上100m9秒58
座標(40.101, 15.1077)に 0.77933776秒で達するのが最長と算出された。
847:１３２人目の素数さん
20/01/21 18:22:24 9Sn3mJld.net
>>794
アホかいな。
正解もう出ててその数値より早いという事は経路の選択も所用時間の計算も両方間違ってるんだよ。
848:１３２人目の素数さん
20/01/21 18:23:30 Y0gh5JcA.net
>>804
コピペのミス
秒数は
OXZ
[1] 10.77933776
849:１３２人目の素数さん
20/01/21 19:02:51 LzNLfIhD.net
>>760
プールを 0<x<10,0<y<10 と座標設定。
(a,b) に向かうとする。ただし、0<b<a<10
考えるべき方法は次の(a)～(c)で、それぞれ必要な時間を最後に記すと
(a)原点から直接(a,b) この時、必要な時間は、sqrt(a^2+b^2)
(b)(x,0)まで行ってそこから(a,b)へ この時、必要な時間は、x/2+sqrt((a-x)^2+b^2)
極値を取るのはx=a±b/√3だから、プラスを取って代入し、a/2+((5/6)√3)b
(c)(0,0)→(10,0)→(10,y)→(a,b)へ この時、必要な時間は、5+y/2+sqrt((10-a)^2+(b-y)^2)
同様に、y=b-(1/√3)√(a^2-20a+100)の時、(√3/2)(10-a)+5+b/2
(a,b)地点によって、最適な方法が変化する。図示は某所に下式を入力して欲しい。
min{√(a^2+b^2),a/2+((5/6)√3)b,(√3/2)(10-a)+5+b/2} 0<a<10,0<b<10
最も時間がかかる場所は、方法(a)と方法(c)で必要な時間が同じで、かつ、x=y上
つまり、sqrt(a^2+b^2)=(√3/2)(10-a)+5+b/2,a=b を解いて
a=b=10(-2+√3+2√(2-√3))=7.6732698797896034292...
必要な時間は上の値の√2倍で10(√3-1)(2-√(2-√3))=10.8516423317474258765...
850:１３２人目の素数さん
20/01/21 19:13:42 Y0gh5JcA.net
>>807
極値を取るのはx=a±b/√3 これで√3がでてくる理由がわかりました。
ありがとうご�
851:ｴいした。
852:１３２人目の素数さん
20/01/21 19:28:25 9Sn3mJld.net
>>802
それはホント？
本問微分可能関数の極値ではなく、誤差はグリッドから真の極小点までの距離に正比例する。
例えば誤差を5桁、にしようと思えばグリッド巾は10^(-5)、メッシュ数は10^10の100億個取らないといけない。
比例定数が幾ばくか助けてくれたとしても本問単純なモンテカルロ法やメッシュ法でそこまでの精度が出るとは思えないんだけど。
853:１３２人目の素数さん
20/01/21 19:31:54 Y0gh5JcA.net
>>809
荒いグリッドで極値を与えるx, yの近似値がでてくるからそれを挟むように次の計算で
グリッドの上限と下限を狭くしていけばいい。
人間ニュートンハフソンン法ｗ
x=seq(40.099,40.101,by=0.0001)
y=seq(15.107,15.109,by=0.0001)
854:イナ ◆/7jUdUKiSM
20/01/21 19:33:27 WFhY3+vZ.net
前>>795
周りが遅いから今は俺が最速なだけ。もっと速い角度がないか探してる。
855:１３２人目の素数さん
20/01/21 19:36:38 9Sn3mJld.net
>>810
ニュートンラフソン使うなら極値をとる点なり、極値そのものを与える方程式なりがわかってるのが大前提で>>802はそうでなく単純なメッシュ法で求めたんでしょ？
そもそも本問極値を求める方程式はただの一次方程式にしかならないんだからニュートンラフソンもへったくれもないよ。
856:１３２人目の素数さん
20/01/21 19:36:48 Y0gh5JcA.net
>>809
［０，１００］を１００分し４０が返ってきたら次は[３９，４１]を１００分して４０．０を得る。
その次は[３９．９,４０．１」を１００分する。この繰り返し。
もとの正方形プールの10.77350269秒はそうやってもとめた。
857:１３２人目の素数さん
20/01/21 19:38:45 9Sn3mJld.net
>>813
なるほど。
全領域をメッシュしたんじゃないのか。
それならできるな。
納得しました。
858:１３２人目の素数さん
20/01/21 19:40:11 Y0gh5JcA.net
>>812
いや、時間を求める関数はわかっているよ。
複素平面で絶対値を計算して速度で割っただけ。
関数化するとこんな感じ。
f=function(x) x/vs + abs(x-z)/vw
あとは、ニュートン法でRに最小値の数値解をださせるだけ。
859:１３２人目の素数さん
20/01/21 19:47:14.09 9Sn3mJld.net
>>800
でもこの問題問題文の文章からして正解は>>784の4.6033388‥だと思う。
やはり数学的に厳密に解釈しようとすれば>>754のようにならざるを得ないし、だとすると正解は>>784になるしかないと思う。
多分そのサイトの解答はa/b>π+1のとき捕獲可能であるの証明に誤りがある(Bの最適な逃走戦略を見つけきれてない)のだと思う。
そのサイト答えがπ+1としか書いてないからわかんないけど。
860:１３２人目の素数さん
20/01/21 20:46:55 i1LLIpbZ.net
他人の解答を見る気はないって言ってるからほっとくしかないんだろうが、
監視員から(8m,8m)にかかる時間を計算すればいい。>>788よりかかるぞ
861:１３２人目の素数さん
20/01/21 20:55:34 Y0gh5JcA.net
>>817
最後は監視員がプールの水を抜いて最速は１０秒という答を出すのだと思うんだんが。
862:１３２人目の素数さん
20/01/21 20:57:41 Y0gh5JcA.net
>>818
１０秒じゃなくて５√２秒（＝7.071秒）だった。
863:１３２人目の素数さん
20/01/21 21:07:46 9Sn3mJld.net
>>807はイナ？
864:イナ ◆/7jUdUKiSM
20/01/21 21:20:07 WFhY3+vZ.net
前>>811
>>760問題。
>>771を再度検証する。
プールの端を直角に折れて斜め45°に飛びこむ場合。
プールの対角線上のじゅうぶん遠いところに、救出に最長時間が必要な地点があると考える。
x秒かかるとすると最初に監視員がいるコーナーからx(m)離れている。
縁を端まで行くと10/2=5(秒)かかる。
コーナーで折れて最長時間が必要な救出地点を斜め45°に見る縁まで、
{x/√2-(10-x/√2)}/2
=x/√2-5(秒)
水に入って(10-x√2/2)√2(秒)泳ぐ。
救出時間で等式を作ると、
5+x/√2-5+(10-x√2/2)√2=x
x/√2+10√2-x=x
10√2=(2-1/√2)x
分母を払って、
20=(2√2-1)x
x=20(2√2+1)/7
=10.93835060……(秒)
＜10.53849……(秒)
60°のときに及ばない。
やっぱり一様分布か。
865:１３２人目の素数さん
20/01/21 21:32:26 m9UBU6An.net
>>777
#A^#B=#A^B=#{f:B→A}
0^0=#Φ^Φ=#{f:Φ→Φ}=1
866:１３２人目の素数さん
20/01/21 21:34:29 m9UBU6An.net
>>775
>偏角が0なら実数という法則を成立させるには0の偏角=0と定義でいいと思うけど。
０の偏角に０が有ればいい�
867:ﾅしょ
868:１３２人目の素数さん
20/01/21 22:00:02.52 LzNLfIhD.net
>>807
答えが異なっていたのでアップしたが、致命的なミスを発見
×：極値を取るのはx=a±b/√3だから、プラスを取って代入し、a/2+((5/6)√3)b
○：極値を取るのはx=a±b/√3だから、マイナスを取って代入し、(1/2)a+((√3)/2)b
これにより、以下も訂正
×：min{√(a^2+b^2),a/2+((5/6)√3)b,(√3/2)(10-a)+5+b/2} 0<a<10,0<b<10
○：min{√(a^2+b^2),(1/2)a+((√3)/2)b,(√3/2)(10-a)+5+b/2} 0<a<10,0<b<10
×：最も時間がかかる場所は、方法(a)と方法(c)で必要な時間が同じで、かつ、x=y上
×：つまり、sqrt(a^2+b^2)=(√3/2)(10-a)+5+b/2,a=b を解いて
○：最も時間がかかる場所は、方法(b)と方法(c)で必要な時間が同じで、かつ、x=y上
○：つまり、(1/2)a+((√3)/2)b=(√3/2)(10-a)+5+b/2, a=b を解いて
○：a=b=5+5/√3=7.88675...、時刻は5+10/√3=10.77350269...
869:１３２人目の素数さん
20/01/21 22:41:06 Y0gh5JcA.net
対角線上でも最短時間のルートが青から赤に突然変わる。
URLﾘﾝｸ(i.imgur.com)
870:イナ ◆/7jUdUKiSM
20/01/21 22:59:49 WFhY3+vZ.net
前>>821
縁から60°方向に飛びこむとき、救出時間で等式を作ると、
5+{x/√2-(10-x/√2)/√3}/2+(10-x/√2)(2/√3)=x
分母を払って、
10√6+x√3-10√2+x+40√2-4x=2x√6
10√6+30√2=x(2√6+3-√3)
x=(10√6+30√2)/(2√6+3-√3)
=10.8516423……(秒)
45°は超えたけどなぁ。これ以上は、もしや手前から飛びこむか。プールサイド2倍速で走ってこけてもなんにもならんからね。
871:イナ ◆/7jUdUKiSM
20/01/21 23:29:27 WFhY3+vZ.net
前>>826
>>751
いい勘してる。さすが俺。
4.6倍かぁ。すごいね。
その速さを引きだしたBもたいしたもんだ。
872:１３２人目の素数さん
20/01/22 01:48:55 xZx9jgfS.net
>>784
π = tanθ -θ
= cot(π/2 -θ) -θ
≒ 1/(π/2 -θ) -(π/2 -θ)/3 - θ
= 1/(π/2 -θ) -(1/2)π +(2/3)(π/2 -θ),
より2次方程式
(2/3)(π/2 -θ)^2 -(3/2)π(π/2 -θ) +1 = 0,
これを解いて
(π/2 - θ) = 0.218991
θ = 1.351805
マクローリン展開
cot(x) = 1/x -(1/3)x -(1/45)x^3 -(2/945)x^5 - ････
873:１３２人目の素数さん
20/01/22 02:08:58 vRiVJkwC.net
1/cos(1.351805)=4.60309‥‥
その近似だと求められてる6桁一致までいかないね。
874:イナ
20/01/22 04:50:53.85 it61/f5D.net
前>>827縁と水中の速度比が2:1だからカットする縁と水中の距離の比が1:2になる最適な地点から最適な角度で飛びこむとかなりのタイムを期待できる。救出時間で等式を作ると、
5+{x/√2-(10-x/√2)/2}/2+(10-x/√2)(√5/2)={x/√2-(x/√2)/2}/2+(x/√2)(√5/2)
分母を払って、
20√2+2x-10√2+x+20√10-2x√5=2x-x+2x√5
30√2=-2x+4x√5
15√2=(2√5-1)x
x=15√2(1+2√5)/19
x=(15√2+30√10)/19
=6.10955438……(秒)
＜10.8516423……(秒)
速すぎる。わからん。
875:１３２人目の素数さん
20/01/22 07:58:37 jM4eJElw.net
対角線上が到達に一番時間がかかることが分かったけど
監視員の近くにいればプールサイドを通らず直接ジャンプして水中を進めばいいんじゃないかと思う。
URLﾘﾝｸ(i.imgur.com)
図でいうと青のOXZ1でなくて緑のOZ3を選択。
どれくらい近くだと直接ジャンプすべきかをプログラムで探索させたけど直接ジャンプの方が時間がかかるようだ。
プールサイドの歩行速度が遅ければ直接ジャンプの方が速いはずと考えて探索させると
水泳速度を１として歩行速度が√２以下なら監視員の近くは直接ジャンプが速いようだ。
√２が正しいのか、どれくらい近ければ直接ジャンプすべきなのかは、また後で考える。
876:１３２人目の素数さん
20/01/22 09:02:44 /I6vaW/w.net
>>826
プールの水を凍らせてスピードスケートにすれば最速だよね。
877:１３２人目の素数さん
20/01/22 13:36:03 crsPene3.net
>>760の問題はどっかのAO入試の問題らしいけどどこのなんだろ？
程よい解き心地の良問だね。
このスレの住人には結構受かりそうにないのがいるなww
878:１３２人目の素数さん
20/01/22 14:15:23 MGz/KyFY.net
自分は立式や最適な飛び込み角度を考えるのに、ホイヘンスの原理やスネルの法則を考えたな
879:イナ ◆/7jUdUKiSM
20/01/22 14:19:35 it61/f5D.net
わかったかも。前>>830最適な角度で飛びこんで10.85を切る。救出時間で等式を作ると、
5+{x/√2-(10-x/√2)(1/t)}/2+(10-x/√2)√(1+t^2)/t)=x
分母を払って、
10t√2+2tx-20√2+x+10√130-2x√(1+t^2)=2tx√2
x=
あとは5秒以内にこっち側の縁から飛びこんで連立。
目標。
x=10.7……(秒)
＜10.8516423……(秒)
880:１３２人目の素数さん
20/01/22 14:20:40 FzGnA9Ra.net
そうだな。
スネルの法則知ってれば最速到達の経路は一瞬で出る。
でもAO入試でスネルの法則よりって書いて許してもらえるか微妙だから法則で求めた領域が正しい事の検証の論述は必要だろうけど、それでもまともに円の通過領域求めたり所要時間最小の角度を微積で求めたりするよりははるかに楽になるね。
881:１３２人目の素数さん
20/01/22 14:21:49 sJxssOgt.net
>>835
残念ながらわかってない。
やり直し。
882:１３２人目の素数さん
20/01/22 14:56:09 jM4eJElw.net
>>833
検索したら東京工業大学。
うかりそうもない計算マニア？は東京大学卒の芸人と聞いております。
883:１３２人目の素数さん
20/01/22 15:31:17 99kbu1Vi.net
>>838
そうなんだ。thx
884:１３２人目の素数さん
20/01/22 15:35:12 Zb3S28FJ.net
>>838
ホントだ。2007年みたいですね。
885:１３２人目の素数さん
20/01/22 17:44:16 jM4eJElw.net
プールの中の対角線上の点だけを考える
URLﾘﾝｸ(i.imgur.com)
緑のOZがroute 1 青のOXZをroute 5 赤のOPYZをroute 7 として（番号は区別さえできればなんでもいい）
対角線上のｘ座標(=y座標同じ)とプールサイドの走行速度（陸上速度）をグラフにしてみた。
URLﾘﾝｸ(i.imgur.com)
対角線上の位置に関わらず走行速度がある値（多分√２）以下では直接プールに飛び込むのが最速。
それ以上になると
プールサイドの１辺の途中からプールにジャンプ(OXZの青ルート）
か
プールサイドの２辺めの途中からプールにジャンプ（OPYZの赤ルート）
になるようだ。
青ルートと赤ルートの境界をグラフにしたのが
URLﾘﾝｸ(i.imgur.com)
以上、本日の観察日記でした。
886:１３２人目の素数さん
20/01/22 17:55:14 l1lTbxJu.net
>>841
> 青ルートと赤ルートの境界をグラフにしたのが
2つのルートの境界線は>>799の通り角から伸びる直線
887:１３２人目の素数さん
20/01/22 18:14:54 jM4eJElw.net
Bが円周を描き続けるのは反則負けだとしても、こういうふうに渦巻曲線を描いて円周に近づいてAのスキを見計らって円周に直行すれば捕まらない気がする。
URLﾘﾝｸ(i.imgur.com)
888:１３２人目の素数さん
20/01/22 19:34:22.59 jM4eJElw.net
>>842
>841の曲線グラフは横軸が対角線上の点の座標(x=yの点の値)
縦軸は陸送速度0.5から4m/秒で描いてみた。
>799はプールの座標で縦横10m
889:１３２人目の素数さん
20/01/22 20:14:35 7n0H2YC6.net
>>843
それBは初期状態よりも不利になってるよね
890:１３２人目の素数さん
20/01/23 01:16:24.21 L8diiD+d.net
>>828 を改良
π = tanθ - θ = cot(π/2 -θ) - θ ≒ 1/(π/2 -θ) -θ
より
π/2 - θ ≒ 1/(π+θ) ≒ 2/(3π) = 0.2122
(1/45)(π/2 -θ)^2 ≒ 0.001
ここまで準備して
π = tanθ -θ
= cot(π/2 -θ) -θ
= 1/(π/2 -θ) -(1/3)(π/2 -θ) -(1/45)(π/2 -θ)^3 -θ
= 1/(π/2 -θ) -(1/2)π +(2/3)(π/2 -θ) -(1/45)(π/2 -θ)^3
≒ 1/(π/2 -θ) -(1/2)π +(2/3 - 0.001)(π/2 -θ),
より2次方程式
(2/3 -0.001)(π/2 -θ)^2 -(3/2)π(π/2 -θ) +1 = 0,
これを解いて
(π/2 - θ) = 0.2189803
θ = 1.35181605
a = 1/cosθ = 4.603323
891:イナ
20/01/23 02:30:32.78 wc6308KN.net
前>>835
こっち側の縁から飛びこんで直角に飛びこむよりもショートカットするときの水中と縁の辺の比を1:tとすると、
x={x(1-t)/√2}/2+x√(1+t^2)/√2
2√2=1-t+2√(1+t^2)
t+2√2-1=2√(1+t^2)
t^2+2(2√2-1)t+9-4√2=4t^2+4
3t^2-2(2√2-1)+4√2-5=0
t={2√2-1+√(9-4√2-12√2+15)}/3
={2√2-1+√(24-16√2)}/3={2√2-1+2√(6-4√2)}/3
={2√2-1+2√(6-2√8)}/3
={2√2-1+2(√4-√2)}/3
={2√2-1+2(2-√2)}/3
={2√2-1+4-2√2)}/3
={2√2+3-2√2)}/3
=1
あれ？ 45°かぁ。
45°より60°のほうが速かったはず。あいだのt=4/7ぐらいでぎりぎり10.7秒台が出るか思たけど。
60°─x=10√6+30√2)/(2√6+3-√3)
=10.8516423……
45°─x=20(2√2+1)/7
=10.93835