루카스 수를 위한 O(log n) 중간·리플 알고리즘

초록

본 논문은 루카스 수열의 개별 항을 O(log n) 시간에 계산할 수 있는 두 가지 새로운 알고리즘, ‘Middle’와 ‘Ripple’를 제시한다. 기존의 피보나치 기반 방법보다 적은 연산으로 루카스 수를 직접 구할 수 있음을 증명하고, 복잡도와 상수 요인 측면에서 선행 연구보다 효율적임을 보인다.

상세 분석

논문은 루카스 수열 Lₙ을 정의하고, 기존 O(n) 반복 방식과 O(log n) 피보나치 기반 방법의 한계를 지적한다. 핵심 아이디어는 “이중화(doubling) 공식”을 활용해 n을 2진수로 분해하면서 필요한 연산을 로그 단계로 축소하는 것이다. 저자는 L₂k = L₂k−(−1)ᵏ·2 라는 기본 식을 도출하고, 이를 변형해 L₂k+2 = L₂k+1 + (−1)ᵏ·2 형태로 만든다. 이 식을 바탕으로 두 개의 변수 LL, LM, LH를 사용해 한 번의 반복에서 두 번의 제곱 연산과 부호 조정만으로 다음 단계 값을 얻는다. ‘Middle’ 알고리즘은 입력 n을 이진수로 변환한 뒤, 각 비트에 따라 p = ±1을 설정하고 위 식을 순차적으로 적용한다. 이 과정에서 필요 없는 전후 조건(pre/post‑condition)을 없애고, 상수 연산 수를 최소화한다는 점이 특징이다. 복잡도 분석에서는 각 비트당 두 번의 제곱, 두 번의 부호 연산, 세 번의 덧셈이 수행되므로 전체 연산량은 O(log n)이며, 실제 구현 시 상수 계수가 기존의 피보나치 기반 O(log n) 알고리즘보다 작다.

‘Ripple’는 동일한 이중화 원리를 재귀 형태로 구현한다. 짝수 n에 대해서는 Lₙ = L_{n/2}² − p·2 를, 홀수 n에 대해서는 Lₙ = L_{⌈n/2⌉}² − L_{⌊n/2⌋}² + p·4 로 정의한다. 여기서 p는 n/2의 짝/홀에 따라 ±1로 결정된다. 재귀 깊이는 로그 수준이며, 각 호출마다 제곱 연산 하나와 부호 연산 하나만 필요하므로 메모리 사용량도 최소화된다.

두 알고리즘 모두 Lₙ을 직접 계산함으로써 Lₙ = F_{n−1}+F_{n+1} 와 같은 피보나치 변환을 거칠 필요가 없으며, √5 를 이용한 변환을 통해 피보나치 수를 구하는 경우보다 약 2배 정도 빠른 상수 요인을 제공한다. 논문은 또한 기존 연구