高校数学の質問スレ Part410
■ このスレッドは過去ログ倉庫に格納されています
【質問者必読！！】
まず>>1-4をよく読んでね
数学@5ch掲示板用 掲示板での数学記号の書き方例と一般的な記号の使用例
http://mathmathmath.dotera.net/
・まずは教科書、参考書、web検索などで調べるようにしましょう。（特に基本的な公式など）
・問題の写し間違いには気をつけましょう。
・長い分母分子を含む分数はきちんと括弧でくくりましょう。
（× x+1/x+2 ； ○((x+1)/(x+2)) ）
・丸文字、顔文字、その他は環境やブラウザによりうまく表示できない場合があります。
どうしても画像を貼る場合はPCから直接見られるところに見やすい画像を貼ってください。
ピクトはPCから見られないことがあるので避けてください。
・質問者は名前を騙られたくない場合、トリップを付けましょう。
（トリップの付け方は 名前(N)に 俺！#oretrip ←適当なトリ)
・質問者は回答者がわかるように問題を書くようにしましょう。
でないと放置されることがあります。
（変に省略するより全文書いた方がいい、また説明なく習慣的でない記号を使わないように）
・質問者は何が分からないのか、どこまで考えたのかを明記しましょう。
それがない場合、放置されることがあります。
（特に、自分でやってみたのに合わないので教えてほしい、みたいなときは必ず書くように）
・回答者も節度ある回答を心がけてください。
・970くらいになったら次スレを立ててください。
※前スレ
高校数学の質問スレ Part409
http://rio2016.5ch.net/test/read.cgi/math/1608682829/ 入試だけでは測れない能力もあるからね。ワンチャンの試験では失敗もあるし。 そんなの極僅かでしょ？
ワンチャンスのテストで実力が出せる奴が本当の実力者
入試で測れないなら前期後期関係ないやん >>305ですけどこの問題除いても3完してるので可能性は残ってそうです >>314
だから「一部」なんでしょ。
統計的にならしちゃえば前期で受かるほうが優秀って
ことになるのは当たり前だし、そういう話ではない。 >>301
>295のアルゴリズムって
1000以下の合成数は√1000=31.68以下の素数 2 3 5 7 11 13 17 19 23 29 31の倍数であるということなので
それでプログラムを組めば
n=1000
pmax=floor(sqrt(1000))
p=(1:pmax)[-outer(2:pmax,2:pmax)][-1]
f=function(x) all(x%%p!=0)
primes=sort(c(p,(2:n)[sapply(2:n,f)]))
length(primes)
> n=1000
> pmax=floor(sqrt(1000))
> p=(1:pmax)[-outer(2:pmax,2:pmax)][-1]
> f=function(x) all(x%%p!=0)
> (primes=sort(c(p,(2:n)[sapply(2:n,f)])))
[1] 2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83
[24] 89 97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199
[47] 211 223 227 229 233 239 241 251 257 263 269 271 277 281 283 293 307 311 313 317 331 337 347
[70] 349 353 359 367 373 379 383 389 397 401 409 419 421 431 433 439 443 449 457 461 463 467 479
[93] 487 491 499 503 509 521 523 541 547 557 563 569 571 577 587 593 599 601 607 613 617 619 631
[116] 641 643 647 653 659 661 673 677 683 691 701 709 719 727 733 739 743 751 757 761 769 773 787
[139] 797 809 811 821 823 827 829 839 853 857 859 863 877 881 883 887 907 911 919 929 937 941 947
[162] 953 967 971 977 983 991 997
> length(primes)
[1] 168
√1000=31.68以下の素数
2 3 5 7 11 13 17 19 23 29 31は
結局のところ、しらみつぶしに列挙しているだけだから、
最初からプログラムで列挙させるのと何も変わらんよ。 百万以下の素数の数は664579以下であることを示せ、という問題すると
そのアルゴリズムでは、
1000までの素数を列挙してプログラムを組むことになるからね。 数理で暗算で列挙できるレベルまでに絞って解答できるように入学試験は作られているけど
暗算で列挙するのが前提なら、最初から計算機で列挙したって同じことに思えるんだなぁ。
関数f(n)はn以下の素数の数を表す関数とする。
例 ： f(1000)=168
この値が2021になる最小のｎの値を求めよ。
fをプログラムして
> y=sapply(1:100000,f)
> min(which(y==2021))
[1] 17579
で終了 >>317
まず素数のリストp_1(=2),p_2(=3),p_3(=5),...を作っておいて、
関数f(N,p_n)を
f(N,p_1)=[N/p_1]=[N/2]
f(N,p_n)=[N/p_n]- Σ(k=1,n-1)f([N/p_n],p_k)
と定義して、f(1000,2),f(1000,3),...を計算するってこと。
計算量が減るんじゃない？ >>319
１から虱潰し探すより、fは増加関数だからニュートンラフソンで計算した方が計算時間が短縮できる。
uniroot(function(n,u=2021) fn(n)-2021, c(1,1e5))$root
fn(17578)
fn(17579)
> uniroot(function(n,u=2021) fn(n)-2021, c(1,1e5))$root
[1] 17578.55
> fn(17578)
[1] 2020
> fn(17579)
[1] 2021 >>320
>素数のリストp_1(=2),p_2(=3),p_3(=5),...を作っておいて
1000までの素数リストは
(1:1000)[-outer(2:1000,2:1000)][-1]
このスクリプトなら、for loopを使わない行列計算だから結果がでるのも高速。
一行で完成するのに、なんで31以下の素数リストを列挙する手間をかける必要があるんだ、と思う。
百万以下の素数の数は664579以下であることを示せ、という問題だと、百万の平方根=1000以下の素数のリストが必要になる。 なんでやめろと言われても人の嫌がることをし続けるんですか？ 自分が他人の嫌がる事をしても他人がそれを止められないのをみて自分の優越性に感じる人種がいるんだって
他人が嫌がる事で自分の欲求を達成してる事に何の躊躇もない時点で完全に人格が破綻しているのを理解する能力がない
もうこの年齢になったら矯正も効かない
無視するしかない もっと切実に目を背けたいことがあるんじゃないか
なんか必死だし >>323
>301でプログラム依頼がきているからね。
代わりにあんたが答えてくれてもいいんだが。
百万以下の素数の数は664579以下であることを示せ、プログラムなしでやってみてくれ。 プログラムの数値解なんて、もう高校数学でも何でもないだろ
専用スレでやってくれよ
無ければ、自分で立てろ >>313
高校数学の範囲を超えるかもしれんがこんな計算をしてみるのも暇つぶしにいいな。
合格可能性はまったく未知でその確率分布は一様分布を仮定する。
ある受験生が一回受験して不合格になった。
次の受験で合格する期待値とその95％信頼区間を求めよ。
類題は前スレのこれ
エロ本自動販売機に何冊かに１冊無修正が紛れ込んでいるという噂があったな。
（問題）
１０冊買ってみたが全部モザイク付きであったとする。噂が正しい確率を求めよ。 プログラム計算は高校数学と分けるべきだな
それとも見てもらえないと思ってんのか？ >>331
それも小学生がウンコとかチンコとかいうのと同じ
精神のレベルが小学生レベルで成長が止まってる >>331
次の受験で合格する期待値とその95％信頼区間を出して励ましてあげればいいのに！ >>334
合格率の事前確率分布を一様分布と設定したとき、
一度不合格になったという条件付き確率（＝合格率）の期待値。 72になった
lambda = 50 + ( sqrt 2501 )
a 0 = 1
a n = floor $ ( * lambda ) $ fromInteger $ a $ n-1
main = do
print $ map ( flip mod 100 ) $ map a [ 0.. 100 ]
----
[1,0,0,99,99,98,98,97,98,56,96,76,48,36,16,80,48,80,40,4,96,80,4,64,40,96,20,72,60,56,56,52,44,60,68,72,28,92,36,36,92,40,92,92,16,56,16,16,80,8,12,48,28,48,28,60,48,16,60,68,28,8,76,48,44,32,72,96,16,16,44,28,56,92,60,84,36,68,12,40,40,8,32,0,64,28,24,88,24,92,56,4,76,84,76,16,56,80,12,8,72] >>336は正しいと思うんだけど
これを高校数学ではどう解くんだろうか なんだかな
プロおじ封じを目論んだのは理解できるが
そのために出題が超高校級になったんじゃ本末転倒 >>338
なんで？
オレのプログラムあってると思うけど？
オレなんか思い違いしてる？ >>339
すまん
コレの答えが51の理由書いてくれん？
超高校級でもいいから 1234567890
24680
3692581470
48260
50
62840
7418529630
86420
9876543210
0
周期性と1の位の対称性と、1つとして同じ数がこないこと、偶数の場合奇数の周期がないので5巡でリセットされること。5または10を軸としての対称性、これらを一意に証明する方法はありますか? 100<λ<101 .
a_n = [λ*a_{n-1}]で、λは無理数、a_{n-1}は整数だから、
a_n < λ*a_{n-1} < 1+a_n . よって (a_n)/λ < a_{n-1} < 1/λ + (a_{n-1})/λ .
よって [ (a_n)/λ ] = a_{n-1}-1 .
よって a_{n+1}=[λ*a_n]=[(100+1/λ)*a_n]
=100*a_n+[(a_n)/λ]=100*a_n + a_{n-1}-1 .
よって mod 100 で a_{n+1}≡a_{n-1}-1 . >>344
なるほど
これなら高校生でも解けるかも >>318
>>327
百万以下の奇数は 50万個以下だから明らかぢゃね？
p(78498) = 999983,
p(78499) = 1000003, >>344
ホントだ
あってる
という事は計算精度がオーバーフローしたのか >>346
さすがに素数と合成数反対なんだと信じたい
それにしたって5までの篩で終わっちゃうけどw >>301
f <- function(N){
+ a=numeric()
+ a[1]=0
+ a[2]=0
+ # a[3]=(a[1]-1)%%100
+ # a[4]=(a[2]-1)%%100
+ for(n in 2:N){
+ a[n+1]=(a[n-1]-1)%%100
+ }
+ a[N]
+ }
> f(100)
[1] 51 >>346
ご指摘の通り。出願ミス。
1000万以下の素数の数は664579以下であることを示せの間違い。 > f(1:200)
[1] 0 0 99 99 98 98 97 97 96 96 95 95 94 94 93 93 92 92 91 91 90 90 89 89 88 88 87 87 86 86 85
[32] 85 84 84 83 83 82 82 81 81 80 80 79 79 78 78 77 77 76 76 75 75 74 74 73 73 72 72 71 71 70 70
[63] 69 69 68 68 67 67 66 66 65 65 64 64 63 63 62 62 61 61 60 60 59 59 58 58 57 57 56 56 55 55 54
[94] 54 53 53 52 52 51 51 50 50 49 49 48 48 47 47 46 46 45 45 44 44 43 43 42 42 41 41 40 40 39 39
[125] 38 38 37 37 36 36 35 35 34 34 33 33 32 32 31 31 30 30 29 29 28 28 27 27 26 26 25 25 24 24 23
[156] 23 22 22 21 21 20 20 19 19 18 18 17 17 16 16 15 15 14 14 13 13 12 12 11 11 10 10 9 9 8 8
[187] 7 7 6 6 5 5 4 4 3 3 2 2 1 1
> f(201:400)
[1] 0 0 99 99 98 98 97 97 96 96 95 95 94 94 93 93 92 92 91 91 90 90 89 89 88 88 87 87 86 86 85
[32] 85 84 84 83 83 82 82 81 81 80 80 79 79 78 78 77 77 76 76 75 75 74 74 73 73 72 72 71 71 70 70
[63] 69 69 68 68 67 67 66 66 65 65 64 64 63 63 62 62 61 61 60 60 59 59 58 58 57 57 56 56 55 55 54
[94] 54 53 53 52 52 51 51 50 50 49 49 48 48 47 47 46 46 45 45 44 44 43 43 42 42 41 41 40 40 39 39
[125] 38 38 37 37 36 36 35 35 34 34 33 33 32 32 31 31 30 30 29 29 28 28 27 27 26 26 25 25 24 24 23
[156] 23 22 22 21 21 20 20 19 19 18 18 17 17 16 16 15 15 14 14 13 13 12 12 11 11 10 10 9 9 8 8
[187] 7 7 6 6 5 5 4 4 3 3 2 2 1 1
200を周期に繰り返すみたいだな。 一般解は
fn <- function(n) (100-floor((n-1)/2))%%100
floorはガウス記号と同じ
%%100は100で割った剰余を返す
でよさげ。 プログラムキチガイが出て行かないのなら
新しいスレを建てるしかないな
手計算による高校数学質問スレ
みんなでこっちに移動すればよい
どうせプログラムキチガイも新しいスレに来るだろうけどね プログラムキチガイって小中学校の質問スレにもいるんだな
初めて知ったわ
まさに害悪 え？なに？
a_{n+1}≡a_{n-1}-1 (mod 100) が示されてなお
ループ回さないと正解にたどり着けなかったの？
即割り算じゃないの？ >>361
設定通りやるのがシミュレーションの基本。
1から100まで順に足して総和を求めよという問題はその手順を踏む。処理しきれなくなったら公式や定理等を利用する。 >>352
100での剰余だから周期100かと思ったら周期200なんだな。 >>347
>という事は計算精度がオーバーフローしたのか
Haskellでもそうなのか？
R言語には多倍数精度計算Rmpfrというパッケージがあるのだが、Wolframの結果と乖離していたから、信用していない。
漸化式 a_0=1, a_{n+1}=[λ* a_n ] (n=0,1,2,…) の計算に使ってみたけど精度不足なので諦めた。 >>329
信頼区間の算出法にはいろいろ流儀がある。
まあ、設問では事前確率分布を一様分布と設定してあるからベイズでやれという意味と解釈される。
エロ本ネタで種々の流儀で信頼区間を算出してグラフ化。
https://i.imgur.com/ijAwNUh.png
Bayesの事前確率分布はデフォルトのJefferey分布のまま。
ネタがネタなので一部にモザイクをかけたｗ >>365
数学の勉強なんかなんもした事ないのにそんな事わからんやろ？
お前の問題がいつも問題になってないと指摘されるのは、その「統計問題で許される事実上の既成事実化された暗黙の了解」を逸脱してるからやろ？
なんも勉強する事ない俺様ワールドで残りの人生一人で生きればいいやん？ >>365
理屈と膏薬はどんなところにもつく、という格言があるが、一番上の確率が負の値というのはどう理屈をつけるんだろ？ >>361
計算機に頼ってると、こんな簡単なこともわからなくなるということを示すいい例かと 私立医を必死にバカにするくせに自分はそれ未満なんだよなこいつ 次のスレでは、1にプログラムネタの書き込み禁止と明記しよう 前>>208
>>272
9000いくらのところに36差あるのをみつけた。
5000〜10000の中では最大。
だれか先に答えたはるからそれだと思う。 >>350
p(664579) = 9 999 991.
p(664580) = 10 000 019. >>322
>一行で完成するのに、なんで31以下の素数リストを列挙する手間をかける必要があるんだ、と思う。
人の手間より計算量で考えるべきでしょう。
√Nまでの素数をリストアップしてから計算するプログラムと
Nまでの素数リストを計算するプログラムとどちらが計算量が
多いかで比較すべき。 そもそもN以下の素数をリストアップするのに√N以下の素数のリストを利用するのはエラトステネスの篩そのもので別段新しくともなんともない話
しかし勉強不足のおばかちゃんは大発見と一人興奮してクソコード垂れ流す >>376
エラトステネスの篩と比べて計算量が多いか少ないかって話でしょ。
しらんけど。 >>377
もちろん変わらん
エラトステネスの篩そのまんま
なんも数学の勉強したことないアホが今まで発見されたことない新アルゴリズム見つけられるわけがない >>378
それは>>322のプログラムのことでしょ？ >>379
何が言いたいのか知らんがこんな便所の落書きで突然今まで知られてなかった素数リストアップの革命的方法なんか生まれるはずもなかろうに
ましてやあのブロおじじゃます無理
小学生でも知ってる話を“再発見”して喜んでるだけ
アホ丸出し これだけ色々言われたら普通の人ならもう書き込まないんだろうけど
害悪おじさんは普通じゃないからまた書き込むんだろうね >>380
なんか話がズレてるね。
誰かほかの人と勘違いしてんじゃないの？
素数のリストアップの方法の話なんかしてないんだけど。 >>383
そうだな
こんなアホな話いつまでもいつまでもダラダラ引きずってスレ汚すのは良くないな
あのアホと関わると人生の損以外にはならんからな
まともな人間が決して関わってはいかない奴 問題の意味は小学生にもわかる問題（解法は小学校の範囲を超える）
厚さが一定で20cm×10cmの楕円形ステーキを2:1に分割したい。
切断線の長さは最低にしたい、どこを切ればよいか？ >>250
30ずつのグループで考えたのはなぜですか >>386
2,3,5の剰余だとギリギリ足りないから(30の剰余で8個も素数候補があると8×33+10=274)生まれる発想じゃない？ >>387
20kだと1,7,11,13,17,19だからmod5で0-4まで揃えられないってことですか
この解答すごいですね イヤそれは30ごとではダメだから次どうしようの話やろ
30は素数小さい方から３つとって2×3×5=30だからというただそれだけ >>376
この1行コードは√N以下の素数のリストを利用していないよ。
(1:1000)[-outer(2:1000,2:1000)][-1]
あんたが、コードも読めないだけ。 >>385
短軸方向に切るという条件を外した問題
厚さが一定で20cm×10cmの楕円形ステーキを長軸に対して60°の角度で切って体積を2:1に分割したときの切断線の長さは？ >>390
outer関数って、色々な言語に使われていると思う。
Wolframにもあるし、pythonだとnumpy.outerとして同機能の関数がある。
知らなきゃ調べりゃいいのに。 相手に不愉快な思いをさせる事に喜びを感じる事が自分の人間性の恥ずべき部分だと認識できない
改めようともしない
改める事は自分の負けを意味するから
小学生の発想 >改める事は自分の負けを意味するから
これしっくり来たわ 普通に考えて
(1:1000)[-outer(2:1000,2:1000)][-1]
の1行プログラムで平方根をとっていると思う方がどうかしているよなぁ。
使っている関数はouterだけなのに。
sqrtの文字も入っていないのにどういう風に理解したのだろうね。
perlもCもpythonもHaskellも平方根計算はsqrt関数を使うはず。エクセルは大文字でSQRT。 もうほっとこう
構うといつまでも続けるし
そうやって他人に迷惑をかけてる自分の行為を止められない事を「自分が優秀で力がある」と認識するタイプなので構えば構うほどコイツにとっては喜び
他人を不愉快な気分にさせる事で自分の力を誇示してきた人生なんやろ
もちろん全てのレスは他人を不愉快にさせる事のためだけに書かれてる
ほっとくしかないよ >>392
１行なのだから読むも読まないもないだろうに。読んでも理解できなかったってことだろ。 >>398
> 普通に考えたら
普通じゃないからキチガイなんです >>377
１０００を超える合成数も計算させているから、計算量は多いよ。
outerを使ったコードも読めないのに
>378みたいに、もちろん変わらんというアホもいるが。
でも１００００程度なら１秒もかからず答が返ってくる。
やってみたら、0.6秒で10000までの素数をリストアップ。
> system.time((1:10000)[-outer(2:10000,2:10000)][-1])
user system elapsed
0.60 0.03 0.66
まあ、省スペースな分だけ計算量が多くてCPUに負荷はかかっていると思う。
10万にするとエラーがでた。
Error: cannot allocate vector of size 74.5 Gb 「pが素数ならp^4+14は素数ではないことを示せ」
という京大の問題だけどなんで5ではなくて14に問題設定したのだろう？
p^2+5でもいい気がするし。 pが奇数のとき偶数でアウト。p=2 のとき 9 でアウト。
それぢゃあ中学校の問題だろ。 >>407
2より大きい素数は奇数で、奇数の四乗が偶数にはならんでしょ。 よく>>406のレベルで数学板に書き込む気になるよな 前>>208
>>385
楕円型ステーキを圧縮して、
半径10の円型を2:1に分ける切断面の長さを考えると、
2θ-2sinθcosθ=2π/3
2θ-sin2θ=2π/3
20sinθ=16.66
θ=56.4068°
16.66ぐらいなんだけど、
妙に数字が並ぶところを見ると50/3なのかな？ ■ このスレッドは過去ログ倉庫に格納されています
