Теорема о делении с остатком. Алгоритм Евклида — различия между версиями

Текущая версия на 13:07, 24 июня 2015

Напомним теорему о делении с остатком:

Теорема о делении с остатком

Для любых целых $a$ и $b$ , $b \not= 0$ , существует единственный набор целых чисел $q$ и $r$ , что $a = bq + r$ и $0 \le r < |b|$ , где $|b|$ — модуль числа $b$ .

Легко доказывается, что для любых целых чисел $a$ и $b$ , $b \not= 0$ деление с остатком возможно и числа $q$ и $r$ определяются однозначно.

Рассмотрим следующий пример:

Пример

Пусть модуль $p = 6$ .

Тогда имеем шесть классов разбиения множества целых чисел по модулю 6:

$K_0 = \left\{\ldots,-18,-12,-6,0,6,12,18,\ldots \right\}, r=0$ ;

$K_1 = \left\{\ldots,-17,-11,-5,1,7,13,19,\ldots \right\}, r=1$ ;

$K_2 = \left\{\ldots,-16,-10,-4,2,8,14,20,\ldots \right\}, r=2$ ;

$K_3 = \left\{\ldots,-15,-9,-3,3,9,15,21,\ldots \right\}, r=3$ ;

$K_4 = \left\{\ldots,-14,-8,-2,4,10,16,22,\ldots \right\}, r=4$ ;

$K_5 = \left\{\ldots,-13,-7,-1,5,11,17,23,\ldots \right\}, r=5$ ,

где через $r$ обозначен остаток от деления целого числа на 6.

В данном примере полная система наименьших неотрицательных вычетов есть множество $\left\{0, 1, 2, 3, 4, 5 \right\}$ ;

полная система наименьших положительных вычетов – множество $\left\{1, 2, 3, 4, 5, 6 \right\}$ ;

полная система наименьших по абсолютной величине вычетов – множество $\left\{-2,-1, 0, 1, 2, 3 \right\}$ ;

приведённая система вычетов – множество $\left\{1,5 \right\}$ , так как $\varphi(6) = 6\,\left(1-\frac{1}{2}\right)\left(1-\frac{1}{3}\right) = 2$ ;

фактор-множество $Z{|{\equiv}_6} = \left\{K_0, K_1, K_2, K_3, K_4, K_5 \right\}$ .

Один из методов выполнения арифметических операций над целыми числами основан на простых положениях теории чисел. Идея этого метода состоит в том, что целые числа представляются в одной из непозиционных систем – в системе остаточных классов. А именно: вместо операций над целыми числами оперируют с остатками от деления этих чисел на заранее выбранные числа – модули $p_1, p_2, \ldots, p_n$ . Чаще всего модули $p_1, p_2, \ldots, p_n$ выбирают из множества простых чисел.

Пусть

$A \equiv {\alpha}_1 \pmod{p_1}$ ,

$A \equiv {\alpha}_2 \pmod{p_2}$ ,

$\ldots$ ,

$A \equiv {\alpha}_n \pmod{p_n}$ .

Так как в кольце целых чисел имеет место теорема о делении с остатком, т. е. $\forall a \in Z, \forall b \in Z, b\not= 0 \, \exist q \in Z, \exist r \in Z (0 \le r <b):\, a=bq+r$ , то кольцо Z, по определению, является евклидовым. Таким образом, в качестве чисел ${\alpha}_i, i=0,\ldots,n$ можно выбрать остатки от деления числа $A$ на $p_i$ соответственно.

Рассмотрим гомоморфное отображение:

$Z \to Z{|{\equiv}_{p_1}} \times Z{|{\equiv}_{p_2}} \times \ldots \times Z{|{\equiv}_{p_n}}$ .

Тогда каждому целому числу $A$ можно поставить в соответствие кортеж $\left\{{\alpha}_1, {\alpha}_2, \ldots, {\alpha}_n \right\}$ наименьших неотрицательных вычетов по одному из соответствующих классов.

Важно отметить, что при преобразовании числа нет потери информации, если выполнено условие $A < p_1 \cdot p_2 \cdot \ldots \cdot p_n$ , поскольку всегда, зная $\left({\alpha}_1, {\alpha}_2, \ldots, {\alpha}_n \right)$ можно восстановить само число $A$ . Поэтому кортеж можно рассматривать как один из способов представления целого числа $A$ – модулярное представление, или представление в системе остаточных классов (СОК).

Для дальнейшего используем расширенный алгоритм Евклида или его аналог – алгоритм нахождения линейного представления наибольшего общего делителя целых чисел: если числа а и b одновременно не равны нулю, то существуют целые числа $x$ и $y$ , такие, что $a \cdot b = a \cdot x + b \cdot y$ .

Действительно, пусть $d$ – наименьшее целое положительное число вида $a \cdot x + b \cdot y$ , например, $d = a \cdot x_0 + b \cdot y_0$ , где числа $x_0$ и $y_0$ не обязательно определены однозначно. Существование числа $d$ следует только из принципа полной упорядоченности. Очевидно, что $d > 0$ . Остаётся показать, что $d = (a,b)$ . Для этого надо проверить выполнение двух условий: а) $d|a$ и $d|b$ б) если $c|a$ и $c|b$ то $c|d$ .

От противного: допустим, что свойство а) не выполняется, для определенности положим, что не выполнено $d|b$ . Тогда по теореме о делении с остатком $a = dq + r$ , $0 \le r < |d|$ , и, следовательно,

$r=b-dq = b-(ax_0+by_0)q = a(-qx_0)+b(1-qy_0)$ ,

что противоречит минимальности $d$ .

Выполнение свойства б) проверяется непосредственно:

$c|a \Rightarrow \left(\exist q_1 \in Z\right)\left(a = q_1\, c\right)$ ;

$c|b \Rightarrow \left(\exist q_2 \in Z\right)\left(a = q_2\, c\right)$ ;

$d=ax_0 + by_0 \Rightarrow d = q_1 c x_0 + q_2 c y_0 \Rightarrow d = \left(q_1 x_0 + q_2 y_0\right) c|d$ .

Алгоритм Евклида для целых чисел

Пусть $a$ и $b$ — целые числа, не равные одновременно нулю, и последовательность чисел

$a > b > r_1 > r_2 > r_3 > r_4 > \cdots >r_n$

определена тем, что каждое $r_k$ — это остаток от деления предпредыдущего числа на предыдущее, а предпоследнее делится на последнее нацело, то есть

$a = bq_0 + r_1$

$b = r_1q_1 + r_2$

$r_1 = r_2q_2 + r_3$

$\cdots$

$r_{k-2} = r_{k-1} q_{k-1} + r_k$

$\cdots$

$r_{n-2} = r_{n-1}q_{n-1}+ r_n$

$r_{n-1} = r_n q_n$

Тогда НОД $(a,b)$ , наибольший общий делитель $a$ и $b$ , равен $r_n$ , последнему ненулевому члену этой последовательности.

Существование таких $r_1, r_2, ...$ , то есть возможность деления с остатком $a$ на $b$ для любого целого $a$ и целого $b\ne 0$ , доказывается индукцией.

Корректность этого алгоритма вытекает из следующих двух утверждений:

Пусть $a = bq + r$ , тогда НОД $(a,b)$ = НОД $(b,r)$ .
НОД $(0,r)$ = $r$ для любого ненулевого $r$ (так как 0 делится на любое целое число, кроме нуля).

Расширенный алгоритм Евклида для целых чисел

Рассмотрим расширенный алгоритм Евклида для нахождения линейного представления наибольшего общего делителя $(a,b) = ax + by$ .

Значения $x$ и $y$ вычисляются в серии шагов, в каждом из которых мы выражаем $a_i$ (вычисленное в процессе работы алгоритма Евклида) в форме $a x_i + b y_i$ . А именно, рассмотрим последовательность

$a_0 = a$ , $a_0 = a x_0 + b y_0$ ,

$a_1 = b$ , $a_1 = a x_1 + b y_1$ ,

$a_2 = a_0 - a_1 q_1$ , $a_2 = a x_2 + b y_2$ ,

$a_3 = a_1 - a_2 q_2$ , $a_3 = a x_3 + b y_3$ ,

$\cdots$

$a_i = a_{i-2} - a_{i-1} q_{i-1}$ , $a_i = a x_i + b y_i$ ,

$\cdots$

$a_k = a_{k-2} - a_{k-1} q_{k-1}$ , $a_k = a x_k + b y_k$ ,

$0 = a_{k-1} - a_k q_k$ , $0 = a x_{k+1} + b y_{k+1}$ ,

В левом столбце алгоритма записана последовательность делений, которая получается в результате работы алгоритма Евклида и которая разрешена относительно остатков. Согласно теореме Ламе (1844 г.) число делений, которое необходимо выполнить для нахождения $(a,b)$ , не превосходит числа цифр в меньшем из чисел $a$ и $b$ , умноженного на 5 (оценка наихудшего случая для алгоритма Евклида). Теорема Ламе доказывается на основе последовательности Фибоначчи.

В правом столбце алгоритма каждый остаток выражен через $a x_i + b y_i$ . Надо вычислить $x_i$ и $y_i$ . Очевидно, что $x_0 = 1, y_0 = 0$ и $x_1 = 0, y_1 = 1$ . Сравнивая обе части на i-м шаге, получим

$a_i = a x_i + b y_i = a_{i-2} - a_{i-1} q_{i-1} = \left(a x_{i-2} + b y_{i-2} \right) - \left(a x_{i-1} + b y_{i-1} \right) q_{i-1} =$

$\left(a x_{i-2} - x_{i-1} q_{i-1} \right) + \left(b y_{i-2} - y_{i-1} q_{i-1} \right)$ ,

откуда получается следующая индуктивная процедура вычисления $x_i$ и $y_i$ :

$q_{i-1} = \left \lfloor \frac{a_{i-2}}{a_{i-1}} \right \rfloor$ ,

$a_i = a_{i-2} - a_{i-1} q_{i-1}$ ,

$x_i = x_{i-2} - x_{i-1} q_{i-1}$ ,

$y_i = y_{i-2} - y_{i-1} q_{i-1}$ .

Пример (Расширенный алгоритм Евклида).

Применим расширенный алгоритм Евклида к числам $a = 342, b = 612$ .

Весь алгоритм представим в виде следующей таблицы.

$\begin{array}{|c|r|r|r|r|r|r|r|} \hline Iteration & q & a_0 & a_1 & x_0 & x_1 & y_0 & y_1 \\ \hline 0 & - & 342 & 612 & 1 & 0 & 0 & 1 \\ \hline 1 & 0 & 612 & 342 & 0 & 1 & 1 & 0 \\ \hline 2 & 1 & 342 & 270 & 1 & -1 & 0 & 1 \\ \hline 3 & 1 & 270 & 72 & -1 & 2 & 1 & -1 \\ \hline 4 & 3 & 72 & 54 & 2 & -7 & -1 & 4 \\ \hline 5 & 1 & 54 & 18 & -7 & 9 & 4 & -5 \\ \hline 6 & 3 & 18 & 0 & 9 & -34 & -5 & 19 \\ \hline \end{array}$

Заметим, что равенство $a_0 = a x_0 + b y_0$ выполняется на каждом шаге итерации. Алгоритм выдаёт

$d = 18, x = 9, y = -5$ ,

и тогда

$18=342 \cdot 9 + 612 \cdot (-5)$ .

@@ Строка 30: / Строка 30: @@
 где через <math>r</math> обозначен остаток от деления целого числа на 6.
-В данном примере полная система наименьших неотрицательных вычетов есть множество <math>\left\{0, 1, 2, 3, 4, 5 \right\}</math>; полная система наименьших положительных вычетов – множество <math>\left\{0, 1, 2, 3, 4, 5, 6 \right\}</math>; полная система наименьших по абсолютной величине вычетов – множество <math>\left\{-2,-1, 0, 1, 2, 3 \right\}</math>; приведённая система вычетов – множество <math>\left\{1,5 \right\}</math>, так как <math>\phi(6) = 6\,\left(1-\frac{1}{2}\right)\left(1-\frac{1}{3}\right) = 2</math> ; фактор-множество <math>Z{|{\equiv}_6} = \left\{K_0, K_1, K_2, K_3, K_4, K_5 \right\}</math>.
+В данном примере полная система наименьших неотрицательных вычетов есть множество <math>\left\{0, 1, 2, 3, 4, 5 \right\}</math>;
+полная система наименьших положительных вычетов – множество <math>\left\{1, 2, 3, 4, 5, 6 \right\}</math>;
+полная система наименьших по абсолютной величине вычетов – множество <math>\left\{-2,-1, 0, 1, 2, 3 \right\}</math>;
+приведённая система вычетов – множество <math>\left\{1,5 \right\}</math>, так как <math>\varphi(6) = 6\,\left(1-\frac{1}{2}\right)\left(1-\frac{1}{3}\right) = 2</math>;
+фактор-множество <math>Z{|{\equiv}_6} = \left\{K_0, K_1, K_2, K_3, K_4, K_5 \right\}</math>.
@@ Строка 49: / Строка 57: @@
 :<math>Z \to Z{|{\equiv}_{p_1}} \times Z{|{\equiv}_{p_2}} \times \ldots \times Z{|{\equiv}_{p_n}} </math> .
-Тогда каждому целому числу А можно поставить в соответствие кортеж <math>\left\{{\alpha}_1, {\alpha}_2, \ldots, {\alpha}_n \right\}</math> наименьших неотрицательных вычетов по одному из соответствующих классов.
+Тогда каждому целому числу <math>A</math> можно поставить в соответствие кортеж <math>\left\{{\alpha}_1, {\alpha}_2, \ldots, {\alpha}_n \right\}</math> наименьших неотрицательных вычетов по одному из соответствующих классов.
 Важно отметить, что при преобразовании числа нет потери информации, если выполнено условие <math>A < p_1 \cdot p_2 \cdot \ldots \cdot p_n</math>, поскольку всегда, зная <math>\left({\alpha}_1, {\alpha}_2, \ldots, {\alpha}_n \right)</math> можно восстановить само число <math>A</math>. Поэтому кортеж  можно рассматривать как один из способов представления целого числа <math>A</math> – модулярное представление, или представление в системе остаточных классов (СОК).
@@ Строка 88: / Строка 96: @@
 : <math>r_{n-2} = r_{n-1}q_{n-1}+ r_n</math>
 : <math>r_{n-1} = r_n q_n</math>
 Тогда '''НОД''' <math>(a,b)</math>, наибольший общий делитель <math>a</math> и <math>b</math>, равен
@@ Строка 130: / Строка 139: @@
 В правом столбце алгоритма каждый остаток выражен через <math>a x_i + b y_i</math>. Надо вычислить <math>x_i</math> и <math>y_i</math>. Очевидно, что <math>x_0 = 1, y_0 = 0</math> и <math>x_1 = 0, y_1 = 1</math>. Сравнивая обе части на ''i''-м шаге, получим
-:<math>a_i = a x_i + b y_i = a_{i-2} - a_{i-1} q_{i-1} = \left(a x_{i-2} + b y_{i-2} \right) - \left(a x_{i-1} + b y_{i-1} \right) q_{i-1} = \left(a x_{i-2} - x_{i-1} q_{i-1} \right) + \left(b y_{i-2} - y_{i-1} q_{i-1} \right)</math>,
+:<math>a_i = a x_i + b y_i = a_{i-2} - a_{i-1} q_{i-1} = \left(a x_{i-2} + b y_{i-2} \right) - \left(a x_{i-1} + b y_{i-1} \right) q_{i-1} =</math>
+:<math> \left(a x_{i-2} - x_{i-1} q_{i-1} \right) + \left(b y_{i-2} - y_{i-1} q_{i-1} \right)</math>,
 откуда получается следующая индуктивная процедура вычисления <math>x_i</math> и <math>y_i</math>:
@@ Строка 150: / Строка 161: @@
-<math> \begin{array}{|c|r|r|r|r|r|r|r|} \hline Iteration & q & A_0 & a_1 & x_0 & x_1 & Y_0 & y_1
+<math> \begin{array}{|c|r|r|r|r|r|r|r|} \hline Iteration & q & a_0 & a_1 & x_0 & x_1 & y_0 & y_1
 \\ \hline 0 & - & 342 & 612 &  1 &   0 &  0 &  1
 \\ \hline 1 & 0 & 612 & 342 &  0 &   1 &  1 &  0
@@ Строка 162: / Строка 173: @@
-Заметим, что равенство <math>a_0 = a x_0 + b y_0</math> выполняется на каждом шаге итерации. Алгоритм выдаёт <math>d = 18, x = 9, y = -5</math> и тогда <math>18=342 \cdot 9 + 612 \cdot (-5)</math>.
+Заметим, что равенство <math>a_0 = a x_0 + b y_0</math> выполняется на каждом шаге итерации. Алгоритм выдаёт
+:<math>d = 18, x = 9, y = -5</math>,
+и тогда
+:<math>18=342 \cdot 9 + 612 \cdot (-5)</math>.

Теорема о делении с остатком. Алгоритм Евклида — различия между версиями

Текущая версия на 13:07, 24 июня 2015

Навигация

Персональные инструменты

Пространства имён

Варианты

Просмотры

Действия

Поиск

Навигация

Инструменты