최대 균형 이분 그래프 완전쌍 찾기 효율적 정확 알고리즘

초록

본 논문은 최대 균형 이분 클리크(균형 이분 완전쌍) 문제에 대해 새로운 상한 전파(Upper Bound Propagation) 기법을 제안하고, 이를 기존의 분기‑한정(B&B) 알고리즘 BBClq에 통합한 ExtBBClq와 부분 집합만 열거하는 ExtUniBBClq 두 가지 정확 알고리즘을 설계한다. 또한 상한값을 이용한 새로운 유효 부등식들을 도입해 MIP 모델을 강화한다. 실험 결과, 제안된 알고리즘과 강화된 MIP 모델이 무작위 그래프와 실제 대규모 데이터셋에서 기존 방법보다 현저히 빠른 해결 시간을 보였다.

상세 분석

본 연구는 최대 균형 이분 클리크 문제(MBBP)의 NP‑hard 특성에도 불구하고, 실용적인 정확 해법을 제공하기 위해 세 가지 핵심 아이디어를 결합한다. 첫 번째는 각 정점이 포함될 수 있는 최대 균형 이분 클리크의 절반 크기에 대한 상한값을 사전 계산하는 Upper Bound Propagation(UBP) 절차이다. UBP는 세 가지 수학적 명제를 기반으로 한다. 명제 1은 정점의 차수가 바로 상한이 된다는 단순하지만 기본적인 관찰이며, 명제 2는 두 정점이 공유하는 이웃의 개수를 이용해 더 강력한 상한 y_v 를 도출한다. 여기서 w_{vu}=|N(v)∩N(u)| 를 계산하고, 최소 y_v 개수의 정점들이 w_{vu}≥y_v 를 만족하도록 하는 최대 y_v 를 이진 탐색으로 찾는다. 명제 3은 정점 v의 이웃 중 상한값이 최소 z_v 이상인 정점이 z_v 개 이상 존재할 때, z_v 가 v가 포함될 수 있는 균형 클리크의 절반 크기에 대한 상한이 된다는 것을 증명한다. 이 세 명제를 순차적으로 적용하고, 명제 3에 기반한 전파 과정을 반복함으로써 ub_v 를 점진적으로 감소시킨다. 전파는 ub_v 가 0 이하로 내려가지 않으며, 실험적으로는 수십 번의 반복만에 수렴한다.

두 번째 기여는 UBP에서 얻은 정밀 상한을 기존 BBClq 알고리즘에 통합한 ExtBBClq이다. BBClq는 후보 집합 C_A, C_B 를 유지하며 재귀적으로 A와 B를 교대로 확장하는 전통적인 B&B 방식이다. ExtBBClq는 분기 시점에 ub_{A∪{v}} 를 확인해 현재 전역 하한(lb)보다 작으면 해당 분기를 즉시 차단한다. 이는 기존 BBClq가 색칠 기반 상한을 사용해 비효율적으로 많은 불필요한 분기를 탐색하던 문제를 크게 완화한다. 또한, UBP는 정점 별 상한을 사전 계산하므로 탐색 중 추가 연산이 거의 필요하지 않아 전체 탐색 비용을 크게 감소시킨다.

세 번째 기여는 UBP 기반 상한을 활용한 새로운 유효 부등식들을 MIP 모델에 추가함으로써 모델을 강화한 점이다. 기존 MIP 모델은 보완 그래프에 대한 변수와 제약을 사용했으나, 상한값을 직접 변수 ub_v 로 도입하고, ub_v ≥ x (x는 목표 절반 크기) 와 같은 선형 부등식을 추가한다. 또한, 명제 2와 3에서 유도된 집합 관계를 이용해 “정점 v가 선택되면, 그 이웃 중 ub_u ≥ x 인 정점이 최소 x 개 존재한다”는 형태의 제약을 삽입한다. 이러한 강화는 CPLEX와 같은 상용 솔버가 초기 탐색 단계에서 더 강력한 컷을 생성하도록 유도하여, 특히 밀도가 높은 대규모 그래프에서 해결 시간을 크게 단축시킨다.

마지막으로, ExtUniBBClq 라는 새로운 탐색 스킴을 제안한다. BBClq가 모든 가능한 균형 클리크를 완전 탐색하는 반면, ExtUniBBClq는 반집합(half‑set)만을 열거한다. 즉, A 혹은 B 중 하나만을 확장하면서 ub_v 로 가능한 최대 절반 크기를 미리 판단하고, 해당 반집합이 목표 크기 이하이면 더 이상 확장을 중단한다. 이 방식은 탐색 트리의 폭을 크게 줄이며, 특히 목표 절반 크기가 작을 때 효율이 극대화된다.

실험에서는 무작위 밀도 0.10.9 범위의 그래프와, 소셜 네트워크, 생물학적 네트워크 등 실제 대규모 데이터(수만수십만 정점)를 대상으로 비교하였다. 결과는 ExtBBClq와 ExtUniBBClq가 기존 BBClq 대비 평균 3배~10배 빠른 해결 시간을 보였으며, 강화된 MIP 모델은 CPLEX가 동일 인스턴스를 해결하는 데 소요되는 시간의 40% 이하로 감소하였다. 특히 정점 수가 30,000을 초과하는 경우에도 UBP가 몇 차례의 전파만으로 수렴해 전처리 비용이 무시할 수준임을 확인하였다.

요약하면, 본 논문은 상한 전파라는 독창적인 전처리 기법을 통해 정점 별 상한을 정밀히 추정하고, 이를 기존 B&B와 MIP 두 갈래의 정확 해법에 효과적으로 통합함으로써 MBBP 해결에 있어 현존 최고의 성능을 달성하였다.