分からない問題はここに書いてね448

さあ、今日も１日がんばろう★☆

前スレ
分からない問題はここに書いてね447
https://rio2016.5ch.net/test/read.cgi/math/1537106483/

削除依頼を出しました

しつもんいいすか

いいよ。
でも完全性定理と数理論理のモデル理論はダメだからね。

なぜですか？

面白くないからに決まってるジャン

画像貼るんでちとまってください

おう、
待つよ！

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

https://i.imgur.com/81mXtaZ.jpg

1枚目の問題の解答が2枚目なのですが、途中に出てくるmaxは正しくはminではないでしょうか

その通り

∫1→∞　1/x^2　dx
について
1→∞[-2/x]
1→∞[0+2]=2
とするんですが違和感がありすぎます
x=1のとき高さが2であるとしかこれは言っていなくて
積分ならば当然面積なので
�凾ｱれね直角三角形の面積
高さが2だから面積が2なんて数学はなくて
当然(底辺×高さ)／２ですよね
しかるに
∫1→∞　1/x^2　dx
は
底辺∞高さ2の直角三角形だから
(∞×2)/2=∞
とこうなるんじゃないのん？
おらが間違えたんならどう間違ったか解説してくれんかのう

そもそも１だし

高さ1底辺∞の直角三角形＝高さ1底辺∞の帯状領域（直方体）
だし
積分で面積を求めた領域≒（x=1付近を除けば）点(0,1)から延びる半直線
だから
違和感持つ方が無理だろ

> 1→∞[0+2]=2
左辺の変な記号が一つ上の行と同じ意味なら右辺は0だなｗ

そかそか

∫1→∞　1/x^2　dx
について
1→∞[-1/x]
1→∞[0+1]=1
とするんですが違和感がありすぎます
x=1のとき高さが1であるとしかこれは言っていなくて
積分ならば当然面積なので
�凾ｱれね直角三角形の面積
高さが1だから面積が1なんて数学はなくて
当然(底辺×高さ)／２ですよね
しかるに
∫1→∞　1/x^2　dx
は
底辺∞-1高さ1の直角三角形だから
((∞-1)×1)/2=∞
とこうなるんじゃないのん？
おらが間違えたんならどう間違ったか解説してくれんかのう

∫1→∞　1/x^2　dx
の面積って1なの？∞なの？どっち？

dx＝1として計算してみるよ
1→∞[-1/x]

x=1　ｙ＝1
x=2　ｙ＝0.5
x=3　ｙ＝0.33
x=4　y=0.25

ここまでの面積
1+0.5+0.33+0.25=2.08

∫1→∞　1/x^2　dx
の面積は1にはどうしても思えなくて∞だろ

ああ、ちゃうな

dx＝1として計算してみるよ
y=1/x^2

x=1　ｙ＝1
x=2　ｙ＝0.25
x=3　ｙ＝0.11
x=4　y=0.0625

ここまでの面積
｜1+0.25+0.11+0.0625｜=1.4225

∫1→∞　1/x^2　dx
の面積は1にはどうしても思えなくて∞だろ

だから無限遠の重力ポテンシャルをゼロにしたなんちゃら宇宙速度
ってそもそも面積∞だから使い物にならんよねえ

さあ、どうやって誤魔化しますかぁ？

三角形の合同条件3つが合同条件になる証明ができない高校生です笑笑
検索しても出なかったので詳しい方お願いできないでしょうか。
メモします

それは例えばヒルベルトの公理系からスタートするのかR^2の座標とっていいのかでも話がだいぶ違うな。
後者でいいん？

>>16
> x=1　ｙ＝1
> x=2　ｙ＝0.25
> x=3　ｙ＝0.11
> x=4　y=0.0625
>
> ここまでの面積
> ｜1+0.25+0.11+0.0625｜=1.4225

お前が足してるのは左上の点が1/x^2のグラフ上にある長方形の面積だから
1/x^2のグラフから大いにはみ出てる（積分を上から評価することはできる
右上の点がグラフ上にあるようにするなら最初の1x1の長方形は範囲外だから足したらダメ
そうするとお前の示したのは 0.4225 < 求める積分 < 1.4225 になるということだけ

そもそも一般に
∫1→∞　1/r^2　dr　＝　1
と
∫0→1　1/r^2　dr　は相似なのに　＝　∞　
↑
なんで1と∞を混在して採用しているわけさ？
ご都合主義なん？

相似なんだから1か∞に統一すべきじゃないのか？

だから、結局これは

∫0→∞　1/r^2　dr　＝　3　or　∞

3なの？∞なの？

∫1→∞　1/x^2　dx　＝　1
とする
∫0→1　1/x^2　dx　＝　∞
だけど
∫1→∞　1/x^2　dx　＝　1
と縦横入れ替わっただけで相似だから
∫0→1　1/x^2　dx　＝　1
とする
それに1×1=1
ゆえに
∫0→∞　1/x^2　dx　＝　1+1+1　=　3

こういうことにはならんのかえ？

ちょっと違うな

∫1→∞　1/x^2　dx　＝　1
とする
∫0→1　1/x^2　dx　＝　∞
だけど
∫1→∞　1/x^2　dx　＝　1
と縦横入れ替わっただけで相似だから 1×1=1の部分を足して
∫0→1　1/x^2　dx　＝　2
とする
ゆえに
∫0→∞　1/x^2　dx　＝　1+2　=　3

こういうことにはならんのかえ？

そもそも1/x^2のグラフはｙ軸に対して対称であって縦横が相似じゃないぞ

老子とプリンストン大学数学科の教授の中で断然トップの人はどっちの方が頭が良いですか？

分子が1、分母がn桁の正整数である有理数全体からなる集合をS_nとする。
S_nの要素のうち、循環節の長さを最小とするものを1つ取り、その長さをm[n]とする。同様に循環節の長さを最大とするものについてその長さをM[n]とする。

（1）m[n]を求めよ。

（2）以下を示せ。
(a) lim[n→∞] m[n]/M[n] = 0
(b) M[n]≦M[n+1]
(c) M[n]<10^n

(10,a) = 1のとき
1/a の循環節の長さ = 10の Z/aZの乗法群での位数。
とくにそれはaより小さいからa<10^nのとき
1/a の循環節の長さ < 10^n。
またa|bのとき
1/a の循環節の長さ≦1/b の循環節の長さ。
pを素数としてa = p^e、vをp進付値mを10の Z/pZの乗法群での位数とするとき
v(10^(mn) −1) = v(10^m−1)+v(n)
により10のZ/aZの乗法群での位数はmp^(e-v(10^m-1))。
特にp = 7のとき10のZ/(p^e)Zの乗法群での位数は6・7^(e-1)。
10^(n-1)<7^e<10^n であるn,eをとるとき1/7^eの循環節の長さは
6・7^(e-1)であり特に
M[n] ≧ 6・7^(e-1) > 6/7 10^(n-1)。

>>26
y=1/x^2
∫0→1　1/x^2　dx　＝　∞

y=1/x^2
yとx入れ替えて
x=1/y^2
y=1/√x
∫1→∞　1/√x　dx　＝　1
1→∞[2√x]=1→∞[∞-1]=∞

あれ？なんでこっちは収束しないん？

>>26
y=1/x^2
∫0→1　1/x^2　dx　＝　∞

y=1/x^2
yとx入れ替えて
x=1/y^2
y=1/√x
∫1→∞　1/√x　dx　＝　1
1→∞[2√x]=1→∞[∞-2]=∞

あれ？なんでこっちは収束しないん？

これの18問ってどうやって解けば良いの？
http://www.ms.u-tokyo.ac.jp/kyoumu/b20170524.pdf

前スレの992
点T(1,t)で円2つが交わるとすれば線分OTの垂直二等分線の第一象限で切り取られた部分が2円の中心間距離l。
l=(t^2+1)^(3/2)/2tはすぐ出てくるのであとは微分してください
おわり

前スレ993

2∫0→t (a^2-2(a-√(a^2-t^2))^2)dtで出てくるやろ

>>32
最近はどうか知らんが
内部の学生有志で解答作ってないの？

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

宿題の答えを聞いているような感じ

>>33

〔前スレ．992〕
　xy平面上に，原点Oでそれぞれx軸，y軸に接する2円があり，この2円は点P(1,p) (p>0) で交わっている。
この2円の中心間の距離の最小値を求めよ。

>>32
(1)が5になった。自信なし。
(2)(1)のAF(X) = Qを満たすXをX0とすると
EG(X0)⊂X0⊂AF(X0) = Q
だから
X0∈{ AF(EG(X0)) = Q}
により
min {|X| ; AF(EG(X0)) = Q} ≦ 5。

>>36
x y の法6の類は交代していくだけなので
sin((xn -x0)π/3) + sin((yn -y0)π/3) = 0

>>38

2円の中心は、線分OPの垂直二等分線とx軸，y軸の交点。
（(pp+1)/2，0)　(0，(pp+1)/(2p))
その距離の2乗は
　L(p)^2 = (pp+1)^3 /(2p)^2 = (27/16) + (1/4)(pp+4)(pp-1/2)^2 ≧ 27/16,
　L(p) ≧ L(1/√2) = (3√3)/4,

二円の直径(半径)は、それぞれ直角三角形の相似で簡単にわかる.
M = 4L^2 = (pp + 1)^2 + ( p + 1/p )^2 = q^2 + 3q + 1/q + 3 (q = pp と置いた)
M'= 2q + 3 - 1/(qq) = (2q^3 + 3qq - 1) = (2q - 1)(q+1)^2 /(qq)
増減表より M は q=1/2 にてミニマム値をとる事がわかる. (条件 q＞０)
よって L_min = (1/2) * √(1/4 + 3/2 + 2 + 3)　= 1/4 √27 = (3√3)/4.

前スレ>>897
サイコロを繰り返し投げ, 出た目が直前の回に出た目の約数でなくなったら終了する。
n回目にサイコロを投げ, かつその目が1である確率p[n]を求め, n回目に終了する確率をp[n]とp[n+1]を用いて表せ。

n回目が1になるのは, 次のような経過の場合である:
6→1, 6→3→1, 6→2→1, 5→1, 4→1, 4→2→1, 3→1, 2→1, 1

∴n回目が1である確率P[n]は,
P[n]={1+5･C(n-1, 1)･3･C(n-1, 2)}/6^n
=(3n²+n-2)/(2･6^n)
を得て,
n-1回目に終了していない確率は, 6･P[n]なので
, n回で終了する確率は,
6(P[n]-P[n+1])=(15n²-n-14)/(2･6^n)
を得る。

n回目が1である確率から, 直ちにn回で終了する確率が求められるところが面白いと感じますね。

ある牧場では１００頭の羊を放すと１５日間で牧草がなくなり、
１２０頭の羊を放すと１０日間で牧草が食べつくされました　
この牧場で８０頭の羊を１０日間放した後、
さらに何頭xかの羊を加えたところ、
加えてから４日間で牧草は食べつくされました
後から加えた羊は何頭ですか
ただし、牧草は１日に一定量a生え、また、
どの羊も１日で同じ量uの牧草を食べるものとします

（ヒント：最初からある草の量をbとおく）

ニュートンのパチモンか

１００

問題文見間違えた >46はなしで

80

>>34

〔前スレ．993〕
　aを正の定数とする。
xyz空間において，円柱 yy + zz ≦ aa と角柱 |x| + |z|≦ a との共通部分をKとする。
(1) Kの体積を求めよ。
(2) Kの表面積を求めよ。

>>49
(1)
ｚ＝一定の平面で切ると、
　|x| ≦ a - |z|,
　|y| ≦ √(aa-zz),
の長方形。

V = 8∫[0，a] (a-z)√(aa-zz) dz = (2π - 8/3)a^3 = 3.61651864a^3

a,b,cは素数で、2≦a≦b≦cかつa+b>cを満たす。
AB=c,BC=a,CA=bである△ABCの面積をS(a,b,c)とする。

（1）有理数pと自然数nを用い、S(a,b,c)=p√nと表したとき、n=1とならないことを示せ。

（2）次の命題の真偽を述べよ。
「どのような素数qについても、a,b,cをうまく選ぶことで、n=qとなるようにできる」

高専２年
行列
(2)は簡単ですが(1)の固有値が求まりません
お願いします

https://i.imgur.com/3dtlh8V.jpg

わからないんですね

>>50 (補足)

∫(a-z)√(aa-zz) dz
　= ∫[(1/2)a^3 -aaz -azz +z^3]/√(aa-zz) dz + (a^3)/2･∫1/√(aa-zz) dz
　= (1/6) (2aa+3az-2zz) √(aa-zz) + (a^3)/2･arcsin(z/a) +c,

「無」に勝るものは何もありませんか？

>>44
増える草の量+最初の草の量-食べる草の量=0
として式を作る。
15a+b-100*15u=0
10a+b-120*10u=0
14a+b-(80*14+x*4)u=0
これを解くとx=80

>>38
勝手にtで置いてたけどpだったかw

面白スレの795で、宝は2つのまま、縦と横のマス数をそれぞれn、n+1と置いたとき、横に沿って探した方が相手より先に見つけやすいことは3,4の場合でそうだったことから容易に想像出来るが、その証明は出来るだろうか？

縦nマス、横n+1マスのn(n+1)マスのうちランダムに選ばれた2マスにそれぞれ宝が眠っている。
縦1列を探し終えたらすぐ右の1列に移って宝を探していく方法をとるP君と、横1行を探し終えたらすぐ下の1行に移って宝を探していく方法をとるQ君が、同時に左上の地点から探索を開始した。
例えば、n=3の時はP君はAEIBFJCGKDHLの順で探す。Q君はABCDEFGHIJKの順で探すことになる。

ABCD
EFGH
I JK L

1つの地点を捜索するのにかかる時間は同じで、相手が1度探し終えた地点を重複して調べることも当然ある。
相手より先に宝を見つけた方を勝者とする。同時の場合は引き分けとする。
どちらの方が有利になるだろうか？

え？3x4なら横からやったほうがいいの？
直観的には同じだけど…

50の(2)ってどうすればいいの？

>>62
49のでした

>>62
49のでした

>>64
切断面は 半分の楕円が4つなので簡単、残りの円柱側面は積分で求める。
S = 2 * (π a (√(2) a) ) + 4 a² ∫ [0, +π] dθ (1- sinθ)
以下略

α，β，γ は α>0，β>0，γ>0，α+β+γ=π を満たすものとする．このとき， sinαsinβsinγ の最大値を求めよ．

最もエレガントな回答を教えてください。
ごちゃごちゃ一つ固定して微分すればすぐ解けますが
対称性から一発で解けたりしませんか？

>>66
面積に直したら、３項の相加相乗の問題に帰着するから一瞬じゃないの？

あんまり一瞬でもないな
適当すぎｗ

>>66
よく知られてるのは log sin x の凸性使うやつだな。

無に勝てるものはありますか？

z/{((z-1)^2)((z-2)^3)}

の各特異点における留数を求めるのって
z=1 だったら
(z-1)^5をかけて4回も微分して極限をとるっていうことしないといけないのってめちゃくちゃ手間がかかると思うんですけど
そうする以外に簡単にもとまる方法ってないですか？

>>66
f(α,β,γ) = sinαsinβsinγ と置く.
領域境界では f = 0 、領域内点では f ＞ 0 .
境界が素直なので f の勾配ベクトルが平面 α + β + γ = π と直交する点を探せばよい.
つまり cosα sinβ sinγ = sinα cosβ sinγ = sinα sinβ cosγ より
tanα = tanβ = tanγ ∴ α = β = γ = π/3
f = (√(3)/2)^3 = (3/8)√3 を得る.

>>71
z/{((z-1)^2)((z-2)^3)}
= {(z-1) + 1}/{((z-1)^2)((z-1 - 1)^3)} (以降 h = z-1 と置く)
= -(1/h + 1/h^2) * (1 + h + h^2 + ...)^3
= -(1/h + 1/h^2) * (1 + 3h + ...)
= -1/h^2 - 4/h - ...
1/h の係数だけ拾えばよい

(z-2 + 2)/{((z-2 + 1)^2)((z-2)^3)}
= (1/h^2 + 2/h^3) { 1 - h + h^2 + ... }^2
= (1/h^2 + 2/h^3) { 1 - 2h + 3h^2 +... }
以下略

>>73
おおおおおおおお
確かに！！！！！！！！
ありがとうございます

>>51
これお願いします
（2）がわかりません

>>69

最もエレガントな解答は log(sin(x)) の凸性使えば一発で出ますが >>69

GM-AM で下準備
　sinα sinβ sinγ ≦ {(sinα + sinβ + sinγ)/3}^3
してから sin(x) の凸性使う
　(sinα + sinβ + sinγ)/3 ≦ sin{(α+β+γ)/3} = sin(60ﾟ) = (√3)/2,
ほうが簡単かもです。

こういうのをゴリ押しで解こうとするたび思うんだが、sinxをexp(ix)で表しても手間は減らないもの？

>>65
ありがとうございます

>>60
コンピュータでシミュレーションしてみた。

n=3のときは (P1st:：P君が先に見つける宝の埋没場所の組み合わせ数）

> t342=treasure(3,4,2)
P1st Q1st even
26 27 13

n=4のときは

> t452=treasure(4,5,2)
P1st Q1st even
84 83 23

常に横に探す方が有利ではないようだ。

Rでのコードはここ

http://tpcg.io/d6OYvn

>>79
ｎを変化させてP,Qが先に見つける宝の配置を計算させてみた。

大きいほうが有利になる。

> t(sapply(1:15,treasure1))
P1st Q1st even
[1,] 0 0 1
[2,] 4 5 6
[3,] 26 27 13
[4,] 84 83 23
[5,] 203 197 35
[6,] 413 398 50
[7,] 751 722 67
[8,] 1259 1210 87
[9,] 1986 1910 109
[10,] 2986 2875 134
[11,] 4320 4165 161
[12,] 6054 5845 191
[13,] 8261 7987 223
[14,] 11019 10668 258
[15,] 14413 13972 295

2点(0,0,0),(2,0,1)を通る直線をl,2点(1,-2,0),(0,-4,-1)を通る直線をmとし、l,mをz軸のまわりに、1回転して得られる曲面をそれぞれα、βとする。

>>81
2平面z=0,z=5とαで囲まれた部分をA,2平面z=0,z=5とβで囲まれた部分をBとするとき、共通部分A∩Bの体積を求めよ

>>82
詳しい解説お願いします。

>>82
>>83
自分の答えは2511π/15となったんですがあっていますか？

約分

>>85約分すればあっていますか？

>>60>>61
Ωの部分集合を事象と言う
Ω自身は全事象と言う

Ω＝｛A，B，C，D，E，F，G，H，I，J，K，L｝となる

各 i （１≦i≦１２） が根元事象である

最初に宝が出るという事象A＝｛宝｝で確率P(A)は

P(A)＝１／１２　となる

最初に探す方向を i
列が変わる時を ｊ として

最初に宝が出るという事象Aと事象Bを考える.

A＝｛（i，j）|　i または j が宝｝
B＝｛（i，j）|　i または j が宝｝

Ω＝｛（i，j）|１≦i≦ｎ，１≦j≦ｎ＋１｝となり

このｎ（ｎ＋１）通りの各要素が根元事象

縦方向に探査する場合

Ω＝｛（i，j）|１≦i≦ｎ，１≦j≦ｎ＋１｝から

#A＝ｎ（ｎ＋１）−ｎ（ｎ−１）＝２ｎ

#Aは事象Aに含まれる要素の個数

横方向に探査する場合

Ω＝｛（i，j）|１≦i≦ｎ＋１，１≦j≦ｎ｝から

#B＝ｎ（ｎ＋１）−ｎ（ｎ−１）＝２ｎ

最初に宝が出る確率は

∴P(A)＝P(B)＝２ｎ／ｎ（ｎ＋１）

>>84
計算ミスしてました
156πです

σをn次の置換とする。R^nからR^nへの写像で、(x_1,...,x_n)を(x_σ(1),...,x_σ(n))にうつすものは連続であることを示して下さい。

直観的に考えたら違う理由が思いつかないから書いたんだけど…
何故違うかもしれないと考えたのかわからないレベルで違う理由が思いつかない

ABCDEFGHIJK
AEIBFJCGKDHL
と並んでる状態で、A-Kのうち２個がランダムで当たり
最初の当たりが左に近いのはどっち？ってことじゃん

>>80では有意差が有るように見えるけど、何故なのかよくわからない

>>87
読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。

>>90
別スレの解説をコピペ

なるほどねえ
確かにQの方が微妙に先に見つける場合が多いな

Pが先に見つけるのは以下の26通り
CE,DE,DI,EF,EG,EH,EI,EJ,EK,EL,FG,FH,FI,FJ,FK,FL,GI,GJ,HI,HJ,IJ,IK,IL,JK,JL,KL

Qが先に見つけるのは以下の27通り
BC,BD,BF,BG,BH,BI,BJ,BK,BL,CD,CF,CG,CH,CJ,CK,CL,DF,DG,DH,DJ,DK,DL,GH,GK,GL,HK,HL

同時に見つけるのは以下の13通り
AB,AC,AD,AE,AF,AG,AH,AI,AJ,AK,AL,BE,CI

>>91
具体的な反例を伴わないのは詭弁ですよ

>>93
既に>80で実証済

>>80

n=2

ABC
DEF

の場合

短軸方向探索Pが先に宝を発見する埋め方：４通り
> print(matrix(LETTERS[t232$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4]
[1,] C D D E
[2,] D E F F

長軸方向探索Qが先に宝を発見する埋め方：５通り
> print(matrix(LETTERS[t232$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5]
[1,] B B B C C
[2,] C E F E F

同時に宝を発見する埋め方：６通り
> print(matrix(LETTERS[t232$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] A A A A A B
[2,] B C D E F D

なんか納得できない結果が出てきてて頭がぐるぐるううううう

そんなの当たり前じゃん(´・ω・`)

等確率にしかならないのに無理やり差異を
見つけようとしているもん

>>96
>95の操作をn=20までやってみた。

> t(sapply(1:20,treasure1))
P1st Q1st even
[1,] 0 0 1
[2,] 4 5 6
[3,] 26 27 13
[4,] 84 83 23
[5,] 203 197 35
[6,] 413 398 50
[7,] 751 722 67
[8,] 1259 1210 87
[9,] 1986 1910 109
[10,] 2986 2875 134
[11,] 4320 4165 161
[12,] 6054 5845 191
[13,] 8261 7987 223
[14,] 11019 10668 258
[15,] 14413 13972 295
[16,] 18533 17988 335
[17,] 23476 22812 377
[18,] 29344 28545 422
[19,] 36246 35295 469
[20,] 44296 43175 519

シミュレーションしても>92の結果に合致。

> x=c(1,1,rep(0,10))
> PQ <- function(){
+ Q=sample(x)
+ z=matrix(Q,ncol=4,byrow=T)
+ P=as.vector(z)
+ c( even=which.max(P) == which.max(Q),
+ p1st=which.max(P) < which.max(Q),
+ q1st=which.max(P) > which.max(Q))
+
+ }
> k=1e6
> re=replicate(k,PQ())
> mean(re['even',]) ; 13/(26+27+13)
[1] 0.197025
[1] 0.1969697
> mean(re['p1st',]) ; 26/(26+27+13)
[1] 0.393803
[1] 0.3939394
> mean(re['q1st',]) ; 27/(26+27+13)
[1] 0.409172
[1] 0.4090909

>>96
宝の埋め方の組み合わせを列挙して分類したら
>95のようになるのは同意？

>>92
この結果面白いね
問題が２つ見つけるまでやって多く獲った方どっち？だったらイーブンだけど、１つ目を先に獲った方が勝ち、とすると差が出る
この場合、２番目を先に見つける確率にもきっと差があるのだろう

>>100
納得できないのは直観的に納得できないだけで、そういうことになるよなぁとはわかっていると思います

今ちょっと考えているのが、遅く見つけたほうが勝ちというルールで行うなら
Ｑ：ABCDEFGHIJKL
Ｐ：AEIBFJCGKDHL
では、Ｐ君の方が勝率は高いということ。
じゃあ、Ｑに対してＰ以上に勝率の高い文字列（検索順序）は存在するはずだけど
それらを具体的に求める方法は？
とか考えてしまう。
で、頭がぐーるぐーるるるるるる

>>101
先に２つの宝を発見した方を勝者とするのでやってみた。

n=2

ABC
DEF

の場合


> t232=treasure2(2,3,2)
P1st Q1st even
5 4 6

短軸方向探索Pが先に２つの宝を発見する埋め方：５通り
> print(matrix(LETTERS[t232$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5]
[1,] A A B B D
[2,] D E D E E

長軸方向探索Qが先に２つの宝を発見する埋め方：４通り
> print(matrix(LETTERS[t232$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4]
[1,] A A B C
[2,] B C C D

同時に２つめの宝を発見する埋め方：６通り
> print(matrix(LETTERS[t232$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] A B C C D E
[2,] F F E F F F

先に２つの宝を発見した方を勝者とするのでやってみた。

n=3

ABCD
EFGH
IJKL

の場合

> t342=treasure2(3,4,2)
P1st Q1st even
27 26 13
> #短軸方向探索Pが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B B C C C D E E E E F
[2,] E F I J K E F I J K I J K K F I J K I
[,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27]
[1,] F F G G G I I J
[2,] J K I J K J K K
> #長軸方向探索Qが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B C C C C C D D D D D
[2,] B C D G H C D G H D E F G H E F G H I
[,20] [,21] [,22] [,23] [,24] [,25] [,26]
[1,] E E F F G H H
[2,] G H G H H I J
> #同時に２つめの宝を発見する埋め方
> print(matrix(LETTERS[t342$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,] A B C D D E F G H H I J K
[2,] L L L J L L L L K L L L L

>>104
R　の　コードは　 http://tpcg.io/RZo4hd

ABCD
EFGH

を一般化するとこうなるかな？

横一列に並んだ ABCDEFGH の中からランダムに2つを選んで宝を隠しておく
ABCDEFGH の中から宝を探す順番は8！通りある
最初の宝を見つけた時点で終わるものとするとき、
8！通りの探し方の中で最も有利な探し方はどのような探し方か

もっとも有利なんてないやろ。
じゃんけんと一緒。
どんな列取ってきてもその先頭文字を末尾に回した探索には負ける。

　A..B..C..D
A■■■□
E■■■■
I ■□■■

全部の中で一番がないというのは>>107の考察通りだと思う。
けど、ある特定の列に対して最も勝率が高いのはどれだろうとは気になる。

けど、先頭文字を末尾に回した奴が一番勝率高くなるのかな。
ABCDEFGに対してなら
BCDEFGAが一番勝率高い気がする

最初に当たり一つ引けばそこでゲーム終了だから
二つ目の当たりとの組み合わせは考慮しなくていい

当たりがどの座標のマスに置かれても
要素の個数は変化しないので
どの方向からの探査によっても確率は変化しない

これ、当たりが１個でも、探索順番によって勝率変わってくるな
やっと構造がなんとなくわかってきた
自分の脳みその弱さが悲しくなってくる

>>109
部屋の数についての帰納法でいけるんじゃね？
主張は
部屋の数が n の時 P:A[1]A[2]…A[n] に引き分けないという条件下で勝つ確率最大なのは Q:A[2]A[3]…A[n]A[1]。
以下Qの探索順をB[i]とする。
n=3では多分成立。
n<k で成立として n=k のとき。
P が Q に勝つのはA[1]とA[2]以外に宝が配置されるときでその確率は (n-2)/C[n,2]。
引き分けるのはA[1]、A[2]に配置されるときで確率1/C[n,2]。
よってQがPに勝つ確率は (C[n] - 1 - (n-2)/(C[n,2] - 1)。
容易にA[1]≠B[1]の場合はコレより確率は大きくならないとわかる。
A[1] = B[1]の場合を考えればよい。
このとき引き分けないという条件下では宝箱はA[1]以外の２つに配置される場合でその場合Qの勝つ条件付き確率の最大値は (C[n-1] - 1 - (n-3)/(C[n-1,2] - 1)。
多分 (C[n] - 1 - (n-2)/(C[n,2] - 1) ≧ (C[n-1] - 1 - (n-3)/(C[n-1,2] - 1)より成立。

>>112
おお、そういう風に片付くのか
帰納法で出来ないかとも考えたけど、自分の頭では無理だったのです
これで自分はスッキリしました！

>>112
一列じゃなくて長方形型 n×n+1の配置じゃないの？
nの次は (n+1)(n+2)では

>>98の結果をみると20までだが
縦nマス、横n+1マスのn(n+1)マスのとき
n=1でイーブン
n=2,3で長軸方向探索が有利
n=4以上で短軸方向探索が有利となっているので
数学的帰納法はn=3で適応できないと思う。

>>114 115
考えてる問題が違う。

>>109

# ABCDEFGに対してなら
# BCDEFGAが一番勝率高い気がする
library(gtools)
n=7
k=2
perm=permutations(n,n)
Q=perm[1,]
np=nrow(perm)
p1st=numeric(np)
for(i in 1:np){
P=perm[i,]
tre=combn(n,k)
nt=ncol(tre)
re=numeric()
for(j in 1:nt){
re[j]=min(which(tre[1,j]==P),which(tre[2,j]==P))-
min(which(tre[1,j]==Q),which(tre[2,j]==Q))
}
p1st[i]=sum(re<0)
}
plot(p1st)
p1st[which.max(p1st)]
(p.max=which(p1st==15))
print(matrix(LETTERS[perm[p.max,]],ncol=7),quote=F)

#

>>117
# ABCDEFGに対してなら
# BCDEFGAが一番勝率高い気がする

一番勝率高い探索順は４通りあった

> print(matrix(LETTERS[perm[p.max,]],ncol=7),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] B C D E F A G
[2,] B C D E F G A
[3,] B C D E G A F
[4,] B C D E G F A

>>118
なるほどね
先回り側がEの次の部屋へ進むってことは当たりはFGだからどっちに進んでも同じか

そうして、
先回り側の順序の最後の2つは決して実行されない
その４つの順序のどれでも最後から3番目のFかGまでで決着が付くから

宝を２個先にみつけた方が勝者とすると
ABCDEFGに対して一番勝率高い探索順は？

これも４通り出てきた。

[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] A C D E F G B
[2,] B C D E F G A
[3,] C A D E F G B
[4,] C B D E F G A

要するに｢相手に『自分が探索済のマス』を探させる｣｢自分は『相手が探索済のマス』を探さない｣の2つを出来るだけ守っていればいい話だから相手の探索方法に対応する最適解の議論はあまり意義がないのではと思う

しかし、n=4から先はずっと短軸探索が有利になるのか。長軸側が逆転することはなさそうだし、n≧4の場合について｢短軸探索が有利である｣は成り立ちそう。これを証明することは出来ないだろうか…

>>124
俺もそっちに興味があるが、証明できる頭脳はない。

宝箱１個なら、なんとか証明できそうな感じだし、そこから拡張すれば宝箱２個でもいけるのかなぁ
整数苦手だからよくわかんない

>>126
そうは問屋が卸さないみたいだよ。

縦４マス、横５マスで宝箱を１から７まで増やしてみると

宝が６個になると短軸有意から長軸有意に逆転した。

処理速度の制約であまり大きな数字で検証できないのだが。

> sapply(1:7,function(k) treasure(4,5,k))
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
短軸有利 9 84 463 1776 5076 11249 19797
長軸有利 9 83 453 1753 5075 11353 20057
同等 2 23 224 1316 5353 16158 37666

>>127
気長にやってみた。

> sapply(1:20,function(k) treasure0(4,5,k))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11]
短軸有利 9 84 463 1776 5076 11249 19797 28057 32243 30095 22749
長軸有利 9 83 453 1753 5075 11353 20057 28400 32528 30250 22803
同等 2 23 224 1316 5353 16158 37666 69513 103189 124411 122408
[,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
短軸有利 13820 6656 2486 695 137 17 1 0 0
長軸有利 13831 6657 2486 695 137 17 1 0 0
同等 98319 64207 33788 14114 4571 1106 188 20 1

４×５の場合
宝：１個　同等
宝：２〜５個　短軸有利
宝：６〜１３個　長軸有利
宝：１４〜２０個　同等

と推移した。

>>127
えー
なかなか簡単にはいかせてもらえないね
考えながら仕事片付けるとしよう

>>129
それに宝箱１のときは、イーブンだろ。

あれーほんとだ
どこかで派手な勘違いをしたままでやってたぽい
そりゃあっちこっちぐるぐるうううになるわ…

あかん、今考える余裕ないｗ
>>131の同じになるってのもとりあえず保留

>>124
宝２個でn=30まで計算させてみた。４以上で短軸有利は不変だった。

> t(sapply(1:30,treasure1))
P1st Q1st even
[1,] 0 0 1
[2,] 4 5 6
[3,] 26 27 13
[4,] 84 83 23
[5,] 203 197 35
[6,] 413 398 50
[7,] 751 722 67
[8,] 1259 1210 87
[9,] 1986 1910 109
[10,] 2986 2875 134
[11,] 4320 4165 161
[12,] 6054 5845 191
[13,] 8261 7987 223
[14,] 11019 10668 258
[15,] 14413 13972 295
[16,] 18533 17988 335
[17,] 23476 22812 377
[18,] 29344 28545 422
[19,] 36246 35295 469
[20,] 44296 43175 519
[21,] 53615 52305 571
[22,] 64329 62810 626
[23,] 76571 74822 683
[24,] 90479 88478 743
[25,] 106198 103922 805
[26,] 123878 121303 870
[27,] 143676 140777 937
[28,] 165754 162505 1007
[29,] 190281 186655 1079
[30,] 217431 213400 1154

正の整数の組(x,y)であって,x!+y!=x^yを満たすようなものを全て求めよ
の解説をして頂けませんか？

答えは2,2 2,3だと思うのですが解答が無くて
よろしくお願いします

>>89
お願いします

>>134
xは偶数しかありえないのでx=2mとおけば

>>136
y≦x-1のとき,
x!+y!=y!(x!/y!)+y!=y!((x!/y!)+1),
(3≦)(x!/y!)+1=x･(x-1)!/y!+1とxは互いに素だから, x!+y!≠x^y.
すなわちx≦y.
3≦xのとき,
x!+y!=x!(1+(y!/x!))は(x-1)(≧2)の倍数.
x-1とxは互いに素であり, x!+y!≠x^y.
すなわちx≦2.
1)x=1のとき, 与式を満足させるyはない.
2)x=2のとき, 2+y!=2^y.
y≧4とすれば,
2+y!=2+24･(y!/4!)>2+3･2^(k-1)>2^k.
すなわちy≦3.
よって求める組は(x,y)=(2,2), (2,3).

できました！

宝箱問題、
もとの 4x3 型の12部屋で宝箱の数を変えてみると
1と8以上で有利不利無し、それ以外は長軸優先有利となるな
初見での印象よりも随分奥深いなこれ

>>138
計算式お願いする

数列の項を並べ替えてできる数列の収束性、極限値は如何？

もう少し正確にいうと、
全単射関数　n : N -> N で 数列 a[ i ] を n で並べ替えた数列b[ i ]を
b[ i ] = a[ n(i) ] で定義する。
b[ i ] の収束性、極限値はどうなるでしょう?

プログラムで計算したので式はなんとも

部屋が
ABCD
EFGH
IJKL

として
宝物10個のときはABが空きなら縦の勝ち、
AEが空きなら横の勝ち

縦勝ちの宝物9個の配置
CDEFGHIJK
CDEFGHIJL
CDEFGHIKL
CDEFGHJKL
CDEFGIJKL
CDEFHIJKL
CDEGHIJKL
CEFGHIJKL
DEFGHIJKL

横勝ち
BCDFGHIJK
BCDFGHIJL
BCDFGHIKL
BCDFGHJKL
BCDFGIJKL
BCDFHIJKL
BCDGHIJKL
BCFGHIJKL
BDFGHIJKL

以下各個数での勝敗の数
treasures 1: p win 5 q win 5 even 2
treasures 2: p win 26 q win 27 even 13
treasures 3: p win 73 q win 76 even 71
treasures 4: p win 133 q win 140 even 222
treasures 5: p win 167 q win 176 even 449
treasures 6: p win 148 q win 153 even 623
treasures 7: p win 91 q win 92 even 609
treasures 8: p win 37 q win 37 even 421
treasures 9: p win 9 q win 9 even 202
treasures 10: p win 1 q win 1 even 64
treasures 11: p win 0 q win 0 even 12
treasures 12: p win 0 q win 0 even 1

>>139
Rでよければこんな感じ

# 宝の数を変化させる
treasure0 <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
return(c(短軸有利=sum(re<0),長軸有利=sum(re>0),同等=sum(re==0)))
}
sapply(1:12,function(k) treasure0(3,4,k))

> sapply(1:12,function(k) treasure0(3,4,k))
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
短軸有利 5 26 73 133 167 148 91 37 9 1 0 0
長軸有利 5 27 76 140 176 153 92 37 9 1 0 0
同等 2 13 71 222 449 623 609 421 202 64 12 1

>>138

>128に書いたけど

４x５だと宝箱を増やすと途中で短軸有利から長軸有利に変わっちゃうので自分でもびっくりした。

> 141
よろしければプログラムコードをアップしていただけませんか？Pythonでしょうか？

>>143
ちょっと整理してました。

NB. n comb n returns all n length set from 0..m-1
comb =: dyad define
if. x=1 do.
(1,~y)$i.y
elseif. x=y do.
(1,y)$i.y
elseif. do.
((y-1) ,/"0 1 (x-1) comb y-1 ), x comb y-1
end.
)

NB. usage: 3 4 game 2
game =: dyad define
p =. ,/ |: x $ i. */x
q =. i. */x
g =. y comb */x
d =. (<./"1)@(g &((i."1 0)~))
r =.(d p)-(d q)
y, (+/ r<0), (+/ r>0), (+/ r=0)
)

NB. run 3 4 games n for n in 1..12
smoutput 'tre p q even'
smoutput 3 4 game "1 0 (1+i. 12)

>>144のコードはマイナー言語J
ここで実際に動かしてみることができます

https://goo.gl/znRTwf

一般に,
U[j=1, n]A_j=ℝ となるn個の集合 A_j (*1) について,
j=1,2,...,n で a_j∈A_j となるような変数 a_j を取り,
lim[a_j→α] f(a_j) =k (*2) が全ての j について言えたならば,
lim[x→α] f(x) =k (*2) が言えますか。
例えば, p∈ℚ, q∈ℝ\ℚ とすると, p と q を合わせれば全実数を取ります。このとき,
lim[p→α] f(p) =lim[q→α] f(q) =k
かつ
lim[x→α] f(x) ≠k
となる f(x) は存在しますか。

(*1)αに十分近い要素も含む
(*2)離散的極限

いいえ

>>145
お手数かけました。
残念ながら自分の知識ではアラビア文字のように理解不能でした。

>>60
一つ質問ですが
スタート地点Aに宝があるとゲームスタートと同時に
同着でゲーム終了になるけど、ポイントAに宝は設置されるのですか？

σをn次の置換とする。
R^nからR^nへの写像で、(x_1,...,x_n)を(x_σ(1),...,x_σ(n))にうつすものは連続であることを示して下さい。

>>147
示してください

>>140
条件収束する級数を考えればa[i]とb[i]の収束性に関係がないことは明らか

>>149
その場合は引き分けで終了。
宝の置き方はランダム。
12C2=66通りに等確率で配置。

>>145
数字を増やしたらサイトの時間制限を超えて結果がでなくて残念。
尚、>142のRはメモリ不足で停止しました。

NB. usage: 5 6 game 2

NB. run 5 6 games n for n in 1..30
smoutput 'tre p q even'
smoutput 5 6 game "1 0 (1+i. 30)

SMアウトプットとか、なんかヤラシイな、おい。

>>155
分散分析でF分布の値の比に　F-ratio というのが出てくるの知ってた？

可換環論ではAss、穴(Ann)、ホモロジー、(チェイン)ホモトピー、……汚い言葉がいっぱい出てくるよ！やったね！

>>140
a[i] → c とする。
e>0 とする。
|a[i] - c| ≧ e である i は有限個。
∴ |b[i] - c| ≧ e である i は有限個。
∴ b[i] → c。

>>146
離散的極限って離散位相での極限？
だったら Aj が disjoint な集合なら
a[1]→α、a[1]∈A[1]、a[2]→α、a[2]∈A[2] 自体が起こりえないやろ？
誘導位相？

>>154
12部屋から6部屋選ぶ組み合わせは924通りしかないのに
20部屋から10部屋だと184756通り、
30部屋から15部屋だと155117520通り、
という感じなのでどうしても時間やメモリを食いますよね

>>60
場合分けなどが面倒くさくて疲れ果てたけど、計算結果は>>133と一致。
P1st(n)-Q1st(n) が（偶奇によらず） (n^2-2n-6)(n-1)/6 になったので、n=2,3でＱが、n≧4でＰが有利。

コードはSagemath。
from sage.calculus.calculus import symbolic_sum
,var m,l,k,a,n
P1 = (symbolic_sum((m-1)*(m)-2*l-1, l,1,m-2)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,2*k) + symbolic_sum((m-1-k)*(m-k)+3*k-2*l-1, l,2*k+1,m-2), k,1,a-2)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,m-2), k,a-1,m-2)
).substitute({a:m/2}).substitute({m:n+1})
P2 = (symbolic_sum((m-1)*(m)-2*l-1, l,1,m-2)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,2*k) + symbolic_sum((m-1-k)*(m-k)+3*k-2*l-1, l,2*k+1,m-2), k,1,a-1)
+ symbolic_sum(symbolic_sum((m-1-k)*(m-k)+k-1-l, l,k,m-2), k,a,m-2)
).substitute({a:(m-1)/2}).substitute({m:n+1})
def P1st(x):
return P1.substitute({n:x}) if mod(x,2) == 1 else P2.substitute({n:x})
Q1 = (symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,l-1), l,0,a)
+ symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,m-l-1) + symbolic_sum((m-k)^2+k-l,k,m-l,l-1), l,a+1,m-1)
+ symbolic_sum(k^2-2*m*k+m^2+k-m,k,0,m-2)
).substitute({a:(m-1)/2}).substitute({m:n})
Q2 = (symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,l-1), l,0,a-1)
+ symbolic_sum(symbolic_sum((m+1-k)*(m-k)+k-2*l-1, k,0,m-l-1) + symbolic_sum((m-k)^2+k-l,k,m-l,l-1), l,a,m-1)
+ symbolic_sum(k^2-2*m*k+m^2+k-m,k,0,m-2)
).substitute({a:m/2}).substitute({m:n})
def Q1st(x):
return Q1.substitute({n:x}) if mod(x,2) == 1 else Q2.substitute({n:x})

P1 == 1/24*(6*n^3 + 20*n^2 - n - 27)*(n - 1) # nが奇数のとき
P2 == 1/4*n^4 + 7/12*n^3 - 7/8*n^2 - 13/12*n + 1 # nが偶数のとき
Q1 == 1/24*(6*n^2 + 10*n - 3)*(n + 1)*(n - 1) # nが奇数のとき
Q2 == 1/24*(6*n^2 - 2*n - 5)*(n + 2)*n # nが偶数のとき

>>161

>133です。労作ありがとうございます。

コードは全く読めないのですが、宝の数を増やしての計算はこのコードで可能なのでしょうか？

４×５の場合で宝を増やすと

宝：１個　同等
宝：２〜５個　短軸有利
宝：６〜１３個　長軸有利
宝：１４〜２０個　同等

に変化したので差分はどんな関数なのだろうかとか、
５×６ではどうなるのか（メモリ不足で実行できませんでした）とか興味があります。

>>161
(n^2-2n-6)(n-1)/6 をグラフ表示してみました。

http://i.imgur.com/Qel6ZAy.png

>>162
>>161は多項式にまでするために、部屋をn x (n+1)、宝を２個と特殊化したものです。

#nloc(m,n,k,l)は縦m、横nの部屋で横優先が部屋(k,l)で初めて宝を発見する場合で
#宝が置かれても縦優先に先を越されない部屋の数。
def nloc(m,n,k,l):
q,r = divmod(n*k+l,m)
return (n-q)*(m-k)+q-1-l + ((k-r) if r > k else 0)

#nwin(m,n,c)は部屋が縦m、横n、宝がc個で横優先が勝つ宝の配置の数
def nwin(m,n,c):
return sum(binomial(nloc(m,n,k,l),c-1) for k in range(m) for l in range(n) if k*(n-1)<l*(m-1))

縦優先は縦横を替える。

>>164
レスありがとうございます。
コードは読めないのですが、
部屋数から宝部屋の組合せを列挙してどちらが縦横どちらが先にみつけるかを探る手続きで必要な計算式をプログラムが絞りだしてくれるという理解でいいのでしょうか？

>>165
いえ、計算式そのものです。数式で書けば
nwin(m,n,c) := Σ[(k,l)∈{0,…,m-1}×{0,…,n-1}, k*(n-1)<l*(m-1)] binomial((n-q)*(m-k)+q-1-l + (k-r)δ(r > k), c-1)、
ただし、n*k+lをmで割った商をq、余りをrとし、δ(P)をPが真なら1、偽なら0である関数とする。

>>159
例えば→√2を考えたい時、qの近づけ方は問題ないんでしょうが、pの近付き方を、p_n→√2になるような有理数列p_n上で考えることは出来ないんでしょうか。
→0なんかも、実際に関数に0を入れるわけではなくギリギリまで近付けるように、p自身が√2を取れないのは、定義できないほどの大問題でしょうか。

>>168
Ωの部分集合を事象と言う
Ω自身は全事象と言う

最初に探す方向を i
行または列が変わる時を ｊ として
P君とQ君のうちどちらが先に宝を見つけるのかという
事象Aと事象Bを考える.

A＝｛（i，j）|　i または j が宝｝
B＝｛（i，j）|　i または j が宝｝

>>168

別スレでは等確率とデタラメ書いてたよなぁ。

https://rio2016.5ch.net/test/read.cgi/math/1540218853/87


読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。

>>168
すでに正解とPCでのカウントの照合が終わっているのに

読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。

別スレだって( ´,_ゝ｀)

>>166
とすると、部屋数が増えたり宝が増えても数式として算出可能なのでしょうか？

>>172
なんだ、このスレでデタラメ書きつづけてたのかよｗ

>>152 , >>158
ご返答ありがとうございます。

>>167
もちろんそういう近づけかたを考えてもいいけど、その近づけかたを離散位相といってはいけない。
流石にこの程度の基本的な単語は正確に意味を確認するようにしないといかん。

>>175
>>152は無視してな
完全に寝ぼけてたわ

>>173
すみません。発言がよくわかりません。
「数式として算出」とは？ >>166は数式ではない？
ここでいう数式とは多項式などのΣのない形のものでしょうか？
「部屋数が増えたり」も、もともと部屋をn×(n+1)などとしていて大きさを変えられますよ？

宝の数が2以外でも（3なら3と）固定されていて部屋の形がn×(n+1)または(n+1)×nなら>>161を少し変えれば
Σのないnについての多項式が得られる、とは言えます。

念のため書いておくと
式>>166は部屋の縦横、宝の数が任意だが、Σがある。
部屋の形がn×(n+1)または(n+1)×nとして適当に場合分けすることにより数式処理ソフトで
Σの計算できるようにした、そのコードが>>161。

宝箱問題。申し訳ないのですが、プログラミングに詳しくないものでさっぱりです…
高校生にもわかるようにどなたか解説していただけませんか？(入試数学の解答のような形式であればありがたいです)

()>>173
可能だけど立式するのは結構な手間
Σを用いた式として立式して sagemath というソフトで簡略化して n の多項式にしたのが>>161
どう考えても面倒なので161さん以外の誰もやっていなかった

宝箱の数をkとして立式することは可能だろうけれども、
更なる面倒さに付き合ってくれる人がいなければここには書かれない
こんなとこでどうでしょうか。

sagemath はスマートフォンなどでも使うことができて、
僕も今初めて使うので適当ですが例えば>>161の最初の
P1の式(Σを含むもの)の簡略化などはiPhoneアプリでも以下のようにして行えました
アプリ起動して「+」ボタンで新しい式を入力するモードにして

var l,a,n,k = var('l','a','n','k')
a=m/2
m=n+1
sum((m-1)*(m)-2*l-1, l,1,m-2)
+ sum(sum((m-1-k)*(m-k)+k-1-l, l,k,2*k) + sum((m-1-k)*(m-k)+3*k-2*l-1, l,2*k+1,m-2), k,1,a-2)
+ sum(sum((m-1-k)*(m-k)+k-1-l, l,k,m-2), k,a-1,m-2)

と入力して「evaluate」ボタンでこの式を評価(簡略化)

https://i.imgur.com/4bpYZLg.jpg

弥勒（僧）とシュリニヴァーサ・ラマヌジャンはどっちの方が賢いですか？

c　のプログラムが無かったようなので　宝二つ、m×(m+1)型（m=1〜69）を作ったので参考にあげておきます。

http://codepad.org/pbmWeZZ5

少し説明を加えておくと、マスに1から順に番号をあたえます。
配列P[c]には、c番目 のマスをPは何番目に調査するか
配列Q[c]には、c番目 のマスをQは何番目に調査するか　を入れておきます。

i番目とj番目のマスに宝があるとき、P[i]とP[j]を比べて小さい方の値で、Pは宝を発見し、
Q[i]とQ[j]を比べて小さい方の値で、Qは宝を発見します。
この値を比べ、PとQどちらが早く発見したかを判定すると言うだけのものです。

縦、横のマスの数の変更や、宝の数の変更も難しくないと思うので、興味がある方はどうぞ。
（本当は、配列は一つで十分なんだけど、可読性や対称性を考えて書いておきました）

xyz空間の球B:x^2+y^2+z^2=1の表面または内部に点Pをとる。
Pを通り方向ベクトル(1,2,0)に平行な直線lとBとの共有点を考えるとき、以下の問いに答えよ。

（1）lとBの共有点の個数を場合を分けて答えよ。

（2）共有点の個数が2個のときを考える。共有点の一方をS、他方をTとする。
P(x,y,z)とするとき、長さの積PS・PTをx,y,zで表せ。Pが表面上にあるときはP=Sとして考えよ。

（3)Bを平面x=u(-1≦u≦1)で切った切断面D_u上を点Pが動く。P(u,y,z)においてy^2+z^2の取りうる最大値Mをuで表せ。
さらにz=0のとき、積分 I_u = ∫[0→M] (PS・PT) dy をuで表せ。

（4）（3）で求めたI_uに対して定積分K = ∫[-1→1] I_u du を求め、さらに比の値K/(4π/3)を求めよ。

n≧3とする。
次の和を求めよ。
Σ[k=1,2,...,n-1] {(n,k)・(n+1,k-1)}

サイコロを振り、出た目に応じて点Pを動かす。最初点Pは(0,0)にある。
点Pが(a,b)にあるとき、偶数の目が出たら(a+1,b+1)に移動させ、奇数の目が出たら(a+1,b-1)に移動させる。
このとき、以下の事象が起こる確率を求めよ。

（1）Pが半直線y=x(x≧1)の上に乗る。
（2）Pが直線y=2x+1の上に乗る。
（3）m,nを整数の定数とし、Pが半直線y=mx+n(x≧1)の上に乗る。必要があればm,nの値に応じて場合分けして答えよ。

>>178
レスありがとうございました。
多項式で与えられたので他のソフトでも>163のように
簡単にグラフ化できました。
そういう意味で数式と書いたつもりでした。

>>183
いつもｃのコードありがとうございます。

このコードだと縦横マスを増やすのは容易でも、宝の数を増やすには for loopを

for(i=1,Pwin=Qwin=Draw=0;i<mn;i++)for(j=i+1;j<mn;j++) for(k=j+1;k<mn;k++) for(l=k+1;l<=mn;l++)

という具合に増やす必要がありますよね？

そんな感じですね。細かいところですが、少し修正を施すと、
for(i=1,Pwin=Qwin=Draw=0;i<mn-2;i++)for(j=i+1;j<mn-1;j++) for(k=j+1;k<mn;k++) for(l=k+1;l<=mn;l++)
で、空回りを回避してます。

もし、このアルゴリズムで、宝の数を一般数化するなら、i,j,k,...の変数を配列にしてループにいれるか、
再帰関数化するか等の方がスマートですが、二つで固定なら、提示したような感じがシンプルですね。

しかし、宝の数可変を前提にプログラムを組むなら、別の方策を取ります。
Qは時刻 c に最初の宝を見つけるので、
・Pの宝の発見時刻が全てcより大きい　→　Qの勝ち
・Pの宝の発見時刻にcを含み、残りは全てcより大きい　→　引き分け
・それ以外　→　Pの勝ち
です。
Qは、時刻cに、マスcを調査するので、マスc+1、c+2、...の中に、P[x]>c　を満たす
マスがいくつあるかをあらかじめカウントし、テーブル化すれば、あとは、
二項係数の積の和だけの、プログラムとなると思います。

>>185

m+1≧n≧1 のとき
Σ[k=1，n] C(n，k) C(m，k-1) = C(m+n，n-1)

∵ (1+x)^n (1+x)^m を展開したときの x^(n-1) の係数だから。

Σ[k=1,n-1] C(n,k) C(n+1,k-1) = C(2n-1,n-1) - C(n+1,2)

>>182

蝉「おまえさ、人としじみのどっちが偉いか知ってるか？」

伊坂幸太郎「グラスホッパー」角川文庫 (2007)

>>189
部屋の数=mn、宝の数trでmnＣtr個の組み合わせを返すサブルーチンが必要になって、ここが処理のボトルネックになるんじゃないかと思うのですが。

数学とはなんでしょうか？
何が数学の本質なんでしょうか？
論理的な体系の構築？ 定理の創出？

>>161の若干の一般化とその導出を備忘録的に書いておきます。

>>60
まず、部屋を探る順番が一般の場合を考える。
部屋がNあり、その集合をＲとする。Ａ君、Ｂ君が探る順番を表わす全単射写像をそれぞれf,gとする：
f,g: Ｒ→{0,1,…,N-1} （順番は0から始まるとする。）
部屋自体の位置はなんら答えに影響しない。
σ=g・f^{-1} と置くと、σは{0,1,…,N-1}の置換。（・は写像の合成）
Ａ君がi番目に探る部屋はＢ君がσ(i)番目に探る部屋ということ。
以下、「Ａ君がi番目に探る部屋」のことを「部屋i」ということにする。

求めたいのは、「Ａ君がＢ君よりも早く宝を見つける宝の配置の数」であるが、宝の数をcとすると、それは
　Σ[σ(i)>i] binomial(#{j| j>i, σ(j)>i}, c-1)　（0≦i,j≦N-1、binomialは二項係数）
である。
なぜか？
「Ａ君が初めて宝を見つける部屋（部屋iとする）」で場合分けしよう。
（つまり部屋0〜i-1には宝がなく、部屋iに宝がある場合）
部屋iはＢ君がσ(i)番目に探る部屋だからσ(i)>iでないと
少なくともＢ君はＡ君よりも前か同時に部屋iで宝を見つける
（Ｂ君はその前に別の部屋で宝を見つけることもある）ことになりＡ君は勝てない。
したがって、σ(i)>iが必要。
残りのc-1個の宝は部屋i+1〜N-1にあるが、宝がある部屋を部屋jとすると、
やはりσ(j)>iでないといけない。逆に全部の宝でそうであればＡ君が勝つ。
よって、残りのc-1個の宝が置かれてもいい部屋の数は#{j| j>i, σ(j)>i}だけあり、
全部そこに置かれる場合はbinomial(#{j| j>i, σ(j)>i}, c-1)通り。
したがって、上記のようになる。

続く

>>194
続き

部屋が縦m、横nで、Ａ君は横1行を探し終えたらすぐ下の1行に移り、
Ｂ君は縦1列を探し終えたらすぐ右の1列に移るという場合を考える。
つまり、m=4,n=3の場合、Ａ君は
０１２３
４５６７
８９1011
Ｂ君は
０３６９
１４７10
２５８11
という順番で探す。
このとき、σ=0,3,6,9,1,4,7,10,2,5,8,11。
一般には、σ(nk+l)=ml+k (0≦k≦m-1, 0≦l≦n-1)。

ここまでをPythonで表すと：
#二項係数。SageMathでは定義ずみ
def binomial(n,r):
from math import factorial as f
return f(n)//f(r)//f(n-r) if r>=0 and n-r>=0 else 0

#置換p、宝c個で勝つ宝の配置の数
def nwinperm(p,c):
N = len(p)
return sum(binomial(len([j for j in range(i+1, N) if i<p[j]]),c-1)
for i in range(N) if i<p[i])

#部屋が縦m、横nのときの置換
def rectperm(m,n):
return [m*l+k for k in range(m) for l in range(n)]

#部屋が縦m、横n、宝がc個で横優先が勝つ宝の配置の数
def nwinrect0(m,n,c):
return nwinperm(rectperm(m,n),c)

続く

>>196
続き

部屋がm×(m+1) (n=m+1) のとき。
(m-1)l>(n-1)k ⇔ 0≦k≦m-2 かつ k+1≦l≦m。
(nk+l)/m = k + (k+l)/m より k+l<mのときq=k,r=k+l、k+l≧mのときq=k+1,r=k+l-m。
r>k (k+l<m)とr≦k (k+l≧m)とに分けるように場合分けをする：
�@0≦k≦[(m-1)/2], k+1≦l≦m-k-2 のとき r>k、
�A[(m+1)/2]≦l≦m-1, m-1-l≦k≦l-1 または �Bl=m, 0≦k≦m-2 のとき r≦k。

m=6のとき
�@�@�@�@�A�B
××�@�@�A�A�B
×××�A�A�A�B
××××�A�A�B
×××××�A�B
×××××××

m=7のとき
�@�@�@�@�@�A�B
××�@�@�@�A�A�B
×××�@�A�A�A�B
××××�A�A�A�B
×××××�A�A�B
××××××�A�B
××××××××

後はΣの計算。>>60に合わせるとＱ君がＡ君の立場でmが>>60でのn。
#以下 SageMathコード
,var m,n,l,k,q,r,c
T2 = (n-1-l) + (n-q)*(m-1-k)
T1 = T2 - (r-k)
#mが奇数の場合：
Q1 = (sum(sum(binomial(T1.subs({n:m+1,q:k,r:k+l}),c-1), l,k+1,m-k-2), k,0,(m-1)/2-1)
+ sum(sum(binomial(T2.subs({n:m+1,q:k+1,r:k+l-m}),c-1), k,m-1-l,l-1), l,(m+1)/2,m-1)
+ sum(binomial(T2.subs({n:m+1,l:m,q:k+1,r:k}),c-1), k,0,m-2)
).subs({m:n,c:2}).simplify_full().factor()
#mが偶数の場合：
Q2 = (sum(sum(binomial(T1.subs({n:m+1,q:k,r:k+l}),c-1), l,k+1,m-k-2), k,0,m/2-2)
+ sum(sum(binomial(T2.subs({n:m+1,q:k+1,r:k+l-m}),c-1), k,m-1-l,l-1), l,m/2,m-1)
+ sum(binomial(T2.subs({n:m+1,l:m,q:k+1,r:k}),c-1), k,0,m-2)
).subs({m:n,c:2}).simplify_full().factor()
def Q1st(x):
return (Q1 if mod(x,2) == 1 else Q2).subs({n:x})

続く

>>197
続き

部屋がm×(m-1) (n=m-1) のとき。
(m-1)l>(n-1)k ⇔ 1≦l≦m-2 かつ 0≦k≦l。
(nk+l)/m = k + (l-k)/m より q=k,r=l-k。
r>k (l>2k)とr≦k (l≦2k)とに分けるように場合分けをする：
�@0≦k≦[(m-3)/2], 2k+1≦l≦m-2 のとき r>k、
�A1≦k≦[(m-3)/2], k≦l≦2k または �B[(m-1)/2]≦k≦m-2, k≦l≦m-2 のとき r≦k。

m=7のとき
�@�@�@�@�@
�A�A�@�@�@
××�A�A�A�@
×××�B�B�B
××××�B�B
×××××�B
××××××

m=8のとき
�@�@�@�@�@�@
�A�A�@�@�@�@
××�A�A�A�@�@
×××�B�B�B�B
××××�B�B�B
×××××�B�B
××××××�B
×××××××

後はΣの計算。>>60に合わせるとＰ君がＡ君の立場でmが>>60でのn+1。
#mが偶数の場合：
P1 = (sum(sum(binomial(T1.subs({n:m-1,q:k,r:l-k}),c-1), l,2*k+1,m-2), k,0,m/2-2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,2*k), k,1,m/2-2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,m-2), k,m/2-1,m-2)
).subs({m:n+1,c:2}).simplify_full().factor()
#mが奇数の場合：
P2 = (sum(sum(binomial(T1.subs({n:m-1,q:k,r:l-k}),c-1), l,2*k+1,m-2), k,0,(m-3)/2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,2*k), k,1,(m-3)/2)
+ sum(sum(binomial(T2.subs({n:m-1,q:k,r:l-k}),c-1), l,k,m-2), k,(m-1)/2,m-2)
).subs({m:n+1,c:2}).simplify_full().factor()
def P1st(x):
return (P1 if mod(x,2) == 1 else P2).subs(n=x)

以上、整理して少し異なったけど>>161の導出でした。

>>198
P1 == 1/24*(6*n^3 + 20*n^2 - n - 27)*(n - 1) # nが奇数のとき
P2 == 1/4*n^4 + 7/12*n^3 - 7/8*n^2 - 13/12*n + 1 # nが偶数のとき
Q1 == 1/24*(6*n^2 + 10*n - 3)*(n + 1)*(n - 1) # nが奇数のとき
Q2 == 1/24*(6*n^2 - 2*n - 5)*(n + 2)*n # nが偶数のとき

それだけ前置きやってｋを含めた式が作れないのですか？

>>199
「ｋを含めた式」って何？
あと>>178見て
宝の数が任意のものならΣが取れないでしょう。

ｋが任意でΣを含む高次方程式プリース

論理的に考えて「最強」は存在しませんか？

私は何人かの方のコンピュータによる解法はすごいと思いました
正直、Pythonはわからないし、sagemathは数学そのものなのでまだ理解できていないので
読めたのはＣだけですが…

久しぶりにまともなスレになった気がします

かわいそうな人がいるな。

三角形ABCにおいて、
AからBCへ下した垂線をAD,
BCの中点をMとする。
BD > CDとすると
BD^2 - CD^2 = 2BC･MD を示せ。

>>206
中線定理やろ

>>206
「垂線の定理」と言うらしい

左辺を因数分解すれば簡単

スイセンか
ナルシストやね

BD-CD=2MD
⇔
BD-MD=CD+MD

こんな感じか

>>206
BD^2 - CD^2 = (BD+CD) (BD-CD)
= BC ( (BM+MD) - (CM - MD) )
= 2BC･MD
(∵ BM=CM )

特に垂線である意味がないし「垂線の定理」って何かの間違いでは？

左辺がAB^2 - AC^2なら意味あるな。

>>210
誰が面白い‥‥

>>213
AB^2 - AC^2 = (BD^2 + h^2) - (CD^2 + h^2) = BD^2 - CD^2 = ... = 2 BC・MD
なるほど

自然数a,b,cは以下の2つの等式を共に満たす。
a+b^2=c^3
a^2-b(b+c)=a+b+c

（1）このような(a,b,c)を一組求めよ。
（2）（1）で求めたもの以外に(a,b,c)の組が存在するなら、全て決定せよ。

>>216
Prelude> [(a,b,c)|a<-[1..100],b<-[1..100],c<-[1..100],a+b^2==c^3,a^2-b*(b+c)==a+b+c]

[(4,2,2)]

Prelude> [(b,c)|b <-[1..1000],c<-[1..1000],(c^3-b^2)*(c^3-b^2-1)==(b+c)*(b+1)]
[(2,2)]

変な質問ですいません
最初にピタゴラスの定理を証明した人って、どういう発想で定理が正しいと考えたのでしょうか？直感でしょうか、経験的によく知られていたのでしょうか？

ピタゴラスとエウクレイデスはどっちの方が賢いですか？

>>216
a=-1，b=-3，c=2

a=0，b=-1，c=1
a=0，b=0，c=0
a=-1，b=1，c=0
a=4，b=2，c=2

>>219
証明はどうだかわからないがピタゴラスが多いついたのはタイルを見て予想したと言われているらしい
直角二等辺三角形を敷き詰めると直角二等辺三角形の場合にはピタゴラスの定理が成り立つことがすぐにわかる
ピタゴラスはそこから直角三角形なら常に成り立つのではないかと考えたということのようだ

リーマン予想を証明したいのですが、まずは何から勉強をした方が良いのでしょうか？

証明した人に聞いてください
例えばAtiyahとかdeBrangeさんなど

風呂場の壁のタイルじゃないかね

http://www.wfg-bluebonnet.com/blog/garden-life/%E3%83%94%E3%82%BF%E3%82%B4%E3%83%A9%E3%82%B9%E3%81%AE%E3%82%BF%E3%82%A4%E3%83%AB.2jpg.jpg

リーマン予想が証明されたとしたら、残りの他の全ての数学の未解決問題を自分一人で解決したい。
そのためにはやはり、数学の全分野だけでなく、物理学とか哲学とか計算機科学の全分野も究めないと無理なレベルでしょうか？

とりあえず、二項定理くらいはわかるようになりましょうよ、ヒマラヤさん

>>221
証明は？

b^2=c^3-a
a=c^3-b^2

a^2-b(b+c)=a+b+c
a^2-b^2-bc=a+b+c
a^2-c^3+a-bc=a+b+c
a^2-c^3-bc-b-c=0
(c^3-b^2)^2-c^3-bc-b-c=0

の整数解を求める

>>142
宝の数を変化させるコードをHaskellに移植してみた。

import Data.List
import Data.List.Split
m = 5 -- 縦マス(短軸）
n = 6 -- 横マス(長軸)
k = 5 -- 宝の数
q = [0..m*n-1]
matQ = chunksOf m q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換

combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
treasure = combinations k q -- 宝の組み合わせ
ip y = minimum $ map(\x -> elemIndices x p!!0) y　--　宝の、配列ｐでのindex列を求めて最小値を返す
iq y = minimum $ map(\x -> elemIndices x q!!0) y

idxp = map ip treasure -- 宝の組み合せで実行して
idxq = map iq treasure

p_q = zipWith (-) idxp idxq -- 差をとって大小判別
p1st = length $ filter (<0) p_q　-- 短軸方向探索ｐが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q

main = do
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw


Prelude> :main
p1st = 54036, q1st = 55469, draw = 33001

無限ホテルのパラドックス読んでてわからないことがあって、新しい宿泊客のために既存の客が部屋を一つづつずらすってあるけど、あれは何でそうなるの?
ネットで調べたけどそれらしい答えが無くて困ってる
無限ホテルが集合論のお話で、ホテルは可算無限集合、無限に居る宿泊客全員も可算無限集合で、どっちも無限としての大きさが合うから部屋は過不足なく用意されるって話だってところまではネットで読んだ
で、Wikipediaには順序数? の計算ルールが書いてあって、1+ωとω+1は違うってあったからこれが部屋移動の理由かと最初は思った
でも無限ホテルって無限人の来客があってもokってあるから、これってω+ωでどこに客をぶちこんでも意味変わらないなと
だからこの予想は違うと今は思ってる
この疑問のしっくり来る(理解できる)解説が見つからなくてずっとモヤモヤしてるので、誰か教えてくれるとありがたいです

--　バグ修正（行と列を間違えていた(._.)

import Data.List
import Data.List.Split
m = 5 -- 縦マス(短軸）
n = 6 -- 横マス(長軸)
k = 5 -- 宝の数
q = [0..m*n-1]
matQ = chunksOf n q
matP = transpose matQ --行列を転置して
p = concat matP -- 配列に変換

combinations :: Int -> [a] -> [[a]]
combinations 0 _ = [ [] ]
combinations n xs = [ y:ys | y:xs' <- tails xs, ys <- combinations (n-1) xs']
treasure = combinations k q -- 宝の組み合わせ
ip y = minimum $ map(\x -> elemIndices x p!!0) y　--　宝の、配列ｐでのindex列を求めて最小値を返す
iq y = minimum $ map(\x -> elemIndices x q!!0) y

idxp = map ip treasure -- 宝の組み合せで実行して
idxq = map iq treasure

p_q = zipWith (-) idxp idxq -- 差をとって大小判別
p1st = length $ filter (<0) p_q　-- 短軸方向探索ｐが先に宝をみつける
q1st = length $ filter (>0) p_q
draw = length $ filter (==0) p_q

main = do
putStrLn $ "p1st = " ++ show p1st ++ ", q1st = " ++ show q1st ++ ", draw = " ++ show draw

>matrix.exe
p1st = 55469, q1st = 54036, draw = 33001

>>231
先頭を開けて1人追加するのは 1+ω = ω
倍の部屋番号へ移して ω 人追加するのは 2ω = ω

リーマン予想とP≠NP予想はどっちの方が証明するのが難しいですか？

ABC
EFG

ｎ＝２の６マスでP君Q君のそれぞれのファーストの
組の総数をお願いします<(_ _)>

P勝ち：EG FG EF BF
Q勝ち：BG CG BC CE CF
引き分け：AB AC AD AE AF AG BE

p win : CE, EF, EG, FG
q win : BC, BF, BG, CF, CG
even : AB, AC, AE, AF, AG, BE
かと思った

質問では DEF が EFG になってるのから俺はちゃんとその通りにやってるのにお前らときたら自由だな…
>>236に至ってはよく見るとABCDEFGの7種使ってるし

>>233
ありがとうございます
連続の質問になって申し訳ないんで付けど、2ωとω+ωってこれは違うものなんですか?

自由ついでに分かりやすいように数字に置き換えてみた
１個目だけじゃなく、２個目の宝を先に見つけることも考えたら
結局、PQで差はないという直感どおりの結果になるな

１２３
４５６

１２　・Ｑ
１３　・Ｑ
１４　・Ｐ
１５　・Ｐ
１６　・・
２３　ＱＱ
２４　・Ｐ
２５　ＱＰ
２６　Ｑ・
３４　ＰＱ
３５　Ｑ・
３６　Ｑ・
４５　ＰＰ
４６　Ｐ・
５６　Ｐ・

別スレでこんなの見つけたんですが、これどこで証明されてるかご存知の方います？
https://rio2016.5ch.net/test/read.cgi/math/1537516085/256

＞オイラー定数をγと置く。nの約数の総和をσ(n)と置く。RHは
＞
＞σ(n)＜(e^γ)*n*log(log n)　(∀n＞5040)
＞
＞と同値であることが知られている。

自己レス

とりあえず元論文はコレらしい

[24] G. Robin, Grandes valeurs de la fonction somme de diviseurs et hypoth`ese de Riemann,
J. Math. Pures Appl. 63 (1984), 187–213.

英語で読めるのないかなぁ？

５×６マスで宝の数を１０まで増やしていくと、

D:\bin>for %i in (1,2,3,4,5,6,7,8,9,10) do treasure 5 6 %i

D:\bin>treasure 5 6 1
p1st = 14, q1st = 14, draw = 2

D:\bin>treasure 5 6 2
p1st = 203, q1st = 197, draw = 35

D:\bin>treasure 5 6 3
p1st = 1801, q1st = 1727, draw = 532

D:\bin>treasure 5 6 4
p1st = 11418, q1st = 11008, draw = 4979

D:\bin>treasure 5 6 5
p1st = 55469, q1st = 54036, draw = 33001

D:\bin>treasure 5 6 6
p1st = 215265, q1st = 211894, draw = 166616

D:\bin>treasure 5 6 7
p1st = 685784, q1st = 680768, draw = 669248

D:\bin>treasure 5 6 8
p1st = 1827737, q1st = 1825076, draw = 2200112

D:\bin>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184

D:\bin>treasure 5 6 10
p1st = 7995426, q1st = 8023257, draw = 14026332



１：同等
１〜８：短軸探索有利
９、１０：長軸探索有利
という結果になった。

Haskellのコードはここ

--exe Fileにコンパイルしてコマンドラインから実行できるように改変（但し、エラー処理皆無）
https://rio2016.5ch.net/test/read.cgi/math/1490734993/209

>>235
$Rscript main.r
P1st Q1st even
3 4 13

６マスで宝を３個にしてみた

$Rscript main.r
P1st Q1st even
3 4 13

P 1st
[,1] [,2] [,3]
[1,] C C D
[2,] D D E
[3,] E F F

Q 1st
[,1] [,2] [,3] [,4]
[1,] B B B C
[2,] C C E E
[3,] E F F F

even
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,] A A A A A A A A A A B B B
[2,] B B B B C C C D D E C D D
[3,] C D E F D E F E F F D E F

Rのスクリプトをここに置いたから数値を変更して実行可能

http://tpcg.io/X3bC0A

>>241
自由ついでに５×６マスで宝が５個、先に全部の宝を見つけた方が勝者とすると

これくらい差が出る

D:\bin>treasure2 5 6 5

p1st = 54036, q1st = 55469, draw = 33001

確率にすると
> treasure2(5,6,5)
P1st Q1st even
54036 55469 33001
> 54036/(54036+55469+33001)
[1] 0.379184
> 55469/(54036+55469+33001)
[1] 0.3892398

なので差があると直観するかどうかは個々人の感性だな。

Haskel だの R だの。そういうのは他でやれよ...。

>>247
使える人間にとっては電卓みたいなもんだよ。
log2の計算にいちいちマクローリン展開して手書き計算しないだろ。

ここまでjuliaが出てこなかった
juliaが流行しているのは自分の周りだけなのかな
（NGに巻き込まれて見えてないだけだったらゴメン）

>>249
どれかを移植して実力を示していただけたらうれしい。

５✕６マスで宝が１５個の時の計算とかまだ誰も出してない。

>>241
先に２個の宝をみつけた方なら

123
456

12 Q
13 Q
14 P
15 P
16 =
23 Q
24 P
25 P
26 =
34 Q
35 =
36 =
45 P
46 =
56 =

にならない？

>>241
２個を先にみつけるじゃなくて

これは１個めの発見はQの方が確率が高くて、２個めに発見はPの方が確率が高いというだけの話だったみたいね。

>>250
juliaが周りで流行ってるだけで自分自身はＣの人（Ｃソース書いてくれた人とは別人）

5x6ますで宝１５個とか、ID:TZqGbv4dにお願いしたらすぐやってくれるんじゃない？
完全に作り直してるし。

>>252
i 番目をどちらが先にみつけるかを計算してみた。

４×５マスに宝が５個あるとき

> treasures(4,5,5)
p1st q1st even
[1,] 1948 9680 3876
[2,] 5488 10016 0
[3,] 7752 7752 0
[4,] 10016 5488 0
[5,] 9680 1948 3876

１個め２個めは短軸方向探索のQが、４個め５個めは長軸方向探索のPが、先にみつける宝の配置の組み合わせが多い。３個めは同じ。

全体としてはイーブンだが、

勝者は１個めを先にみつけた方にするか、全部を先にみつけた方にするかで結果が変わる。

Rのコードはここに置いたので数値を変えて実行可能。

http://tpcg.io/Ph7TUQ

n元集合からk個の元を取り出す順番を考慮して可能な場合を数え上げるとn*(n-1)*....*(n-k+1)通りあるというのはより原始的なものから導かれるものですか？

>>248
「プログラムで、ごり押し計算」
「マクローリン展開して手書き計算」
俺は後者の方が美しく感じるけどな。
実は前者で計算したのに後者を装ってほしいくらい。 (※ 私見です)

>>254
単なる確認なんだけれども、
「i番目をどちらが見つけるか」というのは
先にi個見つけた方を勝ちとするのではなくて

例えばi=2だと
Pが１つ発見、Qが１つ発見⇒2番目を見つけたQの勝ち、ということですか？

ム板でやれ

>>257
>254の計算は各人にとってi番めの計算。
例えばi=2だと
Pが１つ発見、Qが１つ発見だと勝敗は未決で
どちらが発見者にとって２個めを発見したらそれが勝者として数えた。

んで、
ここまで答が出せた

254 名前：１３２人目の素数さん[sage] 投稿日：2018/10/30(火) 13:03:57.24 ID:TZqGbv4d
>>252
i 番目をどちらが先にみつけるかを計算してみた。

４×５マスに宝が５個あるとき

> treasures(4,5,5)
p1st q1st even
[1,] 1948 9680 3876
[2,] 5488 10016 0
[3,] 7752 7752 0
[4,] 10016 5488 0
[5,] 9680 1948 3876

１個め２個めは短軸方向探索のQが、４個め５個めは長軸方向探索のPが、先にみつける宝の配置の組み合わせが多い。３個めは同じ。

全体としてはイーブンだが、

勝者は１個めを先にみつけた方にするか、全部を先にみつけた方にするかで結果が変わる。

Rのコードはここに置いたので数値を変えて実行可能。

http://tpcg.io/Ph7TUQ

>>260
全体を眺めると直感的通り互角。
局所でみると濃淡があるということと理解した。

ゲルト・ファルティングスとアラン・コンヌの知能指数はどれくらいですか？

「真理」というのは存在するのでしょうか？
「真理」の探究は意味があるのでしょうか？

マイケル・アティヤとエドワード・ウィッテンはどっちの方が賢いですか？

>>253
ここまでは算出できたが、宝を１４にしたらエラー終了した。

D:\bin>treasure 5 6 11
p1st = 13346984, q1st = 13395944, draw = 27884372

D:\bin>treasure 5 6 12
p1st = 19312228, q1st = 19372871, draw = 47808126

D:\bin>treasure 5 6 13
p1st = 24301031, q1st = 24358063, draw = 71100756

>>92
■引き分けの組み合わせは勝敗と無関係なので除外

宝が２個以上の時、
スタート地点のAマスと対極にある最終マスのLには
P君もQ君もどちらも決してたどり着くことはできないので
このLマスと組みとなる宝の配置は重複情報で意味を持たない
ので除外する

Pが先に見つけるのは以下の21通り
CE,DE,DI,EF,EG,EH,EI,EJ,EK,FG,FH,FI,FJ,FK,GI,GJ,HI,HJ,IJ,IK,JK,

Qが先に見つけるのは以下の22通り
BC,BD,BF,BG,BH,BI,BJ,BK,CD,CF,CG,CH,CJ,CK,DF,DG,DH,DJ,DK,GH,GK,HK,

となる

>>260
何マスだろうが、宝が何個であろうが
出発点と終点が同じであれば
PQの宝を得られる個数の期待値は同じということだな

>>266
BE　と　CI　の扱いは？

>>267
期待値は宝の数なわけで、元の問題は１個めをみつけるステップの数を比較しているんだと思う。

>>267
宝を先にみつけたら独り占め、同時にみつけたら折半　というルールなら手に入れる宝の数の期待値は同じになるだろうね。

>>268
BI, CIは引き分けで除外なのでは

>>265
メモリを食わないコードを書いてみた
今思ったけど再帰で書いた方が読みやすかったか
https://ideone.com/HaAqJO

>>270
n x n+1 の部屋を縦横に調べる2人の場合は先着する部屋数が等しくなるからそうなるね
p君 ABCDEFGHIJKL
q君 BCDEFGHIJKLA
とかなら殆どq君が独り占め

>>271
https://ideone.com/HaAqJO

ありがとうございます、これを待望しておりました。

P君Q君問題から得られる知見
早い者勝ちなら先回りすることが勝つ秘訣

>>274
>244に数値を挙げたけど宝の数が増えると逆転しちゃう。
個人的にはどこが逆転する境なのか算出方法が知りたいところ。

>>233
質問してばっかりだったので反省して自分で調べてみたんですけど
ω+ω=ω×ω=ω2だってことでした
でも、これは「無限ホテルのω号室の次の部屋からω人の客を泊めた」って事ですよね?
だから無限ホテルの話にあるように1号室→2号室、2号室→4号室、3号室→6号室とずらして、間に入れ込めば2ω=ωになって万事解決って事で合ってますか?

>>276
ω×ωはちがうかった……
これじゃω^2になっちゃう

q1..q2..q3..q4
q5..q6..q7..q8
q9q10q11q12

p1..p4..p7..p10
p2..p5..p8..p11
p3..p6..p9..p12

同じ座標なら数字の小さいほうが勝ち

［q2とq10］ & ［p4とp6］に宝が配置された時は

互いに数字の小さいほうを選んで勝負

q2 vs p4 で　q2の勝ちとなる

この後にq10とp6の探査をしても
情報としての価値はゼロ

>>276
＞だから無限ホテルの話にあるように1号室→2号室、2号室→4号室、3号室→6号室とずらして、間に入れ込めば2ω=ωになって万事解決って事で合ってますか?

そうそう

まだ続いているようなので、>>189の後半で示したようなアイデアで、宝の数可変版の
プログラムを書いてみました。
多倍長を使える処理系を用いればいいのかもしれませんが、実数型で誤魔化しました。
故に大きな数字のところでは誤差があります。

http://codepad.org/VN03aiqT

>>281
同じ方針のものがPythonで>>194-198にある

てか>>194-199に書いてある事がちゃんと読めれば宝の数が何個になっても場合わけ＋多項式で記述できるのはすぐわかる。
読めよ。数学板なんだから。

>>281
いつもありがとうございます。
いやぁ、この出力は圧巻ですね。
Haskell先生もびっくり。

>>280
ありがとうございます
おかげさまですごくしっくりきました

>>282　>>283
失礼しました。
数列を無理矢理分数式化する人や、価値の無い長い文章を投下する人がいるので、
読み飛ばしていました。
宝箱が二つの場合は、多項式での表現が完成していたんですね。
あのようなσやδを含む式を整理する数式処理ツールがあったとは驚きです。
二個で可能だったのだから、もっと多くの場合でも、可能なんでしょうね。

P1 == 1/24*(6*n^3 + 20*n^2 - n - 27)*(n - 1) # nが奇数のとき
P2 == 1/4*n^4 + 7/12*n^3 - 7/8*n^2 - 13/12*n + 1 # nが偶数のとき
Q1 == 1/24*(6*n^2 + 10*n - 3)*(n + 1)*(n - 1) # nが奇数のとき
Q2 == 1/24*(6*n^2 - 2*n - 5)*(n + 2)*n # nが偶数のとき

多項式ってこれだけ？
ｋは変えられないし出力は意味不明だしナニコレ？

>>204の式ならｋ＝５５４２２２，ｎ＝３２２３００９８８とかでも
数秒で出力してくれるよ

>>204の式ならkにどんな整数をいれても正解にならん。n=3でやってみろよ。
でn=3の場合66通り全部書きだして比較してみろよ。
実際書き出してみた正解とひとつも合わない式になんの意味がある？

>>204の式は１１C２＝５５通りで計算してある

>>289
Prelude Data.Ratio> print [(n+1)*(n^2+2*n-1-k)%(n^2*(n+2)-n*k)|let n = 3,k<-([0..14]++[16..30])]
[56 % 45,26 % 21,16 % 13,11 % 9,40 % 33,6 % 5,32 % 27,7 % 6,8 % 7,10 % 9,16 % 15,1 % 1,8 % 9,2 % 3,0 % 1,8 % 3,2 % 1,16 % 9,5 % 3,8 % 5,14 % 9,32 % 21,3 % 2,40 % 27,22 % 15,16 % 11,13 % 9,56 % 39,10 % 7,64 % 45]
Prelude Data.Ratio>
kに0〜30何入れても正解なんぞ出てこんやろ？

>>290
ｋに５００〜８００００だとどうですか？

ｋ＝５５４２９９７４７２１２，ｎ＝３２１２３０１０９８８５５

でも出力できたよ
ためしてごろうじろう

>>291
k>15だとすべて4/3より大きい値しかでないからアウト。何入れてもだめ。

>>292
n = 3〜100までいれて全滅の式にそんな値いれても糞の意味もない。

正確に一致しなくてもどちらが勝者になるかが
わかればいいと思う

ｋ＝５７２３４５７７５４２９９７４７２１２，ｎ＝３２１２３０１０９８８５５でも
出力できたぞ

>>286
>あのようなσやδを含む式を整理する数式処理ツールがあったとは驚きです。
σを処理できないから>>195-196、δを処理できないから>>196-197を人手で行っている
SageMathにやらせているのはn乗の和の公式さえあれば高校生ができる計算

>二個で可能だったのだから、もっと多くの場合でも、可能なんでしょうね。
>>196-197のsubs({m:n+1,c:2})の2を3に変えれば宝が3個の場合の多項式が得られる

>>295
ｋ＝１７４５６６１９２５１，ｎ＝１３２１２３でちゃんと１が出力される

さすが

>>294
あほか？n=3〜100で正しい数値出してない式になんの信憑性がある？
正しい答え出なきゃなんの意味もない。

>>284
>>195-196のPythonをHaskellにすればいい
Haskellにもリスト内包表記があるんだから

>91で
　読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。
と書いたが犠牲者が出ているようだな。

間違えること自体は悪いことじゃないから、間違えたことがわかれば間違えたと書いておくか
そのまま消えてしまうだけで別にかまわないのに。

なにが無駄ってこいつ
>>204
＞計算知能にそのまま入力するだけで通分と約分を
＞自動計算してくれるので試してごろうじろう
＞
＞■Wolfram入力例
＞
＞（ｎ＋１）（ｎ^2＋２ｎ−１−ｋ）／｛ｎ^2（ｎ＋２）−ｎｋ｝，ｋ＝２，ｎ＝３
ってわざわざ全角で書いてコピペで入力できなくできないようにしてくれてる所。
wolfram 日本語版だけは全角でも入力できるけどその他のツールは全滅。
いちいち半角に打ち直さんといかん。
脳みそ1ccしかないんちゃうかと。

>>302
俺は>205の助言の意味が分かったので>204のidを速攻でNGidに登録したよ。

日本製のエディタには全角半角変換できるのがあるよ。

例えば、http://mana.ikuto.com/

f(x)=-x²+ax+bがあり, y=f(x)は点(-2,1)を通る。
x∈[-3,3]で動くとき最大値Mと最小値mを, aについて次の2つ場合分けすることによって与えよ。
(1)a≧??のとき,
x=3でM=??a-??,
x=-3でm=-??a-??
(2)a≦??のとき,
x=-3でM=-??a-??,
x=3でm=??a-??

となっているのですが、これで場合分けは足りているのですか？

いや、そもそも数学の掲示板で数式全角で書いてる時点でアホだよ。
あとで数式コピペしてソフトに貼り付けるなんて普通にするじゃん。
* はさすがに見苦しいから我慢するけど、全部大文字にするのは意味わからん。
しかも
＞計算知能にそのまま入力するだけで通分と約分を
＞自動計算してくれるので試してごろうじろう
といいながらだよ？
アホじゃね？

>>299
御助言にしたがってHaskelｌに移植しました。

import System.Environment

choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]

nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0

nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin n m c

main = do
argList <- getArgs -- m : 縦マス(短軸） n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)

おかげ様でこういうのも瞬時に計算してくれました。

１０×２０マスで宝が１００個

>takara 10 20 100
p1st = 15057759425309840160151925452579572328997602171271937639470, q1st = 15057796557877993527038542474310161591275806044157319150135, draw = 60432921540347294111327092128863840691952977587098698541050

>>305
数学板は例外かもしれないが、マクロウイルスが貼られるのの予防か半角で投稿すると拒絶されることがあるな。
httpを貼ろうとするとはねられるときには全角にすることもあるな。まあ、数文字大文字に留めるけど。

>>304
誰も答えていないしみんな困ってるんだと思うが、すべての場合を調べているわけではない、と考えればいいだけの話。
というかそうとしか捉えられないw

>>306
既約分数で表示してくれ

a,b,cは自然数とする。
このとき、以下の不等式を満たす(a,b,c)が存在するような自然数Nの最大値を求めよ。
N≦a^2+b^2+c^2≦2018

>>306
タイプミスで draw が間違ってますよ

>>312

ご指摘ありがとうございました。

×　draw m n c = choose(m*n,c) - nwin m n c - mwin n m c

○　draw m n c = choose(m*n,c) - nwin m n c - mwin m n c

>>306
ご指摘を受けたのでデバッグしたのを投稿します。

import System.Environment

choose (n,r) = product[1..n] `div` product[1..n-r] `div` product[1..r]

nloc m n k l = do
let q = div (n*k+l) m
r = mod (n*k+l) m
in (n-q)*(m-k) + q-1-l + if r>k then k-r else 0

nwin m n c = sum[choose ((nloc m n k l), c-1) | k<-[0..m-1], l<-[0..n-1], k*(n-1) < l*(m-1)]
mwin m n c = sum[choose ((nloc n m k l), c-1) | k<-[0..n-1], l<-[0..m-1], k*(m-1) < l*(n-1)]
draw m n c = choose(m*n,c) - nwin m n c - mwin m n c

main = do
argList <- getArgs -- m : 縦マス(短軸） n : 横マス(長軸) k : 宝の数
let m = read (argList !! 0)
n = read (argList !! 1)
k = read (argList !! 2)
putStrLn $ "p1st = " ++ show(mwin m n k) ++ ", q1st = " ++ show(nwin m n k) ++ ", draw = " ++ show(draw m n k)

>>304
変な問題だけど、次の2つの場合、すなわち
・x=3のとき最大、x=-3のとき最小 (a≧6のときか？)
・その逆 (a≦-6のときか？)
に分けて??を埋めよという問題なのだろうから、
その2つのときだけ考えて答えれば良いのではないだろうか

「分けて」ってのが変だよね
次の2つの場合について、ならわかるんだけど。

俺の最大の夢は、「「無」になってもう二度と「有」にならない」ことだ。
どうすればこれを実現できるのでしょうか？
自殺をしても無駄なのでしょうか？

高専2年
行列の固有値と対角化
(4)が全然わかりません
よろしくお願いします

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

先に１個めの宝を見つけるには短軸探索と長軸探索とどちらが有利かは宝の数によって変わるのでグラフにしてみた。
縦５横６のとき宝の数を１から３０まで増やして長軸探索が先にみつける確率と短軸探索がさきにみつける確率の差を描いてみた。
http://i.imgur.com/7qGjOJX.png
縦５横６のときだと宝の数は９から２１のときが長軸探索が有利となった。
短軸有利→長軸有利→同等となるようで、再逆転はないもよう。
縦ｍ横ｍ＋１として長軸探索が有利になる宝の数の上限と下限を算出してみた。

[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
[1,] 0 2 2 6 9 13 17 23 29 36 43 52 61 71 82 93 105 118 132 147
[2,] 0 3 7 13 21 31 43 57 73 88 105 118 135 152 166 185 202 220 242 253

グラフにしてみた。
http://i.imgur.com/PiL9xyH.png

>>311
これをお願いします

>>311
a=44,b=9,c=1のとき2018-a^2-b^2-c^2=0

2018-a^2-b^2-c^2，a=44,b=9,c=1

∴N=2018

N＝2018　
(a,b,c)＝(44,9,1)、(43,12,5)

2018-a^2-b^2-c^2，a=41,b=16,c=9

∴N=2018

(44,9,1)
(43,12,5)
(41,16,9)
(35,27,8)
(34,29,11)
(33,23,20)

a=36,b=19,c=19
a=35,b=28,c=3
a=35,b=27,c=8

∴N=2018

>>324
(34,29,11)は違う

計算機実験は大事だと思うけどダンプリストみたいなの延々載せられてもなんかもにょる。

>>317,327

a[n] = a^n、b[n] = ((b+c)^n + (b-c)^n)/2、c[n] = ((b+c)^n - (b-c)^n)/2とおいて
A^n = [[a[n],0,0],[0,b[n],c[n]],[0,c[n],b[n]]]、
[1,0,0]A = a[1,0,0]A、[0,1,1]A = (b+c)[1,0,0]、[0,0,1]A = (b-c)[1,0,0]A。
(1) c[-1]。
(2) a,b+c,b-c。
(3) a^2=4 ⇔ a=±2、
　 (b+c)^2 = 2i、(b-c)^2 = -2i ⇔ (b,c) = (1, i)、(-1, -i)、(i, 1)、(-i, -1)。
(4) b[n] + c[n] = (b+c)^n。