- 1 名前：１３２人目の素数さん mailto:sage [2020/01/27(月) 20:12:01 ID:QSsw4R/8.net]
- 過去ログ置き場(1-16問目)
www3.tokai.or.jp/meta/gokudo-/omoshi-log/
まとめwiki
www6.atwiki.jp/omoshiro2ch/
1 cheese.5ch.net/test/read.cgi/math/970737952/
2 natto.5ch.net/test/read.cgi/math/1004839697/
3～6「datが存在しません。」
7 science2.5ch.net/test/read.cgi/math/1064941085/
8 science3.5ch.net/test/read.cgi/math/1074751156/
9 science3.5ch.net/test/read.cgi/math/1093676103/
10 science4.5ch.net/test/read.cgi/math/1117474512/
11 science4.5ch.net/test/read.cgi/math/1134352879/
12 science6.5ch.net/test/read.cgi/math/1157580000/
13 science6.5ch.net/test/read.cgi/math/1183680000/
14 science6.5ch.net/test/read.cgi/math/1209732803/
15 science6.5ch.net/test/read.cgi/math/1231110000/
16 science6.5ch.net/test/read.cgi/math/1254690000/
17 kamome.5ch.net/test/read.cgi/math/1284253640/
18 kamome.5ch.net/test/read.cgi/math/1307923546/
19 uni.5ch.net/test/read.cgi/math/1320246777/
20 wc2014.5ch.net/test/read.cgi/math/1356149858/
21 wc2014.5ch.net/test/read.cgi/math/1432255115/
22 rio2016.5ch.net/test/read.cgi/math/1464521266/
23 rio2016.5ch.net/test/read.cgi/math/1497416499/
24 rio2016.5ch.net/test/read.cgi/math/1502016223/
25 rio2016.5ch.net/test/read.cgi/math/1502032053/
26 rio2016.5ch.net/test/read.cgi/math/1518967270/
27 rio2016.5ch.net/test/read.cgi/math/1532793672/
28 rio2016.5ch.net/test/read.cgi/math/1540739963/
29 rio2016.5ch.net/test/read.cgi/math/1548267995/
なお、削除依頼は不要です。
※前スレ
面白い問題おしえて～な 30問目
https://rio2016.5ch.net/test/read.cgi/math/1572866819/
- 116 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 11:42:45 ID:MOGD/Do4.net]
- >>111
所要時間の式を偏微分して極値を出すのではないの？
- 117 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 12:12:04 ID:04w+XRU0.net]
- >>114
所要時間のなす関数は最大値を与える点で偏微分不可能です。
理由は二次元の場合と同じく、関数の定義にminが入るから。
明らかに無視できる経路を除いて最短経路になる候補が６個あり、所要時間=min{f1,f2,‥,f6}の形になる。
各々のfiは偏微分可能ですが、求める点はいずれのfiの極値にもなってはいません。
x=y=zに制限してもダメ。
手持ちの解答の方針としては
・まず６個に絞る。
・x=y=zに絞る。
・実質二個になる。
・min{f1,f2}の最大値は？
です。
６個に絞るのはめんどくさいだけ。
x=y=zに絞るところが手持ちの解はあまり綺麗でない。
以下は簡単。
- 118 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 12:19:10 ID:04w+XRU0.net]
- あ、ウソ言った。
・６個に絞る。
・実質２個に絞る。
・x=y=zに絞る
でした。
やってる事は東工大のと同じ。
- 119 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 12:37:52.30 ID:MOGD/Do4.net]
- >>114
個々のfをwolfram使って偏微分しようと思っていたけど無駄なんだな。
確かに自分のプログラムコードでもminを使っている。
- 120 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 12:38:51.80 ID:MOGD/Do4.net]
- >>112
話題は立方体に移っているよ。
- 121 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/03(月) 15:59:52 ID:avp8Qlns.net]
- 前>>112
>>113偏微分。
それだと思う！
入水角度θと監視員が最初にいる地点から対角線上にある救出場所までの距離xという2つの変数がある。
xが一次だから解けたのかもしれない。
- 122 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 16:24:28.58 ID:Bd06CPXX.net]
- >>81
単純化のためp≧q≧rとし
経路a: (0,0,0) -> (x,y,0) -> (p,q,r)
経路b: (0,0,0) -> (0,y,z) -> (p,q,r)
経路c: (0,0,0) -> (x,0,z) -> (p,q,r)
とするとき、経路aの所要時間
t=√(x^2+y^2)/2+√((p-x)^2+(q-y)^2+r^2)
の極小値(∂t/∂x=0,∂t/∂y=0)を計算すると
t=(√(p^2+q^2)+√3 r)/2,
(x/p=y/q=1-r/(√3 √(p^2+q^2))のとき)
で、これは経路a～cで最も小さい。
この所要時間はp=q=rのとき最大となり
t=((√2+√3)/2)p ----(1)
経路d: (0,0,0) -> (10,y,z) -> (p,q,r)
経路e: (0,0,0) -> (x,10,z) -> (p,q,r)
経路f: (0,0,0) -> (x,y,10) -> (p,q,r)
とするとき、経路aの所要時間
t=√((10+y)^2+z^2)/2 + √((10-p)^2+(q-y)^2+(r-z)^2), (y<zのとき)
t=√(y^2+(10+z)^2)/2 + √((10-p)^2+(q-y)^2+(r-z)^2), (y≧zのとき)
の極小値(∂t/∂x=0,∂t/∂y=0)を計算すると
t=(√3 (10-p)+√(q^2+(10+r)^2))/2,
((q-y)/y=(r-z)/(10+z)=(10-p)/(-(10-p)+√3 √((10+q)^2+r^2))のとき)
で、これは経路d～fで最も小さい。
この所要時間はp=q=rのとき最大となり
t=(√3 (10-p)+√(p^2+(10+p)^2))/2 ----(2)
(1)(2)を連立させて
√(p^2+(10+p)^2)=(√2+2√3)p-10√3
これを解くと
p=(5/6)(15-4√6+√(249-96√6))
のとき
t=(5/12)(3+√6)(5√3-4√2+√(83-32√6))
=11.69815627...
- 123 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 17:22:13.38 ID:lGSYI3JC.net]
- >>119
(1)(2)を連立させての意味が直ぐには理解できなかったのでグラフにしてみました。
https://i.imgur.com/GG3h127.jpg
- 124 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 19:21:36 ID:lGSYI3JC.net]
- wolframに
local minimum sqrt(x^2+y^2)/2+sqrt((p-x)^2+(q-y)^2+r^2) where 0<x<10 and 0<y<10
local minimum sqrt((10+y)^2+z^2)/2 + sqrt((10-p)^2+(q-y)^2+(r-z)^2) where 0<y<10 and 0<z<10 and y<z
local minimum sqrt(y^2+(10+z)^2)/2 + sqrt((10-p)^2+(q-y)^2+(r-z)^2) where 0<y<10 and 0<z<10 and y>=z
を入力したけど、どれも上手くいかなかった。
- 125 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 19:29:46 ID:lGSYI3JC.net]
- 所要時間はp=q=rのとき最大 というのが私には明らかでないので
座標をいれたら所要時間を計算する関数sim2を作ってコンピュータに最大値を探索させてみた。
探索を始める初期値によって収束しないこともあるので初期値を乱数発生させて収束したら表示するように設定。
> while(opt$convergence!=0){ # 初期値を乱数発生させて収束するまで繰り返す
+ opt=optim(par=sample(0:10,3),sim2,control = list(fnscale=-1),method='N')
+ }
> opt
$par
[1] 7.436222 7.436221 7.436221
$value
[1] 11.69816
$counts
function gradient
308 NA
$convergence
[1] 0
$message
NULL
コンピュータでの探索値では収束したらp=q=rになった。
- 126 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 19:38:43 ID:Bd06CPXX.net]
- >>122
>所要時間はp=q=rのとき最大 というのが私には明らかでないので
pを固定させてq,rをp≧q≧rの範囲で動かすことを考える。
このとき、所要時間はqまたはrの単調増加関数だから明らか。
- 127 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 21:19:57 ID:lGSYI3JC.net]
- >>123
立方体でなくて直方体のときも所要時間最大の点は
原点と最遠の頂点を結ぶ線上にあるのかな？
- 128 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 21:54:52.14 ID:lGSYI3JC.net]
- 数値を変えて
オリンピックサイズ・プール50m×25mの水の入ったプールの一つの角に監視員を置く。
水深2.5mとする。
この監視員は世界記録で移動するものとする。
水泳100m自由形 46秒91で水中を移動
陸上100m9秒58でプールを囲む面を」移動
この監視員がプールのどこへでも到達しうるには，最短で何秒必要か計算せよ。
をやってみたけど、最遠の頂点が一番時間がかかるという結果になったので面白みがなかった。
ただ、所要時間最大点はこの頂点と原点を結ぶ線上にあるという前提でのプログラムなので結果には自信がない。
- 129 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/03(月) 22:59:09 ID:avp8Qlns.net]
- 微分して極値を与える角度と距離だと思うんだよ。
／∥__`∥￣￣∥;;;;;;
∥∩∩ ∥ □ ∥;;;;;;
（(-_-)∥ ∥;;;;;;
（っ⌒⌒ ｡∥╂─╂
■`(_)_)ц~ ∥╂─╂
＼■υυ■___∥､＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼`前>>118＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼＼
- 130 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 23:19:00 ID:SKsq1rTN.net]
- >>125
> 水深2.5mとする。
この情報いる？
それはともかく、対角までの時間は、
75*0.0958=7.185
で、例えばプールの中心までは
(25-12.5tan(asin(9.58/46.91)))*0.0958+12.5cos(asin(9.58/46.91))*0.4691≒7.88
じゃないの？
> 所要時間最大点はこの頂点と原点を結ぶ線上にあるという前提
そんな根拠はない、というか間違いだろう
ぱっと考えられるのが、対角の2等分線上が考え付くが、それを採用するにも根拠がいる
- 131 名前：１３２人目の素数さん mailto:sage [2020/02/03(月) 23:19:19 ID:to5eQB6u.net]
- 陸上の速度をv、水中の速度をwとし、m=w/√(v^2-w^2)とする。
プールを0<x<a、0<y<bとする。
辺y=0から入水してt秒�
- 132 名前：繧ﾉ到達できる領域はmx+y≦mvt、
辺x=0から入水してt秒後に到達できる領域はmy+x≦mvt、
辺y=bから入水してt秒後に到達できる領域は-y+mx≦-b+mvt、
辺x=aから入水してt秒後に到達できる領域は-x+my≦-a+mvt
である。
方程式
mx+y=mvt‥?、my+x=mvt‥?、
-y+mx=-b+mvt‥?、-x+my≦-a+mvt‥?
において
???を連立して得られるtをt1、
???を連立して得られるtをt2とすれば到達時刻の最大値はmin{t1,t2}である。 [] - [ここ壊れてます]
- 133 名前： 【大吉】 mailto:sage [2020/02/04(火) 00:08:14 ID:+IjSdzOF.net]
- 前>>126
>>54修正。
向こう側の縁からθの角度で飛びこんだ監視員が泳ぐ時間は、
(10-x/√2)/sinθ(秒)
表記ミスがあった。計算が間違ってなければいいんだけど。
- 134 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 03:29:19 ID:W/1szoPy.net]
- >>127
z軸もあるから水深は必要。
- 135 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 03:33:34 ID:W/1szoPy.net]
- >>123
経路 a のt をqで偏微分すると
(q - y)/√((p - x)^2 + (q - y)^2 + r^2)
増加関数と言いるんだろうか？
- 136 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 04:18:31.12 ID:LNHsvcqa.net]
- >>131
そっちじゃなくて、tの極小値のほう
t=(√(p^2+q^2)+√3 r)/2,
これは明らかにqまたはrの増加関数
- 137 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 05:36:20 ID:W/1szoPy.net]
- >>127
立体だと複雑になるので平面で考えて
横２０ｍ縦１０ｍのプールで陸上速度毎秒２ｍ、水中速度毎秒１ｍで１５．５秒で到達できる範囲を描画してみました。
https://i.imgur.com/xZrdUpX.jpg
ご指摘の通り、対角線上に所要到達時間最大点があるというのは間違いであると確認できました。
- 138 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 05:49:23 ID:W/1szoPy.net]
- >>133
すいません、プログラムにバグを発見したので撤回します。m(__)m
- 139 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 06:22:55.67 ID:W/1szoPy.net]
- 気づいたバグを修正して長方形プールで描画しました。
対角線と対角二等分線をあわせて描画しました。
横２０ｍ縦３０ｍのプールで陸上速度毎秒２ｍ、水中速度毎秒１ｍで２６秒で到達できる範囲
https://i.imgur.com/vfPdQff.jpg
横３０ｍ縦２０ｍのプールで陸上速度毎秒２ｍ、水中速度毎秒１ｍで２６秒で到達できる範囲
https://i.imgur.com/WS21BMz.jpg
>127の直感通り、対角の２等分線上に所要時間最頂点が位置するようです。
- 140 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 06:27:09.62 ID:W/1szoPy.net]
- >>132
レスありがとうございます。
立法体なのでp≧q≧rという仮定が許されるということと理解しました。
- 141 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 07:24:13 ID:W/1szoPy.net]
- >81の問題を立方体から直方体に拡張して考えてみた。
オリンピックサイズ・プール50m×25mで水深2.5mの水の入った直方体プールの一つの角に監視員を置く。
この監視員は世界記録で直方体の面上や水中を移動するものとする。
水泳100m自由形 46秒91で水中を移動
陸上100m9秒58でプールを囲む面を」移動
この監視員がプールのどこへでも到達しうるには，最短で何秒必要か計算せよ。
立方体でなくて直方体のときには、所要時間最大の点は原点と最遠の頂点を結ぶ線上にはない、ということを教えていただいたのでプログラムを組み直した。
所要時間最大点の座標
par
[1] 49.980916 24.788643 2.288643
所要時間
$value
[1] 5.552414
という数値がでてきた。
- 142 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 07:48:42 ID:W/1szoPy.net]
- 探索初期値設定により、結果がばらつくけど
多数派意見(？）は
> opt
$par
[1] 49.06521 23.86881 1.36881
$value
[1] 5.855706
$counts
function gradient
256 NA
$convergence
[1] 0
$message
NULL
になった。
確かに、この方が到達時間が長い。
- 143 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 10:56:49 ID:3+QKrfHh.net]
- >>128
??の交点が頂点(a,b)にある角の二等分線上lなのでt1での???の交点もt2での???の交点もl上。
よくよく考えたらt1=t2だった。
- 144 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 11:21:43 ID:3+QKrfHh.net]
- >>139
ウソ書いた。
a,bの大小とt1,t2の大小は一致するでした。
- 145 名前：イナ mailto:sage [2020/02/04(火) 11:44:43.02 ID:+IjSdzOF.net]
- 前>>129問題(前スレ760)
向こう側の縁からθの角度で飛びこんだ監視員が泳ぐ時間は、
(10-x/√2)/sinθ(秒)
突き当たりを直角に曲がって飛びこむまでの時間は、
{10-(10-x/√2)(1+1/cosθ)}(1/2)
=5-(5-x/2√2)(1+1/cosθ)
=-5/cosθ+x/2cosθ√2+x/2√2(秒)
これらと最初のコーナーまでの5秒を足すと救出時間は、
(1/sinθ)(10-x/√2)-(1/2cosθ)(10-x/√2)+(1/2)(10+x/√2)
微分すると、
-1/sinθ√2+1/2cosθ√2+1/2√2=0
-2/sinθ+1/cosθ+1=0
-2cosθ+sinθ+sinθcosθ=0
(2-sinθ)cosθ=sinθ
cosθ=sinθ/(2-sinθ)
sin^2θ+{sinθ/(2-sinθ)}^2=1
sin^2θ(2-sinθ)^2+sin^2θ=(2-sinθ)^2
4sin^2θ-4sin^3θ+sin^4θ+sin^2θ=4-4sin^2θ+sin^4θ
4sin^3θ-9sin^2θ+4=0
sin^3θ-(3sinθ/2)^2+1=0
sin^3θ+1=(3sinθ/2)^2
3sinθ/2-√(sin^3θ+1)=0になるθを探すと、
θ=57.465773447629°のとき、
(3/2)(sin57.465773447629°)-√((sin57.465773447629°)^3+1)=0となり極小値を与える。
すなわち向こう側の縁から57.465773447629°の方向に飛びこむと最速で救出できる。
監視員が最初にいるコーナーの対角コーナーから、泳ぐ経路に引いた垂線が、向こう側の縁となす角は、
90°-57.46773447629°=32.53226552371°
対角線となす角は、
45°-32.53226552371°=12.46773447629°
監視員が10秒過ぎてから救出地点に達するまでに泳いだ時間t(秒)または距離t(m)は、
(10√2-10-t)sin12.465773447629°=t
(10√2-10)sin12.465773447629°=t+tsin12.465773447629°
t=(10√2-10)sin12.465773447629°/(1+sin12.465773447629°)
=0.735371693
到達時間10+t=10.735371693(秒)
- 146 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 12:03:49.77 ID:3+QKrfHh.net]
- xで微分してそれが0になるθ探してどーするん？
微分の意味がまるで分かってない。
結局意味もわからずやり方だけ覚えたらいいと思ってるから一つも前進しない。
- 147 名前：イナ mailto:sage [2020/02/04(火) 12:39:12.69 ID:+IjSdzOF.net]
- 前>>141問題(前スレ760)再考察。
救出する最遠方地点は監視員が最初にいる位置から対角線上x(m)にあると見て、向こう側の縁からθの角度で飛びこんだ監視員が泳ぐ時間は、
(10-x/√2)/sinθ(秒)
突き当たりを直角に曲がって飛びこむまでの時間は、
{10-(10-x/√2)(1+1/cosθ)}(1/2)
=5-(5-x/2√2)(1+1/cosθ)
=-5/cosθ+x/2cosθ√2+x/2√2(秒)
これらと最初のコーナーまでの5秒を足すと救出時間は、
(1/sinθ)(10-x/√2)-(1/2cosθ)(10-x/√2)+(1/2)(10+x/√2)
xで微分し、
-1/sinθ√2+1/2cosθ√2+1/2√2=0とすると、
-2/sinθ+1/cosθ+1=0
-2cosθ+sinθ+sinθcosθ=0
(2-sinθ)cosθ=sinθ
cosθ=sinθ/(2-sinθ)
sin^2θ+{sinθ/(2-sinθ)}^2=1
sin^2θ(2-sinθ)^2+sin^2θ=(2-sinθ)^2
4sin^2θ-4sin^3θ+sin^4θ+sin^2θ=4-4sin^2θ+sin^4θ
4sin^3θ-9sin^2θ+4=0
sin^3θ-(3sinθ/2)^2+1=0
sin^3θ+1=(3sinθ/2)^2
3sinθ/2-√(sin^3θ+1)=0になるθを探すと、
θ=57.465773447629°のとき、
(3/2)(sin57.465773447629°)-√{(sin57.465773447629°)^3+1}=0となり極小値を与える。
すなわち向こう側の縁から57.465773447629°の方向に飛びこむと最速で救出できる。
監視員が最初にいるコーナーの対角コーナーから、泳ぐ経路に引いた垂線が、向こう側の縁となす角は、
90°-57.46773447629°=32.53226552371°
対角線となす角は、
45°-32.53226552371°=12.46773447629°
監視員が10秒過ぎてから救出地点に達するまでに泳いだ時間t(秒)または距離t(m)は、
(10√2-10-t)sin12.465773447629°=t
(10√2-10)sin12.465773447629°=t+tsin12.465773447629°
t=(10√2-10)sin12.465773447629°/(1+sin12.465773447629°)
=0.735371693
到達時間10+t=10.735371693(秒)
- 148 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 12:54:31.79 ID:VWzue31P.net]
- >>143
直前のレス読んでるか？
xで微分してそれが0になるところ求めてどーするん？
それで何で所要時間が最小になるθが見つかるの？
微分というのが何か？
それで何故最小値が求まるのかという当たり前の理屈が分かってないから答えられないんだよ。
何度も解答見直した？
xで微分した。
＝0としてθについて解いた。
あれ？なんでコレで答え見つかるんだっけ？と自分に問い直してみないの？
- 149 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/04(火) 13:38:50 ID:+IjSdzOF.net]
- 前>>143
>>144前スレ760を見たら、なんで答えがみつかるかを説明せよとは問われてない。
ただ最短となる時間を計算せよとある。
だから計算した。そんなけ。AO入試ってなんだ？ と思って調べたら、論文みたいだった。答えはこうじゃないかああじゃないかと思案検討し計算する姿勢が求められてるんじゃないかと思う。
なんでxで微分して答えがみつかるか知りたい気もするし、べつに知りたくない気もする。
入水角度が60°のときも計算した。60°のときは計算しやすいけど最短でということでは角度が甘いと思った。
- 150 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 13:51:00 ID:3+QKrfHh.net]
- >>145
なんで答えが微分でもとまるか書けといわれてないから書かなくていい、分かってなくていいって思ってるからいつまで経ってもデキフるようにだけならないんだよ。
思案検討ってなんで微分したら答えがわかるという事は思案したの？
してないよね？
なーんにも考えてないよね？
なんとなく最小値求める時は微分。
でもθで微分なんてできない。
よーしxで微分してみよう！
おぉできた。
60°っぽいぞ！
きっとみんなの答えより正確なハズだ！
カッコいい！オレ！
‥‥
そういうのは思案とはいわん。
- 151 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 14:42:24.87 ID:VWzue31P.net]
- >>81です。>>119さんの解答がほぼ用意してた解答です。
ひとつだけコメント。
たとえば経路a: (0,0,0) -> (x,y,0) -> (p,q,r)において(x,y,0)についての極小値の出し方なのですが、
これは距離関数d(A,P)のPについての全微分が
d d(A,P)=e(A,P) dP (e(A,P)はAPベクトルと同じ向きの単位ベクトル、以下同じ）
になることを用いると意味がはっきりします。
この時の所要時間Tは(p,q,r)をAとおいて陸上の速度をv、水中の速度をwとして
T = d(O,P)/v + d(A,P)/w
なので
dT = (e(O,P)/v + e(A,P)/w)dP
となります。
これが任意のz=０内のdPについて0になるのはe(OP)/v + e(A,P)/wがxy平面の法線ベクトルと平行になるときで、
すなわちe(OP)/v + e(A,P)/wのxy平面への射影が0になるときです。
これはAxyから平面へおろした垂線の足HがOPの外分点であり、
かつe(A,P)をxy平面へ射影したものの長さがw/v=1/2となるとき、すなわち∠APHが60°となるときです。
よってこの場合PはHからOの方向へPH/√3だけ移動した点なので
f1(p,q,r)=(√(p^2+q^2)-r/√3)/2 + 2r/√3/1 = √(p^2+q^2)/2 + √3/2r
が経路aの極小値です。
経路b,cは文字入れ替えるだけ、経路dについては同様に考えて
f4(p,q,r)=√((10+q)^2+r^2)/2 + √3/2(10-p)
となります。
- 152 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/04(火) 15:38:31 ID:+IjSdzOF.net]
- 前>>145
>>146困ったら微分。
それしかない。
60°のときを考えるのはだれでもする。けどそのまま答えは60°のときとするのは高校生まで。
大人は困ったら微分する。
60°のときじゃない、と思ってθと置いたわけで、苦しんで微分するために置いたんじゃない。
未知数xの方程式を微分して極値を与える角度θを出したら60°になるよりちょっと手前で入水したらいいとわかった。問題ない。
- 153 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 16:10:36 ID:3+QKrfHh.net]
- >>148
ちょっと確認させて欲しい。
> 未知数xの方程式を微分して極値を与える角度θを出したら60°になるよりちょっと手前で入水したらいいとわかった。問題ない。
コレは本気で書いてるのか、それともココで引き下がったらレスバに負けるから間違ってるの承知でむりくり押し通してやろうと考えてるのかどっち？
もしかしてxで微分してもいいと本気で思ってるん？
xで微分しようがθで微分しようが好きな方で微分していいと本気で思ってるの？
- 154 名前：イナ mailto:sage [2020/02/04(火) 17:11:43.29 ID:+IjSdzOF.net]
- 前>>148
>>149どうやって解いたんだ？ と思って解きなおしたら何度やっても解けなくて、計算間違いかなぁと思ってあきらめかけた。
計算間違いじゃなくxで微分して極小値を与える角度θを出したんだとわかった。
一度はやろうとしたxとθの両方で微分するとどうなるか、またθで微分するとどうなるか、ぜひやってみてほしい。
xで微分して極小値を与える角度θを出して救出時間を出したのはまだ俺だけだと思う。今のところ正しいかどうか比べるものがない。なぜかみんな三次元がいいとか言って潜水してしまって、無人島にいる感じ。入水角度θ=60°のときより速いことは調べた。
- 155 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 17:27:27.89 ID:3+QKrfHh.net]
- >>150
だからxで微分しても正しい答えはでないと何度も指摘してるじゃん？
入水角が60°でない経路は最小にはなり得ません。
もし本気で出てる答え5+10/√3より小さい答えが出たと言い張るなら既出の答えの最小到達時間が最大になる点
(5(1+1/√3),5(1+1/√3))
=(7.886751345948, 7.886751345948)
に
5+10/√3 = 10.773502691896
より先に到達できる経路を明示しないとダメ。
わかる？明示？
要するにx=x.xxx‥の地点から入水したら10.77350‥より早く到達できるというx.xxx‥を一つでも見つければいい。
まぁやってごらんなさいな。
- 156 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 17:32:07.28 ID:W/1szoPy.net]
- >>148
びぶんのことはびぶんでやれ、という高木貞治を想い出したよ。
- 157 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 18:06:56 ID:VWzue31P.net]
- 7.886751345948/sin(57.465773447629deg)+7.886751345948(1-1/tan(57.465773447629deg))/2
=
10.7826518083
(10-7.886751345948)/sin(57.465773447629deg)+(10-(10-7.886751345948)(1+1/tan(57.465773447629deg)))/2+5
=
10.7759541902
いずれの経路でも 10.773502691896秒より前に到達できない。
- 158 名前：イナ mailto:sage [2020/02/04(火) 18:49:28.88 ID:+IjSdzOF.net]
- 前>>150
>>61到達時間10+t=10.735371693(秒)
＜10.7735……
入水角度θ(°)、到達時間10+t(秒)、あとは──。
>>151入水地点は、
つきあたりからの距離、
10-(10-x/√2)(1+cosθ/sinθ)に、
θ=57.465773447629°と、
xを代入するとわかる。
xは到達時間、
5+{10-(10-x/√2)(1+cosθ/sinθ)}(1/2)+(10-x/√2)(1/sinθ)=10.735371693にθ=57.465773447629°を代入し、
5+5-(5-x/2√2)(1+cos57.465773447629°/sin57.465773447629°)+(10-x/√2)(1/sin57.465773447629°)=10.735371693
=5+x/2√2+(x/2√2)(0.637910393)-5(0.637910393)+(10-x/√2)(1/sin57.465773447629°)=10.735371693
求めたxを代入すると入水地点もわかるはず。
- 159 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 18:54:24.56 ID:3+QKrfHh.net]
- こいついわれてる事全く理解してない。
真性のバカなんだな。
- 160 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 18:57:33.41 ID:W/1szoPy.net]
- wolframに∂t/∂x=0, ∂t/∂y=0を解いてもらおうと
x/(2 Sqrt[x^2 + y^2]) + (-p + x)/Sqrt[r^2 + (-p + x)^2 + (-q + y)^2]=0
y/(2 Sqrt[x^2 + y^2]) + (-q + y)/Sqrt[r^2 + (-p + x)^2 + (-q + y)^2]=0
を入力すると
r = -(1.73205 sqrt(p^2 + q^2) (p - x))/p, y = (q x)/pと返ってきてx,yについて解いてもらえなかった。
- 161 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 19:18:15 ID:VleZ36bS.net]
- xy平面において、x軸上の正の部分のみ、速度 v、その他の領域は速度 1 で移動できるものとする。
原点にいる人物が、目標地点(cosθ,sinθ) に到達すべく、移動する。
この時、より短時間で目標地点に到達するには、次の戦略αとβ、どちらが有利かを考える。
戦略α：現地点から、直接目標地点の方向へ速度 1 で移動する。
戦略β：x軸に沿って速度 v で移動する。
ε を正の小さな量とする。戦略αあるいはβ取って移動を開始し、εの時間がたった時のそれぞれの到達地点をA,Bとすると
A(εcosθ,εsinθ)、B(vε,0)
目標地点までの距離は、それぞれ、1-ε、√((vε-cosθ)^2+sin^2θ) となるが、さて、どちらが小さいか？
二乗したもの同士の差をとって比べてみると、
(1-ε)^2-((vε-cosθ)^2+sin^2θ) = 1-2ε+ε^2 -v^2ε^2+2vεcosθ-1 = ε(2v cosθ-2)+(1-v^2)ε^2
εは小さな正の量としているので、二次の項を無視すると、cosθ＞1/v で
1-ε＞√((vε-cosθ)^2+sin^2θ) となる。
つまり、目的地との方向のずれがθあるものの、v 倍の速度で移動できるとき、 cosθ＞1/v を
満たすなら、そのコースは直接目的地に向かうより有利である とえる。
この結論は、θとvのみが関与し、他の次元にも適用可。
- 162 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 19:18:48 ID:VleZ36bS.net]
- と同時に、この類いの問題に対し、次の戦略が最速であることを示す。
現在地から目標地点へのベクトル、あるいは、その方向への単位ベクトルをｐ↑、
選択可能ないくつかの速度ベクトルｖ↑が与えられたら、
内積 ｐ↑・ｖ↑ が最大になる速度ベクトルｖ↑ に沿うコースこそ最速コースである。
この戦略に従って、四次元プールの問題を考えるなら、微分は必要なくなる。
（この戦略の背景は、微分の考え方そのものであるが、結論のみを利用するならば、微分は不使用）
目的地を、(p,q,r) ただし、対称性から p≧q≧r として考える。
この方向への単位ベクトルは(p/D,q/D,r/D) 但し、D=√(p^2+q^2+r^2)
直接この方向へ向かう場合、速度ベクトルも(p/D,q/D,r/D)なので、内積は、1
縁を進む場合は、三つの平面の内どれか。p≧q≧r という条件では、平面z=0 上に、最適コースが存在し、
それは、(2p/d,2q/d,0) 但し、d=√(p^2+q^2)
時刻 t まで、移動したとき、(2pt/d,2qt/d,0)に移動しているので、目的地へのベクトルは (p-2pt/d,q-2qt/d,r)
速度ベクトルは(2p/d,2q/d,0)であり、この時、この両者の角度がπ/3だという方程式を解くと、
t=(1/2)d±((√3)/6)r が得られる。マイナスの方を代入して整理すると、残りの距離は((2√3)/3)rで、
トータル (1/2)d-((√3)/6)r+((2√3)/3)r=(1/2){√(p^2+q^2)+(√3)r} の時間がかかる
- 163 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 19:19:23 ID:VleZ36bS.net]
- 以上は、向こう側の「縁」を利用しない場合の最速コースについての議論。
向こう側の縁を利用する場合は、まずは、平面x=10へ下ろした時の足の座標、(10,q,r)へ向かうコースを考える。
立方体の表面しか移動できないので、展開図上で考えることになるが、直角を挟む2辺が10+rとqである直角三角形の
斜辺上にあたるコースを辿りながら、向こう側の平面に到達したときに、(p,q,r)を目指すことになる。
これは、無限に広がるプール、ただし、三つの平面x=0、y=0、z=0上だけは、
速度2で歩けるという条件で、(10+r,q,10-p)を目標にするのと同じ事になる。
こう考えると、先ほどの結果がそのまま使えて、このコースをとった場合のトータル時間は、
(1/2){√((10+r)^2+q^2)+(√3)(10-p)}
最も時間がかかる地点の座標には、(1/2){√(p^2+q^2)+(√3)r}=(1/2){√((10+r)^2+q^2)+(√3)(10-p)}
という条件が加わる。面倒になってきたので、細かいことは省略するが、上の式で、p=q=rとして
方程式を解くと、p=q=r=(5/6)(15-4√6+√(249-96√6))
（これは、>>119さんの結果と一致）
最後端追ったが、以上は、微分を使わない方法である。
- 164 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/04(火) 19
]
- [ここ壊れてます]
- 165 名前：:48:21 ID:+IjSdzOF.net mailto: 前>>154計算のつづき。
5+x/2√2+(x/2√2)(0.637910393)-5(0.637910393)+(10-x/√2)(1/sin57.465773447629°)=10.735371693
x=(5-3.189557196+11.8614066-10.735371693)/(0.838728105-0.353553391-0.225535382)
=11.309854(m)──救出地点までの直線距離
つきあたりから入水地点までの距離は、
10-(10-x/√2)(1+cosθ/sinθ)=10-(10-11.309854/√2)(1+0.637910393)
=6.71971502(m) [] - [ここ壊れてます]
- 166 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 19:52:12 ID:3+QKrfHh.net]
- (7.886751345948,7.886751345948)に10.773502691896秒以内に到達できる地点を探せと言われて7.886751345948の全く出てこない式を立てるのはどういう頭の構造してんの？
- 167 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 20:14:28 ID:W/1szoPy.net]
- >>158
p≧q≧r という条件では、平面z=0 上に、最適コースが存在し、
までは理解できるのですが、
入水する点の座標が
(2p/d,2q/d,0) 但し、d=√(p^2+q^2)
が最適とはどうして分かるのでしょうか？
- 168 名前：イナ mailto:sage [2020/02/04(火) 20:48:03.71 ID:+IjSdzOF.net]
- 前>>160
>>161
救出地点まで遠いほうの縁からの距離は、
11.309854/√2=7.99727446(m)
図を描いて8mぐらいかなぁと思ってたからいい値だと思った。
7.88……だと入水角度も入水地点も変わると思う。
7.88……がどうやって出た値かだよね。
xとθを両方とも微分するか、θで微分して、
x/√2=7.88……ってことなら、あるいはありうるかも。わるい値じゃない。
- 169 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 21:21:08 ID:3+QKrfHh.net]
- >>143
まぁしつこいからマジメにつっこむと
>これらと最初のコーナーまでの5秒を足すと救出時間は、
(1/sinθ)(10-x/√2)-(1/2cosθ)(10-x/√2)+(1/2)(10+x/√2)
xで微分し、
-1/sinθ√2+1/2cosθ√2+1/2√2=0とすると、
xで微分してそれが0になるθとはつまり到達地点(x,x)がどこにあろうと到達時間が一定であるようなθを探している事になる。
そんな地点は存在しないし実際wolfram大先生にグラフ書いてもらってもそんなθは存在してない。
https://www.wolframalpha.com/input/?i=-1%2Fsin%28x%29%2Bcos%28x%29%2B1%2F2%3D0&lang=ja
にもかかわらずどこからかコレが解
θ=57.46773447629
なる謎の数値を導き出す。
そしてこの謎の数値を元にした到達時間の最大値を出して、それが既出の数値より小さいから既出の値は間違ってると騒ぎ立てる。
そしてだったら既出の最大地点
(7.886751345948, 7.886751345948)
に既出の最小値10.773502691896より早く到達できる経路を明示してみろというと、この7.886751345948が全く出てこない式を立式して10.773502691896より小さいと言って得意顔。
バカさの次元の桁が違う。
- 170 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 21:23:45 ID:VleZ36bS.net]
- >>162
>> 入水する点の座標が
>> (2p/d,2q/d,0) 但し、d=√(p^2+q^2)
>> が最適とはどうして分かるのでしょうか？
なるほど、紛らわしい書き方をしてしまったようです。申し訳ありません。
(2p/d,2q/d,0) というのは、入水地点ではなく、速度ベクトルです。
原点から、この方向に、時刻0 から 時刻 t まで移動すると、
(2pt/d,2qt/d,0)
に到ります。この地点から、目的地をみると、(p-2pt/d,q-2qt/d,r)という方向にあります。
このまま、この速度を維持したまま、進んだ方がいいか、戦略をβからαに切り替えた方がよいか、
その判定に用いるのが、
「cosθ＞1/v」
という式です。
この式が不成立になる時刻を求めるための、方程式が
((p-2pt/d,q-2qt/d,r),(2p/d,2q/d,0)) =(1/2)*|(p-2pt/d,q-2qt/d,r)|*|(2p/d,2q/d,0)|
です。（左辺は内積の式であり、右辺は、ベクトルの大きさの積とcos(π/3)で構成されています。）
ここで求まった時刻を、(2pt/d,2qt/d,0) に代入すると、入水地点がわかります。
- 171 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 21:24:31 ID:THlBhxRo.net]
- >>143で救出までに最も長い時間
> 到達時間10+t=10.735371693(秒)
がかかる、と言っている点の座標はどこなん？
まあ、どこだろうが
> θ=57.465773447629°のとき、
の角度で行くより短時間のコースはあるわけだが
- 172 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/04(火) 23:15:55 ID:+IjSdzOF.net]
- 前>>163
>>166救出地点まで遠いほうの縁からの距離は、
11.309854/√2=7.99727446(m)
救出地点を座標でいうと、最初に監視員がいる地点を原点(0,0)、つきあたり方向にy軸をとり、
-xの方向に直角に曲がってy軸から6.71971502mの地点から、
θ=57.7465773447629°の方向に入水して、原点を出てから、10.735371693秒後に、
(x,y)=(-7.99727446,7.99727446)に到達する。
- 173 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 23:34:25 ID:3+QKrfHh.net]
- >>167
> θ=57.7465773447629°の方向に入水して、原点を出てから、10.735371693秒後に、
> (x,y)=(-7.99727446,7.99727446)に到達する。
じゃあその(7.99727446, 7.99727446)の地点に60°で見込む点
(10, 6.841000330368)
から入水して何秒
- 174 名前：かかるかちゃんと計算してみたかね？
その数値は10.735371693より大きいかね？
そういう当たり前の確かめを一つもしないからダメダメなんだよ。 [] - [ここ壊れてます]
- 175 名前：１３２人目の素数さん mailto:sage [2020/02/04(火) 23:47:26 ID:THlBhxRo.net]
- >>167
座標
> (x,y)=(-7.99727446,7.99727446)に到達する。
までの最短時間は
(10+7.99727446-(10-7.99727446)cot(π/3))/2+(10-7.99727446)/sin(π/3)/1≒7.3305
になり、、
(10+7.99727446-(10-7.99727446)cot(θ))/2+(10-7.99727446)/sin(θ)/1
θ=57.7465773447629°
の
> 10.735371693秒
より短いな
- 176 名前： mailto:sage [2020/02/05(水) 00:33:59.16 ID:C9wRmgDi.net]
- 前>>167
>>168第Ⅰ象限には水がないという設定です。
最速になる角度を探したんでほかの角度は60°と90°と45°ぐらい。
入水地点を決めてから角度を決めたんじゃなく、微分して角度が決まってから入水地点を計算した。
- 177 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 01:00:30 ID:gfGkl938.net]
- >>170
こんだけ言われてまだ何言われてるか理解できてないの？
どこまで頭悪いの？
みんなが60°で入水が最速である理由をあれだけ手を変え品を変えいろんな方法で示してたよね？
そのどれ一つとして理解できなかったとしても、そして自分が60°以外の角でより早い経路をみつけたとしても、最低限まず自分が見つけた地点に最速でいける方法がその角度なのか確かめてみろと言ってるんだよ。
なんでそんな簡単なことがわからん？
何よりそんな事まず自分で思いつかないの？
君のそのアポレスがどんだけスレの流れ乱してるからわからんの？
そのアポレスいつまで続けるん？
もう出てけよ。
- 178 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 01:06:11 ID:OkeImVJQ.net]
- 思付直感数学
- 179 名前：イナ ◆/7jUdUKiSM mailto:sage [2020/02/05(水) 01:55:16 ID:C9wRmgDi.net]
- 前>>170
問題見て最初に思いついたのがたしか60°だった。
縁と水中で速さが2:1だから。
その直感は正しいと思ってたけど、微分してθ=57.465773447629°と出て、到達時間を計算した。まだこの段階で半信半疑。
むしろ60°のとき計算したら10秒735切るぐらい速いはずと思って計算したら、
10秒9……って出て、あれ!? ってびっくりした。
θ=57.465773447629°のほうがθ=60°のときよりコンマ2秒速かった。
今は結果を受け入れてる段階。
- 180 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 06:23:22 ID:+pUSmyEU.net]
- >>165
解説ありがとうございました。
最後の方程式をWolframに解いてもらったら
人間技では扱えそうにない答になりました。
Solve[{p - 2 p (x/d), q - 2 q (x/d), r} . {2 (p/d), 2 (q/d), 0} == Norm[{p - 2 p (x/d), q - 2 q (x/d), r}] (Norm[{2 (p/d), 2 (q/d), 0}]/2), x, MaxExtraConditions -> Automatic]
x = (d^3 p^2 (q/d)^2 + d^3 p^2 (p/d)^2 + d^3 q^2 (p/d)^2 + d^3 q^2 (q/d)^2 - sqrt(-d^6 p^2 (r)^2 (q/d)^4 - 2 d^6 p^2 (r)^2 (p/d)^2 (q/d)^2 +
d^6 (-p^2) (r)^2 (p/d)^4 - d^6 q^2 (r)^2 (p/d)^4 - 2 d^6 q^2 (r)^2 (p/d)^2 (q/d)^2 - d^6 q^2 (r)^2 (q/d)^4 + 4 d^4 p^4 (r)^2 (q/d)^2 + 4 d^4 p^4 (r)^2 (p/d)^2 + 8 d^4 p^2 q^2 (r)^2 (p/d)^2 + 8 d^4 p^2 q^2 (r)^2 (q/d)^2 + 4 d^4 q^4 (r)^2 (p/d)^2 + 4 d^4 q^4 (r)^2 (q/d)^2) - 4 d p^4 -
8 d p^2 q^2 - 4 d q^4)/(2 (d^2 p^2 (q/d)^2 + d^2 p^2 (p/d)^2 + d^2 q^2 (p/d)^2 + d^2 q^2 (q/d)^2 - 4 p^4 - 8 p^2 q^2 - 4 q^4))
x = (d^3 p^2 (q/d)^2 + d^3 p^2 (p/d)^2 - d^3 q^2 (p/d)^2 + d^3 q^2 (q/d)^2 + sqrt(-d^6 p^2 (r)^2 (q/d)^4 - 2 d^6 p^2 (r)^2 (p/d)^2 (q/d)^2 +
d^6 (-p^2) (r)^2 (p/d)^4 - d^6 q^2 (r)^2 (p/d)^4 - 2 d^6 q^2 (r)^2 (p/d)^2 (q/d)^2 - d^6 q^2 (r)^2 (q/d)^4 + 4 d^4 p^4 (r)^2 (q/d)^2 + 4 d^4 p^4 (r)^2 (p/d)^2 + 8 d^4 p^2 q^2 (r)^2 (p/d)^2 + 8 d^4 p^2 q^2 (r)^2 (q/d)^2 + 4 d^4 q^4 (r)^2 (p/d)^2 + 4 d^4 q^4 (r)^2 (q/d)^2) - 4 d p^4 -
8 d p^2 q^2 - 4 d q^4)/(2 (d^2 p^2 (q/d)^2 + d^2 p^2 (p/d)^2 + d^2 q^2 (p/d)^2 + d^2 q^2 (q/d)^2 - 4 p^4 - 8 p^2 q^2 - 4 q^4))
- 181 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 08:43:38.73 ID:t1CV2afM.net]
- >>174
なぜ FullSimplify しない？
X=の最初の式を%とすると
FullSimplify[%, d > 0 && p > 0 && q > 0 && r > 0]
1/6 d (3 + (Sqrt[3] r)/Sqrt[p^2 + q^2])
- 182 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 09:26:37.46 ID:VrbXRcrj.net]
- >>174
165です。これは自戒を含めてのコメントになりますが、あの方程式は、手で簡単に計算できます。
お試しください。
- 183 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 09:41:51 ID:PzHdrrq1.net]
- >>175
ありがとうございます。
その機能をはじめて知りました。
- 184 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 14:10:45 ID:VrbXRcrj.net]
- >>174
「お試しください」と書きましたが、実際にお示しします。
あの戦略からの要請、二つのベクトル、P-Vt と V のなす角度がπ/3であるという方程式は
(P - V t).V=(1/2)*|(P -V t)|*|V|
と書けます。ピリオドはベクトルの内積、絶対値記号はノルムを表す記号としてます。
>>165では、無理矢理成分表示で、式を表していたため、見苦しくなりましたが、最初からこう書けばよかったですね。
|V|=2、P.V=p*(2p/d)+q*(2q/d)+r*0=2d、P.P=p^2+q^2+r^2=d^2+r^2 に注意して変形すると
P.V-t*V.V = |P -V t|
2d-4t = √(P.P-2t*P.V+4t^2)
16t^2-16td+4d^2=d^2+r^2-4td+4t^2
12t^2-12td+3d^2-r^2=0
t=(1/12){6d±√(36d^2-12(3d^2-r^2))}=(1/12){6d±(2√3)r}
と、言う具合に、簡単に t を求めることができます。
- 185 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 14:25:04 ID:t1CV2afM.net]
- 二次元平面上に無限に続く、1オームの抵抗で作られた正方形の格子において、
ナイトの動き(桂馬飛び)の位置にある2つのノード間の抵抗は
4/π-1/2 オームであることを示せ。
（Google入社試験 - 難易度を下げるために一部簡単化）
- 186 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 14:36:16 ID:298bnSpu.net]
- >>179
コレは電気抵抗の知識なくても解けるの？
Googleの試験だからそこは知らなくても推定しろなのかな？
とりあえずググったら長さに比例して断面積に反比例するというのしか見つからない。
https://kenkou888.com/category21/dousen_teikou.html
- 187 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 14:38:21 ID:298bnSpu.net]
- あれ？
格子点と格子点を結ぶように1Ωの抵抗が繋がってるという意味？
もしかして？
- 188 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 14:55:22 ID:t1CV2afM.net]
- >>181
そうです。
>>179
の補足ですが、1オームの二次元無限格子の隣接ノード間の抵抗は
対称性の意味を知っていれば中学生で出せます。
より一般的には、任意の二つのノード間の抵抗は
有理数＋有理数×1/πであらわされることを示してください。
- 189 名前：１３２人目の素数さん [2020/02/05(水) 15:02:45 ID:t1CV2afM.net]
- >>180
前提となる物理知識は、中学生レベルのオームの法則とキルヒホッフの法則のみです。
- 190 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 15:28:22 ID:298bnSpu.net]
- つまりijにおける電位をe[i,j]として(0,0)から-1A、(2,1)に+1A流入してるとして
e[i+1,j]+e[i-1,j]+e[i,j+1]+e[i,j-1]-4e[i,j]
=δi0δj0-δi2δj1
のときのe[2,1]-e[0,0]かな？
留数定理の香りがする。
- 191 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 22:11:30 ID:+pUSmyEU.net]
- >>178
どうもありがとうございました。
d=√(p^2+q^2)の情報なしでwolframに入力したので複雑な答で表示されたのだと理解しました。
- 192 名前：１３２人目の素数さん mailto:sage [2020/02/05(水) 23:39:40.39 ID:t1CV2afM.net]
- >>184
ヒント
ローラン展開による母関数
E(z,w)=Σ[i,j:整数] e(i,j) z^i w^j
- 193 名前：イナ mailto:sage [2020/02/06(木) 04:43:23.90 ID:Mv+y98sK.net]
- 前>>173だれか入水角度60°の決め打ちじゃなくて、微分してみたって人いないかなぁ。
- 194 名前：１３２人目の素数さん mailto:sage [2020/02/06(木) 06:18:07.30 ID:Ya801udz.net]
- >>187
前スレで
https://rio2016.5ch.net/test/read.cgi/math/1572866819/807
が偏微分で極値を出している。
プログラムでの数値解と合致した。
立方体の方の計算�
- 195 名前：ﾉうつったら。
オリンピップールの直方体の方が計算のしがいがあると思う。 [] - [ここ壊れてます]
- 196 名前：１３２人目の素数さん mailto:sage [2020/02/06(木) 09:32:28.14 ID:tNI6h0TT.net]
- >>188
前スレの807を書いた者だが、極値は二つ出たが、807では採用する方を誤ってしまった。
訂正内容を824に記してあるので、807を見る場合は、824もセットで見て欲しい。
- 197 名前：１３２人目の素数さん mailto:sage [2020/02/06(木) 10:42:10.03 ID:5WVjoOPr.net]
- >>187
偏微分以外は全部決め打ちと思ってる時点でもうこのスレでレスできるレベルに到達してない。
- 198 名前：１３２人目の素数さん mailto:sage [2020/02/06(木) 22:27:39.52 ID:eS4p1xAB.net]
- > だれか入水角度60°の決め打ちじゃなくて、微分してみたって人いないかなぁ。
イナ以外で60°という角度を使っている人は、思い付きだけで使っているわけでなく、
書くまでもなく計算したり、スネルの法則等の定理を用いて60度を導出しているんだからな
- 199 名前：１３２人目の素数さん mailto:sage [2020/02/06(Thu) 23:33:31 ID:Ya801udz.net]
- タクシー料金の改訂
# 京浜地区
# 旧運賃（小型)
F1=740 # 初乗運賃 Fair
D1=2000 # 初乗り距離 initial Distance
C1=90 # 加算運賃 Charge by distance
B1=288 # 加算距離 charge By distancce
# 新運賃
F2=500
D2=1200
C2=100
B2=264
https://travel.watch.impress.co.jp/img/trw/docs/1224/505/02_o.jpg
距離と新旧運賃および差額をグラフにしてみた。
運賃改定率が8.88％と記載されているのだがどうやって計算するんだろう？
- 200 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 01:43:45 ID:YN6u30Ej.net]
- >>187
横に１０ｍ走って縦に方向を変えてプールサイドからθの角度で座標(p,q）に向かって飛び込む時の所要時間は
10/2+((p-10)/tan(θ)+q)/2+sqrt((10-p)^2+((10-p)/tan(θ))^2)
角度を決めたら縦方向の走行距離が決まってしまう。
これを微分すればいい
D[5 + (q + (-10 + p) Cot[θ])/2 + Sqrt[(10 - p)^2 + (10 - p)^2 Cot[θ]^2], θ] をWolfram先生にお願いすると
導関数は((10 - p) Csc[θ]^2)/2 - ((10 - p)^2 Cot[θ] Csc[θ]^2)/Sqrt[(10 - p)^2 + (10 - p)^2 Cot[θ]^2]
んでもって
solve ((10 - p) Csc[θ]^2)/2 - ((10 - p)^2 Cot[θ] Csc[θ]^2)/Sqrt[(10 - p)^2 + (10 - p)^2 Cot[θ]^2]==0 for θ
導関数が０になるθを求めてもらうと
θ = π/3 θ = -π/3
マイナスだとプールに飛び込めないから、θ = π/3
目的の座標に関わりなく６０°と算出されました。
- 201 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 02:37:34 ID:JwTQ0wHH.net]
- >>192
8.88%をだしたいのなら、例えば、
100m利用、200m利用、...、6600m利用、6700m利用
の料金の合計を、新旧で比較すると、8.8位のアップになる。
1200m以下だと、新運賃は240円安い
1800m位から、逆転し、その後、じわじわ差が大きくなり、
4200m位から、240円位高くなる。
100mから4200m位をまんべんなく利用する人がいたとすると、この改定により、
利用額の増減はほとんど無いという解釈も可能。
距離が大きくなれば、値上げの効果がどんどん大きくなる。
最終的には、1m辺りの加算運賃の比 90/288 ： 100/264 = 33:40
なので、21.212121...%の上昇に近づく。
それが、6.7km辺りでは、8.8%だというだけ。
つまり、8.88位になるよう、最小距離100mと最大距離を6.7kmを恣意的に選んだだけ。
文頭の説明には説得力は全く無い。
恐らく、距離別利用割合のデータに基づいて、新旧の料金比較したのだろう。
この情報が無ければ、8.88%等の数値は出せないと思われる。
- 202 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 03:19:17 ID:9IJwzjmO.net]
- >>179
でけたかも。
まず(0,0)以外で漸化式
4e(i,j)=e(i+1,j) + e(i-1,j) + e(i,j+1) + e(i,j-1)
を満たす列を探す。
e(i,j)=∫[|x|,|y|<π] (1-cos((x+y)i)cos((x-y)j))/(1-cosxcosy)dxdy
がこの条件を満たす。
また|i|,|j|→∞で0に行く。
そこで点i,jに電荷はe[i,j]-e[2-i,1-j]となる。(多分解は一意、ノーチェック)
e[0,0]=0, e[2,1]=32π-4π^2
であるから電位差は64π-8π^2。
e[1,0]=4π^2だから原点から隣接する4点に計16π^2の電流が流れる。
よって求める抵抗値は(32π-4π^2)/16π^2=4/π-1/2である。
またe[i,i]が
e[i,i]=∫(1-cosix))/(1-cos(x)cos(y))dxdy
であるが、yについて先に積分すると
e[i,i]=π∫(1-cosix))/|sin(x)|dx
となり、この値はπの有理数倍になる。
コレと漸化式によりe[i,j]はπとπ^2の有理係数の線形結合である。□
e[i,i]の計算が全く思いつかなかった。
e[i,j]の母関数って作れるのかな？
- 203 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 08:18:59 ID:YN6u30Ej.net]
- >>194
レスありがとうございます。
距離と新旧運賃と差額のグラフのアップロードを忘れておりました。
https://i.imgur.com/bElGqRG.jpg
与えられたデータだけからは平均値上げ率は算出できない思っていたのが確認できました。
ある距離までの乗客数が同じと仮定したときの平均の値上げ率をグラフにすると
https://i.imgur.com/2PkwY3n.jpg
> which.min((crs-0.0888)^2)
[1] 6869
> pir(6869)
[1] 0.08882413
6．9キロくらいの平均で8.9％の値上げ率になりました。
計算したひとはこういう数字を使ったのでしょう。
- 204 名前：イナ mailto:sage [2020/02/07(金) 08:40:02.32 ID:VtLCtPNo.net]
- 前>>187
>>193
角度を決めたら泳ぐ距離が決まる。θで微分するか、θとxの両方で微分するかってとこですか。
- 205 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 10:21:34.08 ID:YN6u30Ej.net]
- >>197
https://i.imgur.com/reg8DOz.jpg
Oから出発してAを経て角度θで入水してS(p,q)に泳ぐとする
AJの長さをxとすると
tan(θ)=(10-p)/(q-x)だから
x=q-(10-p)/tan(θ)
となり、
所要時間の計算からxは消去できて
10/2+((p-10)/tan(θ)+q)/2+sqrt((10-p)^2+((10-p)/tan(θ))^2)
となる。
この極値を与えるθはp,qによらないのは>193に書いた通り。
- 206 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 15:11:11.89 ID:YN6u30Ej.net]
- >>197
>角度を決めたら泳ぐ距離が決まる。
違う、角度を決めたら走る距離も泳ぐ距離も決まる
- 207 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 15:12:19.07 ID:YN6u30Ej.net]
- 走る距離 10+((p-10)/tan(θ)+q)
泳ぐ距離 sqrt((10-p)^2+((10-p)/tan(θ))^2)
- 208 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 20:35:27.63 ID:CyUpE86n.net]
- >>179の解の一意性の証明ができないなぁ。
昔これエレガントな解答を求むかなんかで
4a[ij]=a[i+1j]+a[i-1j]+a[ij+1]+a[ij-1]
をみたす有界な列は定数に限る事を示せ
の形で出題されて2chにえらいエレガントな解答が上がって数セミに載ったっていう事件があったけど、あれどんな証明でしたっけ？
誰か覚えてます？
- 209 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 20:48:28.93 ID:oSzq3jEL.net]
- >>195
正解です。よく特殊解を探せましたね。
その特殊解を(i',j')個すらして符号を変えて重ね合わせて正規化すれば、
2点(i',j')--(0,0)間に1Aの電流を流した時の電位の式が出ます。
想定していた解答は、2点(2,1)--(0,0)間に1Aの電流を流した時の電位の式
e[i+1,j]+e[i-1,j]+e[i,j+1]+e[i,j-1]-4e[i,j]=δi0δj0-δi2δj1
にexp(√-1 (ix+jy))をかけてi,jで和を取ると
(exp(-√-1 x)+exp(√-1 x)+exp(-√-1 y)+exp(√-1 y)-4)E(x,y)=1-exp(√-1 (2x+y))
(ここで E(x,y)=Σ[i,j:整数] e[i,j]
- 210 名前：exp(√-1 (ix+jy)) と置く)
より
E(x,y)=(1-exp(√-1 (2x+y)))/(2cosx+2cosy-4)
これをフーリエ級数の公式(留数定理)
e[i,j]=(1/(2π)^2)∫[0,2π]∫[0,2π]E(x,y)exp(-√-1 (ix+jy))dxdy
を用いて逆変換すると、(2,1)--(0,0) 間の電位は
e[2,1]-e[0,0]=(1/(2π)^2)∫[0,2π]∫[0,2π](exp(-√-1 (2x+y))-1)(1-exp(√-1 (2x+y)))/(2cosx+2cosy-4) dxdy
=(1/(2π)^2)∫[0,2π]∫[0,2π] (1-cos(2x+y))/(2-cosx-cosy) dxdy
=4/π-1/2
一般に(0,0)--(i,j)間の抵抗値は
(1/(2π)^2)∫[0,2π]∫[0,2π] (1-cos(ix+jy))/(2-cosx-cosy) dxdy
=(1/(8π^2))∫[0,2π]∫[0,2π] (1-cos(i(x+y)+j(x-y)))/(1-cosxcosy) dxdy
=(1/(8π^2))∫[0,2π]∫[0,2π] (1-cos((i+j)x)cos((i-j)y))/(1-cosxcosy) dxdy
=(1/(2π))∫[0,π] (1-cos((i+j)x)((1-|sinx|)/cosx)^|i-j|)/|sinx| dx
(0,0)--(i,i)間の抵抗値は
(1/(2π))∫[0,π] (1-cos((i+j)x))/|sinx| dx
=(2/π)(1+1/3+1/5+1/7+...+1/(|2i|-1)) [] - [ここ壊れてます]
- 211 名前：１３２人目の素数さん mailto:sage [2020/02/07(金) 21:19:32.21 ID:CyUpE86n.net]
- >>200
その解は
f(x,y)=Σe[kl]exp(ikx+ily)
が収束すると仮定して(仕事率の有限性から二乗は収束する)みたすべき関数方程式で見つけました。
見つかっちゃえば解答はコレが解だでいいハズなんですが、級数の収束性とかが自明でないのでコレのみが解なのか示せてなくて気持ち悪い。
まぁ入社試験ではそこまで求められないんだろけど。
>>201のエレガントな解答求むのやつは確か与式が等号でなくて不等号だったかな？
しかしそこから等号の有界な非自明解の存在が必要性で出てきて矛盾を導くという流れだったような。
検索しても出てこないなぁ？
- 212 名前：１３２人目の素数さん mailto:sage [2020/02/08(土) 00:08:01.33 ID:aj0WebTe.net]
- >>202
補足
特殊解(原点に8π^2の電流を注入した解)
∫[0,2π]∫[0,2π] (1-cos((i+j)x)cos((i-j)y))/(1-cosxcosy) dxdy
=4π∫[0,π] (1-cos((i+j)x)((1-|sinx|)/cosx)^|i-j|)/|sinx| dx
は|i|=|j|で 4πΣ[n=1,|i|] 1/(2n-1) になって |i|=|j|→∞で発散します。
一方、2点間(0,0)--(2,1)で符号を変えて重ね合わせた解((0,0)--(2,1)間に8π^2の電流を流した解)
∫[0,2π]∫[0,2π](cos((i+j)x)cos((i-j)y)-cos((i+j-3)x)cos((i-j-1)y))/(1-cosxcosy)dxdy
=4π∫[0,π](cos((i+j)x)((1-|sinx|)/cosx)^|i-j|-cos((i+j-3)x)((1-|sinx|)/cosx)^|i-j-1|)/|sinx|dx
はリーマン・ルベーグの補題より|i|,|j|→∞で0に収束します。
ちなみに、3次元の無限格子ではこのような発散は起こりません。
- 213 名前：１３２人目の素数さん mailto:sage [2020/02/08(土) 07:44:31.25 ID:MW5Whxwa.net]
- 隔離中のクルーズ船では
船内の換気が共通らしいから13日後に発症した奴がいるとその近くの部屋のやつはプラスで14日隔離しないといけない
それが今の船内の状況という。
こんな問題を考えてみた（答は自分でもまだ持ってません）。
計算には追加の設定がいるかもしれません。
両隣のどちらかが感染したら14日延長、どの部屋も1日で感染する確率pは1%
部屋の配置は長方形（つまり始まりも終わりもなし）。
発症するか、隔離期間が終われば下船できる。全員定員１の個室として客と乗務員を合わせた人数nは3000人。
クルーズ船から全員下船できる日数の期待値は？
- 214 名前：１３２人目の素数さん mailto:sage [2020/02/08(土) 08:08:03 ID:MW5Whxwa.net]
- >>190
正解が出てから誤答を連投する芸人をどう納得させるかというゲームだと思って俺は楽しんでる。
- 215 名前：１３２人目の素数さん mailto:sage [2020/02/08(土) 08:11:36 ID:MW5Whxwa.net]
- 部屋の配置は長方形（つまり始まりも終わりもなし）。
どの部屋にも両隣があると言う意味。
長方形に意味なし、円形配置でも同じ。
- 216 名前：１３２人目の素数さん mailto:sage [2020/02/08(土) 11:46:54 ID:MU+ZFKMw.net]
- もし 4a[i,j] = a[i,j+1]+a[i,j-1]+a[i+1,j]+a[i-1,j] =: 4(Δa)[i,j] が全ての(i,j)で成り立つ状況なら、
閉領域[-n,n]^2に属する格子点(i,j)全体におけるa[i,j]の最大値、最小値は、
どちらも辺∂([-n,n]^2)上でとらなければならないから、一意性はこれが鍵になったりするのかなあ
状況的にはリウヴィルの定理に似てるから、その証明と同じ手法が使えたりはしないだろうか
|
 |
