안정 매칭 보로노이 다이어그램: 복잡도와 효율적 알고리즘

초록

본 논문은 평면상의 n개의 점 사이트와 각 사이트에 할당된 면적(식욕) 제한을 고려한 “안정 매칭 보로노이 다이어그램”의 조합론적 복잡도와 구성 알고리즘을 연구한다. 최악의 경우 면과 변의 수는 Θ(n²)에 달함을 보이며, O(n²+ε) 상한을 증명한다. 또한, 정확한 대수적 계산이 불가능한 기하학적 기본 연산을 포함하는 O(n³ log n + n² f(n)) 시간 알고리즘을 제시하고, 다각형 볼록 거리 함수에 대해서는 해당 기본 연산을 정확히 수행할 수 있음을 보인다.

상세 분석

이 논문은 기존 보로노이 다이어그램과 안정 매칭 이론을 결합해, 각 사이트가 사전에 지정된 면적(앱피트)만큼만 점들을 할당받는 새로운 평면 분할 구조를 정의한다. 핵심 아이디어는 “원 성장” 과정을 3차원에서 높이‑시간 축으로 해석해, 각 사이트를 정점으로 하는 원뿔을 동일한 속도로 확장시키고, 원뿔이 자신의 할당 면적에 도달하면 성장을 멈추는 방식이다. 이 시각화는 기존 보로노이 다이어그램이 무한히 뻗는 원뿔들의 하부 외피와 달리, 제한된 높이의 원뿔 집합의 하부 외피가 바로 안정 매칭 보로노이 다이어그램임을 보여준다.

조합론적 분석에서는 먼저 각 사이트 s에 대해 “경계 원” B_s 를 정의한다. B_s는 s를 중심으로 그 사이트가 할당받는 모든 점을 포함하는 최소 원이며, 논문은 모든 비공백 면의 합집합이 바로 이러한 경계 원들의 합임을 증명한다(Lemma 3). 이를 기반으로, 직선 경계는 두 사이트의 수직 이등분선에, 곡선 경계는 해당 사이트의 경계 원 위에 놓인다는 성질을 도출한다(Lemma 4). 이러한 구조적 제약을 이용해, 각 경계 원이 O(n)개의 교차점을 가질 수 있음을 보이고, 전체 면·변·정점 수는 O(n²+ε) (임의의 ε>0) 로 상한을 잡는다. 또한, 동일한 앱피트를 갖는 n개의 사이트를 적절히 배치하면 Θ(n²)개의 면과 변이 발생함을 구성 예시를 통해 제시해, 상한이 거의 최적임을 확인한다.

알고리즘적 측면에서는, 이 다이어그램을 정확히 계산하려면 “원뿔 높이와 면적 사이의 비선형 관계”를 해결해야 하는데, 이는 일반적인 대수 연산만으로는 표현할 수 없음을 증명한다(Section 4). 따라서 저자들은 “geometric primitive”이라 명명한 연산을 별도 서브루틴으로 정의하고, 이를 수치적으로 근사하거나, 다각형 볼록 거리 함수와 같이 제한된 거리 메트릭에서는 정확히 계산할 수 있는 방법을 제시한다. 최종 알고리즘은 각 단계에서 현재 가장 작은 원뿔 높이를 찾고, 해당 원뿔을 완성시키는 과정을 n번 반복하며, 전체 복잡도는 O(n³ log n + n² f(n))이다. 여기서 f(n)은 위에서 정의한 기본 연산의 실행 시간이다.

마지막으로, 다각형 볼록 거리 함수(예: L₁, L_∞ 등)에서는 원뿔의 단면이 다각형이 되므로, 경계 원을 다각형 교차 연산으로 정확히 구할 수 있어, 알고리즘을 완전하게 구현할 수 있음을 보인다. 이는 유클리드 거리에서는 불가능하지만, 실용적인 응용(예: 도시 계획, 정치적 구역 설정)에서 중요한 의미를 가진다.