18/10/31 01:39:59.62 VK521Oc+.net
>>286
>あのようなσやδを含む式を整理する数式処理ツールがあったとは驚きです。
σを処理できないから>>195-196、δを処理できないから>>196-197を人手で行っている
SageMathにやらせているのはn乗の和の公式さえあれば高校生ができる計算
>二個で可能だったのだから、もっと多くの場合でも、可能なんでしょうね。
>>196-197のsubs({m:n+1,c:2})の2を3に変えれば宝が3個の場合の多項式が得られる
309:１３２人目の素数さん
18/10/31 01:50:39.23 JttzkDdq.net
>>295
ｋ＝１７４５６６１９２５１，ｎ＝１３２１２３でちゃんと１が出力される
さすが
310:１３２人目の素数さん
18/10/31 01:51:22.71 Demuw4Zw.net
>>294
あほか？n=3～100で正しい数値出してない式になんの信憑性がある？
正しい答え出なきゃなんの意味もない。
311:１３２人目の素数さん
18/10/31 01:52:28.70 VK521Oc+.net
>>284
>>195-196のPythonをHaskellにすればいい
Haskellにもリスト内包表記があるんだから
312:１３２人目の素数さん
18/10/31 06:08:55.55 2LxBlHwr.net
>91で
読んだ人の時間を無駄遣いさせるような明らかな誤答は慎めよ。
と書いたが犠牲者が出ているようだな。
313:１３２人目の素数さん
18/10/31 10:23:25.66 k/QZWhBY.net
間違えること自体は悪いことじゃないから、間違えたことがわかれば間違えたと書いておくか
そのまま消えてしまうだけで別にかまわないのに。
314:１３２人目の素数さん
18/10/31 10:29:51.18 PPhF82WW.net
なにが無駄ってこいつ
>>204
＞計算知能にそのまま入力するだけで通分と約分を
＞自動計算してくれるので試してごろうじろう
＞
＞■Wolfram入力例
＞
＞（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝，ｋ＝２，ｎ＝３
ってわざわざ全角で書いてコピペで入力できなくできないようにしてくれてる所。
wolfram 日本語版だけは全角でも入力できるけどその他のツールは全滅。
いちいち半角に打ち直さんといかん。
脳みそ1ccしかないんちゃうかと。
315:１３２人目の素数さん
18/10/31 13:58:47.90 6U/VyaCA.net
>>302
俺は>205の助言の意味が分かったので>204のidを速攻�
316:ﾅNGidに登録したよ。 日本製のエディタには全角半角変換できるのがあるよ。 例えば、http://mana.ikuto.com/
317:１３２人目の素数さん
18/10/31 14:39:11.34 D1u5pYAL.net
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-??
となっているのですが、これで場合分けは足りているのですか？
318:１３２人目の素数さん
18/10/31 14:39:54.65 /RfK3tjD.net
いや、そもそも数学の掲示板で数式全角で書いてる時点でアホだよ。
あとで数式コピペしてソフトに貼り付けるなんて普通にするじゃん。
* はさすがに見苦しいから我慢するけど、全部大文字にするのは意味わからん。
しかも
＞計算知能にそのまま入力するだけで通分と約分を
＞自動計算してくれるので試してごろうじろう
といいながらだよ？
アホじゃね？
319:１３２人目の素数さん
18/10/31 14:59:23.31 6U/VyaCA.net
>>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
320:１３２人目の素数さん
18/10/31 15:03:15.11 6U/VyaCA.net
>>305
数学板は例外かもしれないが、マクロウイルスが貼られるのの予防か半角で投稿すると拒絶されることがあるな。
httpを貼ろうとするとはねられるときには全角にすることもあるな。まあ、数文字大文字に留めるけど。
321:１３２人目の素数さん
18/10/31 17:14:52.73 42bMLcC4.net
>>304
誰も答えていないしみんな困ってるんだと思うが、すべての場合を調べているわけではない、と考えればいいだけの話。
というかそうとしか捉えられないw
322:１３２人目の素数さん
18/10/31 17:39:16.62 JttzkDdq.net
>>306
既約分数で表示してくれ
323:１３２人目の素数さん
18/10/31 18:31:09.78 JttzkDdq.net
>>305
P(A)をP(B)で割ることによって
P君の勝つ数とQ君の勝つ数が導ける
P(A)／P(B)＝（P君の勝つ数）／（Q君の勝つ数）
｛ｎ（ｎ＋２）－ｋ－１｝／｛ｎ^2（ｎ＋１）－ｋｎ｝
P(A)／P(B)＝――――――――――
｛ｎ（ｎ＋２）－ｋ｝／｛ｎ（ｎ＋１）^2－ｋ（ｎ＋１）｝
＝（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝
∵［ｎ≧２，ｎ（ｎ＋１）－１＞ｋ≧１］
∵の範囲でｎとｋの数値をいろいろと変えることにより
様々な勝率が導ける
計算知能にそのまま入力するだけで約分を
自動計算してくれるので試してごろうじろう
■Wolfram入力例
（ｎ＋１）（ｎ^2＋２ｎ－１－ｋ）／｛ｎ^2（ｎ＋２）－ｎｋ｝，ｋ＝２，ｎ＝３
(n+1)(n^2+2n-1-k)/｛n^2(n+2)-nk｝,k=2,n=3
スタート地点のAマス以外のすべてのマスに
宝がある状態であるｋ＝ｎ（ｎ＋１）－１の時、
必ずP(A)／P(B)＝１になる
ｋ＝ｎ（ｎ＋１）－１の時にP(A)／P(B)≠１となるｎを
見つけることができれば反例になる
見つけてごろうじろう
324:１３２人目の素数さん
18/10/31 18:52:26.89 3bIZfida.net
a,b,cは自然数とする。
このとき、以下の不等式を満たす(a,b,c)が存在するような自然数Nの最大値を求めよ。
N≦a^2+b^2+c^2≦2018
325:１３２人目の素数さん
18/10/31 18:54:45.40 Xi/4xckY.net
>>306
タイプミスで draw が間違ってますよ
326:１３２人目の素数さん
18/10/31 19:12:43.10 6U/VyaCA.net
>>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
327:１３２人目の素数さん
18/10/31 19:15:52.62 6U/VyaCA.net
>>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)
328:１３２人目の素数さん
18/10/31 19:41:15.07 Xi/4xckY.net
>>304
変な問題だけど、次の2つの場合、すなわち
・x=3のとき最大、x=-3のとき最小 (a≧6のときか？)
・その逆 (a≦-6のときか？)
に分けて??を埋めよという問題なのだろうから、
その2つのときだけ考えて答えれば良いのではないだろうか
「分けて」ってのが変だよね
次の2つの場合について、ならわかるんだけど。
329:１３２人目の素数さん
18/10/31 20:07:35.11 oZeu8G8O.net
俺の最大の夢は、「「無」になってもう二度と「有」にならない」ことだ。
どうすればこれを実現できるのでしょうか？
自殺をしても無駄なのでしょうか？
330:１３２人目の素数さん
18/10/31 21:37:19.52 EEWI02Z3.net
高専2年
行列の固有値と対角化
(4)が全然わかりません
よろしくお願いします
URLﾘﾝｸ(i.imgur.com)
331:１３２人目の素数さん
18/10/31 21:49:44.27 6U/VyaCA.net
先に１個めの宝を見つけるには短軸探索と長軸探索とどちらが有利かは宝の数によって変わるのでグラフにしてみた。
縦５横６のとき宝の数を１から３０まで増やして長軸探索が先にみつける確率と短軸探索がさきにみつける確率の差を描いてみた。
URLﾘﾝｸ(i.imgur.com)
縦５横６のときだと宝の数は９から２１のときが長軸探索が有利となった。
短軸有利→長軸有利→同等となるようで、再逆転はないもよう。
縦ｍ横ｍ＋１として長軸探索が有利になる宝の数の上限と下限を算出してみた。
[,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
グラフにしてみた。
URLﾘﾝｸ(i.imgur.com)
332:１３２人目の素数さん
18/10/31 21:58:51.26 3bIZfida.net
>>311
これをお願いします
333:１３２人目の素数さん
18/10/31 22:29:57.00 ldFHIXo+.net
>>311
直感の概算 (a,b,c)＝(40,20,4) N＝2016
微調整 (a,b,c)＝(44,9,1) N＝2018
なんか問題を勘違いしてるかな？
334:１３２人目の素数さん
18/10/31 22:31:04.22 JttzkDdq.net
>>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
335:１３２人目の素数さん
18/11/01 00:10:20.85 AEjEpZy5.net
N＝2018
(a,b,c)＝(44,9,1)、(43,12,5)
336:１３２人目の素数さん
18/11/01 00:29:38.31 b1wO9L0a.net
2018-a^2-b^2-c^2，a=41,b=16,c=9
∴N=2018
337:１３２人目の素数さん
18/11/01 00:31:21.78 AEjEpZy5.net
(44,9,1)
(43,12,5)
(41,16,9)
(35,27,8)
(34,29,11)
(33,23,20)
338:１３２人目の素数さん
18/11/01 00:53:58.60 b1wO9L0a.net
a=36,b=19,c=19
a=35,b=28,c=3
a=35,b=27,c=8
∴N=2018
>>324
(34,29,11)は違う
339:１３２人目の素数さん
18/11/01 02:54:39.51 vtjUzc7H.net
計算機実験は大事だと思うけどダンプ
340:リストみたいなの延々載せられてもなんかもにょる。
341:１３２人目の素数さん
18/11/01 05:19:46.53 xVnRbBm5.net
>>317
17.27 正則行列Ａ = { [a,0,0] [0,b,c] [0,c,b] } について，次の問に答えよ。(九大*)
(1) 行列Ａの逆行列Ａ^(-1) の (2,3) 成分を求めよ。
(2) Ａの固有値を求めよ。
(3) Ａ^2 = { [4,0,0] [0,0,2i] [0,2i,0] } を満たす a,b,c の値を求めよ。iは虚数単位。
(4) nを自然数とし，Ａ^n (i,j)は行列Ａ^n の(i,j)成分を表わすものとする。
そのとき、Ａ^n (2,2) + Ａ^n (3,2) を n, b, c を用いて表わせ。
342:１３２人目の素数さん
18/11/01 05:35:54.50 GatmQtrC.net
>>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。
343:１３２人目の素数さん
18/11/01 06:05:53.60 xVnRbBm5.net
>>317 >>327
(1)
det(Ａ) = a(bb-cc),
Ａ^(-1) = { [1/a，0，0] [0，b/(bb-cc)，-c/(bb-cc)] [0，-c/(bb-cc)，b/(bb-cc)] }
(2)
det(Ａ-λE) = det{ [a-λ，0，0] [0，b-λ，c] [0，c，b-λ] }
= (a-λ)(b-c-λ)(b+c-λ)
∴ λ = a，b±c,
(3)
Ａ^2 = { [a^2，0，0] [0，bb+cc，2bc] [0，2bc，bb+cc] }
∴ a = ±2，(b,c) = (0,±(1+i)) (±(1+i),0)
(4)
Ａ^n = { [a^n，0，0] [0，f_n，g_n] [ 0，g_n，f_n] }
ただし、f_n = {(b+c)^n + (b-c)^n}/2， g_n = {(b+c)^n - (b-c)^n}/2,
(f_n)^2 - (g_n)^2 = (bb-cc)^n,
あとは自分で考えて
344:１３２人目の素数さん
18/11/01 07:42:27.89 xCdOvDq8.net
>>311
Nの最大値は２０１８
顰蹙のプログラム解
Prelude> [(a,b,c)|a<-[1..45],b<-[a..45],c<-[b..45], a^2+b^2+c^2==2018]
[(1,9,44),(3,28,35),(5,12,43),(8,27,35),(9,16,41),(19,19,36),(20,23,33)]
345:１３２人目の素数さん
18/11/01 10:00:52.96 ZI9FoIBR.net
>>326
そこから規則性が見いだせれば理論はあとからついてきたりすることもあるからね。
コラッツの問題みたいに未決のままのもあるけど。
346:１３２人目の素数さん
18/11/01 12:19:48.42 +Vmpp6Zg.net
ここでコード書いてるやつは規則見出して解くなんて気持ちサラサラないやろ？
プログラム書いて遊んでるだけ。
数学的な解出てもガン無視してるし。
347:１３２人目の素数さん
18/11/01 13:15:12.78 xCdOvDq8.net
処理速度が不十分なインタープリタでのコードをコンパイラのコードに移植してくれるのはとても勉強になるので嬉しいね。
>312のような指摘もとてもありがたい。
348:１３２人目の素数さん
18/11/01 13:50:05.89 X0yV8qdr.net
遊ぶなら自分一人でやってればいいのにね。
こんなんできた～ってひけらかしたいんだろ？
349:１３２人目の素数さん
18/11/01 14:33:43.50 DGlwDrwF.net
荒らしよりは意味がある
350:１３２人目の素数さん
18/11/01 15:18:50.87 02Lyc5pT.net
PCでのシミュレーション解を越えた解析解が出たら
それを検証して解析解をPCでの計算に応用。
おかげで>142から>314に進化できた。
プログラミングのトレーニング課題を与えてくれた方に深謝。
引き分けのバグ指摘にも感謝。
数理
351:展開が勉強になるようにコードの議論も俺には嬉しい。 このスレではじめてHaskellの存在を知った初心者なので>299のような適格なアドバイスは嬉しいね。
352:１３２人目の素数さん
18/11/01 16:03:45.77 yg4Nrziz.net
nを自然数、aを実数とするとき、
x^2-(4n-3)x+a/(n^2+n+1) ≦ 0
を満たす整数xが存在するためにn,aが満たすべき条件を述べよ。
353:１３２人目の素数さん
18/11/01 18:18:17.69 AEjEpZy5.net
ａのｂ乗×ｃのｄ乗＝ａｂｃｄ
ａｂｃｄに当てはまる数字は？
※答は１通りしかないようです。
354:１３２人目の素数さん
18/11/01 18:25:46.20 yg4Nrziz.net
>>338
aとcで割れば？
細かい条件は自分でやって
355:１３２人目の素数さん
18/11/01 19:16:38.68 xCdOvDq8.net
>>338
１を許すと沢山ある(1,1,1,1),(1,1,2,1),(1,1,2,2),(1,1,3,1),(1,1,4,1),(1,1,5,1),(1,1,6,1).....けど
(2,2,2,2)が答？
356:１３２人目の素数さん
18/11/01 19:37:28.66 AEjEpZy5.net
>>340
ごめん、ａｂｃｄ は４桁の整数
357:１３２人目の素数さん
18/11/01 19:47:15.96 OC3wBzdi.net
Prelude Data.Ratio> [x | a<-[0..9],b<-[0..9],c<-[0..9],d<-[0..9],let x = 1000*a + 100*b+10*c + d, x == a^b*c^d]
[2592]
358:１３２人目の素数さん
18/11/01 19:47:27.74 AEjEpZy5.net
a^b＋c^d＝1000a+100b+10c+d
359:１３２人目の素数さん
18/11/01 19:52:36.37 AEjEpZy5.net
>>343 訂正
× a^b＋c^d＝1000a+100b+10c+d
○ (a^b)*(c^d)＝1000a+100b+10c+d
360:１３２人目の素数さん
18/11/01 19:59:05.61 xCdOvDq8.net
>>344
これまた顰蹙のダンプリストｗ
Prelude> [(a,b,c,d)|a<-[1..10],b<-[1..10],c<-[1..10],d<-[1..10],a^b*c^d==1000*a+100*b+10*c+d]
[(2,5,9,2)]
361:１３２人目の素数さん
18/11/01 20:01:51.76 xCdOvDq8.net
>>345
失礼しました
>342のコードが正しい
362:１３２人目の素数さん
18/11/01 20:04:00.31 mQEkML9R.net
>>327
＞(4) nを自然数とし，Ａ^n (i,j)は行列Ａ^n の(i,j)成分を表わすものとする。
＞ そのとき、Ａ^n (2,2) + Ａ^n (3,2) を n, b, c を用いて表わせ。
元の質問者の方向きに解き方の解説
行列のn乗の計算は
A を A’ = P ^-1 A P (A’ は対角行列) と対角化して
A’^n = (P ^-1 A P)^n
⇔ A’^n = P ^-1 A^n P
⇔ P A’^n P^-1 = A^n
ここでA’ は対角行列なので
A’^n は各要素をn乗するだけという流れ
問い (1)～(3) は対角化の仕方を調べているうちにわかると思うので略
363:１３２人目の素数さん
18/11/01 22:26:21.80 tizy9POX.net
(1/x)*ln(1+x)>1+ln(2/(x+2)),
x>0
のときの証明方法を教えて下さい
364:１３２人目の素数さん
18/11/01 22:43:00.52 dkftLkCy.net
とりあえず微分すれば何とかなりそう
365:１３２人目の素数さん
18/11/02 00:07:51.28 +UTP9GLJ.net
>>295
４マス３行（３ターン）と３マス４列（４ターン）で一つの宝と出くわす
確率は同じにならない
■３マス４ターンで少なくとも一つの宝と出くわす確率は
#A＝３^4－２^4＝６５なので
P(A)＝６５／８１
■４マス３ターンで少なくとも一つの宝と出くわす確率は
#B＝４^3－３^3＝３７なので
P(B)＝３７／６４
∴P(A)＞(B)
∵P(A)＝６５／８１＝０．８０２
∵P(B)＝３７／６４＝０．５７８
366:１３２人目の素数さん
18/11/02 00:13:02.61 vfB9uvei.net
今日のNGIDがこんなに早い時間にw
367:１３２人目の素数さん
18/11/02 00:31:24.41 +UTP9GLJ.net
ジョーカー１１枚とハートのエース１枚が入った１２枚の
トランプカードをよくシャッフルする
この山札から１ターン３枚を４回ですべて引くのと
１ターン４枚を３回ですべて引く場合も同じ
368:１３２人目の素数さん
18/11/02 00:41:36.07 G7GSas0t.net
この人確率の問題好きなんだろうね。
しょっちゅう確率の問題に手を出してる。
しかし一度たりとも正解の数値と合ってる式出した事ない。
まぁ本人自分の出した答えが間違ってる事すら理解出�
369:�てないのである意味で幸せなのかもしれない。 苦労して立式して合うはずの答えが何故か合わないあの苦々しさに耐えないで済むんだから。
370:１３２人目の素数さん
18/11/02 01:31:44.33 +UTP9GLJ.net
>>352
１ターン３枚を４－１回で引く時に
ハートのエース１枚が出る確率は
P(A)＝１－（３／４）（２／３）（１／２）＝３／４
１ターン４枚を３－１回で引く時に
ハートのエース１枚が出る確率は
P(B)＝１－（２／３）（１／２）＝２／３
∴P(A)＞(B)
371:１３２人目の素数さん
18/11/02 01:36:19.38 +UTP9GLJ.net
３ｘ４の合計１２マスに宝を一つだけ設置した時に
３列ｘ４ターンと４行ｘ３ターンの探査で同じ確立になるという
計算式をお願いします<(_ _)>
■■■■
■■□■
■■■■
372:１３２人目の素数さん
18/11/02 02:15:59.43 YYpR1gsw.net
>>337
だれかこれをお願いします。
nが自然数なので2次不等式を解いてもあまり上手くいきそうにありません
373:１３２人目の素数さん
18/11/02 02:28:39.71 UxWLcMBZ.net
>>356
x = 2n-1で成立
⇔ 4*n^4－2*n^3－4*n+2 ≧ a
374:１３２人目の素数さん
18/11/02 02:38:59.54 YYpR1gsw.net
>>357
ありがとうございます
1つの例ではなくて必要十分な形で占めしていただけませんか
375:１３２人目の素数さん
18/11/02 03:04:32.60 vfB9uvei.net
それで必要十分条件だろ
x^2-(4n-3)x+a/(n^2+n+1) ≦ 0 ... (1)
左辺は x=2n-3/2 のとき最小、
整数の範囲では x=2n-1 または x=2n-2 で最小値
(a - 4n^4+2n^3+4n^2-2n) / (n^2+n+1) となる
この式の分母は正なので分子が0以下なら(1)を満たす
376:１３２人目の素数さん
18/11/02 03:16:04.56 im1SI6w9.net
>>143 >>194 ～ >>198 >>203 >>282 >>299
立山秀利「入門者のPython」講談社BlueBacks (2018/Sep)
398p．1404円
URLﾘﾝｸ(bookclub.kodansha.co.jp)
【執筆時に使用した環境】
・Microsoft Windows 8.1 および 10
・Python version 3.6
・Anaconda 5.2 for Windows
・Spyder 3.2.8
上記以外の環境でご利用の場合、本書の解説どおりに操作を行えない可能性があります。予めご了承ください。
本書に掲載されている情報は、2018年8月時点のものです。実際にご利用になる際には変更されている場合があります。
【サポートページ】
URLﾘﾝｸ(tatehide.com)
377:１３２人目の素数さん
18/11/02 05:54:31.51 wc7sV/cw.net
>>348
左辺 - 右辺のグラフを書いてみた。
URLﾘﾝｸ(i.imgur.com)
378:１３２人目の素数さん
18/11/02 11:47:02.16 Y7Tkqu2S.net
赤いビックリマーク以後の行がよくわかりません
4^k+1を4×4^kと見なすことで
成り立つと仮定された不等式を援用して新たな不等式を考えているらしいことはわかりますが
どう計算したら24k-5>0になるのかがわかりません
4^k-（8k+1）
URLﾘﾝｸ(i.imgur.com)
379:１３２人目の素数さん
18/11/02 12:07:43.45 ACrozris.net
k≧3
24k≧72
24k-5≧72-5=67>0
380:１３２人目の素数さん
18/11/02 12:49:21.13 cBeA3Am5.net
...やっとわかりました
0より大きい24k-5よりもさらに両辺の差は大きいのでもちろんそれは0より大きく、よって不等号の正しさが証明されていたのですね
二回両辺の差を考えようとしてどうしても24k-5が作り出せず混乱していました
ありがとうございました
381:１３２人目の素数さん
18/11/02 13:33:04.31 b
382:Y2r18eX.net
383:１３２人目の素数さん
18/11/02 14:55:26.18 GY5yQIwK.net
解析的整数論専攻で有名な教授って誰がいますか？雪江明彦氏以外で知ってる方いたら教えてください
384:１３２人目の素数さん
18/11/02 15:01:48.97 MJ+cRGf4.net
微分方程式が解けませんでした。どなたかお願いします。
y'+(x^2)y=1
385:１３２人目の素数さん
18/11/02 15:39:32.17 wc7sV/cw.net
>>367
URLﾘﾝｸ(m.wolframalpha.com)
386:１３２人目の素数さん
18/11/02 15:41:23.13 wc7sV/cw.net
>>368
こっちだな
URLﾘﾝｸ(m.wolframalpha.com)
387:１３２人目の素数さん
18/11/02 16:35:27.47 d4cqGK7t.net
3×4＝12マス、宝1個のみ
□■■■ ■■■■
■■■■ ■■■■ PとQが同時に見つける
■■■■ ■■■□
■■■■ ■■■■ ■■■■ ■■■■ ■■■■
□■■■ ■■■■ ■□■■ ■■■■ ■■■■ Pが先に見つける
■■■■ □■■■ ■■■■ ■□■■ ■■□■
■□■■ ■■□■ ■■■□ ■■■■ ■■■■
■■■■ ■■■■ ■■■■ ■■□■ ■■■□ Qが先に見つける
■■■■ ■■■■ ■■■■ ■■■■ ■■■■
388:１３２人目の素数さん
18/11/02 20:30:16.62 wc7sV/cw.net
ABCD
EFGH
IJKL
のように命名すると
縦:m 横: 宝:k での配置を列挙するコードは既出。数値を変えて実行可。
URLﾘﾝｸ(tpcg.io)
389:１３２人目の素数さん
18/11/02 22:15:38.09 +UTP9GLJ.net
>>371
配置の列挙は確率ではありませんよ
宝が一つの時、縦探査のP君が決して取れない宝は２マス
□□□■
□□□■
□□□□
宝が一つの時、横探査のQ君が決して取れない宝は３マス
□□□□
□□□□
■■■□
決して宝を取れないマスが一マス多いQ君が
P君と同じ確立になるのはなぜ？
390:１３２人目の素数さん
18/11/02 22:20:40.10 +UTP9GLJ.net
P君とQ君が決して取れない宝がある列と行のマス数ね
391:１３２人目の素数さん
18/11/02 22:23:49.11 +UTP9GLJ.net
この状態で計算式を作ると
P(A)＝１－（３／４）（２／３）（１／２）＝３／４
P(B)＝１－（２／３）（１／２）＝２／３
∴P(A)＞(B)
>>354と同じ
392:１３２人目の素数さん
18/11/02 22:52:20.36 G7GSas0t.net
>>372
>>370
393:１３２人目の素数さん
18/11/02 23:55:39.98 CO4fuCl5.net
>>371
P（短軸探索）が先、Q（長軸探索）、同時 の配置を表示するスクリプトを書いてみた。
数値を変えて実行できる。
ｍ：短軸 ｎ：長軸 ｋ：宝の数
URLﾘﾝｸ(tpcg.io)
P1st Q1st even
26 27 13
P (= long axis searcher) finds first.
■ ■ □ ■
□ ■ ■ ■
■ ■ ■ ■
■ ■ ■ □
□ ■ ■ ■
■ ■ ■ ■
■ ■ ■ □
■ ■ ■ ■
□ ■ ■ ■
以下略、
394:１３２人目の素数さん
18/11/03 00:14:01.33 zSMa/Wom.net
>>376
５ｘ６で宝が２個のとき
P1st Q1st even
203 197 35
引き分けになる配置は３５通り、３例ほど挙げるとこんな感じ
とても手作業で列挙する気にはならん。
■ ■ ■ ■ □ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
□ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ □ ■ ■ ■
■ ■ ■ ■ ■ ■
■ □ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ □ ■ ■
■ ■ ■ ■ ■ ■
■ ■ ■ ■ ■ ■
■ □ ■ ■ ■ ■
395:１３２人目の素数さん
18/11/03 00:19:05.37 /E6xXixt.net
>>367
y(x) = u(x)e^(-xxx/3)
を与式に入れると
du/dx = e^(xxx/3),
u(x) = u(0) + ∫[0,x] e^(ttt/3) dt,
y(x) = e^(-xxx/3) {y(0) + ∫[0,x] e^(ttt/3) dt },
かな
396:１３２人目の素数さん
18/11/03 02:01:28.43 GJogDojw.net
y = ∫0→x xy dx
この方程式が解けません
教えて下さい
397:１３２人目の素数さん
18/11/03 02:16:31.16 6u03sBH6.net
y=0
398:１３２人目の素数さん
18/11/03 02:26:53.34 GJogDojw.net
y = ce^(x^2/2)
であってますか？
399:１３２人目の素数さん
18/11/03 02:34:04.86 Sa4Jrve0.net
微分する
初期条件y(0)=0
400:１３２人目の素数さん
18/11/03 13:29:52.32 PHFdSdeY.net
確率1/3のくじを１回ひくのと確率1/9のくじを３回ひくのでは、
当たりをひく確率は同じですか？
複数回ひく場合でも前にひいたくじがなくなる訳ではなく
毎回同じ確率で抽選されるという仮定の場合です
401:１３２人目の素数さん
18/11/03 14:05:15.33 Ng1V9LvS.net
1/9^3で3回当たり
24/9^3で2回当たり
192/9^3で1回当たり
8^3/9^3で全て外れ
「一回でも当たる確率」は、(1+24+192)/9^3 = 217/9^3 = 1-(8/9)^3 < 243/9^3 = 1/3
なので、当選確率1/3のクジを一回引くのより小さい
しかし、「当たる回数の期待値」は
(3*1+2*24+1*192)/9^3=(3+48+192)/9^3=243/9^3=1/3
なので、当選確率1/3のクジを一回引くのと同じ
402:１３２人目の素数さん
18/11/03 15:37:38.65 zSMa/Wom.net
>>372
離散量の確率は場合の数をいかに効率的にカウントするかによるね。
手作業だと漏れがでるからプログラムの利用は必須
>377参照。
投稿前に自分でシミュレーション検証して投稿すれば、
>302のように こいつ 呼ばわりされなくて済むんだけどね。
自分で算出した値が別の言語の算出結果と一致したと投稿されるとシミュレーションの正しさが確認できていいね。
俺が鈍足のRコードのをだすと高速のｃが投稿されたり、解析解が投稿されて数理とプログラム論理の勉強になって嬉しいね。
403:１３２人目の素数さん
18/11/03 15:45:00.36 tJ6POH4O.net
>>302
sed使えよｗ
404:１３２人目の素数さん
18/11/03 15:49:27.19 zSMa/Wom.net
>>383
顰蹙のシミュレーション検証
１００万回シミュレーションして頻度をだしてみた
確率1/3のくじを１回ひく
> mean(replicate(k,sample(x,1))
[1] 0.333435
確率1/9のくじを３回ひくのでは、)
> mean(replicate(k,sum(sample(y,3))))
[1] 0.333176
405:１３２人目の素数さん
18/11/03 15:54:56.85 Sa4Jrve0.net
１見学者からのお願いだけど、NGリスト入りしてそうなネタに関連する話と
普通の話は出来れば
406:１３２人目の素数さん
18/11/03 15:55:58.31 Sa4Jrve0.net
>>388続き
普通の話とは出来ればわけておいてほしいな、と思う
407:１３２人目の素数さん
18/11/03 18:19:29.22 Ha92ty6K.net
ここの国では硬貨は7種類流通しています
この7種類の硬貨を使って1円～70円の70通りの支払いができます
ただし一度に使用できる硬貨は3枚以下(同じ硬貨2度使いは可)です
7種類の硬貨はそれぞれ何円だったのでしょうか？
408:１３２人目の素数さん
18/11/03 19:16:17.42 Ha92ty6K.net
とりあえず分かったこと
最低額は1円、最高額は24円以上
409:１３２人目の素数さん
18/11/03 19:28:25.33 dzhfVDG/.net
「2度使いは可」とは３度はダメという意味か？
410:１３２人目の素数さん
18/11/03 19:38:13.92 Ha92ty6K.net
たぶん3度使いも可
411:だと思う 不可だと無理なんじゃね？ 知らんけど
412:１３２人目の素数さん
18/11/03 19:40:09.02 fqIvCrY7.net
とりあえず解を一組でも見つけりゃいいんだよね？
413:１３２人目の素数さん
18/11/03 20:09:37.68 nqMGpkef.net
>>385
３ｘ４の１２マスで宝が一つだけの時、
P君とQ君は互いに最終列と最終行の宝は
取ることができない
□□□■
□□□■
□□□□
□□□□
□□□□
■■■□
つまり、P君の探査範囲は縦３マスｘ３列
Q君の探査範囲は横４マスｘ２行になる
それぞれの探査範囲内でP君とQ君が
少なくとも一つの宝を見つけるという
事象Aと事象Bを考える
P(A)／P(B)＝（P君の勝つ数）／（Q君の勝つ数）
∴P(A)／P(B)＝｛１－｛ｎ／（ｎ＋１）｝^ｎ｝／｛１－｛（ｎ－１）／ｎ｝^（ｎ－１）｝
ｎ＝３のとき、P(A)／P(B)＝３３３／３２０
互いの最終列と最終行にある宝の取れないマスが一つ多い
Q君よりもP君のほうが僅かに確率が上がることが
如実に示される
■Wolfram入力例
｛１－｛ｎ／（ｎ＋１）｝^ｎ｝／｛１－｛（ｎ－１）／ｎ｝^（ｎ－１）｝，ｎ＝３
｛1-｛n/(n+1)｝^n｝/｛1-｛(n-1)/n｝^(n-1)｝,n=3
414:１３２人目の素数さん
18/11/03 20:35:17.71 GvDy13c1.net
たった66通りしかないのによくずっと正解を避け続けられるものだな…
415:１３２人目の素数さん
18/11/03 20:40:39.80 Ha92ty6K.net
>>394
もちろんいいけど、解が2組以上もあるとは思えないというのが最初の直感
こちとらヒントを見ても、さっぱり分からん
416:１３２人目の素数さん
18/11/03 21:14:48.47 OVkXWZOI.net
>>390
1, 4, 5, 15, 18, 27, 34 かな？
417:１３２人目の素数さん
18/11/03 21:34:57.47 vHgKtUFX.net
>>398
素晴らしい❗
Prelude Data.List> filter (<=70) $ let x = [0,1,4,5,15,18,27,34] in nub$sort$[a+b+c | a<-x,b<-x,c<-x]
[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70]
Prelude Data.List> length $ filter (<=70) $ let x = [0,1,4,5,15,18,27,34] in nub$sort$[a+b+c | a<-x,b<-x,c<-x]
71
418:１３２人目の素数さん
18/11/03 21:47:25.87 Ha92ty6K.net
やっと(頭の中での)検算が終わったw
(1,４,５)は全くの盲点だったわ
419:１３２人目の素数さん
18/11/03 22:19:10.88 ZvavWJQe.net
>>395は数学的な意味の確率を理解したいわけでも、計算できるようになりたいわけでもないんだろう。
まぁ確率なんか計算できなくても社会生活困るわけでもないし。
なんとなく割り算して、それで楽しいならそれでいいのかなと思う。
420:１３２人目の素数さん
18/11/03 23:19:31.84 Ha92ty6K.net
トランプの山からA君とB君が交互に1枚ずつ引いて
先にジョーカーを引いたほうが勝ちとする
(引いたカードは山に戻さない)
①トランプ52枚 + ジョーカー1枚
②トランプ52枚 + ジョーカー2枚
先攻勝率は①②で同じ 27/53
後攻勝率は①②で同じ 26/53
不思議だと思わない？(計算めんどい)
421:１３２人目の素数さん
18/11/03 23:37:36.80 nqMGpkef.net
>>350
これはクリプテックスの確率だった
>>374よりも精度を上げることができた
422:１３２人目の素数さん
18/11/03 23:59:56.38 AlD/TWrH.net
pを素数とする。
-p(a+b)+p^2ab+1=0
pa+b=pb
を満たす整数の組(a,b)が存在するかどうか判定せよ。
423:１３２人目の素数さん
18/11/04 00:44:45.11 1wP06nNi.net
(pa-1)(pb-1)≠0
424:１３２人目の素数さん
18/11/04 02:50:46.13 Bn7LN70u.net
y = ∫0→x ∫0→x y dx
y(0)=0
この方程式はどのように解けばいいですか？
答えもお願いします
425:１３２人目の素数さん
18/11/04 03:02:00.16 Bn7LN70u.net
間違え
426:ました こうです y = ∫0→x ∫0→x y dxdx y(0)=0
427:１３２人目の素数さん
18/11/04 03:23:06.58 Bn7LN70u.net
y = ∫0→x ∫0→x y dxdx
y(0)=0
dy(0)/dx=0
初期条件抜けてました
428:１３２人目の素数さん
18/11/04 04:21:30.68 tTiGqsss.net
>>390 >>398
1枚だけで可能 … 1,4,5,15,18,27,34
2枚で可能 … 2,6,8～10,16,19,20,22,23,28,30～33,35,36,38,39,42,45,49,52,54,61,68
残り … 3,7,11～14,17,21,24～26,29,37,40,41,43,44,46～48,50,51,53,55～60,62～67,69,70
429:１３２人目の素数さん
18/11/04 04:43:57.43 tTiGqsss.net
>>407
与式 → y(0) = 0,
xで微分すると
y '(x) = ∫[0,x] y(t) dt → y '(0) = 0,
もう一度xで微分すると
y "(x) = y(x),
これより
y(x) = a･e^x + b･e^(-x)
y(0) = a+b = 0,
y '(0) = a-b = 0,
∴ a=b=0
y(x) = 0
>>408
初期条件を付記する必要はありません。(式から出ます。)
430:１３２人目の素数さん
18/11/04 05:40:42.34 tTiGqsss.net
>>367 >>378
|x| が小さいところでは、マクローリン展開より
y(x) = e^(-xxx/3){y(0) + ∫[0,x] e^(ttt/3) dt}
= e^(-xxx/3) y(0) + x + Σ[k=1,∞] (-1)^k /{4･7････(3k+1)} x^(3k+1)
431:１３２人目の素数さん
18/11/04 06:39:01.88 ol3G48+j.net
>>398
これはどうやって出されたのですか？
組み合わせは1,198,774,720通りなので総当たりは断念しました。
432:１３２人目の素数さん
18/11/04 08:13:32.59 VDxltAIF.net
>>399
表示可能な数字は
[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,72,73,76,79,81,83,86,88,95,102]
で71は無理だな。
433:１３２人目の素数さん
18/11/04 08:16:50.86 zevnpesP.net
>>412
私は 398投稿者ではありませんが、次のプログラムで答が出せます。
単純に、C[70,7]通りとはせず、ある程度変数の範囲を押さえ込むことがこつかも。
void next(int depth,int *x){
if(depth<7){
int i,is=x[depth-1]+1,ib=x[depth-1]*3+1,ie=ib<70?ib:70;
for(i=is;i<=ie;i++){x[depth]=i;next(depth+1,x);}
}else{
int i,j,k,mem[211],flag;
for(i=1;i<211;i++)mem[i]=0;
for(i=0;i<7;i++){mem[x[i]]++;mem[2*x[i]]++;mem[3*x[i]]++;}
for(i=0;i<6;i++)for(j=i+1;j<7;j++){mem[x[i]+x[j]]++;mem[2*x[i]+x[j]]++;mem[x[i]+2*x[j]]++;}
for(i=0;i<5;i++)for(j=i+1;j<6;j++)for(k=j+1;k<7;k++){mem[x[i]+x[j]+x[k]]++;}
for(i=1,flag=1;i<71;i++)if(mem[i]==0)flag=0;
if(flag==1)printf("%d,%d,%d,%d,%d,%d,%d\n",x[0],x[1],x[2],x[3],x[4],x[5],x[6]);
}
}
int main(){int x[7];x[0]=1;next(1,x);return 0;}
URLﾘﾝｸ(codepad.org)
434:１３２人目の素数さん
18/11/04 08:36:53.93 VDxltAIF.net
>>414
いつものｃでの高速解をありがとうございます。
他の組合せがないことがわかって助かりました。
435:１３２人目の素数さん
18/11/04 10:06:08.50 Bn7LN70u.net
>>410
ありがとうございました
436:１３２人目の素数さん
18/11/04 11:37:15.42 Bn7LN70u.net
y "(x) = y(x)^(-2)
何度もすみません
これはどんなふうに解けばいいですか？
437:１３２人目の素数さん
18/11/04 14:03:43.89 Bl3oi46w.net
y"-1/y^2=0 から 0=2y'(y"-1/y^2)=2y'y"-2y'/y^2=(y'^2+2/y)' として
エネルギー積分 y'^2+2/y=C を求め、変数分離形の公式通りに
∫
438:dy/√(C-2/y)=∫ dt=t とする
439:１３２人目の素数さん
18/11/04 15:46:42.35 76EMArHB.net
1 3 7 13 22 34
を一つの数式で表すとどうなりますか？
440:１３２人目の素数さん
18/11/04 16:30:22.70 RDBCTf5Y.net
(n-1)^2 +n
441:１３２人目の素数さん
18/11/04 16:48:26.93 xWdCQ8DY.net
α,β,γ,δはαδ-βγ=1, |α|≦|β|≦|γ|≦|δ|を満たす複素数である。
（1）この4つの複素数全てについて、その実部と虚部が共に整数となる例を1組挙げよ。
（2）（1）において、4つの複素数のいずれも0でないことはあるか。
442:１３２人目の素数さん
18/11/04 17:06:09.24 puykSUeo.net
1,1,1,2
443:１３２人目の素数さん
18/11/04 17:13:46.36 xWdCQ8DY.net
>>421
（3）少なくとも1つが実数でなく、その絶対値が1でないようなa,b,c,dはあるか。
444:１３２人目の素数さん
18/11/04 17:20:21.20 puykSUeo.net
1+I,1+2i,1001(1-2i),2503(1+i)
445:１３２人目の素数さん
18/11/04 17:22:03.54 RDBCTf5Y.net
>>420
間違えた
446:１３２人目の素数さん
18/11/04 17:32:54.95 76EMArHB.net
5移行で一致しない
1 3 7 13 22 34 50 70 95 125 161 203
を一つの数式で表すとどうなりますか？
447:１３２人目の素数さん
18/11/04 17:45:36.46 76EMArHB.net
5ごとに1 3 7 13をループさせるようです
448:１３２人目の素数さん
18/11/04 17:59:07.56 tTiGqsss.net
>>419 >>426
a_n = Σ(k=1，n+1) [kk/4]
生成関数は x/{(1+x)(1-x)^4}
1, 3, 7, 13, 22, 34, 50, 70, 95, 125, 161, 203, 252, 308, 372, 444, 525, 615, 715, …
URLﾘﾝｸ(oeis.org)
449:１３２人目の素数さん
18/11/04 18:00:55.69 3UkFCnqw.net
数列に関してですが
Sn=1+3･2+5･2^2・・・+(2n-1)･2^n-1の問題で
Snを掛ける２してSn-2Snで引くまではわかったんですが、
Sn-2Snで引いた-S=1+2･2+2･2^2+・・・ 2･2^2n-1-(2n-1)･2^nから
2+2^2+2^3・・・+2^n-(2n-1)･2^-1←この一番最後にある「-1」がなぜ付いてるのかがわからないんですがなぜこうなってるんでしょうか？
お願いします
450:１３２人目の素数さん
18/11/04 18:07:35.99 76EMArHB.net
>>428
おお、これは助かりまする
1 3 7 13 22 34 50 70 95 125 161 203
252 308 372 444 525 615
事前準備と見事に一致
451:１３２人目の素数さん
18/11/04 18:12:56.01 tTiGqsss.net
>>419 >>426
2つの数式で書けば
a_{2m} = m(m+1)(4m+5)/6,
a_{2m+1} = (m+1)(m+2)(4m+3)/6
蛇足だけど。
URLﾘﾝｸ(oeis.org)
452:１３２人目の素数さん
18/11/04 18:23:47.81 tTiGqsss.net
>>419 >>426
1つの数式で書けば
a_n = Σ(k=1，n+1) [kk/4] = [ (n+1)(n+3)(2n+1)/24 ],
URLﾘﾝｸ(oeis.org)
453:１３２人目の素数さん
18/11/04 18:31:04.97 SFrSfyf1.net
>>429
括弧の位置おかしくね？
454:１３２人目の素数さん
18/11/04 18:51:07.55 tTiGqsss.net
>>429
r≠1 とする。
S_n = Σ[k=1，n] (2k-1)･r^(k-1),
r･S_n = Σ[k=1，n] (2k-1)･r^k = 1 + Σ[k=1，n+1] (2k-3)･r^(k-1),
辺々引くと
S_n - r･S_n = -1 + 2Σ[k=1,n] r^(k-1) - (2n-1)･r^n
= -1 + 2(r^n - 1)/(r-1) - (2n-1)･r^n,
以下略
455:１３２人目の素数さん
18/11/04 18:53:43.36 3UkFCnqw.net
>>433
2+2^2+2^3・・・+2^n-(2n-1)･2^n-1←でしたすいません
456:１３２人目の素数さん
18/11/04 20:03:23.29 Bn7LN70u.net
>>417なんですが
両辺を二乗して
(y "(x))^2 - 1 = 0とし
y "(x) = ±1
を特性方程式を作って後は解くだけでよいでしょうか？
457:１３２人目の素数さん
18/11/04 20:04:40.09 Bn7LN70u.net
間違えました
458:１３２人目の素数さん
18/11/04 21:01:11.73 lEgi/Aj9.net
さすがに微分方程式の本で簡単なものを買った方がいいように思える
459:１３２人目の素数さん
18/11/04 21:01:22.36 RDBCTf5Y.net
1たす1は2
101たす102は3
8たす1は1
では2たす2は？
460:１３２人目の素数さん
18/11/04 22:02:17.92 Bn7LN70u.net
>>438
いろんなサイト見てるんですけど>>417のような問題が無くて解けないんです
よかったらヒント
461:下さい
462:１３２人目の素数さん
18/11/04 22:11:15.82 LooLpWav.net
今紙ないから確認できないけど>>417は両辺y’かけて積分したらいけそうな気がする。
463:１３２人目の素数さん
18/11/04 22:20:38.55 76EMArHB.net
>>98
完全追尾型多項式が完成しました
宝の個数を２で固定します
P1st =｛12n^4+28n^3-42n^2-52n-3(-1)^n+51｝/48
Q1st =｛12n^4+20n^3-18n^2-20n-3(-1)^n+3｝/48
■Wolframに入力すると既約分数表示になるので御注意
P1st／Q1st
=｛8(n-1)｛(n-2)n-6｝/｛2n(n+2)(6n^2-2n-5)-3(-1)^n+3｝｝+1
464:１３２人目の素数さん
18/11/04 22:46:11.36 JrDn1ZDl.net
>>442
>>161に正しい答えがあるよ
わかりにくいなら最後の辺りだけ見て
P1st は nが奇数の時P1偶数のときP2
Qも同様
465:１３２人目の素数さん
18/11/04 22:55:51.47 76EMArHB.net
>>161は無能
実際に多項式になっていない
466:１３２人目の素数さん
18/11/04 23:08:20.12 lTCeMsqQ.net
>>414
C[70,7]通りのRのスクリプトを書いてみた。
正確にはC[70,8]*C[8,3]=528659651520通リｗｗｗ
他の言語に移植する人いるかなぁ？
is.1_70 <- function(x){
total=NULL
for(i in x){
for(j in x){
for(k in x){
ijk=i+j+k
if(!(ijk %in% total)) total=append(total,ijk)
}
}
}
all(1:70 %in% total)
}
M=69
for(a in 0:M){
for(b in a:M){
for(c in b:M){
for(d in c:M){
for(e in d:M){
for(f in e:M){
for(g in f:M){
for(h in g:M){
y=c(a,b,c,d,e,f,g,h)
if(is.1_70(y)) print(y)
}
}
}
}
}
}
}
}
467:１３２人目の素数さん
18/11/05 00:24:00.08 QvXJrUC9.net
>>445
Haskellに移植。
とりあえずコンパイルエラーは出なかった。
朝までに計算が終わるかどうかは不明。
import Data.List
m = 69
sub x = do
let ijk = filter (<=70).nub $ sort [i+j+k| i<-x,j<-x,k<-x]
all (\y -> elem y ijk ) [0..70]
main = do
print $ [(b,c,d,e,f,g,h)| b<-[0..m],c<-[b..m],d<-[c..m],e<-[d..m],f<-[e..m],g<-[f..m],h<-[g..m],sub [0,b,c,d,e,f,g]]
468:１３２人目の素数さん
18/11/05 00:51:40.42 Un0fMQvD.net
やや速度改善。
import Data.List
firstUnavailable x = let y = 0:x in head $([1..71] ¥¥)$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
next x = [n:x|n<-[head x+1..firstUnavailable x]]
xss = iterate (¥xs->concat [next x|x<-xs]) [[1]]
isGood x = let y = 0:x in (==70)$length $intersect [1..70]$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
main = do
print [x|x<-(xss !! 6),isGood x]
469:１３２人目の素数さん
18/11/05 00:55:04.97 fIzIE6qz.net
>>446
-- 最終行にhが抜けてたので修正。
import Data.List
m = 69
sub x = do
let ijk = filter (<=70).nub $ sort [i+j+k| i<-x,j<-x,k<-x]
all (\y -> elem y ijk ) [0..70]
main = do
print $ [(b,c,d,e,f,g,h)| b<-[0..m],c<-[b..m],d<-[c..m],e<-[d..m],f<-[e..m],g<-[f..m],h<-[g..m],sub [0,b,c,d,e,f,g,h]
470:１３２人目の素数さん
18/11/05 01:07:42.61 fIzIE6qz.net
>>447
いつもありがとうございます。
お見事に算出されました。
*Main Data.List> :main
[[34,27,18,15,5,4,1]]
471:１３２人目の素数さん
18/11/05 01:09:16.43 Un0fMQvD.net
これ以上かくと多分うざいのでラスト。やや速度改善。
import Data.List
firstUnavailable x = let y = 0:x in head $([1..] ¥¥)$nub$sort$[a+b+c|a<-y,b<-y,c<-y]
next x = [n:x|n<-[head x+1..firstUnavailable x]]
xss = iterate (¥xs->concat [next x|x<-xs]) [[1]]
isGood x = (>70) $ firstUnavailable x
main = print [x|x<-(xss !! 6),isGood x]
自宅のパソコンだとghc -O2 で22秒でおわったけどcodepadだとTimeoutした。
472:１３２人目の素数さん
18/11/05 01:10:35.80 KyuWjb44.net
1000枚の1円玉の中に1枚だけ両面とも表の1円玉がある。
この中から1枚だけ選んで10回投げたところ、10回連続で表が出た。
このとき、この選ばれた1円玉が両面表である確率は
普通の1円玉である確率より高い？低い？
473:１３２人目の素数さん
18/11/05 01:10:44.03 Un0fMQvD.net
>>449
ghciでやったんだ。流石にその勇気はなかったｗｗ
474:１３２人目の素数さん
18/11/05 01:16:37.48 /nWeWNpo.net
最大価値のコインの価値が1枚35円以上だと、残りの
6種のコイン2枚以下で34円まで表せないといけない
逆に6種を2枚以下で34円まで表せるなら1枚35円の
コインを追加した7種が3枚までで70円まで表せる。
よってまず6種2枚までで1～34円が全て表せるかを調べて、
それが無理ならコインの価値は最大34円までと限定できる
この先も上から攻めていけば多少探索範囲を限定できると思うが定かではない
475:１３２人目の素数さん
18/11/05 01:28:56.96 /nWeWNpo.net
>>453
＞最大価値のコインの価値が1枚35円以上だと、残りの
＞6種のコイン2枚以下で34円まで表せないといけない
この部分はアプリオリではないか
正しいような気はするが少なくとも数行では証明できなさそう
476:１３２人目の素数さん
18/11/05 02:04:50.67 Un0fMQvD.net
>>451
P(本物|10連続表)
=P(本物＆10連続表)/P(10連続表)
=999/1000・1/1024/P(10連続表)
=999/1024/P(10連続表)/1000
P(偽物|10連続表)
=P(偽物＆10連続表)/P(10連続表)
=1/1000・1/1/P(10連続表)
=1/P(10連続表)/1000
∴ P(本物|10連続表)＜P(偽物|10連続表)
477:１３２人目の素数さん
18/11/05 02:32:14.74 OlP2HpBB.net
>>441
y"(x) = y(x)^(-2)
以下より(x)を省略
(y'^2)' = 2y'y"
(y^(-1))' = -y(x)^(-2) * y'
ここで
y" = y^(-2)の両辺にy'をかけて
y"y' = y(x)^(-2) * y' となり
1/2 * (y'^2)' = - (y^(-1))'
(y'^2)' = (-2y^(-1))'
故に y'^2 = -2y^(-1)
ここまではあってますか？
ここから先が解けるかどうかわかりませんがもう少し考えてみます
478:１３２人目の素数さん
18/11/05 02:39:55.15 KyuWjb44.net
>>455
なるほど、スッキリした
① P(本物＆10連続表)＝(999/1000)*(1/1024)
② P(偽物＆10連続表)＝(1/1000)*(1)
①：②＝999：1024
P(本物|10連続表)＝①/(①+②)＝999/2023
P(偽物|10連続表)＝②/(①+②)＝1024/2023
479:１３２人目の素数さん
18/11/05 03:05:38.40 KyuWjb44.net
1024枚以下なら結論は変わらずで
1025枚以上から結論が逆になるんだな
480:１３２人目の素数さん
18/11/05 03:27:38.95 /nWeWNpo.net
1000人に1人の予言者を探すなら10回では足らず20回は当て物させて確かめないとだめってことだな
481:１３２人目の素数さん
18/11/05 03:54:13.72 B1F8UTQM.net
>>417 >>456
故に (y ')^2 = c - 2/y,
ですね。cは積分定数です。
そこから先は xをyの関数と見て
x = ∫(1/y ') dy
= ∫√{y/(cy-2)} dy
= (1/c)√{y(cy-2)} + {2/c^(3/2)}･log[c√y + √{c(cy-2)} ] + c'
のような式になり、逆関数を求めるのは難しい。
・用途
クーロン散乱・ラザフォード散乱で正面衝突する場合(θ=0)とかに使えるかなぁ。
482:１３２人目の素数さん
18/11/05 04:19:20.08 /nWeWNpo.net
辺々 y' かけて wolfram alpha に
y'' * y' = 1/ y' とか y'' * y'^2 = 1
と入力したら答えでるね
それによると
y'(x) = v(x) とおけば
v' * v^2 = 1
積分して
v^3 / 3 = x + c
以下略
483:１３２人目の素数さん
18/11/05 07:01:50.73 GGbUoN9W.net
>>459
コインを1000回投げて10回以上連続して表がでる確率は？
484:１３２人目の素数さん
18/11/05 07:57:00.13 fIzIE6qz.net
>>462
1000回投げて10回以上連続して表
[1] 0.3854498
1000回投げてちょうど10回連続して表
[1] 0.1700181
1000回投げて20回以上連続して表
[1] 0.000468149
1000回投げてちょうど20回連続して表
[1] 0.0002342865
URLﾘﾝｸ(tpcg.io)
これを参考にプログラム
URLﾘﾝｸ(detail.chiebukuro.yahoo.co.jp)
485:/question_detail/q1149349046
486:１３２人目の素数さん
18/11/05 08:00:43.20 fIzIE6qz.net
コインをN回投げてK回以上連続して表がでる確率を多項式で表現できるのかどうかは知らないので悪しからず。
確率誤答の達人が全角文字で組んでくれるかもｗｗｗ
487:１３２人目の素数さん
18/11/05 15:38:02.87 wfCkOOVj.net
>>462
A[n]：○と×併せてn個を一列に並べた文字列で、○の１０連以上が無いものの数。
B[n]：○と×併せてn個を一列に並べた文字列で、○の１０連以上を含むものの数。A[n]+B[n]=2^n
A[n;k]：A[n]の中で、最後に○がk個連続しているもの。A[n]=A[n;0]+...+A[n;9]
P[n]：A[n;k](k=0～9)とB[n]を並べた、11成分の縦ベクトル
P[1]={1,1,0,0,0,0,0,0,0,0,0}^t
X={{1,1,1,1,1,1,1,1,1,1,0},
{1,0,0,0,0,0,0,0,0,0,0},...,
{0,0,0,0,0,0,0,0,1,0,0},
{0,0,0,0,0,0,0,0,0,1,2}}
P[1000]=X^(999) P[1] の第11成分を2^1000で割ったものが、求める確率
41301272734778977984946818232089531229879543376756574850136155867680807079676964
05909423852137579237591446526939613263507523948827986893531646240157193872907615
64116695521478307244714549348159061083607249922721310512099499789154886902065157
8128373092635280064104398841562373328900104830268510093352961/2^1000
≒0.38544975241248163591...
488:１３２人目の素数さん
18/11/05 16:24:17.34 B1F8UTQM.net
>>456 は c=0 の場合であり、
y = - (9/2)^(1/3)･|x-c '|^(2/3)
と解ける。
489:１３２人目の素数さん
18/11/05 17:35:37.11 KyuWjb44.net
トランプのA～10の10枚とジョーカー1枚の
合計11枚が机の上に裏向きに置いてある。
ランダムに1枚ずつ引いていった場合の、得られた数字の総和の期待値を求めよ。
ただし、ジョーカーを引いた時点で終了するものとし、
Aは数字扱いではなく、最終的に得られた数字の総和が２倍になるものとする。
490:１３２人目の素数さん
18/11/05 17:41:14.94 Qohbqnrn.net
方程式y'=xyln(y')を解け
491:１３２人目の素数さん
18/11/05 17:58:10.16 Qohbqnrn.net
外心をOとする△ABCの頂点の内部に点Kをとり、Oに関してKと対称な点をLとする。
△ABK=△BCLとなる点Kの位置を求めよ。
492:１３２人目の素数さん
18/11/05 18:19:45.10 NWPSgxHY.net
>>467
シミュレーションしたら４５くらいになった。
> summary(re)
Min. 1st Qu. Median Mean 3rd Qu. Max.
41.80 44.29 45.01 45.03 45.81 48.46
493:１３２人目の素数さん
18/11/05 18:28:56.21 NWPSgxHY.net
>>467
１００万回のシミュレーションで総和の分布はこんな感じになった。
URLﾘﾝｸ(i.imgur.com)
494:１３２人目の素数さん
18/11/05 19:09:42.88 Pcec+Aw3.net
アレ？
33になる？
2~10のうちジョーカーより左にくるものの期待値は44/2=22。
確率1/2で2倍になるから33。
シミュ合ってる？
495:１３２人目の素数さん
18/11/05 19:13:15.66 Pcec+Aw3.net
あ、違う。
54/2の3/2倍で41だ。
496:１３２人目の素数さん
18/11/05 19:14:05.08 Pcec+Aw3.net
40.5 orz
497:１３２人目の素数さん
18/11/05 19:17:38.57 UErFb0f/.net
>>472-474
ジョーカーが早めに出る場合はAで2倍になる確率が低い
498:１３２人目の素数さん
18/11/05 19:46:33.82 1RAsBANL.net
Aがあったら全部２倍するんだろ？
Aが出るまでの値だけを２倍？
499:１３２人目の素数さん
18/11/05 19:52:13.15 Qohbqnrn.net
袋の中に大量の赤球と青球があり、どちらを取り出す確率も常に1/2である。
袋の中から球を1個取り出すことをn回繰り返すとき、同じ色がk回連続して取り出され、かつk+1回以上は連続して取り出されない確率をP(n,k)とおく。
以下の問に答えよ。
（1）P(n,k)/P(n+1,k+1)をnとkの式で表せ。
（2）lim[n→∞] Σ[k=1 to n] kP(n,k) を求めよ。
500:１３２人目の素数さん
18/11/05 19:53:16.37 jroNIL0U.net
なぜこの曲線C'の0≦θ≦tでの長さがこうやって求められるのか可能な限り優しく教えてください。
自分は媒介変数表示のまま、√（x'^2 +y'^2）を積分して回答しました。
URLﾘﾝｸ(i.imgur.com)
501:１３２人目の素数さん
18/11/05 19:59:05.80 KyuWjb44.net
>>476
そう、Aを引いたら全部2倍
502:１３２人目の素数さん
18/11/05 19:59:09.74 1RAsBANL.net
>>472
f <- function(x){
i=1
y=numeric()
while(x[i]!=11){ #11:ジョーカーでなければ
y[i]=x[i] # yに保存
i=i+1
}
if(1 %in% y) return(2*(sum(y)-1)) # 1があれば総和から1引いて2倍
else return(sum(y))
}
# simulation ,sample関数で1から11をランダムに並べ替え変え
re=replicate(10^6,f(sample(11))) #100万回fを繰り返す
summary(re)
hist(re,col='lightblue',xlab='sum',main='')
503:１３２人目の素数さん
18/11/05 20:03:58.80 OlP2HpBB.net
>>460>>466
詳しい解説ありがとうございます
微分方程式のサイトを調べてる感じだと後は何とかなりそうな気がします
504:１３２人目の素数さん
18/11/05 20:05:46.39 Xov6oqbO.net
>>479
全部ってのがJoker引くまでの全部なのかAをひくまでの全部なのかを聞かれてるんだよ。
23A47J
の場合
(2+3)×2 + 4 + 7
なのか
(2+3+4+7)×2
なのかどっちにも取れるんだよ。
自分が書いてる文章でホントにちゃんと意味が伝わるのか考えながら書かないと。
505:１３２人目の素数さん
18/11/05 20:08:49.31 UErFb0f/.net
>>482
もとの問題文に明確に書かれているし、
>>476の質問者も単に「全部」と「Aが出るまでの値だけ」の2択で質問しているから紛らわしさはない
506:１３２人目の素数さん
18/11/05 20:13:38.56 zOJn+Nri.net
>>478
図形的な解釈はハイレベル理系数学などに出てるから本屋で見ろ
数式で示すなら
x = rcosθ
y = rsinθ
にパラメータ表示の公式を適用して整理すればいい（ r はθの関数であることに注意）
507:１３２人目の素数さん
18/11/05 20:38:17.38 NWPSgxHY.net
>>482
23A47Jの場合、A=1 J=11
>480のfで
> f(c(2,3,1,4,7,11))
[1] 32
でいいんだよね？
それでよければシミュレーションでいいと思うんだが。
508:１３２人目の素数さん
18/11/05 20:40:54.60 KyuWjb44.net
>>483
言いたかったことを代弁してくれて、ありがとう
509:１３２人目の素数さん
18/11/05 20:53:56.37 jroNIL0U.net
>>484
ありがとうございます
媒介変数表示の時の図形的意味はなんとなく分かっているので
これで一応理解できたと思います
∫√(r^2+r'^2)から直接図形的意味を説明できませんでしょうか？
とても知りたいです
510:１３２人目の素数さん
18/11/05 20:54:36.84 jroNIL0U.net
媒介変数表示の形に戻せば図形的意味は理解できますが、
そういう変形を使わず余弦定理とかを使って証明する方法があれば知りたいということです
511:１３２人目の素数さん
18/11/05 21:02:52.61 Pcec+Aw3.net
>>467
> トランプのA～10の10枚とジョーカー1枚の
> 合計11枚が机の上に裏向きに置いてある。
> ランダムに1枚ずつ引いていった場合の、得られた数字の総和の期待値を求めよ。
> ただし、ジョーカーを引いた時点で終了するものとし、
> Aは数字扱いではなく、最終的に得られた数字の総和が２倍になるものとする。
この文章のどこにエースを引くまでの全部と読み取れる要素があるねん？
むしろ最終的に得られた数の二倍ってJ引くまでのトータルを二倍としか読めない気がするけど。
コレがエースを引くまでのみを二倍の意味ならどう考えても言葉の選択間違ってるやろ？
512:１３２人目の素数さん
18/11/05 21:09:40.27 NWPSgxHY.net
>>485
１１個の順列を列挙して>480のｆを適用して平均値をだそうと思ったのだが、
11! = 39916800 なので PC処理が終わらないので断念。
513:１３２人目の素数さん
18/11/05 22:28:30.12 hlCe+j6H.net
>>442
なかなかやるじゃないですか(´・ω・`)
514:１３２人目の素数さん
18/11/05 22:40:31.92 /AGPvHl/.net
f(x) = x^3に対して f'(x) = 3x^2 なのでf(x)は1階連続微分可能です。
fの逆関数をgとすると
g(x) = x^(1/3)
となって、導関数は
g'(x) = (1/3) x^(-2/3)
となります。するとx=0でg'(x)は無限大に発散して連続でないように思えます。
g(x)は1階連続微分可能ではないということでよいのでしょうか？
もし1階連続微分可能でないとすると、f(x)とg(x)のグラフは回転・反転
させただけで滑らかさは全く変わらないことと不整合なように思えるのですが
どのように考えればよいのでしょうか？
515:１３２人目の素数さん
18/11/05 23:07:22.92 qQ0oc4f2.net
いつからgとg'の定義域が一致すると錯覚していた？
516:１３２人目の素数さん
18/11/05 23:20:02.96 jroNIL0U.net
あーそうか三乗根を自乗するからg'は非負なのか
517:１３２人目の素数さん
18/11/05 23:22:56.68 oUnORFjS.net
>>487
∫√(r^2+(dr/dθ)^2) dθ で、dθが正になるように積分範囲を決定すると
= ∫√((rdθ)^2 + (dr)^2)
だから、三平方の定理で、直角を挟む辺の長さが(rdθ)、(dr)の直角三角形を
考えると斜辺の長さは √((rdθ)^2 + (dr)^2)
この斜辺の長さを足しあわせたものが曲線の長さになる。
非常に物理的な大雑把で直観的な説明。
ハイレベル理系数学持ってないけど、多分似たような説明だと思うので絵だけでも本屋で見てください。
518:１３２人目の素数さん
18/11/06 00:09:04.89 Jhril6/D.net
袋の中に大量の赤球と青球があり、どちらを取り出す確率も常に1/2である。
袋の中から球を1個取り出すことをn回繰り返すとき、同じ色がk回連続して取り出され、かつk+1回以上は連続して取り出されない確率をP(n,k)とおく。
以下の問に答えよ。
（1）P(n,k)/P(n+1,k+1)をnとkの式で表せ。
（2）lim[n→∞] Σ[k=1 to n] kP(n,k) を求めよ。
519:１３２人目の素数さん
18/11/06 00:53:08.60 0/M2gc6l.net
>>472-474
×確率1/2で2倍になるから、54/2の3/2倍で40.5 (54の3/4倍で40.5)
○確率2/3で2倍になるから、54/2の5/3倍で45 (54の5/6倍で45)
任意の数字をＮとすると、Ｎ,Ａ,Ｊを引く順番と期待値は
Ｎ→Ａ→Ｊ (1/6)*(2Ｎ)
Ｎ→Ｊ (1/6)*Ｎ
Ａ→Ｎ→Ｊ (1/6)*(2Ｎ)
Ａ→Ｊ (1/6)*0
Ｊ (2/6)*0
合計すると (5/6)*Ｎ
求める期待値は (5/6)*(2+3+4+5+6+7+8+9+10)＝(5/6)*54＝45
520:１３２人目の素数さん
18/11/06 01:04:46.33 LVSol2sI.net
>>497
なる！
521:１３２人目の素数さん
18/11/06 04:11:21.55 Py2gjw7X.net
>>496
(2)
Q(n,k) を少なくとも一回同じ色がk回連続引く確率とする。
m = [log √n]、l = [n/m] とおく。1～n回のコイントスのなかから連続する m 回のコイントスを重複しないように l 回に分けることが出来る。
T[1]～T[l] をそのような m 回のコイントスとしT[i]がすべて表になる事象をX[i]とするとき
P(X[i])=2^(-[log√n]) > 2^(-log√n) > 1/√n
である。
すべての i でX[i]が起こらない事象をYとするとき
P(Y) < (1-1/√n)^l < (1-1/√n)^(n/m) < exp (n/m) log(1-1/√n) < exp(-(√n/m))。
よって
P(n 回中 [log √n] 回連続表がでる)>(1-exp(-(√n/m)))。
∴Σ[k=1 to n] kP(n,k)
＞Σ[k=[log √n] to n] [log √n]P(n,k)
＞[log √n](1-exp(-(√n/m)))→∞。
522:１３２人目の素数さん
18/11/06 06:54:30.95 FZJllfOU.net
数列{a_n}は
a_1 = 1
a_(3n+1) = a_(2n+1)
a_(3n-1) = a_(2n-1)
a_(3n) = -a_n,
を満たす。この時、 lim(n→∞) (1/n)Σ(k=1，n) a_k
523:を求めよ。 (面白スレ２８ より) このスレも残り半分になりました。
524:１３２人目の素数さん
18/11/06 08:34:55.38 D5qaO8Cz.net
R^2の部分集合で単連結であるが可縮でないものは存在しますか？
525:１３２人目の素数さん
18/11/06 12:58:40.79 UM6as+XG.net
これの解き方がわかりません
考え方を教えてください
URLﾘﾝｸ(i.imgur.com)
526:１３２人目の素数さん
18/11/06 13:15:09.58 cDO4b4Dm.net
>>502
Prelude> [(a4,a3,a2,a1,x2,x1)|a4<-[1..9],a3<-[0..9],a2<-[0..9],a1<-[0..9],x2<-[1..9],x1<-[0..9],
(a4*1000+a3*100+a2*10+a1)*(x2*10+x1)==119868,(a4*1000+a3*100)*(x2*10+x1)==117600]
[(1,4,2,7,8,4),(2,8,5,4,4,2),(4,2,8,1,2,8)]
a=4281 x=28
527:１３２人目の素数さん
18/11/06 13:18:22.55 Nfu+AqXq.net
>>502
aの下二桁とxを掛けると2268
aの上二桁とxを掛けると1176
それぞれ素因数分解する
それらを見比べるとaの下二桁、上二桁にそれぞれ必ず含まれる因数がわかる
そこに残りの因数をどれだけ移せるかを考える
4281と28で合ってる？
528:１３２人目の素数さん
18/11/06 13:30:13.05 cDO4b4Dm.net
>>502
差を取って素因数分解すれば2*2*3*3*3*3*7になるので
あとは自分で考える
529:１３２人目の素数さん
18/11/06 14:15:47.28 1vm1mAgK.net
陰関数の分野なのですが分からないのでお願いします
URLﾘﾝｸ(i.imgur.com)
530:１３２人目の素数さん
18/11/06 16:16:20.93 mYPVHHox.net
>>460
x = ∫(1/y ') dy
= ∫√{y/(cy-2)} dy
= (1/c)√{y(cy-2)} + {2/c^(3/2)}･log[c√y + √{c(cy-2)} ] + c'
ここの2行目から3行目の変形わかりませんでした
出来るもんだと思い込んでただけで全くわかりません
どうやって導き出したんでしょうか？
531:１３２人目の素数さん
18/11/06 17:16:31.42 D5qaO8Cz.net
>>501
リーマンの写像定理から自明でした
532:１３２人目の素数さん
18/11/06 17:55:03.44 KdeHy8c/.net
半径3の円Pの外側に接している半径1の円QはPを一周するといくら回転しますか
533:１３２人目の素数さん
18/11/06 18:06:14.34 0/M2gc6l.net
単純に3回転ではないっていうことなのか
534:１３２人目の素数さん
18/11/06 18:13:13.05 Hp5bh8qb.net
>>501
んなわけない。
開集合じゃなきゃどうすんの？
{x^2+y^2≦1} ∪ [1,2] ∪ {(x-4)^2 + y^2≦4}
とか。
いくらでも複雑な例作れるよ？
んな簡単なわけない。
535:１３２人目の素数さん
18/11/06 18:15:51.91 Jhril6/D.net
次の性質(A)(B)をともに持つ2つの無理数a,bを求めよ。
(A)a^bは自然数
(B)任意の有理数pに対して、a^pは無理数
536:１３２人目の素数さん
18/11/06 18:16:17.29 D5qaO8Cz.net
>>511
申し訳ありません
途中から開集合で考えていたので整合性がとれていませんでした
開集合でなくてもいいならトポロジストの正弦曲線からも作れたりしますよね
537:１３２人目の素数さん
18/11/06 18:21:38.21 hE/3xu/H.net
>>512
a=e、b=log2
538:１３２人目の素数さん
18/11/06 18:23:10.78 0/M2gc6l.net
なるほど、円の中心が動く距離を円周で割るわけね
2π*(3+1)/2π*(1)＝4回転
539:１３２人目の素数さん
18/11/06 18:35:51.54 HFC0B7nW.net
なぜ6π/2πでは間違いなのですか？
540:１３２人目の素数さん
18/11/06 18:58:34.46 wFMVmkUM.net
エストラテネスの篩で限界桁ってどの辺ですか？4桁辺りですか？
541:１３２人目の素数さん
18/11/06 19:13:08.28 CIeCcXbf.net
>>517
限界などない
無限にいける
542:１３２人目の素数さん
18/11/06 19:15:24.37 4DhksDMu.net
4回転だと滑ってない？
543:１３２人目の素数さん
18/11/06 19:44:40.09 LENI/FKH.net
>>506
URLﾘﾝｸ(i.imgur.com)
544:１３２人目の素数さん
18/11/06 19:54:43.26 0/M2gc6l.net
円の中心が動く距離は円周の長さと同じ
545:１３２人目の素数さん
18/11/06 20:24:28.95 FZJllfOU.net
>>507
y = (2/c) cosh(t)^2,
cy - 2 = 2 sinh(t)^2,
546:dy = (4/c) sinh(t) cosh(t) dt, とおいてみる。
547:１３２人目の素数さん
18/11/06 20:46:17.63 jOazYBXJ.net
.
∧__∧
( ´･ω･)∧∧l||l
/⌒ ,つ⌒ヽ ) <>>353
(___ ( __)
"''"" "'ﾞ''` 'ﾞ ﾞﾟ' ''' '' ''' ﾟ`
548:１３２人目の素数さん
18/11/06 20:55:26.30 mYPVHHox.net
>>522
ありがとうございます
coshとか使うの初めて
549:１３２人目の素数さん
18/11/06 21:07:27.04 x61GGXRC.net
>>492どうなりましたか？
550:１３２人目の素数さん
18/11/06 21:09:29.41 0/M2gc6l.net
半径が同じ円なら、外接円は２回転するのか
でも、固定円の円周をちょんぎって直線にしたら1回転になる？
不思議と言えば不思議だな
551:１３２人目の素数さん
18/11/06 21:41:21.11 VbDhzeiW.net
n次元球面上でf(x1,x2,...,xn)=Σ[i=1,n]|xi|^pの極値ってどう求めればいいのでしょうか
552:１３２人目の素数さん
18/11/06 21:48:29.41 Kup5u5BK.net
>>527
とりあえずxi≧0に限定してti = √xiとおいて
f(x) = Σti^(p/2)、Σti = 1
なので凸不等式使えばいいのでは？
553:１３２人目の素数さん
18/11/06 21:49:05.38 Kup5u5BK.net
>>528
訂正 ti = xi^2です。
554:１３２人目の素数さん
18/11/06 21:51:46.89 d1IEKaLp.net
>>526
ぐるぐるした線なら円以上にたくさん回転するよ
こういうの
URLﾘﾝｸ(i.imgur.com)
555:１３２人目の素数さん
18/11/06 21:55:25.11 VbDhzeiW.net
>>528
各成分の正負が一致してる場合はその方法やラグランジュなりで解けるのですが
一致しない場合、各成分で正負が違う場合の求め方が何とも…
556:１３２人目の素数さん
18/11/06 23:11:55.70 JyIr9Vjq.net
>>497
1/6+2/6=1/2の確率でNを引かずにおわるかなぁ？
557:１３２人目の素数さん
18/11/06 23:33:07.59 08uZxk9P.net
>>531
成分が負の場合を扱うなら 負^p を扱う事になるのでそこをどうするか決めないと答え出ない希ガス。
pが整数の時しか考えないとか。
558:１３２人目の素数さん
18/11/07 00:56:52.41 /CQ+FCaa.net
>>352
そのＮは9枚じゃなくて1枚ね。
最後に9枚分の期待値を単純合計してもいいという理屈は
自分でもあまり良く分かってない。(問題>>467)
559:１３２人目の素数さん
18/11/07 01:27:50.13 /CQ+FCaa.net
Ｎ,Ａ,Ｊの合計3枚でも1/2で終わるのが疑問ということなのかな？
いきなり終了のＪが1/3もあるから、そんなに直感に反しないと思うんだけど
560:１３２人目の素数さん
18/11/07 01:45:36.54 igCuCTm9.net
では>>467の出てる答えを清書。
i:2～10に対し確率変数X[i]を
X[i] = 2i (i A J)
i (i J A)
2i (A i J)
0 (A J i)
0 (J i A)
0 (J A i)
とおく。
E(X[i]) = 5/6iである。
よって
E(得点) = Σ E(Xi) = Σ 5/6i = 5/6(2+3+…+10) = 5/6×54 = 45。
561:１３２人目の素数さん
18/11/07 04:44:41.88 LkOhmL9N.net
>>460>>522
x = ∫(1/y ') dy
= ∫√{y/(cy-2)} dy
y = (2/c) cosh(t)^2,
cy - 2 = 2 sinh(t)^2,
dy = (4/c) sinh(t) cosh(t) dt
とおくと
x = ∫{4/c^(3/2)}･cosh(t)^2 dt
= ∫{2/c^(3/2)}･(cosh(2t) + 1) dt
= {1/c^(3/2)}･sinh(2t) + {2/c^(3/2)}t + c'
ここで
y = (2/c) cosh(t)^2
= (2/c)･{(e^(2t) + e^(-2t) + 2)}/4
e^(2t) = s とすると
y = (2/c)･{(s + 1/s + 2)}/4
ｓでそろえると
s^2 - 2(cy-1)s + 1 = 0
{s - ((cy-1))}^2 = 0
s = cy-1
e^(2t) = cy-1
両辺にlogすると
t = (1/2)log(cy-1)
ゆえに
x = (1/c)√{y(cy-2)} + {1/c^(3/2)}･log(cy-1) + c'
となりましたが答えが合いませんでした
どこで間違えたのでしょうか？
562:１３２人目の素数さん
18/11/07 06:20:04.19 a52hrceZ.net
>>514
(B)はどうやって証明する？
563:１３２人目の素数さん
18/11/07 06:35:30.56 p6NUZQ5G.net
>>534-535
[,1] [,2] [,3]
[1,] A J N
[2,] A N J
[3,] J A N
[4,] J N A
[5,] N A J
[6,] N J A
この各行が同様に確からしく起こるってことでいいんだな。
564:１３２人目の素数さん
18/11/07 07:11:39.88 Jai49TZi.net
URLﾘﾝｸ(imgur.com)
OM ≧ OS が成り立つのはなぜでしょうか？
565:１３２人目の素数さん
18/11/07 07:19:10.07 cuI8wqvL.net
>>538
URLﾘﾝｸ(ja.wikipedia.org)
566:１３２人目の素数さん
18/11/07 10:08:30.39 5PMwby1T.net
>>537
s^2 - 2(cy-1)s + 1 = 0
まで正しいが、次から違っている。
{s - (cy-1)}^2 = cy(cy-2),
∴ s = (cy-1) ± √{cy(cy-2)},
e^t = √s = [ √(cy) ± √(cy-2) ] / √2,
t = log[ √(cy) ± √(cy-2)] - (1/2)log(2),
567:１３２人目の素数さん
18/11/07 10:31:01.79 lbG9qHeU.net
>>526
100円玉２個でやってみた
568:１３２人目の素数さん
18/11/07 11:03:47.05 y8xSy0+/.net
1 + 2 + 3 + 4 + … = −1/12
569:１３２人目の素数さん
18/11/07 12:23:18.07 5+J1KYD8.net
>>543
助かった ありがとう
570:１３２人目の素数さん
18/11/07 13:23:20.75 Jai49TZi.net
初等幾何の問題です。
OM ≧ OS が成り立つのはなぜでしょうか？
URLﾘﾝｸ(imgur.com)
571:１３２人目の素数さん
18/11/07 13:46:40.15 vh7jv6Dh.net
>>546
Mのx座標≧Sのx座標、Mのy座標≧Sのy座標だから。
572:１３２人目の素数さん
18/11/07 13:52:53.80 /U8eQOBM.net
>>544
そんなもん有難がってもいい事ないぞ
573:１３２人目の素数さん
18/11/07 14:02:51.99 a52hrceZ.net
次の性質を持つ数aを虚実数と呼ぶ。
・a^2≦0かつa^2>0
虚実数は、実数kと虚数単位iでは表せないことを示せ。
574:１３２人目の素数さん
18/11/07 14:05:07.47 VCjsN+3P.net
存在すらせんな
575:１３２人目の素数さん
18/11/07 14:06:39.25 4Mx2PGQZ.net
この自作問題連投ガイジ、中学生レベルの数学力すらないな
576:１３２人目の素数さん
18/11/07 14:06:52.87 Jai49TZi.net
>>547
ありがとうございました。
577:１３２人目の素数さん
18/11/07 15:48:06.77 a52hrceZ.net
一辺の長さが2の正三角形△ABCが平面z=0の円x^2+y^2=4/3に内接しており、A(4/3,0,0)である。
動点Pについて、f(P)=PA^2+PB^2+PC^2とする。
（1）kをk≧√3である実数とする。動点Pが平面z=0上を動くとき、f(P)=kとなる点Pの軌跡で囲まれる面積を求めよ。
（2）動点Pが空間を動くとき、f(P)=kとなる点Pの軌跡で囲まれる面積を求めよ。f(P)=kとなる点Pの軌跡で囲まれる体積を求めよ。
578:１３２人目の素数さん
18/11/07 15:53:51.49 a52hrceZ.net
>>553
（2）は鮮やかに解いてください
積分してもいいですが
579:１３２人目の素数さん
18/11/07 15:56:40.28 XKKQe3SL.net
さっさと>>492とけや無能ども
580:１３２人目の素数さん
18/11/07 16:05:20.14 PN+gm2kl.net
>>521
タイヤが道路と設置した長さだけ自動車は移動していると言われれば納得できるな。
581:１３２人目の素数さん
18/11/07 16:07:02.47 PN+gm2kl.net
>>556
設置→接地
582:１３２人目の素数さん
18/11/07 16:24:56.36 2p4APZKr.net
>>553
|p-a|^2+ |p-b|^2+ |p-c|^2
=3|p-g|^2-2(|g-a|^2+ |g-b|^2+ |g-c|^2)
583:１３２人目の素数さん
18/11/07 18:40:52.43 63cdf+8Y.net
これの解き方がわかりません
考え方を教えてください
URLﾘﾝｸ(i.imgur.com)
584:１３２人目の素数さん
18/11/07 19:02:50.93 oBXZsNC8.net
>>559
上で回答ついてるじゃんか
585:１３２人目の素数さん
18/11/07 19:26:52.16 PN+gm2kl.net
>>559
差をとると2268、これはaの下二桁とｘを掛けた値。
素因数分解すれば2*2*3*3*3*3*7になる。
aの下二桁をｙとすると一方が22以下なら他方が2桁にならないから23以上。
2268の平方根は47.62352なのでx,yのいずれかは47以下である。
2 2 3 3 3 3 7の積がこの範囲にあるのは27 28 36 42の4つ
その�
586:ﾆきの他方の数は84 81 63 54 あとは自分で考える。
587:１３２人目の素数さん
18/11/07 19:31:20.54 PN+gm2kl.net
>>561
組み合わせを考えるのが面倒だから
2268の約数で23以上、47以下はプログラムにやらせた。
> x=2268/(23:47)
> y=(23:47)[x-floor(x)==0]
> z=2268/y
> rbind(y,z)
[,1] [,2] [,3] [,4]
y 27 28 36 42
z 84 81 63 54
588:１３２人目の素数さん
18/11/07 19:40:28.76 5ORoIydt.net
この回答ヤバすぎだろ
589:１３２人目の素数さん
18/11/07 19:40:48.14 PN+gm2kl.net
>>561
ついでだから続きも書いておく。
xの候補が27 28 36 42 84 81 63 54に絞られたので
119868を割り切るのは
28 42 84の３つ
そのときの商は 4281 2854 1427
でこれがaの候補。
最大は4281でそのときのｘは28
590:１３２人目の素数さん
18/11/07 19:48:34.69 PN+gm2kl.net
Rだと
> x=2268/(23:47)
> y=(23:47)[x-floor(x)==0]
> z=2268/y
> b=c(y,z)
> c=119868/b
> d=b[c-floor(c)==0]
> (a=max(119868/d))
[1] 4281
> 119868/a
[1] 28
>
Haskellだと1行ですんだ。
Prelude> [(a,x)|a4<-[1..9],a3<-[0..9],a2<-[0..9],a1<-[0..9],x2<-[1..9],x1<-[0..9], let a=(a4*1000+a3*100+a2*10+a1),let x= x2*10+x1,a*x==119868,(a4*1000+a3*100)*(x2*10+x1)==117600]
[(1427,84),(2854,42),(4281,28)]
591:１３２人目の素数さん
18/11/07 19:53:55.28 PN+gm2kl.net
こっちの方が可読性がいいかな。
Prelude> [(a,x)|b<-[10..99],c<-[0..99],x<-[10..99],let a=100*b+c, a*x==119868,100*b*x==117600]
[(1427,84),(2854,42),(4281,28)]
592:１３２人目の素数さん
18/11/07 20:39:01.72 LkOhmL9N.net
>>542
ありがとうございます
式を変形して綺麗にしてあるのもようやく理解出来ました
わからないところは>>460の x = (1/c)√{y(cy-2)} + {2/c^(3/2)}･log[c√y + √{c(cy-2)} ] + c'
この式の log[c√y + √{c(cy-2)} ] の部分で+の部分が±ではない理由は何でしょうか？