メインコンテンツへスキップ
教材一覧に戻る
数検1級

数値解析・初等整数論 — ニュートン法・合同式と二次対策

最終確認日:

この章で学ぶこと

仕上げとして、数値解析と初等整数論の基礎、そして二次(数理技能)で得点するための戦略を学びます。

  • ニュートン法による方程式の近似解
  • ユークリッドの互除法と最大公約数
  • 合同式と中国剰余定理
  • 二次(数理技能)対策の進め方

ポイント: 数値解析は「解析的に解けない問題を反復計算で近づける」技法、整数論は「割り算の余り」を主役にした分野です。どちらも二次の選択問題で得点源になりやすい実用的なテーマです。

1. ニュートン法

方程式 f(x)=0f(x) = 0 の解を、初期値 x0x_0 から出発して接線で近づける反復法がニュートン法です。

xn+1=xnf(xn)f(xn)x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)}

例題: f(x)=x22f(x) = x^2 - 2 にニュートン法を適用し、x0=1.5x_0 = 1.5 から 2\sqrt{2} の近似値を1ステップ求めよ。

f(x)=2xf'(x) = 2xx0=1.5x_0 = 1.5f(1.5)=2.252=0.25f(1.5) = 2.25 - 2 = 0.25f(1.5)=3f'(1.5) = 3

x1=1.50.253=1.50.0833=1.4166x_1 = 1.5 - \frac{0.25}{3} = 1.5 - 0.0833\ldots = 1.4166\ldots

検算: 真の値 2=1.41421\sqrt{2} = 1.41421\ldots。1ステップで 1.51.5 から 1.41671.4167 へ大きく近づいた(x12=1.416722.007x_1^2 = 1.4167^2 \approx 2.00722 に近い)。ニュートン法は2次収束するので、さらに数ステップで高精度になる。正しい。

大事: ニュートン法の更新式は「現在の点での接線が xx 軸と交わる点」を次の近似値とするものです。f(xn)=0f'(x_n) = 0 になる点の近くでは破綻するので注意が必要です。

2. ユークリッドの互除法

2つの整数の最大公約数 gcd\gcd は、「大きい数を小さい数で割った余りに置き換える」操作を繰り返すユークリッドの互除法で求まります。

例題: gcd(252,105)\gcd(252, 105) を求めよ。

252=105×2+42252 = 105 \times 2 + 42 105=42×2+21105 = 42 \times 2 + 21 42=21×2+042 = 21 \times 2 + 0

余りが 00 になる直前の除数 2121 が最大公約数。よって gcd(252,105)=21\gcd(252, 105) = 21

検算: 252=21×12252 = 21 \times 12105=21×5105 = 21 \times 5121255 は互いに素なので、2121 が最大公約数で正しい。

ポイント: 互除法を逆にたどると、gcd(a,b)=ax+by\gcd(a,b) = ax + by となる整数 x,yx, y が求まります(ベズーの等式)。これは1次不定方程式 ax+by=cax + by = c の解法や、合同式での逆数計算に使えます。

3. 合同式

整数 a,ba, bmm で割った余りが等しいとき ab(modm)a \equiv b \pmod{m} と書き、合同といいます。合同式はたし算・ひき算・かけ算について通常の等式のように計算できます。

例題: 71007^{100}55 で割った余りを求めよ。

72(mod5)7 \equiv 2 \pmod 5 なので 71002100(mod5)7^{100} \equiv 2^{100} \pmod 522 のべきの余りは周期 44 で循環する(212, 224, 233, 2412^1\equiv2,\ 2^2\equiv4,\ 2^3\equiv3,\ 2^4\equiv1)。100=4×25100 = 4 \times 25 なので

2100=(24)25125=1(mod5)2^{100} = (2^4)^{25} \equiv 1^{25} = 1 \pmod 5

よって余りは 11

検算: 24=161(mod5)2^4 = 16 \equiv 1 \pmod 52525 乗しても 11。フェルマーの小定理(pp が素数で aapp の倍数でないとき ap11(modp)a^{p-1} \equiv 1 \pmod p)からも、241(mod5)2^{4} \equiv 1 \pmod 5 なので 2100=(24)2512^{100} = (2^4)^{25}\equiv1。一致する。

大事: 中国剰余定理は、互いに素な m1,m2,m_1, m_2, \ldots について連立合同式 xa1(modm1), xa2(modm2),x \equiv a_1 \pmod{m_1},\ x \equiv a_2 \pmod{m_2},\ldots がただ1つの解を (modm1m2)\pmod{m_1 m_2 \cdots} でもつことを保証します。たとえば「33 で割ると 22 余り、55 で割ると 33 余る数」は (mod15)\pmod{15}x8x \equiv 8 と一意に定まります(8=32+2, 8=51+38 = 3\cdot2+2,\ 8 = 5\cdot1+3)。

4. 二次(数理技能)対策の進め方

1級の二次は「必須2問 + 選択5問から2問」の記述式です。得点するためのコツを整理します。

  • 得意分野で選択問題を選ぶ: 解析(微積分・微分方程式)、線形代数(固有値・対角化)、確率統計のうち、自信のある分野を選ぶ。
  • 定理の仮定を確認してから使う: 「連続だから最大値をとる」「微分可能だから平均値の定理が使える」など、根拠を明示すると記述の評価が高い。
  • 計算は必ず検算する: 固有値なら「和=トレース・積=行列式」、積分なら「微分して戻す」など、独立な方法で確かめる。
  • 途中式をていねいに書く: 二次は部分点が大きい。最終答だけでなく、立式・変形の流れを論理的に書く。

大事: 1級は範囲が広いので、全分野を浅く より 数分野を確実に が戦略です。本教材なら「解析(第2〜5章)+ 線形代数(第7〜8章)」を主軸に固め、確率統計・複素解析・整数論から得意なものを選択問題用に1〜2分野準備するとよいでしょう。

どう問われるか

  • 一次では「ニュートン法の1〜2ステップ」「最大公約数」「合同式によるべき乗の余り」などの計算が出ます。
  • 二次では「漸化式・近似計算」「不定方程式・合同式の応用」「これまでの全分野の総合問題」が記述式で問われます。

よくまちがえるところ

数値解析と整数論は、誤差の種類合同式で許される操作の理解が勝負どころです。

まちがい直し方
打ち切り誤差と丸め誤差を混同する桁数を増やせば打ち切り誤差も消えると考える打ち切りは有限項・有限回で止めたことによる誤差、丸めは有限桁でしか表せないことによる誤差。原因が別
ニュートン法で導関数が 0 になる点を避けないf(x) = x³ - x に初期値 1/√3 を使うそこで導関数が 0 になり計算が破綻する。初期値を変える
近い数どうしの引き算(桁落ち)を放置する大きい x で √(x+1) - √x をそのまま計算する有理化して 1/(√(x+1) + √x) に直してから計算する
合同式で自由に割り算する2x ≡ 2 (mod 4) から x ≡ 1 と結論する法と互いに素な数でしか割れない。この解は x ≡ 1 と x ≡ 3 の2つ

いちばん多いのは誤差の混同です。打ち切り誤差は、無限級数を有限項で切ったり、反復計算を有限回で止めたりすることから生じます。項数や反復回数を増やせば小さくできます。丸め誤差は、計算機が有限桁でしか数を表せないことから生じ、こちらは桁数(精度)を上げないと小さくなりません。打ち切り誤差は計算の設計で、丸め誤差は計算機の表現能力で決まる、と原因のちがいで整理してください。

合同式の割り算も要注意です。a ≡ b (mod m) の両辺を c で割れるのは、c と m が互いに素なときに限ります。上の例では c = 2 と m = 4 の最大公約数が 2 なので割れません。互いに素でない場合は、法のほうも割って 2x ≡ 2 (mod 4) を x ≡ 1 (mod 2) と読みかえます。

誤差の種類を使い分ける

誤差の名前生じる原因減らし方
打ち切り誤差無限の操作を有限で止めるテイラー展開を3項で切る項数・反復回数を増やす
丸め誤差有限桁でしか表せない1/3 を 0.3333 と表す精度(桁数)を上げる
桁落ち近い数どうしを引く大きい x での √(x+1) - √x式を変形して引き算を避ける
情報落ち大きさの違う数を足す非常に大きい数に非常に小さい数を足す小さい数から順に足す

ニュートン法の性質も整理しておきます。この方法は、現在の点での接線が x 軸と交わる点を次の近似値とするもので、更新式は x(次) = x(今) - f(x(今))/f'(x(今)) です。解の近くでは2次収束するため、正しい桁数がおよそ倍々に増えていきます。

段階x の値2乗した値
初期値1.52.25
1回目1.41666...約 2.0069
2回目1.414215...約 2.0000006

f(x) = x² - 2 の場合、2回の反復で √2 = 1.4142135... にきわめて近づきます。ただし導関数が 0 に近い点では更新の分母が小さくなって発散するので、初期値の選び方が重要です。

大事: 合同式でべき乗の余りを求めるときは、周期を見つけるのが早道です。法 5 における 2 のべきは 2、4、3、1 と周期 4 で循環するので、指数を 4 で割った余りだけを見ればよいことになります。フェルマーの小定理(p が素数で a が p の倍数でないとき a の (p-1) 乗が 1 と合同)を使うと、この周期の上限がすぐ分かります。

自分でチェック

  • 打ち切り誤差と丸め誤差を原因のちがいで説明できる
  • 桁落ちを避けるために式を変形できる
  • ニュートン法の更新式を書き、破綻する条件を言える
  • 合同式で割り算が許される条件を言える
  • べき乗の余りを周期やフェルマーの小定理で求められる

まとめ

  • ニュートン法 xn+1=xnf(xn)f(xn)x_{n+1} = x_n - \dfrac{f(x_n)}{f'(x_n)} は接線で近づける反復法
  • ユークリッドの互除法で gcd\gcd、逆にたどってベズーの等式
  • 合同式はべき乗の余りの周期性・フェルマーの小定理・中国剰余定理
  • 二次は得意分野を選び、仮定の確認・検算・途中式を大切に

これで数検1級の教材はひと通り終わりです。解析と線形代数を主軸に、各章の例題を検算しながら繰り返し解き、一問一答や問題集で計算力を確かめましょう。

※「数検」「実用数学技能検定」は公益財団法人日本数学検定協会の登録商標です。本教材は非公式の学習教材であり、合格を保証するものではありません。

この教材きょうざいやくちましたか?

このしょうれんしゅう

ほかのぶんれんしゅうもチェック