Как вычислить символ якоби
Перейти к содержимому

Как вычислить символ якоби

  • автор:

Вычисление символа Якоби

Не получается написать правильную реализацию алгоритма вычисления символа Якоби J(Q/P). Сначала сам алгоритм: Вход: Q , P — целые числа. Выход: значение символа Якоби. 1) s = 0, u = Q, v = P. 2) Вычисляем r — наименьший положительный остаток при делении u на v . Вычисляем целое k >= 0 и нечетное t , такие, что r = t * 2^k . Вычисляем `s = s + k * (v^2 — 1)/8 + (t — 1)*(v — 1)/4 (mod 2) 3) Если t = 1 , то символ Якоби равен (-1)^s . Конец. 4) Если t >= 3 , то u = v , v = t , переходим на шаг 2. Реализация:

#include int jacobi(int q, int p) < int s = 0, u = q, v = p; int r, k, t; do< // Вычисляем r - наименьший положительный остаток при делении u на v r = u % v; // Вычисляем целое k >= 0 и нечетное t: r = t * 2^k k = t = 0; while(r % 2 == 0) < k++; // Показатель степени двойки в числе r r >>= 1; // Делим r на 2 > t = r; // В t находится результат деления r на 2^k s = (s + k * (v*v - 1)/8 + (t - 1)*(v - 1)/4) % 2; if(t == 1) return (s) ? 1 : -1; // Новая итерация u = v; v = t; >while(t >= 3); > int main()

Здесь J(-104, 997) должен быть равен -1, а моя реализация выдает 1. Еще тесты из хелпа Maple: jacobi(12, 3) = 0 , jacobi(28, 21) = 0 , jacobi(6, 11) = -1 , jacobi(226, 135) = 1 , jacobi(26, 35) = -1 , jacobi(-286, 4272943) = 1 , jacobi(888, 1999) = -1 . Дополнение: в тернарном выражении я ошибся. Если s == 1 (true), то надо возвращать -1, потому что показатель степени нечетный. Может быть, математики что-то дополнят?

Якоби символ

Символ Якоби — теоретико-числовая функция двух аргументов, введённая К. Якоби в 1837 году. Является квадратичным характером в кольце вычетов.

Символ Якоби обобщает символ Лежандра на все нечётные числа, большие единицы. Символ Кронекера — Якоби, в свою очередь, обобщает символ Якоби на все целые числа, но в практических задачах символ Якоби играет гораздо более важную роль, чем символ Кронекера — Якоби.

Определение

P=p_1p_2\cdots p_n

Пусть P — нечётное, большее единицы число и — его разложение на простые множители (среди p1. pn могут быть равные). Тогда для произвольного целого числа a символ Якоби определяется равенством:

\left(\frac<P></p>
<p>\right) = \left(\frac\right)\left(\frac\right)\cdots\left(\frac\right) » width=»» height=»» />,</p><div class='code-block code-block-3' style='margin: 8px 0; clear: both;'>
<!-- 3joomlaumnik -->
<script src=

\left(\frac<p_i></p>
<p>где \right)» width=»» height=»» /> — символы Лежандра.</p>
<p><img decoding=

По определению считаем, что \right)=1″ width=»» height=»» /> для всех a.

Свойства

Значения символа Якоби для аргументов от 1 до 100

  • Мультипликативность: \left(\frac<ab>\right) = \left(\frac\right) \left(\frac\right)» width=»» height=»» />
<ul>
<li>В частности, <img decoding=, то  \left(\frac\right) = \left(\frac\right)
  • \left(\frac<1>\right) = 1″ width=»» height=»» /></li>
<li><img decoding= равен знаку перестановки приведённой системы вычетов по модулю P , которая задается как умножение элементов этой группы на a (где a обязательно взаимно просто с P ).

Важные замечания

О вычислении

Символ Якоби практически никода не вычисляют по определению. Чаще всего для вычисления используют свойства символа Якоби, главным образом — квадратичный закон взаимности.

Более того, несмотря на то, что символ Якоби является обобщением символа Лежандра и определяется через него, чаще именно символ Лежандра вычисляют с помощью символа Якоби, а не наоборот. См. Пример

О связи с квадратичными сравнениями

В отличие от символа Лежандра, символ Якоби нельзя напрямую использовать для проверки разрешимости квадратичного сравнения. То есть, если задано сравнение

x^2 \equiv a \mod<n>,» width=»» height=»» /></td>
<td style=(1)

то равенство единице символа Якоби \left(\frac<n>\right)» width=»» height=»» /> вовсе не означает, что данное сравнение разрешимо. Например, <img decoding= не имеет решений (можно проверить перебором).

\left(\frac<n></p><div class='code-block code-block-7' style='margin: 8px 0; clear: both;'>
<!-- 7joomlaumnik -->
<script src=

Но если \right)=-1″ width=»» height=»» />, то сравнение (1) не имеет решений.

Особенность, используемая в тестах простоты

В общем случае неверно, что для символа Якоби выполняется то же условие, что и для символа Лежандра:

\left(\frac\right) \equiv a^> \mod.» width=»» height=»» /></td>
<td style=(2)

\left(\frac</p><div class='code-block code-block-8' style='margin: 8px 0; clear: both;'>
<!-- 8joomlaumnik -->
<script src=

\right)= \left(\frac\right)\cdot \left(\frac\right)= \left(\frac\right)\cdot \left(\frac\right)=(-1)^>\cdot1=-1″ width=»» height=»» />

7^<\frac<15-1></p>
<p>При этом > \equiv 7^7\equiv13. \mod» width=»» height=»» /></p><div class='code-block code-block-9' style='margin: 8px 0; clear: both;'>
<!-- 9joomlaumnik -->
<script src=

Числа a, взаимно простые с P, для которых не выполнено условие (2), называются Эйлеровыми свидетелями непростоты числа P (поскольку для простого P условие (2) выполнено). Если P — составное число, то такое число a, для которого условие (2) выполнено, называют лжецом для теста Эйлера. Доказано, что для любого составного P есть не более P/2 лжецов, различных по модулю P.

Данное свойство используется в вероятностном тесте Соловея-Штрассена на простоту. В этом алгоритме выбираются случайные числа a и для них проверяется условие (2). Если находится свидетель непростоты, то доказано, что число P – составное. В противном случае говорят, что P — простое с некоторой вероятностью.

Применение

a^<\frac<p-1></p><div class='code-block code-block-10' style='margin: 8px 0; clear: both;'>
<!-- 10joomlaumnik -->
<script src=

Главным образом, символ Якоби используется для быстрого вычисления символа Лежандра. Символ Лежандра, в свою очередь, необходим для проверки разрешимости квадратичного сравнения по модулю простого числа. Но считать его по определению (то есть вычислять >\mod

» width=»» height=»» />) — достаточно долгая по времени процедура. С помощью алгоритма быстрого возведения в степень это делается за O((logp) 3 ) битовых операций (если не использовать быстрое умножение и деление). А вычисление символа Якоби требует только O((logp) 2 ) битовых операций.

Символ Якоби используется в некоторых тестах на простоту, например, в (N+1) – методах и, как уже было сказано, в тесте Соловея — Штрассена.

Алгоритм

Основная идея

Ключевое используемое при вычислении свойство символа Якоби — квадратичный закон взаимности. Благодаря нему алгоритм похож на алгоритм Евклида нахождения наибольшего общего делителя двух чисел, в котором тоже аргументы на каждом шаге меняются местами. Аналогично алгоритму Евклида, при перестановке аргументов больший заменяется на остаток от деления на меньший. Это возможно благодаря периодичности символа Якоби. Однако, поскольку символ Якоби определён только при условии нечётности второго аргумента, то до перестановки выделяется чётная часть первого аргумента.

Формальное описание

Входные данные: a — целое число, b — натуральное, нечётное, больше единицы..

\left(\frac<b></p>
<p><b>Выходные данные:</b> \right)» width=»» height=»» /> — символ Якоби</p>
<pre><b>1</b> <i>(проверка взаимной простоты).</i> Если НОД (<i>a</i>, <i>b</i>) ≠1, выход из алгоритма с ответом <b>0</b>. <b>2</b> <i>(инициализация).</i> <i>r</i>:=1 <b>3</b> <i>(переход к положительным числам).</i> <b>Если</b> <i>a <b>то</b> <i>a:=-a</i> <b>Если</b> <i>b</i> mod 4 = 3 <b>то</b> <i>r:=-r</i> <b>Конец если</b> <b>4</b> <i>(избавление от чётности).</i> <i>t</i>:=0 <b>Цикл</b> ПОКА <i>a</i> – чётное <i>t:=t+1</i> <i>a:=a/2</i> <b>Конец цикла</b> <b>Если</b> <i>t</i> – нечётное, <b>то</b> <b>Если</b> <i>b</i> mod 8 = 3 или 5, <b>то</b> <i>r</i>:=-<i>r</i>. <b>Конец если</b> <b>5</b> <i>(квадратичный закон взаимности).</i> <b>Если</b> <i>a</i> mod 4 = <i>b</i> mod 4 = 3, <b>то</b> <i>r</i>:=-<i>r</i>. c:=a; a:=b mod c; b:=c. <b>6</b> <i>(выход из алгоритма?).</i> Если <i>a</i>≠0, то идти на шаг 4, иначе выйти из алгоритма с ответом <b>r</b>.</pre>
<h4>Комментарии к алгоритму</h4>
<p>В алгоритме везде берётся наименьший положительный вычет (то есть остаток от деления).</p><div class='code-block code-block-13' style='margin: 8px 0; clear: both;'>
<!-- 13joomlaumnik -->
<script src=

\left(\frac<2></p>
<p>На четвёртом шаге используется мультипликативность символа Якоби, а затем вычисляется символ Якоби \right) =(-1)^» width=»» height=»» />. Чтобы избежать лишнего возведения в степень, проверяется только остаток от деления <i>b</i> на 8.</p>
<p>На пятом шаге тоже вместо возведения в степень используется проверка остатков от деления.</p><div class='code-block code-block-14' style='margin: 8px 0; clear: both;'>
<!-- 14joomlaumnik -->
<script src=

O(\log\cdot\log<b></p>
<p>Сложность алгоритма равна )» width=»» height=»» /> битовых операций.</p><div class='code-block code-block-15' style='margin: 8px 0; clear: both;'>
<!-- 15joomlaumnik -->
<script src=

Пример вычисления

Вычисление символа Лежандра с помощью символа Якоби:

\left(\frac\right) = -\left(\frac\right) = -\left(\frac\right) =-\left(\frac\right) = - \left(\frac\right) = - \left(\frac\right) =  = - \left(\frac<2>\right) \left(\frac\right) =- \left(\frac\right) = — \left(\frac\right) = — \left(\frac\right) = 1″ width=»» height=»» /></p>
<h3>Список литературы</h3>
<ul>
<li><i>Василенко О.Н.</i> Теоретико-числовые алгоритмы в криптографии. — Москва: МЦМНО, 2003. — С. 328. — ISBN 5-94057-103-4</li>
</ul>
<ul>
<li><i>Виноградов И. М.</i>Основы теории чисел. — Москва: ГИТТЛ, 1952.</li>
</ul>
<ul>
<li><i>Bach E., Shallit J.</i> Algorithmic Number Theory.Vol. I.. — Massachusetts: MIT Press, 1996. — ISBN 0-262-02405-5</li>
</ul>
<table style= Характеры Квадратичные характеры Символ Лежандра • Символ Якоби • Символ Кронекера — Якоби Характеры степенных вычетов Характер кубического вычета • Характер биквадратичного вычета • Символ степенного вычета просмотр • обсуждение

Алгоритм Соловея-Штрассена

Роберт Соловей и Фолькер Штрассен разработали алгоритм вероятностного тестирования простоты числа, который использует символ Якоби. Определяет числа как составные или вероятно простые. Распознает числа Кармайкла как составные.
Итак, для начала необходимо ввести нужные понятия.
Квадратичный вычет. Если число p — простое и 0 < a < p, то число a является квадратичным вычетом по модулю p, если существуют значения x такие, что
x2 = a (mod p).
Для того, чтобы число a было квадратичным вычетом по модулю n, оно должно быть квадратичным вычетом по модулю всех простых делителей n. Например, если n = 7, то квадратичные вычеты равны 1, 2 и 4.
12 = 1 = 1 mod 7,
22 = 4 = 4 mod 7,
32 = 9 = 2 mod 7,
42 = 16 = 2 mod 7,
52 = 25 = 1 mod 7,
62 = 36 = 1 mod 7.
И наоборот, в следующих уравнениях не существует значений x, которые их удовлетворяют.
x2 = 3 mod 7,
x2 = 5 mod 7,
x2 = 6 mod 7.
Итак, числа 3, 5 и 6 являются квадратичными невычетами по модулю 7.
Если число p — нечетное, то существует ровно (p – 1)/2 квадратичных вычетов по модулю p и столько же квадратичных невычетов по модулю p. Если n — произведение двух простых чисел p и q, то существует ровно (p – 1)(q – 1)/4 квадратичных вычетов по модулю n.
Связь между простыми числами и квадратичными вычетами устанавливается с помощью символов Лежандра и Якоби.
Символ Лежандра, который обозначается как L(a, p) — это функция, определенная, если a — любое целое число, а p — простое число, превышающее 2. Символ Лежандра может принимать значения 0, 1 и –1.
L(a, p) = 0, если a делится на p.
L(a, p) = 1, если a — квадратичный вычет по модулю p,
L(a, p) = –1, если a — квадратичный невычет по модулю p.
Сжато, эти факты записываются так:
L(a, p) = a^((p – 1)/2) mod p.

Алгоритм вычисления символа Лежандра.

1. Если a = 1, то L(a, p) = 1.
2. Если число a четное, то L(a, p) = L(a/2, p)*((-1)^((p^2-1)/8)).
3. Если число a — нечетное и a != 1, то L(a, p) = L(p mod a, a)*((–1)^((a–1)*(p–1)/4)).
Символ Якоби, который обозначается как J(a, n) — это обощение символа Лежандра на составные модули. Это функция, определенная для всех целых чисел a и нечетных целых чисел n. Символ Якоби может принимать значения 0, 1 и –1.
Символ Якоби можно задать следующим образом.
1. Символ Якоби определен только для нечетных чисел n.
2. J(0, n) = 0.
3. Если n – простое число, то J(0, n) = 0, если a делится на n.
4. Если n – простое число, то J(0, n) = 1, если a — квадратичный вычет по модулю n.
5. Если n – простое число, то J(0, n) = –1, если a — квадратичный невычет по модулю n.
6. Если n – составное число, то J(a, n) = J(a, p1)*. *J(a, pm), где p1. pm — разложение n на простые множители.

Алгоритм вычисления символа Якоби.

1. J(1, n) = 1.
2. J(a*b, n) = J(a, n)*J(b, n).
3. J(2, n) = 1, если (n^2 – 1)/8 является четным, и –1 в противном случае.
4. J(a, n) = J((a mod m), n).
5. J(a, b1*b2) = J(a, b1)J(a, b2).
6. Если gcd(a, b) = 1 и, кроме того, числа a и b являются нечетными, то
6.1. J(a, b) = J(b, a), если (a – 1)*(b – 1)/4 является четным числом.
6.2. J(a, b) = –J(b, a), если (a – 1)*(b – 1)/4 является нечетным числом.
Если n — простое число, то символ Якоби эквивалентен символу Лежандра.
Символ Якоби нельзя использовать для проверки, является ли число a квадратичным вычетом по модулю n (кроме случая, когда число n — простое). Если J(a, n) = 1 и n — составное число, то число a не всегда является квадратичным вычетом:
J(7, 143) = J(7, 11) * J(7, 13) = (–1)*(–1) = 1,
хотя не существует целых чисел x таких, что x2  7 (mod 143).

Алгоритм Соловея–Штрассена

1. Выберите случайное число a, меньшее p.
2. Если gcd(a, p) != 1, то число p – составное и тест можно не продолжать.
3. Вычислите j = a^((p–1)/2) mod p.
4. Вычислите символ Якоби J(a, p).
5. Если j != J(a, p), то число p точно не является простым.
6. Если j = J(a, p), то вероятность того, що число p не является простым, не превышает 50%.
Число a, которое не указывает явно, что число p не простое, называется свидетелем. Если число p — составное, то вероятность того, що случайное число является свидетелем, составляет не менее 50%. Вероятность того, что составное число пройдет t испытаний, равняется 1/(2^t).

Символ Якоби и его свойства

Свойства символа Якоби прямо вытекают из соответствующих свойств символа Лежандра. Их доказательство оставляется читателю в качестве самостоятельного упражнения.

[math]a_1\equiv a \pmod n\Rightarrow\left(\cfrac\right)=\left(\cfrac\right)[/math]
[math]\left(\cfrac\right)=\left(\cfrac\right)\left(\cfrac\right)[/math]
НОД [math](a,n)=1\Rightarrow\left(\cfrac\right)=\left(\cfrac\right)[/math]
[math]\left(\cfrac<1>\right)=1[/math]
[math]\left(\cfrac<-1>\right)=(-1)^>[/math]

Рассмотрим нечетные [math]n[/math] и [math]m[/math] :

[math]0\equiv(n-1)(m-1)\pmod 4\Rightarrow n-1+m-1=nm-1\pmod 4\Rightarrow \cfrac+~\cfrac\equiv~\cfrac\pmod 4\Rightarrow\cfrac+\cdots+\cfrac\equiv\cfrac\pmod 2[/math]

[math]\left(\cfrac<2>\right)=(-1)^>[/math]

Аналогично предыдущему докажем, что

Рассмотрим нечетные [math]n[/math] и [math]m[/math] :

[math]0\equiv(n^2-1)(m^2-1)\pmod 16\Rightarrow n^2-1+m^2-1\equiv n^2m^2-1\pmod 16\Rightarrow \cfrac+\cfrac\equiv\cfrac\pmod 2\Rightarrow\cfrac+\cdots+\cfrac\equiv\cfrac<(p_1p_2\cdots p_s-1)^2>\pmod 2[/math]

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *