3차원 곡선 L‑다항식 고속 복원 및 야코비안 연산 구현

초록

본 논문은 평면 사차식 형태의 차수 3 곡선에 대해, 기존 Harvey‑Sutherland 알고리즘이 제공하는 모듈러 $p$ 결과를 정수 계수의 $L$‑다항식으로 복원하는 새로운 Las Vegas 알고리즘을 제시한다. 야코비안 군 연산을 위한 완전 구현(naïve, typical, hybrid)을 구축하고, 평균 $O(p^{1/4+o(1)})$, 최악 $O(p^{1/2+o(1)})$ 시간 복원을 보이며, 하이퍼엘립틱 차수 3 곡선에도 적용 가능함을 보인다.

상세 분석

이 연구는 두 단계로 구성된 Harvey의 평균 다항식 시간 알고리즘을 실용적으로 완성한다. 첫 번째 단계는 이미 구현된 CHS23 코드를 이용해 $\bmod p$ 로 감소된 $L_p(T)$ 를 $O(p^{1/2}\log^2p)$ 시간에 얻는다. 두 번째 단계, 즉 정수 계수 $L_p(T)$ 로 “리프팅” 하는 과정은 기존에 구현된 방법이 없었으며, 저자는 이를 야코비안 군 연산을 활용한 Las Vegas 알고리즘으로 해결한다. 핵심 아이디어는 $L_p(T)\bmod p$ 로부터 곡선의 Jacobian $\operatorname{Jac}(C)(\mathbb{F}_p)$ 에서 무작위 점들을 생성하고, baby‑step‑giant‑step (BSGS) 기법을 적용해 Frobenius 고유값을 추정한 뒤, Hasse‑Weil 구속을 이용해 정수 다항식으로 복원하는 것이다.

야코비안 연산 구현은 크게 세 부분으로 나뉜다. 첫째, “naïve” 알고리즘은 모든 종류의 디바이저에 대해 일반적인 Mumford 표현 없이도 덧셈·스칼라곱을 수행한다. 여기서는 9점이 포함된 삼차곡선 $K$ 와 5점이 포함된 이차곡선 $Q$ 를 찾아 교차점을 계산하고, 선형대수 기반 보간과 이상적인(ideal) 보간을 상황에 맞게 전환한다. 둘째, “typical” 알고리즘은 FOR04에서 제시된 전형적인 디바이저(세 점이 비공선, 서로 다른 $x$ 좌표 등)를 대상으로 Mumford $(u,v)$ 쌍을 이용해 $O(\log p)$ 시간에 연산한다. 셋째, “hybrid” 구현은 전형적인 경우에 FOR04 방식을, 그렇지 않은 경우에 naïve 방식을 자동 전환한다. 이 전환은 사전 처리 단계에서 찾은 접선 $l_g$ 와 무한점 $D_\infty$ 로부터 정의된 표준 형태를 이용해 효율적으로 이루어진다.

또한, BSGS 에서 필요로 하는 “완전 해시”를 위해 저자는 디바이저가 선형 등가인 경우를 정확히 구분하는 정리(Prop 2.2)를 제시하고, 해시값을 다중집합이 아닌 경우에는 추가 점 $P$ 로 대체한다. 이는 군 원소를 유일하게 식별함으로써 충돌을 방지하고, 평균 $O(p^{1/4})$ 시간 복원을 가능하게 하는 핵심적인 설계이다.

복잡도 분석에서는, 최악 경우 $O(p^{1/2+o(1)})$ 를 보장하고, “무작위 곡선·무작위 소수”에 대한 휴리스틱 가정 하에 평균 $O(p^{1/4+o(1)})$ 를 증명한다. 전체 알고리즘을 $B$ 이하의 모든 좋은 소수에 적용하면, 기대 시간은 $O(B^{5/4+o(1)})$ (휴리스틱) 혹은 $O(B^{3/2+o(1)})$ (무조건)이며, 기존 구현들의 $O(B^2)$ 와 비교해 실질적인 속도 향상을 입증한다. 하이퍼엘립틱 차수 3 곡선에 대해서도 Cantor 알고리즘과 결합해 동일한 복원 절차가 적용 가능함을 보여준다.

실험 결과는 표 2·표 3에 요약되어 있는데, 특히 $p\approx10^9$ 수준에서도 평균 $0.3$초 내외의 복원 시간을 기록하며, 기존 Tuitman·Costa 방식보다 1~2자리 빠른 성능을 보인다. 구현 코드는 공개 저장소(