Обратное по модулю число

Обратное по модулю целого a — это такое целое число x, что произведение ax сравнимо с 1 по модулю m^[1]. В стандартных обозначениях модульной арифметики эта эквивалентность записывается как:

ax\equiv 1{\pmod {m}},

что является сокращённым способом записи утверждения, что m делит (без остатка) величину ax − 1, или, выражаясь другим способом, остаток от деления ax на целое m равен 1. Если a имеет обратный по модулю m, то имеется бесконечное количество решений этой эквивалентности, которые образуют класс вычетов для этого модуля. Более того, любое целое, которое эквивалентно a (то есть из класса эквивалентности a) будет иметь любой элемент класса эквивалентности x в качестве обратного элемента. Используя обозначения ${\overline {w}}$ для класса эквивалентности, содержащего w, утверждение выше может быть записано следующим образом: обратный элемент по модулю класса эквивалентности ${\overline {a}}$ есть класс эквивалентности ${\overline {x}}$ , такой что

{\overline {a}}\cdot _{m}{\overline {x}}={\overline {1}},

где символ $\cdot _{m}$ означает умножение классов эквивалентности по модулю m^[2]. Такой вид записи представляет аналог обычной концепции обратного числа в множестве рациональных или вещественных чисел, если заменить числа классами эквивалентности и должным образом определения бинарных операций.

Фундаментальное использование этой операции — решение линейной эквивалентности вида:

ax\equiv b{\pmod {m}}.

Нахождение модульного обратного имеет практическое приложение в области криптографии, например, криптосистема с открытым ключом и алгоритм RSA ^[3]^[4]^[5]. Преимущество для реализации этих приложений в том, что существует очень быстрый алгоритм (расширенный алгоритм Евклида), который может быть использован для вычисления модульных обратных.

Модульная арифметика

Говорят, что для данного положительного числа m два целых числа a и b сравнимы по модулю m, если m делит их разность. Это бинарное отношение обозначается как

a\equiv b{\pmod {m}}.

Это является отношением эквивалентности на множестве целых чисел $\mathbb {Z}$ , и классы эквивалентности называются классами вычетов по модулю m. Пусть ${\overline {a}}$ означает класс эквивалентности, содержащий целое число a^[6], тогда

{\overline {a}}=\{b\in \mathbb {Z} \mid a\equiv b{\pmod {m}}\}.

Линейное сравнение — это сравнение по модулю вида

ax\equiv b{\pmod {m}}.

В отличие от линейных уравнений в целых числах, линейное сравнение может иметь нуль, одно или несколько решений. Если x является решением линейного сравнения, то любой элемент из ${\overline {x}}$ также является решением, так что, если говорится о количестве решений линейного сравнения, подразумевается количество различных классов эквивалентности, которые содержат эти решения.

Если d является наибольшим общим делителем чисел a и m, то линейное сравнение $ax\equiv b{\pmod {m}}$ имеет решения тогда и только тогда, когда d делит b. Если d делит b, то существует ровно d решений^[7].

Модульное обратное целого числа a по модулю m — это решение линейного сравнения

ax\equiv 1{\pmod {m}}.

Ранее говорилось, что решение существует тогда и только тогда, когда НОД (наибольший общий делитель) (a, m) = 1, то есть a и m должны быть взаимно простыми числами. Более того, если это условие выполняется, существует ровно одно решение, то есть, если решение существует, модульное обратное единственно^[8].

Если $ax\equiv 1{\pmod {m}}$ имеет решение, оно обозначается часто следующим образом

x\equiv a^{-1}{\pmod {m}},

однако это можно считать злоупотреблением^[англ.], поскольку может неверно интерпретироваться как обычное обратное числа $a$ (которое, в отличие от модульного обратного, не является целым за исключением случаев, когда a равно 1 или —1). Обозначение было бы приемлемым, если бы a интерпретировался как обозначение для класса вычетов ${\overline {a}}$ , поскольку обратный элемент класса вычетов снова является классом вычетов с операцией умножения, определённой в следующем разделе.

Целые по модулю m

Отношение эквивалентности по модулю m разбивает множество целых чисел на m классов эквивалентности. Операции сложения и умножения можно определить на этих объектах следующим образом: для сложения или умножения классов эквивалентности сначала любым способом выбираются представители этих классов, затем осуществляется обычная операция над отобранными целыми числами и, наконец, результат операции лежит в классе вычетов, являющемся результатом операции над классами. В символической форме, с $+_{m}$ и $\cdot _{m}$ представляющими операции над классами вычетов, эти определения можно записать как

{\overline {a}}+_{m}{\overline {b}}={\overline {a+b}}

и

{\overline {a}}\cdot _{m}{\overline {b}}={\overline {ab}}.

Эти операции вполне определены (что означает, что конечный результат не зависит от выбора представителей).

Классы вычетов по m с этими двумя операциями образуют кольцо, называемое кольцом целых чисел по модулю m. Имеется несколько обозначений, используемых для этих алгебраических объектов, наиболее часто используется $\mathbb {Z} /m\mathbb {Z}$ или $\mathbb {Z} /m$ , однако некоторые элементарные учебники и приложения используют упрощённое обозначение $\mathbb {Z} _{m}$ , если не возникает противоречие с другими алгебраическими объектами.

Классы вычетов целых чисел по модулю m были традиционно известны как классы остатков по модулю m, что отражает факт, что все элементы класса эквивалентности имеют один и тот же остаток от деления на m. Любое множество из m целых, выбранных из разных классов вычетов по модулю m называется полной системой вычетов по модулю m^[9]. Деление столбиком показывает, что множество целых чисел {0, 1, 2, ..., m − 1} образуют полную систему вычетов по модулю m, известную как наименьшая система вычетов по модулю m. При работе с арифметическими задачами иногда удобнее работать с полной системой вычетов и использовать язык эквивалентности, в то время как в других случаях более удобен взгляд с точки зрения классов эквивалентности кольца $\mathbb {Z} /m\mathbb {Z}$ ^[10].

Мультипликативная группа кольца вычетов

Не всякий элемент полной системы вычетов по модулю m имеет обратный элемент, например, нуль обратного не имеет. После удаления элементов полной системы вычетов, что не взаимно просты с m, оставшиеся элементы, которые называются приведённой системой вычетов, все имеют обратные. Число элементов в приведённой системе вычетов равно $\phi (m)$ , где $\phi$ есть функция Эйлера, то есть количество положительных целых чисел меньших m, которые взаимно просты с m.

В кольце с единицей общего вида не всякий элемент имеет обратный, а те, что имеют, называются обратимыми. Поскольку произведение обратимых элементов обратимо, обратимые элементы кольца образуют группу, группу обратимых элементов кольца, и эта группа часто обозначается как $R^{\times }$ , если R является названием кольца. Группа обратимых элементов кольца целых чисел по модулю m называется мультипликативной группой целых по модулю m и она изоморфна приведённой системе вычетов. В частности, её порядок (размер) равен $\phi (m)$ .

Когда m является простым, скажем p, то $\phi (p)=p-1$ и все ненулевые элементы $\mathbb {Z} /p\mathbb {Z}$ имеют обратные, а тогда $\mathbb {Z} /p\mathbb {Z}$ является конечным полем. В этом случае мультипликативная группа целых чисел по модулю p образует циклическую группу порядка p − 1.

Пример

Для иллюстрации вышеприведённых определений рассмотрим пример чисел по модулю 10.

Два числа эквивалентны по 10 тогда и только тогда, когда их разность делится на 10, например

32\equiv 12{\pmod {10}}

поскольку 10 делит 32 − 12 = 20,

111\equiv 1{\pmod {10}}

поскольку 10 делит 111 − 1 = 110.

Некоторые из классов вычетов по этому модулю:

{\overline {0}}=\{\cdots ,-20,-10,0,10,20,\cdots \}

{\overline {1}}=\{\cdots ,-19,-9,1,11,21,\cdots \}

{\overline {5}}=\{\cdots ,-15,-5,5,15,25,\cdots \}

{\overline {9}}=\{\cdots ,-11,-1,9,19,29,\cdots \}.

Линейное сравнение 4x ≡ 5 (mod 10) не имеет решений, поскольку целые, конгруэнтные 5 (то есть, входящие в ${\overline {5}}$ ) все нечётные, в то время как 4x все чётные. Однако линейное сравнение $4x\equiv 6{\pmod {10}}$ имеет два решения, а именно, $x=4$ и $x=9$ . НОД(4, 10) = 2 и 2 не делит 5, но делит 6.

Поскольку НОД(3, 10) = 1, линейное сравнение $3x\equiv 1{\pmod {10}}$ будет иметь решения, то есть, модульное обратное числа 3 по модулю 10 будет существовать. 7 удовлетворяет этому уравнению (21 − 1 = 20). Однако и другие целые удовлетворяют этому уравнению, например 17 и −3 (поскольку 3(17) − 1 = 50 и 3(−3) − 1 = −10). В частности, любое целое число из ${\overline {7}}$ будет удовлетворять уравнению, поскольку эти числа имеют вид 7 + 10r для некоторого r и ясно, что

3(7+10r)-1=21+30r-1=20+30r=10(2+3r),

делится на 10. Это уравнение имеет только один класс вычетов в качестве решений. Решение в этом случае может быть получено просто перебором возможных классов, но для получения решения для больших модулей нужны систематические алгоритмы и они будут предоставлены в следующих разделах.

Произведение классов вычетов ${\overline {5}}$ и ${\overline {8}}$ может быть получено путём выбора элемента из ${\overline {5}}$ , скажем 25, и элемента из ${\overline {8}}$ , скажем −2, и их произведение (25)(−2) = −50 находится в классе эквивалентности ${\overline {0}}$ . Таким образом, ${\overline {5}}\cdot _{10}{\overline {8}}={\overline {0}}$ . Сложение определяется аналогично. Десять классов вычетов вместе с этими операциями сложения и умножения классов вычетов образует кольцо целых чисел по модулю 10, то есть $\mathbb {Z} /10\mathbb {Z}$ .

Полной системой вычетов по модулю 10 может быть множество {10, −9, 2, 13, 24, −15, 26, 37, 8, 9}, где каждое целое лежит в своём классе вычетов по модулю 10. Минимальной системой вычетов по модулю 10 служит {0, 1, 2, …, 9}. Приведённой системой вычетов по модулю 10 может быть {1, 3, 7, 9}. Произведение любых двух классов вычетов, представленных этими числами, снова является один из этих четырёх классов. Из этого следует, что эти четыре класса вычетов образуют группу, в этом случае циклическую группу порядка 4, имеющую в качестве генератора 3 или 7. Представленные классы вычетов образуют группу обратимых элементов кольца $\mathbb {Z} /10\mathbb {Z}$ . Эти классы вычетов в точности те, что имеют обратные по модулю.

Вычисление

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

Модульное обратное для числа a по модулю m можно найти с помощью расширенного алгоритма Евклида.

Алгоритм Евклида определяет наибольший общий делитель (НОД) двух целых чисел, скажем a и m. Если a имеет обратное по модулю m число, этот НОД должен быть равен 1. Несколько последних равенств, получаемых в процессе работы алгоритма, могут быть решены для нахождения НОД. Затем, используя метод «обратной подстановки», может быть получено выражение, связывающее исходные параметры и НОД. Другими словами, могут быть найдены целые x и y, удовлетворяющие равенство Безу,

ax+my={\text{НОД }}(a,m)=1.

Перепишем это следующим образом

ax-1=(-y)m,

то есть,

ax\equiv 1{\pmod {m}},

и модульное обратное числа a вычислено. Более эффективной версией является расширенный алгоритм Евклида, который с помощью дополнительных равенств сокращает два прохода алгоритма (обратную подстановку можно понимать как прохождение алгоритма в обратном порядке) до одного.

В нотации O большое этот алгоритм работает за время $O(\log(m)^{2})$ в предположении, что $|a|<m$ , и считается очень быстрым. Он обычно более эффективен чем альтернативный экспоненциальный алгоритм.

Использование теоремы Эйлера

В качестве альтернативы расширенному алгоритму Евклида для вычисления модульного обратного элемента можно рассматривать использование теоремы Эйлера^[11].

Согласно теореме Эйлера, если a взаимно просто с m, то есть НОД(a, m) = 1, то

a^{\phi (m)}\equiv 1{\pmod {m}},

где $\phi$ — функция Эйлера. Это следует из факта, что a принадлежит мультипликативной группе $(\mathbb {Z} /m\mathbb {Z} )^{\times }$ тогда и только тогда, когда a взаимно просто с m. Таким образом, модульное обратное можно найти напрямую:

a^{\phi (m)-1}\equiv a^{-1}{\pmod {m}}.

В специальном случае, когда m простое, $\phi (m)=m-1$ и модульное обратное задаётся равенством

a^{-1}\equiv a^{m-2}{\pmod {m}}.

Этот метод, как правило, медленнее расширенного алгоритма Евклида, но иногда используется, если реализация для модульного возведения в степень уже доступна. Недостатки этого метода:

Значение $\phi (m)$ должно быть известно, а наиболее эффективный метод вычисления требует m разложений. Разложение хорошо известно как вычислительно сложная задача. Однако вычислить $\phi (m)$ можно непосредственно, если известно разложение на простые множители числа m.
Относительно высокая стоимость возведения в степень. Хотя это можно реализовать более эффективно с помощью модульного возведения, при больших значениях m наиболее эффективен метод алгоритма Монтгомери. Этот алгоритм сам по себе требует вычисления обратного по модулю m. Без алгоритма Монтгомери стандартный бинарный алгоритм возведения, который требует деления по модулю m на каждом шаге, является медленной операцией при большом m.

Замечательное преимущество этой техники в том, что нет условных ветвлений, которые зависят от значения a, а потому значение a, которое может быть важным секретом в криптосистемах с открытым ключом, может быть защищено от атак по сторонним каналам. По этой причине стандартная реализация Curve25519^?! использует технику вычисления обратного элемента.

Вычисление нескольких обратных

Можно вычислить обратные для нескольких чисел $a_{i}$ по общему модулю m с помощью одного прохода алгоритма Евклида и трёх умножений на каждое дополнительное входное число^[12]. Основной идеей является образование всех $a_{i}$ , обращение его, а затем умножение на на $a_{j}$ для всех $j\neq i$ , чтобы остался только $a_{i}^{-1}$ .

Алгоритм (вся арифметика осуществляется по модулю m):

Вычисляем префиксные произведения ${\textstyle b_{i}=\prod _{j=1}^{i}a_{j}=a_{i}b_{i-1}}$ для всех $i\leqslant n$ .
Вычисляем $b_{n}^{-1}$ с помощью любого доступного алгоритма.
Для i от n до 2 вычисляем
- $a_{i}^{-1}=b_{i}^{-1}b_{i-1}$ и
- $b_{i-1}^{-1}=b_{i}^{-1}a_{i}$ .
Наконец, $a_{1}^{-1}=b_{1}^{-1}$ .

Можно осуществить умножение в виде древесной структуры, а не линейной, чтобы обеспечить параллельность вычислений.

Приложения

Поиск модульного обратного имеет много приложений в алгоритмах, опирающихся на теорию модульной арифметики. Например, в криптографии использование модульной арифметики позволяет осуществлять некоторые операции быстрее и с меньшими требованиями к памяти, в то время как другие операции становится выполнить труднее^[13]. Оба этих свойства могут использоваться во благо. В частности, в алгоритме RSA шифрование и дешифрование сообщения делается с помощью пары чисел, обратных друг другу по некоторому тщательно выбранному модулю. Один из этих чисел делается открытым и он может быть использован в быстрой процедуре шифрования, в то время как другое число используется в процедуре расшифровки и не разглашается. Определение скрытого ключа по открытому считается вычислительно невыполнимой задачей, что и защищает от несанкционированного доступа^[14].

В качестве другого использования в другом контексте рассмотрим задачу точного деления в компьютерах, когда у вас есть список нечётных чисел размером в компьютерное слово, каждое из которых делится на k и вы хотите все их разделить на k. Одно из решений следующее:

Используем расширенный алгоритм Евклида для вычисления $k^{-1}$ , модульного обратного $k{\pmod {2^{w}}}$ , где w равно числу бит в слове. Это обратное существует, поскольку все числа нечётные, а рассматриваются вычеты по модулю, не имеющему нечётных делителей.
Каждое число в списке умножаем на $k^{-1}$ и выбираем наименее значащие биты результата (то есть отбрасываем все биты за пределами компьютерного слова).

Во многих машинах, в частности тех, у которых нет аппаратной поддержки деления, деление является более медленной операцией, чем умножение, так что этот подход может повлечь существенное увеличение скорости. Первый шаг относительно медленный, но его нужно сделать всего один раз.

Модульные обратные числа используются для получения решения системы линейных сравнений, что гарантируется китайской теоремой об остатках.

Например, система

X\equiv 4{\pmod {5}}

X\equiv 4{\pmod {7}}

X\equiv 6{\pmod {11}}

имеет общее решение, поскольку 5, 7 и 11 попарно взаимно просты. Решение даётся формулой

X=t_{1}(7\times 11)\times 4+t_{2}(5\times 11)\times 4+t_{3}(5\times 7)\times 6

где

t_{1}=3

модульное обратное 7 × 11 (mod 5),

t_{2}=6

модульное обратное 5 × 11 (mod 7)

t_{3}=6

модульное обратное 5 × 7 (mod 11).

Тогда,

X=3\times (7\times 11)\times 4+6\times (5\times 11)\times 4+6\times (5\times 7)\times 6=3504

и в приведённом виде

X\equiv 3504\equiv 39{\pmod {385}}

поскольку 385 является НОК чисел 5, 7 и 11.

Модульное обратное фигурирует на видном месте в определении сумм Клоостермана.

См. также

Инверсный конгруэнтный метод — алгоритм генерации псевдослучайных чисел, использующий обратные по модулю числа.
Восстановление рационального числа по его модулю^[англ.]

Примечания

↑ Rosen, 1993, с. 132.
↑ Schumacher, 1996, с. 88.
↑ Stinson, 1995, с. 124–128.
↑ Trappe, Washington, 2006, с. 164−169.
↑ Moriarty, K.; Kaliski, B.; Jonsson, J.; Rusch, A. PKCS #1: RSA Cryptography Specifications Version 2.2 (неопр.). Internet Engineering Task Force RFC 8017. Internet Engineering Task Force (2016). Дата обращения: 21 января 2017.
↑ Other notations are often used, including [a] and [a]_m.
↑ Ireland, Rosen, 1990, с. 32.
↑ Shoup, 2005, с. 15 Theorem 2.4.
↑ Rosen, 1993, с. 121.
↑ Ireland, Rosen, 1990, с. 31.
↑ Koshy, 2007, с. 346.
↑ Brent, Zimmermann, 2010, с. 67–68.
↑ Trappe, Washington, 2006, с. 167.
↑ Trappe, Washington, 2006, с. 165.

Литература

Douglas R. Stinson. Cryptography / Theory and Practice. — CRC Press, 1995. — ISBN 0-8493-8521-0.
Victor Shoup. A Computational Introduction to Number Theory and Algebra. — Cambridge University Press, 2005. — ISBN 9780521851541.
Thomas Koshy. Elementary number theory with applications. — 2nd edition. — Academic Press, 2007. — ISBN 978-0-12-372487-8.
Richard P. Brent, Paul Zimmermann. §2.5.1 Several inversions at once // Elementary number theory with applications. Modern Computer Arithmetic. — Cambridge University Press, 2010. — Т. 18. — (Cambridge Monographs on Computational and Applied Mathematics). — ISBN 978-0-521-19469-3.

Kenneth Ireland, Michael Rosen. A Classical Introduction to Modern Number Theory. — 2nd. — Springer-Verlag, 1990. — ISBN 0-387-97329-X.
- К. Айерлэнд, М. Роузен. Классическое введение в современную теорию чисел / Перевод с английского С. П. Демушкина. — Москва: «Мир», 1987.

Kenneth H. Rosen. Elementary Number Theory and its Applications. — 3rd. — Addison-Wesley, 1993. — ISBN 978-0-201-57889-8.
Carol Schumacher. Chapter Zero: Fundamental Notions of Abstract Mathematics. — Addison-Wesley, 1996. — ISBN 0-201-82653-4.
Wade Trappe, Lawrence C. Washington. Introduction to Cryptography with Coding Theory. — 2nd. — Prentice-Hall, 2006. — ISBN 978-0-13-186239-5.

Ссылки

Weisstein, Eric W. Modular Inverse (англ.) на сайте Wolfram MathWorld.
Guevara Vasquez, Fernando provides a solved example of solving the modulo multiplicative inverse using Euclid’s Algorithm

[_1ef089f524c52bc8-1] Rosen, 1993, с. 132.

[_2270fa8c9b595bed-2] Schumacher, 1996, с. 88.

[_7ed180ea0128122a-3] Stinson, 1995, с. 124–128.

[_282b4ebc56a9d1c8-4] Trappe, Washington, 2006, с. 164−169.

[ietf-5] Moriarty, K.; Kaliski, B.; Jonsson, J.; Rusch, A. PKCS #1: RSA Cryptography Specifications Version 2.2 (неопр.). Internet Engineering Task Force RFC 8017. Internet Engineering Task Force (2016). Дата обращения: 21 января 2017.

[Other-6] Other notations are often used, including [a] and [a]_m.

[_8d183a1fe53e2f85-7] Ireland, Rosen, 1990, с. 32.

[_95869c38ce99ce87-8] Shoup, 2005, с. 15 Theorem 2.4.

[_1ef08af524c52d9c-9] Rosen, 1993, с. 121.

[_8d183a1fe53e2f86-10] Ireland, Rosen, 1990, с. 31.

[_a19ae7ccbf88cbc3-11] Koshy, 2007, с. 346.

[_ff36d609541757ab-12] Brent, Zimmermann, 2010, с. 67–68.

[_cfa2d7c797a63955-13] Trappe, Washington, 2006, с. 167.

[_cfa2d7c797a63957-14] Trappe, Washington, 2006, с. 165.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]