ログインしてさらにmixiを楽しもう

コメントを投稿して情報交換!
更新通知を受け取って、最新情報をゲット!

独学ノート(土筆の子)コミュの中国の剰余定理

  • mixiチェック
  • このエントリーをはてなブックマークに追加
m_1,m_2,...m_nを、どの2つをとっても互いに素
(最大公約数が1であること)な整数としたとき、
1次連立合同方程式

x≡b_1 (mod m_1) (1)
x≡b_2 (mod m_2) (2)
......
x≡b_n (mod m_n) (n)

は、積M=m_1m_2...m_nを法としてただ一つの解を
持つ。
-------------------------------------

3 〜 5 世紀ごろ成立したとされる中国の算術書 『孫子算経』 には 「3 で割れば 2 余り,5 で割れば 3 余り, 7 で割れば 2 余るような数は何か」 という問いと解法が書かれている。中国の剰余定理はこの問いを他の整数についても適用できるように一般化したもの。

x≡2(mod.3)、x≡3 (mod.5)、x≡2 (mod.7)
excelで105までの表を作って求めると答えは、23

定理の趣旨は、(3,5)=(5,7)=(3,7)=1のとき、
3*5*7=105の中に、一つの解をもつ。 とのこと。

証明:
n=2 のときを考える。
   (1)より、x=b_1+m_1y (y∈Ζ(整数))
これを(2):x≡b_2 (mod m_2)に代入し整理すると、
m_1y≡b_2-b_1 (mod m_2)
(m_1, m_2)=1 なので、この方程式は、法m_2でただ一つの解を持つ。
y≡y_0 (mod m_2) つまり、
   y=y_0+m_2q (q∈Ζ) ⇒ x=b_1+m_1(y_0+m_2q)
=b_1+m_1y_0+m_1m_2q
≡b_1+m_1y_0 (mod m_1m_2)

x=(b_1+m_1y_0)+m_1m_2q を、
x≡b_3 (mod m_3)に代入して整理
    m_1m_2q≡b_3-(b_1+m_1y_0) (mod m_3)
(m_1m_2, m_3)=1 なので、qは、法m_3でただ一つの解を持つ。
q≡y_1 (mod m_3) つまり、
q=y_1+m_3r (r∈Ζ)
⇒ x=b_1+m_1{y_0+m_2(y_1+m_3r)}
=b_1+m_1y_0+m_1m_2y_1+m_1m_2m_3r
≡b_1+m_1(y_0+m_2y_1) (mod m_1,m_2,m_3)

x=b_1+m_1y_0+m_1m_2y_1+m_1m_2m_3r
-----
x=b_1+m_1(y_0+m_2(y_1+m_3(y_2+m_4(.....)))) (mod m_1,m_2,...m_n)


別証(ガウスの方法)
M=m_1m_2...m_n, M_i=M/m_i とする
   仮定から、(m_i, M_i)=1 (i=1,2,...n)がわかる。
M_if_i≡1 (mod m_i)
 を満たす整数f_i(i=1,2,...,n) が求められる。
   t_i=M_if_iとすると、
   M_iはm_j (j≠i)の倍数なので、
   t_i≡1 (mod m_i), t_i≡0 (mod m_j) (i≠j)
となることがわかる。これから、
   
   t_1b_1+t_2b_2+...+t_nb_n≡b_i (mod m_i) (i=1,...,n)
なので、
   x=t_1b_1+t_2b_2+...+t_nb_n が一つの解。

   もし、
   y≡b_1 (mod m_1) (1)'
y≡b_2 (mod m_2)  (2)'   
......
y≡b_n (mod m_n) (3)'

であれば、対応する式を引いて
   x≡y (mod m_1) (1)''
x≡y (mod m_2)  (2)''   
......
x≡y (mod m_n) (3)''
を得る。これは、x-yが互いに素な m_1,m_2,...m_nの
   倍数を意味する。
   したがって、x-yはそれらの最小公倍数Mで割れるので、
   x≡y (mod M) でMを法として一致する。

例:3で割れば2余り、5で割れば3余り、7で割れば5が余る数を求めよ。

解答
x=b_1+m_1(y_0+m_2(y_1+m_3(y_2+m_4(.....)))) (mod m_1,m_2,...m_n)
=2+3*(y_0+5*(y_1+7*r)) (mod 3*5*7)
下の式からy_0=y=2, y_1=z=4
x=2+3(2+5*(4+7r))=2+6+3*20+3*5*7r=68+105r

-----------------------------------

x≡2 (mod 3)
   x=2+m_1y=2+3y
x≡3 (mod 5) 2+3y≡3
  3y ≡3-2 =1 (mod 5) つまり 3y ≡6 (mod 5)より
y≡2 (mod 5)
y=2+5z
x≡5 (mod 7)
x=2+m_1y=2+3y=2+3(2+5z )=8+15z≡5 (mod 7)   
15z≡-3 ≡4 (mod 7) つまり、14z+z≡4 (mod 7)より
 z≡4(mod 7)
 z=4+m_3r (r∈Ζ) =4+7r
-----------------------------------------

別解:
方程式で表すと、x≡2 (mod 3)
x≡3 (mod 5) 
x≡5 (mod 7)
となる。定理の証明の記号を用いて、
  M=3*5*7=105, M_1=5*7=35, M_2=3*7=21, M_3=3*5=15
したがって、
  35f_1=1 (mod 3) ⇒ f_1=2 ⇒ t_1=70
21f_2=1 (mod 5) ⇒ f_2=1 ⇒ t_2=21
15f_3=1 (mod 7) ⇒ f_3=1 ⇒ t_3=15
これから
  x=t_1b_1+t_2b_2+t_3b_3
  x=70*2+21*3+15*5=1400+63+75=278 (mod 105)=68 (mod 105)
⇒ x=68+105k (k∈Ζ) 

参照:草場公邦著「ガロアと方程式」

コメント(5)

ガウスの方法のおもしろいところは、
成り立つような解を想定し、それが唯一であることをあとで証明する
ところ。これを演繹的と言うのかもしれない。

M_if_i≡1 (mod m_i)
 を満たす整数f_i(i=1,2,...,n) が求められる。

ところがずーっと判らなかったところだが、
これは、フェルマーの小定理 a^(p-1)≡1 (mod p)に似ていて、
a^p≡a (mod p)になることが次に来る。

RAS暗号でも、
M^(n(p-1)(q-1)+1)=M^(ED)=(M^E)^D mod[n(=pq)]
公開鍵、n, E 秘密鍵、p,q, (D)
で合同式の基本的巡回の性質が使われている。







Table[ChineseRemainder[{a, b, c}, {3, 5, 7}], {a, 0, 2, 1}, {b, 0, 4,
1}, {c, 0, 6, 1}]

{{0, 15, 30, 45, 60, 75, 90},
{21, 36, 51, 66, 81, 96, 6},
{42, 57, 72, 87, 102, 12, 27},
{63, 78, 93, 3, 18, 33, 48},
{84, 99, 9, 24, 39, 54, 69}},

{{70, 85, 100, 10, 25, 40, 55},
{91, 1, 16, 31, 46, 61, 76},
{7, 22, 37, 52, 67, 82, 97},
{28, 43, 58, 73, 88, 103,13},
{49, 64, 79, 94, 4, 19, 34}},

{{35, 50, 65, 80, 95, 5, 20},
{56, 71, 86, 101, 11, 26, 41},
{77, 92, 2, 17, 32, 47, 62},
{98, 8, 23, 38, 53, 68, 83},
{14, 29, 44, 59, 74, 89, 104}}}
Mathematicaで、105算をすべて求めてみた。
たとえば、3,5,7で割って2,4,6余る最少の数は104
この、中国剰余のガウスの回答は、基底(イデアルとでもいうのであろうか)
を取っていくことが、次のラグランジの補間式に似ている。

http://mathworld.wolfram.com/LagrangeInterpolatingPolynomial.html

それにしても、昨晩は、何十年かでやっとわかったすっきり感に
感動した。
http://wankora.blog31.fc2.com/blog-entry-1790.html

このページが面白い。
中国剰余の定理、ガウスの方法は、
抽象代数の、イデアルに導いてくれるようだ。

ログインすると、みんなのコメントがもっと見れるよ

mixiユーザー
ログインしてコメントしよう!

独学ノート(土筆の子) 更新情報

独学ノート(土筆の子)のメンバーはこんなコミュニティにも参加しています

星印の数は、共通して参加しているメンバーが多いほど増えます。