[논문 리뷰] An O(M(n) log n) algorithm for the Jacobi symbol

Richard P. Brent, Paul Zimmermann|arXiv (Cornell University)|2010. 04. 13.

Algorithms and Data Compression참고 문헌 2인용 수 4

한 줄 요약

이 논문은 연속 분수나 가우스 항등식에 의존하지 않고 이진 재귀 GCD 방법을 사용하여 두 개의 $n$-비트 정수의 야코비 승자를 계산하는 $O(M(n)\text{log }n)$ 알고리즘을 제시한다. 이 방법은 몫 수열을 효율적으로 처리함으로써 비제곱 시간 복잡도를 달성하며, 약 10,000자리 십진수보다 큰 입력에서 GMP의 제곱 시간 구현보다 빠른 것으로 나타났다.

ABSTRACT

The best known algorithm to compute the Jacobi symbol of two n-bit integers runs in time O(M(n) log n), using Schönhage's fast continued fraction algorithm combined with an identity due to Gauss. We give a different O(M(n) log n) algorithm based on the binary recursive gcd algorithm of Stehlé and Zimmermann. Our implementation - which to our knowledge is the first to run in time O(M(n) log n) - is faster than GMP's quadratic implementation for inputs larger than about 10000 decimal digits.

연구 동기 및 목표

• 두 개의 $n$-비트 정수의 야코비 승자를 계산하기 위한 비제곱 시간 알고리즘을 개발하기.
• 슐로닝게의 연속 분수 알고리즘 또는 가우스 항등식에 의존하지 않도록 하여 더 단순한 대안을 제공하기.
• 오직 이진 GCD 기법만을 사용하여 $O(M(n)\log n)$ 시간 복잡도를 달성하기.
• 알고리즘을 구현하고 벤치마크를 수행하여 큰 입력에서 기존의 제곱 시간 구현보다 성능 향상을 보여주기.
• 고성능 수론 계산을 위한 증명 가능하게 효율적이고 실용적인 방법을 제공하기.

제안 방법

• 스테플레와 지머만의 이진 재귀 GCD 알고리즘을 변형하여 야코비 승자를 직접 계산하도록 적용한다.
• 양수 몫을 사용하는 이진 나눗셈을 통해 양수성을 유지하고 종료를 보장한다.
• 양수 몫을 사용함에도 불구하고 종료를 보장하기 위해 수정된 정지 기준을 도입한다.
• 항등식을 이용한 야코비 승자의 기호적 변환을 수행한다: $(bc|a)=(b|a)(c|a)$, $(2|a)=(-1)^{(a^2-1)/8}$, 및 홀수 $a,b$에 대해 $(b|a)=(-1)^{(a-1)(b-1)/4}(a|b)$.
• 반복을 '좋음', '나쁨', 또는 '나쁜'으로 분류하고, '나쁜' 반복 수열을 하나의 '무해한' 단계로 통합하여 최악의 경우 복잡도를 $O(M(n)\log n)$으로 감소시킨다.

실험 결과

연구 질문

• RQ1 연속 분수나 가우스 항등식에 의존하지 않고 오직 이진 GCD 기법만을 사용하여 야코비 승자를 $O(M(n)\log n)$ 시간 내에 계산할 수 있는가?
• RQ2 야코비 승자를 위한 이진 GCD 기반 알고리즘의 최악의 경우 시간 복잡도는 무엇이며, 이를 $O(M(n)\log n)$로 줄일 수 있는가?
• RQ3 새로운 알고리즘의 성능은 큰 입력에서 기존의 GMP의 제곱 시간 알고리즘과 비교해 어떻게 되는가?
• RQ4 제안된 알고리즘에서 평균적으로 한 반복당 비트 수축 비율은 얼마이며, 표준 이진 GCD와 비교해 어떻게 되는가?
• RQ5 알고리즘이 실용적으로 효율적으로 구현될 수 있으며, 어느 입력 크기에서 기존 방법을 능가하는가?

주요 결과

• 제안된 알고리즘은 $O(M(n)\log n)$ 시간 복잡도를 달성하여 최고의 이론적 한계와 일치한다.
• 이 알고리즘은 슐로닝게의 연속 분수 알고리즘 또는 가우스 항등식을 사용하지 않고도 이 한계를 달성한 최초의 알고리즘이다.
• 구현 결과, 약 10,300자리 십진수(34,240비트)보다 큰 입력에서 GMP의 제곱 시간 구현을 능가한다.
• 100만 개의 64비트 단어 입력에 대해 약 83.1초와 4750만 번의 반복이 소요되었으며, 이론적 비트 수축 비율은 1회 반복당 0.651993비트였다.
• 실험 결과, 약 53.33%의 반복이 '좋음', 26.67%는 '나쁨', 20%는 '무해한' 것으로 나타나 히ュ리스틱 분포로 수렴함을 보였다.
• $O(n^2)$ GMP 구현과 $O(M(n)\log n)$ 구현 간의 교차점은 약 535단어(약 34,240비트 또는 10,300자리 십진수)에서 발생한다.