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

1 名前：１３２人目の素数さん mailto:sage [2021/02/17(水) 00:42:07.42 ID:pOGUunX7.net] さあ、今日も１日がんばろう★☆ 前スレ 分からない問題はここに書いてね 465 rio2016.5ch.net/test/read.cgi/math/1608546793/ (使用済です: 478) 404 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 07:46:37.54 ID:ut6hP2sf.net] >>373 賞品が賞金1万円として太郎が勝ったら次郎と折半することにすると 談合することで太郎と次郎の賞金獲得金額の期待値は2000円から約2340円に上がるということになるなぁ。 シミュレーションプログラムを組んで100万回ジャンケンさせてみた。 > mean(replicate(1e6,sim(5))) [1] 0.469078 イナ氏の解の通りの値（近似値）が返ってきたのでバグはなさそう。 オマケ Rのコードはここ https://egg.5ch.net/test/read.cgi/hosp/1612996282/830 405 名前：１３２人目の素数さん [2021/03/27(土) 09:20:40.03 ID:NO5KFqjv.net] >>396 正解！ 406 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 10:24:57.05 ID:uBGH9Pl9.net] >>390 これお願いします 407 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 10:44:52.83 ID:Uaw5muz7.net] >>397 期待値も応召義務も分かってないお前は出禁だぞ 408 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 14:56:57.35 ID:uBGH9Pl9.net] xyz空間において、x軸からの距離が1以内かつ、y軸からの距離が1以内かつ、z軸からの距離が1以内である領域Dを考える。 Dに含まれる正四面体の中で一辺の長さが最大であるものを求めよ。 409 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 15:22:36.54 ID:i1fChCKG.net] Dは球 xx+yy+zz ≦ 3/2 に含まれ、等号成立は 　(x,y,z) = (±1/√2, ±1/√2, ±1/√2) のとき。 　(1/√2, 1/√2, 1/√2),　(1/√2, -1/√2, -1/√2),　(-1/√2, 1/√2, -1/√2),　(-1/√2, -1/√2, 1/√2) および 　(1/√2, 1/√2, -1/√2),　(1/√2, -1/√2, 1/√2),　(-1/√2, 1/√2, 1/√2),　(-1/√2, -1/√2, -1/√2) は正四面体をなし、一辺の長さは２ 410 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 16:02:13.35 ID:0sVMdYrW.net] 円というのは多角形の究極の姿なんだよ。 言うなれば無限多角形。 この「無限」さ故に円周率が無限になるんだよ。 411 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 16:15:56.41 ID:/54LmiBv.net] 言葉遊びで考えたつもり 412 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 16:41:31.82 ID:q+bVIr97.net] 重心を(0,0)に、垂心を(1,0)に持つ三角形の、内心からなる領域を求めよ って問題をやってみたい気がするね そもそも(3,1)がそれに含まれるとは思えない 413 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 17:05:54.29 ID:j6Ec7ZOp.net] >>376 負の二項分布の期待値の公式から　n/p + n/(1-p)　でいいのかな？ 414 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 17:15:37.45 ID:j6Ec7ZOp.net] >>406 加算するのはだめだな。 415 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 17:46:47.98 ID:duPI2AD9.net] >>384 nとｐに具体的な数値を入れて100万回シミュレーションで検証。 （実は、シミュレーションの検証） n=10 p=0.25 > mean(replicate(1e6,sim(n,p))) [1] 13.30754 > E(n,p) [1] 13.30616 シミュレーションのコード（R言語） https://egg.5ch.net/test/read.cgi/hosp/1612996282/831 Wolframに入れてみたけど、超幾何関数とかでてきた。 Sum[(n + k) Binomial[n + k - 1, n - 1] p^(n - 1) (1 - p)^k p, {k, 0, n - 1}] + (n + k) Binomial[n + k - 1, k] p^k (1 - p)^(n - 1) (1 - p) 416 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 18:08:32.25 ID:Uaw5muz7.net] プロおじは引っ込んでろ 期待値の勉強でもしてろ 417 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 18:27:08.88 ID:uBGH9Pl9.net] >>408 無駄な解答ご苦労さまです 418 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 18:36:23.83 ID:duPI2AD9.net] >>408 ｐの値を変化させて100回の表もしくは裏がでるまでの試行回数をグラフ化。 当然、p=0.5のときが最大の左右対称のグラフになった。 https://i.imgur.com/bBFQmJl.png 419 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 18:39:56.73 ID:duPI2AD9.net] >>410 シミュレーションプログラムを作るのが楽しいんだなぁ。 >359のシミュレーションは予想外の値が返ってきたのだが、数理解と近似していてホッとした。 420 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 18:44:11.25 ID:X8dc1afi.net] それを他人が見て楽しいと思えない事を想像できないのが発達障害 421 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 19:01:16.31 ID:duPI2AD9.net] >>390 三角形ABCの座標をA(a1,a2) B(b1,b3),C(c1,c2)とすると 6連立方程式 a1+b1+c1=0 a2+b2+c2=0 (a1*(a2*(b1-c1)-b1*b2+c1*c2)+(b2-c2)*(a2^2-a2*(b2+c2)+b1*c1+b2*c2))/(a1*(c2-b2)+a2*(b1-c1)-b1*c2+b2*c1)=1 (a1^2*(b1-c1)+a1*(a2*b2-a2*c2-b1^2+c1^2)+a2*(c1*c2-b1*b2)+(b1-c1)*(b1*c1+b2*c2))/(a1*(b2-c2)+a2*(c1-b1)+b1*c2-b2*c1)=0 (a1*sqrt((b1-c1)^2+(b2-c2)^2)+b1*sqrt((c1-a1)^2+(c2-a2)^2)+c1*sqrt((a1-b1)^2+(a2-b2)^2))/(sqrt((a1-b1)^2+(a2-b2)^2)+sqrt((b1-c1)^2+(b2-c2)^2)+sqrt((c1-a1)^2+(c2-a2)^2))=3 (a2*sqrt((b1-c1)^2+(b2-c2)^2)+b2*sqrt((c1-a1)^2+(c2-a2)^2)+c2*sqrt((a1-b1)^2+(a2-b2)^2))/(sqrt((a1-b1)^2+(a2-b2)^2)+sqrt((b1-c1)^2+(b2-c2)^2)+sqrt((c1-a1)^2+(c2-a2)^2))=1 を解けばいいみたいだな。 422 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 19:05:45.21 ID:duPI2AD9.net] >>413 プログラムで見取り図や３D動画を作れないような人は楽しめないだろうなぁ。 423 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 19:17:13.42 ID:X8dc1afi.net] >>415 そうやって全部自分のいい方にしか解釈できない もちろん心の中の理性は自分の方が間違っていふ事を百も理解してる しかし自分が”負ける”という自分にとって耐えられない事態に陥る事を許さない人間は客観的事実を捻じ曲げて解釈してでも自分にとって都合の良い解釈を採用する その事を恥ずかしいと思う“心の制御装置”はもうとっくに壊れてしまっている 424 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 20:06:10.13 ID:i1fChCKG.net] >>393 ｐ∫[-π,π] cos(mt) cos(nt) dt 　= (1/2)∫[-π,π] {cos((m+n)t) + cos((m-n)t)} dt 　= π (δ_{-m,n} + δ_{m,n}), 425 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 21:01:08.17 ID:Uaw5muz7.net] >>415 誰にも相手にされてないどころかゴミ扱いのにレスし続けて楽しいか？社会や家族はおろか5chですら必要とされてないなんて哀れだな。 426 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 21:04:43.57 ID:Vv3CtIpR.net] シミュレーション向きの問題です 不等式 y≦x^2-4 で表されるxy平面上の領域に含まれる線分の中で、最長のものを求めよ。 427 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 21:11:08.55 ID:X8dc1afi.net] ∞ 428 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 21:14:52.69 ID:YQhhnaZL.net] 存在しないってことになるんじゃ？ 429 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 21:36:36.26 ID:X8dc1afi.net] 存在しないな 430 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 21:38:50.13 ID:pE5L1fWz.net] 何がどうシミュレーション向きなのかさっぱり不明だが そもそもここはクイズスレでも自作問題を出題するスレでもない 431 名前：１３２人目の素数さん mailto:sage [2021/03/27(土) 22:12:28.94 ID:/54LmiBv.net] 終わらないシミュレーションで厄介払いできる の意味じゃね？ 432 名前：１３２人目の素数さん [2021/03/27(土) 23:36:44.95 ID:nTP63oRo.net] n次の巡回行列とn次元のベクトルの積はΘ(n*log(n))で計算できることを示せ。 433 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 00:00:18.70 ID:sQHcAkjP.net] n回の掛け算とn-1回の足し算なんだから当たり前じゃないの？ 434 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 02:01:15.88 ID:Eu8CzLjp.net] シミュレーション向きかどうか分からない問題です。 不等式 　x^2 - 4 ≦ y ≦ 0 で表されるxy平面上の領域に含まれる線分の中で、最長のものを求めよ。 435 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 02:38:02.88 ID:FT1ycNw8.net] シミュレーションとか必要ない √(71/4+2√2)≒4.53634512848683 >>427 436 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 02:44:33.29 ID:sQHcAkjP.net] f(x)=x^2-4としてmax{√((x-2)^2+f(x)^2)}の最大値が線分の長さの最大値 これはx=1/√2-1の時最大をとる よって求める線分は(2,0)と(1/√2-1,f(1/√2-1))を結ぶ線分とそのy軸対称の２つである https://www.wolframalpha.com/input/?i=maximize+sqrt%28%28x-2%29%5E2%2B%28x%5E2-4%29%5E2%29+in+-2+%3C%3D+x+%3C%3D+2&lang=ja 437 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 02:49:27.20 ID:sQHcAkjP.net] >>426 勘違いorz 438 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 04:19:47.88 ID:T+i43wxS.net] >>405 の解として (1/2, 0) からの距離が 1/2 未満の領域 という予想を置いておく 439 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 04:33:57.49 ID:T+i43wxS.net] 言い換えると、 正三角形を除く三角形の内心は、重心-垂心間の線分を直径とする円の内側に存在する ということになる 証明はまだ 440 名前：１３２人目の素数さん [2021/03/28(日) 06:05:12.70 ID:/jK5jGei.net] >>426 当たり前なのは、O(n^2)で計算できることです。 441 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 06:44:30.96 ID:Jul26fm0.net] >>425 n次の巡回行列とn次元のベクトルの積は、n次元ベクトルの巡回畳み込みに相当する n次元ベクトルの巡回畳み込みは、畳み込み定理(Convolution theorem)により 高速フーリエ変換(FFT)を使ってn*log(n)のオーダーで計算できる 442 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 08:57:46.20 ID:iiFYE3Wp.net] >>429 線分の一方の端点が(2,0)になるのは感覚的には分かるけど、証明できる？ 443 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 09:03:41.21 ID:iiFYE3Wp.net] ごめん自己解決 444 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 09:45:19.86 ID:sQHcAkjP.net] >>434 wikiにはnが2べきならnlon(n)って明言してるけどそれ以外でもnlog(n)でできるんですか？ 445 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 10:05:04.30 ID:/Dlncm5P.net] シミュレーション向きの問題です 不等式 x^2-4≦y≦0 で表されるxy平面上の領域に含まれる線分の中で、最長のものを求めよ。 446 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 10:09:51.08 ID:/Dlncm5P.net] 数値計算向けの問題です 方程式 x^3-x=3 の各実数解の小数点以下2桁目の数字を求めよ。 447 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 10:13:29.34 ID:/Dlncm5P.net] 数値計算向けの問題です 整数2nCnの桁数をnで表せ。 注：2nCn=C[2n,n] 448 名前：１３２人目の素数さん [2021/03/28(日) 10:50:33.87 ID:/jK5jGei.net] >>437 できます。 449 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 12:20:07.69 ID:/Dlncm5P.net] nを正の整数の定数とする。 xyz空間の立方体領域 D: -n≦x≦n,-n≦y≦n,-n≦z≦n を考える。 Dに含まれる格子点で、z=xyかつy=xzを満たすものの個数をnで表せ。 450 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 12:39:45.75 ID:sQHcAkjP.net] >>431 どうやってやるんですか？ とりあえずwikiに載ってるアルゴリズムはnが自明でない分解n=mlを持つとき周期mでの変換と周期lでの変換に還元して‥として行くとあるのでnが小さい素因子の積になってないと大して高速化できない、というか全然高速化できない 451 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 14:25:46.01 ID:Eu8CzLjp.net] >>442 条件式を辺々掛けて 　yz = xxyz, 　(x+1)(x-1)yz = 0, 　x=±1 または y=0 または z=0,　　(←整域) 後の2つは、条件式から y=z=0 　(1,y,y)　(-1,y,-y)　(x,0,0) ∴　6n+1 個 452 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 14:45:46.48 ID:Eu8CzLjp.net] >>439 　f(x) = x^3 - x - 3　とおく。 ニュートン法で 　x　→　x - f(x)/f '(x) = x - (x^3 -x-3)/(3xx-1), 　5/3　→　5/3 + 1/198　→　5/3 + 1/198 - 991/57316842　→　1.67169988 実数解は１つしかねぇ 453 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 15:22:16.93 ID:Eu8CzLjp.net] 　y = x^3 - x - 3 　　= (x-5/3)^3 + 5(x-5/3)^2 + (22/3)(x-5/3) - 1/27, x=5/3 で接線を曳く. 　y ≒ (22/3)(x-5/3) - 1/27, y=0　とおくと　x-5/3 ≒ 1/198 　x ≒ 5/3 + 1/198 = 1.671717… 454 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 17:10:01.34 ID:Eu8CzLjp.net] どうでもいいことだが 　X = (3/2 - 180/121)^{1/3} + (3/2 + 180/121)^{1/3} = 1.671698593 は 　X^3 - [1 - (1/242)^2]^{1/3}･X - 3 = 0, の実数解 455 名前：１３２人目の素数さん [2021/03/28(日) 17:58:38.66 ID:/jK5jGei.net] >>443 ゼロで埋めて、nが2の累乗になるようにすればいいのではないでしょうか？ 456 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 18:43:31.48 ID:Eu8CzLjp.net] >>435 　a = (1 + 1/√2)/2 = 0.8535534 　b = (5/2+√2)/2 = 1.95710678 　L = √(71/4 + 2√2) = 4.5363451285 とおく。 問題の領域は 直径がLである3円 　(x-a)^2 + (y+b)^2 ≦ (L/2)^2, 　(x+a)^2 + (y+b)^2 ≦ (L/2)^2, 　x^2 + (y+b-a√3)^2 ≦ (L/2)^2, により二重に被覆されている。doubly covered. どの点も 2つ以上の円に含まれる。 どの2点も同じ円に含まれるから 距離 ≦ L. 457 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 19:49:01.62 ID:fpb6SHPW.net] >>414 この方程式をプログラムで数値解を出して 重心が(0,0),垂心が(1,0)となる三角形をAの座標を乱数発生させて描いてみた。 https://i.imgur.com/RLXA6g8.png Gが重心、oが垂心、Iが内心 458 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 20:01:29.09 ID:gC1h70qb.net] >>448 無理でしょ？ (a,b,c,d,e)と(x,y,z,u,v)を周期5で考えて畳み込んだものと(a,b,c,d,e,0,0,0)と(x,y,z,u,v,0,0,0)を周期8で畳み込んだものは一致しないでしょ？ 459 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 20:10:43.19 ID:fpb6SHPW.net] >>450 乱数発生を1000回繰り返して内心となる座標を重ね合わせてみた。 https://i.imgur.com/njz826H.png 内心の分布は円形に収まりそう。(3,1)が内心となることはなさそうだな。 460 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 20:14:27.21 ID:fpb6SHPW.net] >>452 シミュレーションの結果は>432の予想を支持する。 プログラミングが楽しめた。 461 名前：１３２人目の素数さん [2021/03/28(日) 21:00:42.43 ID:/jK5jGei.net] >>451 0で埋めると巡回行列じゃなくなってしまいますね。 462 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 21:05:40.37 ID:iiFYE3Wp.net] >>452 おー分かりやすい 何個か円の外にはみ出てる点あるけどこれは何でなん？ 463 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 21:11:51.83 ID:gC1h70qb.net] >>454 少なくともFFTを使って高速化できるのはnに何の要件もなければ無理なんじゃないかな？ 少なくともwikiではnが小さい素因子をたくさんかけた形の高速化法しか載ってない 一般の場合でもできるなら方法そのものを載せるのは無理でも論文へのリファレンスがないのは考えられないし 出題ミスかな？ FFTでできないから一般にも無理とは言えないけど 464 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 21:20:31.38 ID:gC1h70qb.net] O:外心 I内心 H:垂心 OI^2/OH^2 https://www.wolframalpha.com/input/?i=+plot+%28%28sin%28x%2By%29%2Bsin%28x%29cos%282y%29%2Bsin%28y%29cos%282x%29%29%5E2%2B%28sin%28x%29sin%282y%29-sin%28y%29sin%282x%29%29%5E2%29%2F%28%281%2Bcos%282x%29%2Bcos%282y%29%29%5E2%2B%28sin%282x%29-sin%282y%29%29%5E2%29&lang=ja 465 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 21:31:11.30 ID:fpb6SHPW.net] >>455 多分、三角形を形成しないような乱数が選択されたのではと推測。 466 名前：１３２人目の素数さん mailto:sage [2021/03/28(日) 22:57:13.38 ID:9hClS6gi.net] >>458 プロおじは出禁だぞ。 467 名前：１３２人目の素数さん [2021/03/29(月) 00:35:57.63 ID:GgCLqWW4.net] 今回は役に立ってるからな おまえさんにゃ分が悪い 468 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 00:41:01.25 ID:PU1lRMXx.net] ミスリードしてるからあかんやろ 469 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 00:58:14.08 ID:JEuzD2OD.net] プロおじは期待値から勉強し直してきなさい。 470 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 01:22:09.05 ID:JXTeJTxs.net] 訂正 https://www.wolframalpha.com/input/?i=+plot+%28%28sin%28x%2By%29%2Bsin%28x%29cos%282y%29%2Bsin%28y%29cos%282x%29%29%5E2%2B%28sin%28x%29sin%282y%29-sin%28y%29sin%282x%29%29%5E2%29%2F%28%281%2Bcos%282x%29%2Bcos%282y%29%29%5E2%2B%28sin%282x%29-sin%282y%29%29%5E2%29%2F%28sin%28x%29%2Bsin%28y%29%2Bsin%28x%2By%29%29%5E2&lang=ja 471 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 02:37:12.22 ID:rt8CA3eO.net] 延々とプロおじ追い出そうとしてる奴いるけど正直うっとうしい レス内容が気に入らないなら見なきゃいいだけなのに 472 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 02:43:06.46 ID:aDVxIJBu.net] ハサミは使いよう等と言う 出禁にするのは勿体無い 473 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 02:48:42.75 ID:/r3M6nmW.net] 公園で遊んでたらオナニーしてる人がいて困ってる状況ですし コテつけてるならNGすればいいけど、そうしない時点で「プロおじ=公衆オナニー見せつけて喜んでるキチガイ」なんだもの 474 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 02:51:59.03 ID:9N6BVMjD.net] >>466 >>466 475 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 04:38:45.90 ID:JxT5eGiE.net] 五者択一の問題に連続４問正解したら合格の試験がある。 １問解答するのに１０００円を徴収される。 太郎君は１０万円を準備して試験に臨み、問題文は読まずにランダムに解答することにした。 (1)太郎君の合格する確率はいくらか？ (2)太郎君の合格確率を1/2にするにはいくら準備すればよいか？ 476 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 06:46:22.46 ID:Q46gHliY.net] wolframalpha>>(越えられない壁)>>プロおじ >>464 左翼は、よく左翼のふりをして外国贔屓の外国右翼に利用される｡ 下手にプロおじを静観しない方がスレの為｡ 477 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 13:29:04.38 ID:pDr3G3SZ.net] >>468 （１）は　表のでる確率が1/5のコインを100回投げて4回以上連続して表がでる確率と等しいかな？ 478 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 13:58:58.60 ID:lBxffmcv.net] プロおじは社会や家庭だけでなく5chにも居場所なんかないからな。 479 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 14:23:56.95 ID:5dpYUdde.net] >>445-447 更にどうでもいいことだが 　X = 5/3 + 1/198, は 　X^3 - (１ + 3/198^2)X - (３ - 2/198^3) = 0, の実数解 480 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 14:58:15.46 ID:5dpYUdde.net] 更に更にどうでもいいことだが 　X = 5/3, は 　X^3 - (１ - 1/45)X - 3 = 0, 　X^3 - X - (３ - 1/27) = 0, の実数解 481 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 15:02:28.77 ID:pDr3G3SZ.net] >>468 (１)　 0.1175311 (2) 542000円 482 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 15:14:24.79 ID:pDr3G3SZ.net] >>445 更にどうでもいいことだがｗ グラフ化 https://i.imgur.com/rQVl1Yv.png Newton-Raphson法での数値解 > f <- function(x) x^3 - x - 3 > curve(f(x),-5,5) > abline(h=0,lty=3) > uniroot(f,c(-5,5),tol=1e-24)$root [1] 1.67169988166 483 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 16:36:11.44 ID:c+KCNM8F.net] シミュレーション向きの問題です 以下の連立不等式で表される3次元空間の立体 x^2+y^2≦1 y^2+z^2≦1 z^2+x^2≧1 を図示せよ。 484 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 17:18:02.69 ID:5dpYUdde.net] >>468 (1)　0.1176679987 (2)　542000円 n問解答後に、最後のk問が正解の確率を Q_k(n) とする。 　(n-k が不正解で、n-k+1〜n が正解) 漸化式 　Q0(n) = (1-p){Q0(n-1) + Q1(n-1) + Q2(n-1) + Q3(n-1)}, 　Q1(n) = p･Q0(n-1), 　Q2(n) = p･Q1(n-1), 　Q3(n) = p･Q2(n-1), n問目で合格する確率は 　P(n) = p･Q3(n-1), 特性方程式 (1/(1-p))･t^4 - t^3 - p･t^2 - p^2･t - p^3 = 0, 特性値 (p=1/5) 　α = 0.998713391320282 541問までで合格　0.49989971143 542問までで合格　0.50054314480 485 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 17:54:16.00 ID:JXTeJTxs.net] 単にn回目までに合格する確率S(n)求めるだけならk=4として S(n)=S(n-1)+qp^k(1-S(n-k-1)) の方が早いがな 特殊解S(n)=1もすぐ見つかるし この漸化式から求まるP(n)の特性方程式は x^(k+1)-x^k=p^(k+1)-p^k もうこのネタ何十回見たやろ 486 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 19:31:54.37 ID:5dpYUdde.net] >>435 この領域を平行線で挟んだときの幅を考える。 平行線の傾きが正のときは点 (-2,0) を通り、 傾きが負のときは点 (2,0) を通る。 ∴　どちらにしても　幅 ≦ L, ∴　この領域内の線分の長さ ≦ L, 参考 　支持関数 (support function) 　凸領域 487 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 20:43:18.88 ID:L8k5fESM.net] >470の確率と一致するはずと思っての答が>474でしたが、>477と一致しないので プログラムを見直したらインデックスがひとつずれておりました。 >474は撤回して以下に修正 > P(100,1/5,4) [1] 0.1176679986993025 （2）の方は> P(541:543) [1] 0.4998997114295112 0.5005431448015157 0.5011857503265590 なので542000円のまま。 シミュレーションプログラムを組んで検証してみよう。 488 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 21:00:28.87 ID:L8k5fESM.net] 100万回のシミュレーションだと太郎君の合格確率は > mean(y[,2]) [1] 0.117868 と出た。 ちなみに太郎君が合格したときに手元に残るお金の平均値は > (100-mean(z))/10 # 合格したときの残金 [1] 4.91万円になった。案外手元に残るもんだな。 オマケ（Rのコードはここ） https://egg.5ch.net/test/read.cgi/hosp/1612996282/849 オマケのおまけ：総解答数と合格確率のグラフ https://i.imgur.com/NBTEzhb.png 489 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 21:04:17.06 ID:L8k5fESM.net] >>476 罵倒厨が３D見取り図を作る練習問題にいいかもな。 粘土でつくるかもしれんが。 490 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 21:17:05.17 ID:c+KCNM8F.net] シミュレーション向きの問題です （1）A=C[20212022,2021],B=C[20211011,2020]とするとき、KA=LBを満たす正の奇数K,Lが存在することを示せ。 （2）Aを4で割った余りを求めよ。 491 名前：１３２人目の素数さん mailto:sage [2021/03/29(月) 21:21:53.39 ID:lBxffmcv.net] 誰にも相手にされてないプロおじ虚しくないのか？ 492 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 01:11:01.96 ID:9OPtgdh7.net] https://www.wolframalpha.com/input/?i=binomialcoeff%2820212022%2C2021%29%2Fbinomialcoeff%2820211011%2C2020%29&lang=ja 493 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 01:11:15.54 ID:gtiIL6DW.net] >>432 0＜α＜π/2, -α＜β＜αとして、A(cos2α,sin2α), B(cos2β,sin2β), C(cos2α,-sin2α)をとる。 任意の三角形はこの形式に相似変換できる。 外心O=0, 重心G=((2cos2α+cos2β)/3,(sin2β)/3) 垂心H=(2cos2α+cos2β,sin2β) G-Hの中点M=(2(2cos2α+cos2β)/3,2(sin2β)/3) 内心I=(cos(β+α)+cos(β-α)-1,sin(β+α)+sin(β-α))=(2cosαcosβ-1,2cosαsinβ) 線分MHとMIの距離の二乗の差を考える。 MH^2-MI^2 =(((2cos2α+cos2β)/3)^2+((sin2β)/3)^2)-((2cosαcosβ-1-2(2cos2α+cos2β)/3)^2+(2cosαsinβ-2(sin2β)/3)^2) =4cosα(cosβ-cosα)((2cosα-cosβ)^2+1-(cosβ)^2)/3 �@ 0＜cosα＜cosβ＜1 なので �@≧0、よって MH≧MI 等号は(2cosα-cosβ)^2+1-(cosβ)^2=0のときのみ。cosβ=1,cosα=1/2なのでこのとき△ABCは正三角形。 よって、 ＞正三角形を除く三角形の内心は、重心-垂心間の線分を直径とする円の内側に存在する □ 494 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 03:32:52.72 ID:1nCYsZUW.net] >>456 FFTは任意のnで行えるし実装もされている（rader's fft algorithmという）。 wikipediaの情報が誤っているのは執筆者の学が足りないだけ。 495 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 03:43:46.05 ID:9OPtgdh7.net] 外接円の半径を1、外心O、垂心H、内心Iとする OI^2=3-2cosA-2cosB-2cosC OH^2=3+2cos2A+2cos2B+2cos2C OI.OH=2cosA+2cosB+2cosC-cos(A-B)-cos(B-C)-cos(C-A) 3IG.IH =(OH-3OI).(OH-OI) =2(cos2A+2cos2B+2cos2C) -14(cosA+cosB+cosC) +4(cos(A-B)+cos(B-C)+cos(C-A)) +12 =4(cosA+cosB+cosC)^2-10(cosA+cosB+cosC)+6 ≦0 (∵ cosA+cosB+cosC≦3/2) 496 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 03:53:39.25 ID:9OPtgdh7.net] 訂正 1≦cosA+cosB+cosC≦3/2 497 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 06:17:06.69 ID:hjMIVqxL.net] >>483 >468で太郎が合格したときに手元に残る金の期待値の方がシミュレーション向きだと思う。 498 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 07:22:13.42 ID:Yc3/2/EC.net] >>486 力作の投稿ありがとうございます。 三角形の五心の座標を算出する関数をおもちゃ箱から取り出して 　正三角形を除く三角形の内心は、重心-垂心間の線分を直径とする円の内側に存在する を体感してみました。 https://i.imgur.com/nGyW1AK.gif もとを辿れば、 >414の連立方程式をプログラムで数値解を出そうと探索させてみたけどして最適解がみつからず、 プログラムが間違っているのかと訝っていましたが、題意を満たす三角形は存在しないことが証明されて感動しました。 499 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 08:46:53.71 ID:Jgfj8cZF.net] >>487 そりゃ離散フーリエ変換自体は任意のnで定義されてるでしょ？ 問題はその計算量 これもwikiだのみだけど 畳み込み積分の計算をゼロ埋めなしのFFTで行う場合、計算時間はNの性質に強く依存する。最悪の場合、N-1 が素数 N2 により N-1= N2 と表され、また N2–1が素数 N3 により N2–1 = N3 と表され、以下同様に続いていく場合である。このような場合、レーダーのアルゴリズムの再帰が連続することになり、O(N²)の計算時間がかかる可能性がある。このような性質をもつNは ソフィー・ジェルマン素数と呼ばれ、上記の数列は一次の Cunningham(ビル-カニンガム)チェーンと呼ばれる。し� 500 名前：ｩしながら、これまでの研究ではカニンガムチェーンの成長はlog2(N)よりも遅いことが分かっているため、レーダーのアルゴリズムの再帰によりかかる計算時間はO(N²)よりかは速いと思われる。幸いにも、畳み込み計算にゼロ埋めを用いたFFTを使えば計算時間はO(N log N)のオーダーになることが保証されている。 とあるけど？ ゼロ埋め使えないし [] [ここ壊れてます] 501 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 11:04:46.46 ID:ux9g6nBw.net] >>492 イヤ、訂正 これはなんかwikiの方がおかしい気もする アルゴリズムの内容も読んでみたらnlog(n)で計算できる気もする 502 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 12:31:15.21 ID:1nCYsZUW.net] >>493 なぜゼロつめができないと決めつける？ 非巡回畳み込みはゼロつめで高々2倍のnで2^kの計算に持ち込めて定数を無視すると O(2n*log(2n))=O(n*log(n))で計算できる。 巡回畳み込みは非巡回畳み込みをn回足せばよいだけだからこれもn*log(n)オーダで計算できる。 したがって任意の非巡回、巡回畳み込みはn*log(n)オーダで計算でき、 任意のnの離散フーリエ変換もこのオーダーで計算できる。 またレーダーの方法以外にも素数nの離散フーリエ変換をO(n*log(n))で計算するアルゴリズムは少々技巧的だが存在する。 503 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 12:41:51.51 ID:85330UQl.net] >>494 イヤ、実際できないでしょ？ [[a,b,c],[b,c,a],[c,a,b]]と[p,q,r] のサイズが2冪になるように0詰めしたら [[a,b,c,0],[b,c,0,a],[c,0,a,b],[0,a,b,c]]と[p,q,r,0] になってこの積は高速に計算できるかもしれないけど、それは元の行列とベクトルの積とは一致しない オーディオ機器とかへの応用ででNのサイズに特に意味がないなら好きなだけ0詰めすればいいけど、今数学の問題で「行列とベクトルの積をFFTの理論を応用して高速に行いたい」というテーマなんだから0詰めなんてできないじゃん？ 504 名前：１３２人目の素数さん mailto:sage [2021/03/30(火) 12:49:56.29 ID:1nCYsZUW.net] >>495 訂正:高々2倍→高々4倍 (a,b,c)×(p,q,r)は (a,b,c,0,0,0,0,0)×(p,q,r,0,0,0,0,0)を計算して、 巡回してない部分を足し合わせることで計算可能

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef