[논문 리뷰] Computing $\pi(x)$ Analytically

David J. Platt|arXiv (Cornell University)|2012. 03. 26.

Analytic Number Theory Research참고 문헌 10인용 수 9

한 줄 요약

이 논문은 리만의 명시적 공식과 간격 산술을 사용하여 x = 10^24까지 π(x), 즉 소수 카운팅 함수를 조건 없이 엄밀하게 분석적 알고리즘으로 계산하는 방법을 제시한다. 멕스웰 변환 쌍과 오차 함수 스무딩을 활용함으로써 기존 방법보다 개선되었으며, 검증된 오차 한계와 간격 산술을 통해 수치적 안정성과 정확성을 확보하여 π(10^24)의 조건 없는 첫 번째 계산을 달성한다.

ABSTRACT

We describe a rigorous implementation of the Lagarias and Odlyzko Analytic Method to evaluate the prime counting function and its use to compute unconditionally the number of primes less than $10^{24}$.

연구 동기 및 목표

• 기존의 조합적 걸러내기 기법으로는 도달할 수 없는 범위까지 π(x)를 엄밀하고 조건 없는 방법으로 계산하기 위한 방법 개발.
• 고정밀 분석적 수론 계산에서 수치적 불안정성과 부동소수점 오차 누적 문제를 극복하기 위한 방법.
• 간격 산술과 명시적 오차 한계를 사용하여 라자리우스-오들리즈코 분석적 방법의 변종을 구현하고 검증하기.
• 리만 가설을 가정하지 않고도 π(10^24)를 조건 없이 계산하여 계산적 소수 이론 분야에서 오랫동안 남아 있던 문제를 해결하기.
• 미래의 고정밀 π(x) 평가를 위한 재현 가능하고 검증 가능하며 수치적으로 안정된 프레임워크 제공.

제안 방법

• 페론의 공식을 적용하여 π*(x)를 로그 ζ(s)의 복소 적분과 오차 함수에서 유도된 스무딩 커널 함수 φ(t)의 곱으로 표현한다.
• 갈웨이가 정의한 멕스웰 변환 쌍을 사용: φ(t) = ½ erfc( log(t/x)/(√2 λ) ) 및 φ̂(s) = x^s / s ⋅ exp(λ²s²/2), 수렴성과 국소화 사이의 균형을 확보한다.
• 전체 계산 과정에서 부동소수점 오차를 엄밀하게 경계하기 위해 간격 산술을 사용하여 히우리스틱 오차 추정을 대체한다.
• 비자명한 제타 함수 근의 합을 잘라내는 데 발생하는 오차에 대해 리만-폰 망골트 함수 N(t)과 尾 추정을 사용하여 명시적 경계를 적용한다.
• 제로 합 잘라내기의 정밀한 제어를 가능하게 하기 위해 레마 A.1을 사용하여 알려진 점근 전개와 오차 항을 이용해 N(t)를 경계한다.
• 비자명한 제로가 [T₁, T₂] 내부와 외부에 있을 경우 각각 E₁과 E₂의 별도 오차 경계를 도입하며, 스티엘제스 적분과 지수 감쇠 추정을 사용한다.

실험 결과

연구 질문

• RQ1 라자리우스-오들리즈코 분석적 방법을 간격 산술과 엄밀한 오차 경계를 사용하여 조건 없는 방식으로 π(x)를 대규모로 계산할 수 있는가?
• RQ2 고정밀 분석적 수론 계산에서 수치적 불안정성과 부동소수점 오차 누적 문제를 어떻게 제어할 수 있는가?
• RQ3 오차를 최소화하면서도 계산 가능성을 유지하기 위해 스무딩 파라미터 λ와 잘라내기 임계값 T₁, T₂ 사이의 최적 균형은 무엇인가?
• RQ4 기존에 알려진 근과 검증되지 않은 근 위치를 모두 고려한 명시적 오차 경계를 사용하여 비자명한 제타 함수 근의 합을 잘라내는 것이 가능한가?
• RQ5 리만 가설을 가정하지 않고도 분석적 방법과 검증된 수치 기법을 사용하여 π(10^24)를 조건 없이 계산할 수 있는가?

주요 결과

• 논문은 리만 가설을 가정하지 않고도 π(10^24)를 성공적으로 조건 없이 계산하였으며, 이는 이 규모에서 처음으로 이루어진 계산이다.
• 구현은 간격 산술을 사용하여 부동소수점 오차를 엄밀하게 제어하여 하드웨어나 컴파일러 이상 상황에서도 정확성을 보장한다.
• 시간 복잡도는 O(x^{1/2+ε})로, 대규모 x에 대해 조합적 걸러내기 기법의 O(x^{2/3} log^{-2}x)보다 크게 향상되었다.
• 리만-폰 망골트 함수 N(t)를 사용하여 제로 합 잘라내기의 명시적 오차 경계를 유도하였으며, 레마 A.1을 통해 제로 분포에 대한 정밀한 제어가 가능했다.
• 기각된 제로에서 온 오차 기여는 E₁과 E₂로 경계되었으며, E₂는 실수부가 ½이 아니어도 커널의 지수 감쇠 덕분에 무시할 만큼 작다.
• 최종 결과는 여러 독립적 검증을 통해 확인되었으며, 더 작은 척도에서의 π(x) 값과의 교차 검증 및 소수 정리와의 일관성으로 검증되었다.