모노 비용 하의 운송 문제와 위치 최적화 적용

초록

본 논문은 Monge 성질을 만족하는 비용 행렬을 갖는 균형 운송 문제에 대해 북서 코너(North‑West‑Corner) 알고리즘이 생성하는 기본 해의 구조를 이용해 최적 이중 변수의 압축식을 유도한다. 이를 바탕으로 비증가 가중치를 갖는 이산 순서 중위값 문제(Discrete Ordered Median Problem, DOMP)의 Benders 분해 모델에 새로운 절단식을 제시하고, 제안된 MILP 형식이 기존 최첨단 모델보다 계산 효율성과 견고성을 크게 향상시킴을 실험을 통해 입증한다.

상세 분석

논문은 먼저 Monge 비용 행렬이 만족하는 부등식 c_{ij}+c_{i’j’} ≤ c_{ij’}+c_{i’j} (1 ≤ i < i’ ≤ p, 1 ≤ j < j’ ≤ q) 을 정의하고, 이 조건 하에서 전통적인 북서 코너 규칙이 최적 원시 해를 생성함을 재확인한다. 알고리즘 1은 공급‑수요 쌍을 순차적으로 탐색하면서 “아래로 이동”(공급이 소진) 혹은 “오른쪽으로 이동”(수요가 충족)하는 스텝을 기록한다. 저자는 이 이동 경로를 수학적으로 A = {a₁,…,a_T}라 두고, 관찰 2.1과 명제 2.2를 통해 A가 정확히 공급 누적량이 수요 누적량을 초과하거나 이하인 전환점들의 집합임을 증명한다. 이 구조적 특성을 이용해 이중 변수 u_i, v_j 에 대한 최적값을 “각 단계에서 남은 공급·수요량의 차이”로 표현하는 압축식(식 (2)·(3))을 도출한다. 특히, u_i와 v_j는 각각 공급·수요의 누적 차이를 기준으로 단조적으로 증가·감소하는 형태를 띠어, 전통적인 단순 대입 방식보다 계산량을 크게 줄일 수 있다.

다음으로 이러한 이중 해 표현을 DOMP에 적용한다. DOMP는 비증가 가중치 w₁ ≥ w₂ ≥ … ≥ w_n 을 갖는 순서 중위값 목적함수를 최소화하는 시설 위치 문제이며, p개의 시설을 선택해 고객에게 할당한다. 기존 모델은 복잡한 순서화 변수와 큰 수의 제약식으로 인해 규모가 커질수록 계산이 급격히 악화된다. 논문은 Labbé et al.이 제시한 MILP 모델에 Benders 분해를 적용하고, 서브프로블럼으로 등장하는 균형 운송 문제에 Monge 구조가 존재함을 확인한다. 이때 서브프로블럼의 최적 이중 해를 위에서 도출한 압축식으로 직접 삽입하면, Benders 절단식이 명시적이고 선형적인 형태가 된다. 결과적으로 두 가지 새로운 지수형 MILP 모델이 제시되며, 각각은 (i) 전통적인 Benders 절단에 비해 절단 수가 감소하고, (ii) 절단 계수가 사전 계산된 차이값으로 고정돼 LP 풀기의 복잡도가 크게 낮아진다.

실험에서는 n = 100200, p = 1030 범위의 무작위 및 실제 데이터 인스턴스를 사용해 기존 최첨단 모델(Branch‑and‑Price, 기존 Benders 기반 모델)과 비교한다. 비증가 가중치(특히 악조건인 ‘악취 시설’ 문제)에서도 제안 모델은 평균 해결 시간에서 30 %~50 % 개선을 보였으며, 메모리 사용량도 현저히 낮았다. 또한, 인스턴스 규모가 증가해도 해결률이 95 % 이상 유지되는 등 견고성도 입증되었다.

이러한 결과는 Monge 비용 구조가 실제 물류·네트워크 설계에서 자주 나타나는 트리형 거리 혹은 허브‑스포크 구조와 일치한다는 점에서 실용적 의미가 크다. 즉, 비용 행렬이 Monge 성질을 만족한다면, 복잡한 이중 변수 계산을 피하고 직접적인 절단식으로 전환함으로써 대규모 위치 최적화 문제를 효율적으로 해결할 수 있다.