Символ Якоби - Jacobi symbol
k п | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 1 | ||||||||||||||||
3 | 0 | 1 | −1 | ||||||||||||||
5 | 0 | 1 | −1 | −1 | 1 | ||||||||||||
7 | 0 | 1 | 1 | −1 | 1 | −1 | −1 | ||||||||||
9 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | ||||||||
11 | 0 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | ||||||
13 | 0 | 1 | −1 | 1 | 1 | −1 | −1 | −1 | −1 | 1 | 1 | −1 | 1 | ||||
15 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | −1 | 1 | 0 | 0 | −1 | 0 | −1 | −1 | ||
17 | 0 | 1 | 1 | −1 | 1 | −1 | −1 | −1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | 1 | 1 |
В Символ Якоби является обобщением Символ Лежандра. Представлен Якоби в 1837 г.,[1] это представляет теоретический интерес в модульная арифметика и другие отрасли теория чисел, но его основное применение - в вычислительная теория чисел, особенно проверка на простоту и целочисленная факторизация; это, в свою очередь, важно в криптография.
Определение
Для любого целого числа а и любое положительное нечетное целое число п, символ Якоби (а/п) определяется как продукт Лежандровые символы соответствующие простым факторам п:
куда
разложение на простые множители п.
Символ Лежандра (а/п) определено для всех целых чисел а и все нечетные простые числа п к
Следуя обычному соглашению для пустого продукта, (а/1) = 1.
Когда нижний аргумент - нечетное простое число, символ Якоби равен символу Лежандра.
Таблица значений
Ниже приводится таблица значений символа Якоби. (k/п) с п ≤ 59, k ≤ 30, п странный.
k п | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
3 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 |
5 | 1 | −1 | −1 | 1 | 0 | 1 | −1 | −1 | 1 | 0 | 1 | −1 | −1 | 1 | 0 | 1 | −1 | −1 | 1 | 0 | 1 | −1 | −1 | 1 | 0 | 1 | −1 | −1 | 1 | 0 |
7 | 1 | 1 | −1 | 1 | −1 | −1 | 0 | 1 | 1 | −1 | 1 | −1 | −1 | 0 | 1 | 1 | −1 | 1 | −1 | −1 | 0 | 1 | 1 | −1 | 1 | −1 | −1 | 0 | 1 | 1 |
9 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 |
11 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | 0 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | 0 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 |
13 | 1 | −1 | 1 | 1 | −1 | −1 | −1 | −1 | 1 | 1 | −1 | 1 | 0 | 1 | −1 | 1 | 1 | −1 | −1 | −1 | −1 | 1 | 1 | −1 | 1 | 0 | 1 | −1 | 1 | 1 |
15 | 1 | 1 | 0 | 1 | 0 | 0 | −1 | 1 | 0 | 0 | −1 | 0 | −1 | −1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | −1 | 1 | 0 | 0 | −1 | 0 | −1 | −1 | 0 |
17 | 1 | 1 | −1 | 1 | −1 | −1 | −1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | 1 | 1 | 0 | 1 | 1 | −1 | 1 | −1 | −1 | −1 | 1 | 1 | −1 | −1 | −1 | 1 |
19 | 1 | −1 | −1 | 1 | 1 | 1 | 1 | −1 | 1 | −1 | 1 | −1 | −1 | −1 | −1 | 1 | 1 | −1 | 0 | 1 | −1 | −1 | 1 | 1 | 1 | 1 | −1 | 1 | −1 | 1 |
21 | 1 | −1 | 0 | 1 | 1 | 0 | 0 | −1 | 0 | −1 | −1 | 0 | −1 | 0 | 0 | 1 | 1 | 0 | −1 | 1 | 0 | 1 | −1 | 0 | 1 | 1 | 0 | 0 | −1 | 0 |
23 | 1 | 1 | 1 | 1 | −1 | 1 | −1 | 1 | 1 | −1 | −1 | 1 | 1 | −1 | −1 | 1 | −1 | 1 | −1 | −1 | −1 | −1 | 0 | 1 | 1 | 1 | 1 | −1 | 1 | −1 |
25 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 0 |
27 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 | 1 | −1 | 0 |
29 | 1 | −1 | −1 | 1 | 1 | 1 | 1 | −1 | 1 | −1 | −1 | −1 | 1 | −1 | −1 | 1 | −1 | −1 | −1 | 1 | −1 | 1 | 1 | 1 | 1 | −1 | −1 | 1 | 0 | 1 |
31 | 1 | 1 | −1 | 1 | 1 | −1 | 1 | 1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | −1 | 1 | −1 | −1 | 1 | −1 | −1 |
33 | 1 | 1 | 0 | 1 | −1 | 0 | −1 | 1 | 0 | −1 | 0 | 0 | −1 | −1 | 0 | 1 | 1 | 0 | −1 | −1 | 0 | 0 | −1 | 0 | 1 | −1 | 0 | −1 | 1 | 0 |
35 | 1 | −1 | 1 | 1 | 0 | −1 | 0 | −1 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | −1 | −1 | 0 | 0 | −1 | −1 | −1 | 0 | −1 | 1 | 0 | 1 | 0 |
37 | 1 | −1 | 1 | 1 | −1 | −1 | 1 | −1 | 1 | 1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | −1 | −1 | −1 | 1 | −1 | −1 | −1 | 1 | 1 | 1 | 1 | −1 | 1 |
39 | 1 | 1 | 0 | 1 | 1 | 0 | −1 | 1 | 0 | 1 | 1 | 0 | 0 | −1 | 0 | 1 | −1 | 0 | −1 | 1 | 0 | 1 | −1 | 0 | 1 | 0 | 0 | −1 | −1 | 0 |
41 | 1 | 1 | −1 | 1 | 1 | −1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | −1 | −1 | 1 | −1 | 1 | −1 | 1 | 1 | −1 | 1 | −1 | 1 | −1 | −1 | −1 | −1 | −1 |
43 | 1 | −1 | −1 | 1 | −1 | 1 | −1 | −1 | 1 | 1 | 1 | −1 | 1 | 1 | 1 | 1 | 1 | −1 | −1 | −1 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | −1 | −1 |
45 | 1 | −1 | 0 | 1 | 0 | 0 | −1 | −1 | 0 | 0 | 1 | 0 | −1 | 1 | 0 | 1 | −1 | 0 | 1 | 0 | 0 | −1 | −1 | 0 | 0 | 1 | 0 | −1 | 1 | 0 |
47 | 1 | 1 | 1 | 1 | −1 | 1 | 1 | 1 | 1 | −1 | −1 | 1 | −1 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | 1 | −1 | −1 | 1 | 1 | −1 | 1 | 1 | −1 | −1 |
49 | 1 | 1 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 1 | 1 | 0 | 1 | 1 | 1 | 1 | 1 | 1 | 0 | 1 | 1 |
51 | 1 | −1 | 0 | 1 | 1 | 0 | −1 | −1 | 0 | −1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 1 | 0 | −1 | 1 | 0 | 1 | −1 | 0 | −1 | 1 | 0 |
53 | 1 | −1 | −1 | 1 | −1 | 1 | 1 | −1 | 1 | 1 | 1 | −1 | 1 | −1 | 1 | 1 | 1 | −1 | −1 | −1 | −1 | −1 | −1 | 1 | 1 | −1 | −1 | 1 | 1 | −1 |
55 | 1 | 1 | −1 | 1 | 0 | −1 | 1 | 1 | 1 | 0 | 0 | −1 | 1 | 1 | 0 | 1 | 1 | 1 | −1 | 0 | −1 | 0 | −1 | −1 | 0 | 1 | −1 | 1 | −1 | 0 |
57 | 1 | 1 | 0 | 1 | −1 | 0 | 1 | 1 | 0 | −1 | −1 | 0 | −1 | 1 | 0 | 1 | −1 | 0 | 0 | −1 | 0 | −1 | −1 | 0 | 1 | −1 | 0 | 1 | 1 | 0 |
59 | 1 | −1 | 1 | 1 | 1 | −1 | 1 | −1 | 1 | −1 | −1 | 1 | −1 | −1 | 1 | 1 | 1 | −1 | 1 | 1 | 1 | 1 | −1 | −1 | 1 | 1 | 1 | 1 | 1 | −1 |
Характеристики
Следующие ниже факты, даже законы взаимности, являются прямым выводом из определения символа Якоби и соответствующих свойств символа Лежандра.[2]
Символ Якоби определяется только тогда, когда верхний аргумент («числитель») является целым числом, а нижний аргумент («знаменатель») является положительным нечетным целым числом.
- 1. Если п является (нечетным) простым числом, то символ Якоби (а/п) равно (и записывается так же) соответствующему символу Лежандра.
- 2. Если а ≡ б (мод п), тогда
- 3.
Если фиксирован верхний или нижний аргумент, символ Якоби является полностью мультипликативная функция в оставшемся аргументе:
- 4.
- 5.
В закон квадратичной взаимности: если м и п - нечетные положительные взаимно простые целые числа, то
- 6.
и его добавки
- 7.
- 8.
Объединение свойств 4 и 8 дает:
- 9.
Как символ Лежандра:
- Если (а/п) = −1, тогда а является квадратичным невычетом по модулю п.
- Если а это квадратичный вычет по модулю п и gcd (а,п) = 1, то (а/п) = 1.
Но, в отличие от символа Лежандра:
- Если (а/п) = 1, тогда а может быть или не быть квадратичным остатком по модулю п.
Это потому, что для а быть квадратичным вычетом по модулю п, это должен быть квадратичный вычет по модулю каждый главный фактор п. Однако символ Якоби равен единице, если, например, а является невычетом по модулю ровно двух простых делителей числа п.
Хотя символ Якоби нельзя единообразно интерпретировать в терминах квадратов и неквадратов, его можно единообразно интерпретировать как знак перестановки посредством Лемма Золотарева.
Символ Якоби (а/п) это Dirichlet персонаж по модулю п.
Вычисление символа Якоби
Приведенные выше формулы приводят к эффективному О (бревно а бревно б)[3] алгоритм вычисления символа Якоби, аналогичный Евклидов алгоритм для нахождения НОД двух чисел. (Это не должно вызывать удивления в свете правила 2.)
- Уменьшите «числитель» по модулю «знаменатель», используя правило 2.
- Извлеките любой четный «числитель», используя правило 9.
- Если «числитель» равен 1, правила 3 и 4 дают результат 1. Если «числитель» и «знаменатель» не взаимно просты, правило 3 дает результат 0.
- В противном случае «числитель» и «знаменатель» теперь являются нечетными положительными взаимно простыми целыми числами, поэтому мы можем перевернуть символ, используя правило 6, а затем вернуться к шагу 1.
Реализация в Lua
функция Якоби(п, k) утверждать(k > 0 и k % 2 == 1) п = п % k т = 1 пока п ~= 0 делать пока п % 2 == 0 делать п = п / 2 р = k % 8 если р == 3 или же р == 5 тогда т = -т конец конец п, k = k, п если п % 4 == 3 и k % 4 == 3 тогда т = -т конец п = п % k конец если k == 1 тогда возвращаться т еще возвращаться 0 конецконец
Пример расчетов
Символ Лежандра (а/п) определено только для нечетных простых чисел п. Он подчиняется тем же правилам, что и символ Якоби (т. Е. Взаимность и дополнительные формулы для (−1/п) и (2/п) и мультипликативность «числителя».)
Проблема: учитывая, что 9907 - простое число, вычислите (1001/9907).
Использование символа Лежандра
Использование символа Якоби
Разница между этими двумя вычислениями заключается в том, что при использовании символа Лежандра «числитель» должен быть разложен на простые степени, прежде чем символ будет перевернут. Это делает вычисление с использованием символа Лежандра значительно медленнее, чем вычисление с использованием символа Якоби, поскольку нет известного алгоритма полиномиального времени для факторизации целых чисел.[4] Фактически, именно поэтому Якоби ввел этот символ.
Проверка на первичность
Есть еще одна причина различия символов Якоби и Лежандра. Если Критерий Эйлера формула используется по модулю составного числа, результат может быть или не быть значением символа Якоби, а на самом деле может даже не быть -1 или 1. Например,
Итак, если неизвестно, п простое или составное, мы можем выбрать случайное число а, вычислить символ Якоби (а/п) и сравните ее с формулой Эйлера; если они отличаются по модулю п, тогда п составной; если они имеют одинаковый остаток по модулю п для множества различных значений а, тогда п "вероятно, премьер".
Это основа для вероятностного Тест на простоту Соловея – Штрассена и усовершенствования, такие как Тест на простоту Baillie-PSW и Тест на простоту Миллера – Рабина.
В качестве косвенного использования его можно использовать как процедуру обнаружения ошибок во время выполнения Тест на простоту Лукаса-Лемера что даже на современном компьютерном оборудовании может занять недели при обработке Числа Мерсенна над (наибольшее известное простое число Мерсенна на декабрь 2018 г.). В номинальных случаях символ Якоби:
Это также верно для окончательного остатка и, следовательно, может использоваться как проверка вероятной действительности. Однако, если в оборудовании возникает ошибка, существует 50% -ная вероятность того, что вместо этого результат станет 0 или 1 и не изменится с последующими условиями (если не возникнет другая ошибка и она не вернется к -1).
Смотрите также
- Символ Кронекера, обобщение символа Якоби на все целые числа.
- Символ остатка мощности, обобщение символа Якоби на вычеты более высоких степеней.
Примечания
- ^ Якоби, К. Дж. Дж. (1837). "Über die Kreisteilung und ihre Anwendung auf die Zahlentheorie". Bericht Ak. Wiss. Берлин: 127–136.
- ^ Ирландия и Розен стр. 56–57 или Леммермейер стр. 10
- ^ Коэн, стр. 29–31.
- ^ В числовое поле сито, самый быстрый из известных алгоритмов, требует
Рекомендации
- Коэн, Анри (1993). Курс вычислительной алгебраической теории чисел. Берлин: Springer. ISBN 3-540-55640-0.
- Ирландия, Кеннет; Розен, Майкл (1990). Классическое введение в современную теорию чисел (второе издание). Нью-Йорк: Springer. ISBN 0-387-97329-X.
- Леммермейер, Франц (2000). Законы взаимности: от Эйлера до Эйзенштейна. Берлин: Springer. ISBN 3-540-66957-4.
- Шаллит, Джеффри (Декабрь 1990 г.). «О худшем случае трех алгоритмов вычисления символа Якоби». Журнал символических вычислений. 10 (6): 593–61. Дои:10.1016 / S0747-7171 (08) 80160-5. Технический отчет по информатике PCS-TR89-140.
- Валле, Брижит; Леме, Шарли (октябрь 1998 г.). Анализ среднего случая трех алгоритмов вычисления символа Якоби (Технический отчет). CiteSeerX 10.1.1.32.3425.
- Эйкенберри, Шона Мейер; Соренсон, Джонатан П. (октябрь 1998 г.). «Эффективные алгоритмы вычисления символа Якоби» (PDF). Журнал символических вычислений. 26 (4): 509–523. CiteSeerX 10.1.1.44.2423. Дои:10.1006 / jsco.1998.0226.
внешняя ссылка
- Вычислить символ Якоби показывает этапы расчета.