평면 델로네 삼각분할의 형식적 증명

초록

본 논문은 평면 상에서 Delaunay 삼각분할을 얻기 위한 알고리즘의 정밀한 형식 증명을 제시한다. 초기 삼각분할에서 비 Delaunay 가장자리를 반복적으로 플립하여 Delaunay 성질을 만족할 때까지 진행한다. 하이퍼맵을 이용한 조합적 모델링과 좌표 임베딩을 결합해 삼각분할·임베딩·플립 연산의 불변성을 보이고, 유한 집합 위의 비순환 관계가 잘 정의됨을 이용해 종료성을 증명한다.

상세 분석

이 논문은 평면 Delaunay 삼각분할 알고리즘을 형식적으로 검증하기 위해 두 가지 핵심 기술을 도입한다. 첫 번째는 하이퍼맵(hypermap)이라는 조합적 구조를 활용해 삼각분할을 순수히 이산적인 데이터 구조로 모델링한다는 점이다. 하이퍼맵은 darts, permutations σ₁, σ₂, σ₃ 로 정의되며, 각각 정점, 변, 면에 대응한다. 논문에서는 기존에 개발해 온 하이퍼맵 사양 프레임워크를 확장해, 각 dart에 2차원 좌표 (x, y)를 부착함으로써 기하학적 임베딩을 정의한다. 이때 좌표 부착은 “consistent embedding” 조건을 만족해야 하는데, 이는 동일한 정점에 연결된 모든 dart이 동일한 좌표를 공유하고, 인접한 면이 서로 겹치지 않으며, 모든 변이 직선 구간으로 표현된다는 것을 의미한다.

두 번째 핵심은 Delaunay 성질을 판별하고 위반된 변을 교정하는 플립 연산이다. 논문은 “legal Delaunay edge”와 “illegal Delaunay edge”를 정확히 정의한다. 법적 변은 그 양쪽 삼각형의 외접원( circumcircle )이 서로 겹치지 않는 경우이며, 불법 변은 외접원 내부에 반대쪽 정점이 존재하는 경우이다. 플립 연산은 불법 변을 제거하고, 그 변을 대각선으로 교체해 새로운 두 삼각형을 만든다. 중요한 것은 이 플립이 하이퍼맵 구조를 보존한다는 증명이다. 즉, σ₁, σ₂, σ₃ 의 순열 관계가 플립 전후에 동일하게 유지되며, 임베딩 조건도 좌표 재배치를 통해 그대로 만족한다.

알고리즘의 종료성을 보이기 위해 저자는 “non‑cyclic relation on a finite set is well‑founded” 라는 일반적인 수학적 원리를 적용한다. 구체적으로, 각 플립은 전체 변 집합에 대해 어떤 정수값(예: 전체 외접원 반지름의 합 또는 변의 길이 사전순)을 감소시킨다. 이 감소 관계가 순환을 포함하지 않으며, 변의 수가 유한하므로 반드시 최소값에 도달한다. 따라서 플립 과정은 유한 단계 내에 종료하고, 최종 상태는 모든 변이 legal Delaunay edge 를 만족하므로 Delaunay 삼각분할이 된다.

논문은 또한 형식 검증 도구(Coq, Isabelle 등)와의 연계 가능성을 논의한다. 하이퍼맵 사양을 타입 이론 기반으로 기술하고, 플립 연산과 종료성 증명을 자동화된 증명 보조기구로 검증할 수 있음을 시사한다. 이는 기존의 실험적 Delaunay 구현이 갖는 “정확성에 대한 의심”을 해소하고, 컴퓨터 과학 이론과 계산기하학 사이의 교량을 놓는 중요한 기여라 할 수 있다.