베트위니스 중심성의 효율적 계산을 위한 병렬·대수 알고리즘과 복잡도 하한

본 논문은 네트워크 분석에서 가장 널리 쓰이는 정점 중요도 지표인 베트위니스 중심성(Betweenness Centrality, BC)의 효율적 계산을 목표로, 두 가지 새로운 접근법과 그에 대한 복잡도 하한을 제시한다. 서론에서는 BC가 사회·통신·생물학·교통 등 다양한 분야에서 어떻게 활용되는지를 설명하고, 기존에 널리 쓰이는 Brandes 알고리즘이 “경로‑비교 기반” 모델에 속함을 밝힌다. 이어서 저자들은 이 모델에 대해 Ω(nm) 이하의 시간으로는 정확한 BC를 구할 수 없다는 하한을 증명한다. 이 증명은 Karger·Koller·Phillips가 제시한 Ω(n³) APSP 하한을 변형해, 3‑partite 그래프 구조와 가중치 조작을 이용해 경로 비교가 없을 경우 베트위니스 값이 변함없이 유지되는 모순을 도출함으로써 이루어진다. 결과적으로, Brandes 알고리즘의 O(nm) 시간 복잡도는 이 모델 내에서 최적에 가깝다는 결론을 얻는다. 다음 섹션에서는 두 가지 새로운 알고리즘을 제시한다. 첫 번째는 “대수적 방법”으로, 행렬 연산을 활용해 BC를 계산한다. 무가중 그래프에 대해서는 Seidel의 O(n^ω·log Diam) APSP 알고리즘을 변형해, 인접 행렬 A를 반복 제곱하면서 각 거리 레벨 l에서 처음 등장하는 비제로 원소를 최단 경로 개수 행렬 Λ에 기록한다. 이 과정은 그래프의 지름(Diam)만큼 반복하면 되므로 전체 복잡도는 O(n^ω·Diam(G))가 된다. 가중 그래프(정수 가중치 1…M)에서는 거리 곱(min‑plus product) 연산을 이용해 O(M·n^ω·Diam(G)) 시간에 거리와 최단 경로 개수를 동시에 얻는다. 여기서 ω는 현재 알려진 행렬 곱셈 지수(≈2.376)이다. 두 번째 단계는 “역방향 패스”로, Brandes 정리의 재귀식을 행렬 형태로 전개한다. 거리 레벨 l을 역순으로 순회하면서 Δ_{l‑1}=((D_l+Δ_l)÷Λ)·A 를 계산한다. D_l은 거리 l에 해당하는 0‑1 마스크 행렬, Δ_l은 레벨 l에서의 의존도 행렬이며, ÷는 원소별 나눗셈을 의미한다. 이 연산은 각 레벨마다 O(n^ω) 시간이 소요되고, 레벨 수는 Diam(G) 이하이므로 전체 역방향 패스도 O(n^ω·Diam(G))에 수행된다. 최종적으로 Δ 행렬을 모두 합산하면 모든 정점의 BC 값을 O(n^2) 추가 시간에 얻을 수 있다. 병렬 구현 부분에서는 OR‑CRCW PRAM 모델을 가정한다. 무가중 그래프에 대해 Ullman·Yannakakis의 무작위 병렬 BFS를 사용해 전방 패스를 O(n) 시간에 O(m log n) 프로세서로 수행한다. 가중 그래프(정수 가중치 1…M)에서는 O(n log² n log M) 시간에 O(m) 프로세서를 사용한다. 후방 패스는 기본적으로 O(n²) 시간이 필요하지만, 정점 차수가 제한된 경우 O(n log m) 시간에 O(m) 프로세서로 최적화한다. 이러한 설계는 실제 대규모 소셜·통신 네트워크에서 직렬 알고리즘보다 훨씬 높은 스루풋을 기대하게 한다. 마지막으로 저자들은 “단일 정점의 베트위니스 계산은 전체 정점 집합에 대한 계산보다 적어도 동등히 어렵다”는 가설을 제시하고, APSP와의 복합적 관계를 통해 이 가설이 사실일 가능성을 논증한다. 전체적으로 이 논문은 (1) 경로‑비교 기반 알고리즘에 대한 Ω(nm) 하한, (2) 행렬 곱셈을 이용한 O(n^ω·Diam)·O(M·n^ω·Diam) 시간 알고리즘, (3) 무작위 병렬 SSSP 기반의 실용적 병렬 구현이라는 세 축에서 베트위니스 중심성 계산의 이론적·실용적 한계를 크게 확장한다.