다수의 점을 포함하는 원 찾기의 알고리즘적 접근

초록

본 논문은 평면상의 n개의 점 집합에서, 두 점을 선택했을 때 그 두 점을 포함하는 모든 원이 최소한 cn개의 점을 포함하도록 하는 쌍을 효율적으로 찾는 여러 알고리즘을 제시한다. 무작위 O(n log n) 알고리즘, O(n²) 정확도 향상 알고리즘, 그리고 볼록 위치와 지오데식 디스크 상황에 대한 선형 시간 알고리즘을 포함한다. 또한 두 색상 점 집합과 단순 다각형 내부 점에 대한 일반화도 다룬다.

상세 분석

논문은 먼저 기존 이론적 결과를 정리하고, “C_S(p,q)”와 “˜C_S(p,q)”라는 두 점 쌍에 대한 최소 포함 점 수를 정의한다. 여기서 C_S(p,q)는 p와 q를 포함하는 원 중 가장 작은 점 개수를, ˜C_S(p,q)는 원 안·밖 양쪽에서 최소 점 수를 고려한 대칭값이다. 저자는 이 값을 효율적으로 계산하기 위해 두 점을 잇는 이등분선(bisector) 위에 존재하는 원의 중심들을 정렬하고, 인접 구간 사이의 무게 차이를 ±1로 업데이트하는 O(n log n) 알고리즘(Algorithm 1)을 제시한다. 이 과정은 각 쌍에 대해 O(n log n) 시간으로 C_S와 ˜C_S를 구할 수 있게 한다.

다음으로 무작위 선택 전략을 도입한다. Ramos와 Viana의 비례식 S_k ≤ 3(k+1)n – 3(k+1)(k+2)를 이용해, α∈(0,1)일 때 전체 쌍 중 적어도 6α·(n choose 2)개의 쌍이 ˜C_S(p,q) ≥ (½ – √(1+2α)/12)·n 를 만족함을 증명한다. 따라서 무작위로 쌍을 선택하고 위의 O(n log n) 검증을 반복하면, 기대 시간 O(n log n) 안에 원하는 쌍을 찾을 수 있다. 성공 확률을 1–1/n로 높이려면 O(n log² n) 시간이 필요하다.

두 번째 알고리즘은 모든 쌍을 전부 검토하면서 ˜C_S(p,q)를 계산하는 O(n²) 방법이다. 여기서는 각 점을 기준으로 다른 점들과의 이등분선을 한 번만 정렬하고, 스위핑을 통해 모든 쌍의 최소 무게를 동시에 업데이트한다. 이 과정에서 얻어지는 최적 상수는 ½ – 1/√12 ≈ 0.224로, 기존 이론적 하한인 n/4·7≈0.25에 근접한다. 또한 이 알고리즘을 확장해 C_S(p,q)를 최대화하는 쌍을 O(n² log n) 시간에 찾을 수 있다.

볼록 위치(convex position)에서는 더 강력한 결과가 가능하다. 저자는 볼록 껍질의 순환 구조를 이용해, 인접한 두 꼭짓점 사이의 쌍을 검사하면 언제나 원을 통과하는 모든 점이 최소 n/3개임을 보인다. 이 검사는 단순히 한 번의 선형 스캔으로 수행되므로 전체 시간 복잡도는 O(n)이다. 동일한 아이디어를 직경 디스크(두 점을 지름으로 하는 원) 문제에도 적용해, 직경 디스크가 최소 n/3개의 점을 포함하도록 하는 쌍을 선형 시간에 찾는다.

색상 구분이 있는 경우(red‑blue)에도 알고리즘을 확장한다. 전체 점을 빨강·파랑으로 절반씩 나누고, 빨강‑파랑 쌍만을 대상으로 위의 무작위·전수 검증 절차를 적용한다. 이때도 동일한 상수 c가 보장되며, 색상 제약이 추가된 상황에서도 O(n log n) 기대 시간 혹은 O(n²) 최악 시간 복잡도를 유지한다.

마지막으로 단순 다각형 내부 점 집합에 대한 일반화도 제시한다. 여기서는 유클리드 거리 대신 지오데식 거리를 사용하고, 이등분선이 곡선(선분·쌍곡선)으로 변하는 점을 고려한다. 저자는 기존 유클리드 알고리즘의 핵심 아이디어—이등분선 위의 중심 정렬과 무게 업데이트—를 지오데식 환경에 맞게 변형함으로써, 동일한 시간 복잡도와 상수를 얻을 수 있음을 보인다. 특히, 다각형 경계와의 교차 처리, 그리고 지오데식 원의 존재 여부를 판단하는 절차가 추가되지만, 전체 구조는 크게 변하지 않는다.

전체적으로 이 논문은 “어떤 두 점을 선택하더라도 그 쌍을 포함하는 모든 원이 충분히 많은 점을 포함한다”는 존재론적 결과를 실제 알고리즘으로 구현하는 데 성공했으며, 무작위·전수·특수 위치·색상·지오데식 등 다양한 상황을 포괄한다는 점에서 이론과 실용성을 동시에 만족한다.