[Baekjoon] 2887번 행성 터널

Updated:

문제

때는 2040년, 이민혁은 우주에 자신만의 왕국을 만들었다. 왕국은 N개의 행성으로 이루어져 있다. 민혁이는 이 행성을 효율적으로 지배하기 위해서 행성을 연결하는 터널을 만들려고 한다.

행성은 3차원 좌표위의 한 점으로 생각하면 된다. 두 행성 A(xA, yA, zA)와 B(xB, yB, zB)를 터널로 연결할 때 드는 비용은 min( xA-xB , yA-yB , zA-zB )이다.

민혁이는 터널을 총 N-1개 건설해서 모든 행성이 서로 연결되게 하려고 한다. 이때, 모든 행성을 터널로 연결하는데 필요한 최소 비용을 구하는 프로그램을 작성하시오.

입력

첫째 줄에 행성의 개수 N이 주어진다. (1 ≤ N ≤ 100,000) 다음 N개 줄에는 각 행성의 x, y, z좌표가 주어진다. 좌표는 -109보다 크거나 같고, 109보다 작거나 같은 정수이다. 한 위치에 행성이 두 개 이상 있는 경우는 없다.

출력

첫째 줄에 모든 행성을 터널로 연결하는데 필요한 최소 비용을 출력한다.

예제

Example 1:

Input: 
5
11 -15 -15
14 -5 -15
-1 -1 -5
10 -4 -1
19 -4 19
Output: 
4

조건

시간 제한 : 1초

메모리 제한 : 128 MB

풀이과정

내 풀이1 (메모리 초과)

sort()와 heapq를 사용해보았지만 메모리 초과가 나왔다. 크루스칼 알고리즘을 적용하는데 있어서 다른 방법이 필요하다.

import heapq

def find_parent(parent, x):
    if parent[x] != x:
        parent[x] = find_parent(parent, parent[x])
    return parent[x]

def union_parent(parent, a, b):
    a = find_parent(parent, a)
    b = find_parent(parent, b)
    if a > b:
        parent[a] = b
    else:
        parent[b] = a
        
# n = 행성의 개수
n = int(input())

# 행성 입력
graph = []
for _ in range(n):
    x, y, z = map(int, input().split())
    graph.append((x, y, z))

# 간선 정리
edges = []
for i in range(n):
    for j in range(i+1, n):
        cost = min(abs(graph[i][0]-graph[j][0]), abs(graph[i][1]-graph[j][1]), abs(graph[i][2]-graph[j][2]))
        heapq.heappush(edges, (cost, i, j))                  

    
# 루트 노드 정의
parent = [0] * (n)
for i in range(n):
    parent[i] = i
    
result = 0
# 크루스칼 알고리즘 적용
while edges:
    cost, a, b = heapq.heappop(edges)
    if find_parent(parent, a) != find_parent(parent, b):
        union_parent(parent, a, b)
        result += cost
        
print(result)

내 풀이 2

책 풀이의 내용을 참고해서 코드를 짜보았다.

import sys
input = sys.stdin.readline
def find_parent(parent, x):
    if parent[x] != x:
        parent[x] = find_parent(parent, parent[x])
    return parent[x]

def union_parent(parent, a, b):
    a = find_parent(parent, a)
    b = find_parent(parent, b)
    if a > b:
        parent[a] = b
    else:
        parent[b] = a
        
# n = 행성의 개수
n = int(input())

# 행성 입력
graph = []
for i in range(n):
    x, y, z = map(int, input().split())
    graph.append((x, y, z, i))

# 간선 정리
edges = []

# x로 정렬
x_graph = sorted(graph, key = lambda x : x[0])
for i in range(n - 1):
    edges.append((x_graph[i+1][0]-x_graph[i][0], x_graph[i][3], x_graph[i+1][3]))
# y로 정렬
y_graph = sorted(graph, key = lambda x : x[1])
for i in range(n - 1):
    edges.append((y_graph[i+1][1]-y_graph[i][1], y_graph[i][3], y_graph[i+1][3]))
# z로 정렬
z_graph = sorted(graph, key = lambda x : x[2])
for i in range(n - 1):
    edges.append((z_graph[i+1][2]-z_graph[i][2], z_graph[i][3], z_graph[i+1][3]))
    
# 루트 노드 정의
parent = [0] * n
for i in range(n):
    parent[i] = i
    
result = 0

edges.sort()
# 크루스칼 알고리즘 적용
for edge in edges:
    cost, a, b = edge
    if find_parent(parent, a) != find_parent(parent, b):
        union_parent(parent, a, b)
        result += cost
        
print(result)

책 풀이

이 문제는 $N - 1$개의 터널을 설치해서 모든 행성이 연결되도록 요구하므로, 최소 신장 트리를 만드는 문제로 이해할 수 있다. 일단 크루스칼 알고리즘의 시간 복잡도는 간선의 개수가 E일 때 $O(ElogE)$이다. 이 문제에서는 기본적으로 임의의 두 노드 사이에 터널을 연결 할 수 있다고 가정하므로, 간선의 개수는 $N(N-1)/2$개가 될 것이다. $N$이 최대 100,000이라는 입력 조건을 감안해보면 이는 매우 큰 수가 될 수 있으므로, 모든 두 행성 간의 거리를 확인하는 방법으로는 문제를 해결할 수 없다.

하지만, 터널의 비용이 $min( X_A-X_B , \, Y_A-Y_B , \, Z_A-Z_B )$라고 정의되어 있다. 이러한 특징을 이용하면 고려할 간선의 개수를 줄일 수 있다. 입력을 받은 뒤에 x축, y축, z축을 기준으로 각각 정렬을 수행한다. 예를 들어 문제에서 나온 것과 같이 모든 행성의 좌표가 (11, -15, -15), (14, -5, -15), (-1, -1, -5), (10, -4, -1), (19, -4, 19)라고 해보자. 이때 x축만 고려해서 정렬을 수행하면 -1,10,11,14,19가 된다. 결과적으로 각 행성의 x축에서의 거리는 차례대로 11, 1, 3, 5가 되는 것이다.

aaa

결과적으로 x축에 대해서는 4개의 간선만 고려하면 되는 것이다. 여기서 알아두어야 할 점은, 만약에 y축과 z축은 무시하고 오직 x축만 존재한다고 했을 때, 이러한 4개의 간선만 이용해도 항상 최소 신장 트리를 만들 수 있다는 점이다. 즉, 이러한 방법을 이용하면 최소 신장 트리를 만들지 못하는 경우는 존재하지 않는다. 결과적으로 x축, y축, z축에 대하여 정렬 이후에 각각 $N-1$개의 간선만 고려해도 최적의 솔루션(Optimal Solution)을 찾을 수 있다는 아이디어를 떠올릴 수 있으면 된다.

따라서 문제 풀이를 위해 고려한 총 간선의 개수는 $3 \times (N-1)$개가 되고, 이를 이용해 크루스칼 알고리즘을 수행하면, 제한 시간 안에 해결할 수 있다.

# 특정 원소가 속한 집합을 찾기
def find_parent(parent, x):
    # 루트 노드가 아니라면, 루트 노드르 찾을 때까지 재귀적으로 호출
    if parent[x] != x:
        parent[x] = find_parent(parent, parent[x])
    return parent[x]

# 두 원소가 속한 집합을 합치기
def union_parent(parent, a, b):
    a = find_parent(parent, a)
    b = find_parent(parent, b)
    if a < b:
        parent[b] = a
    else:
        parent[a] = b
        
# 노드의 개수 입력받기
n = int(input())
parent = [0] * (n + 1) # 부모 테이블 초기화

# 모든 간선을 담을 리스트와 최종 비용을 담을 변수
edges = []
result = 0

# 부모 테이블상에서, 부모를 자기 자신으로 초기화
for i in range(1, n + 1):
    parent[i] = i
    
x = []
y = []
z = []

# 모든 노드에 대한 좌표 값 입력받기
for i in range(1, n + 1):
    data = list(map(int, input().split()))
    x.append((data[0], i))
    y.append((data[1], i))
    z.append((data[2], i))
    
x.sort()
y.sort()
z.sort()

# 인접한 노드들로부터 간선 정보를 추출하여 처리
for i in range(n - 1):
    # 비용순으로 정렬하기 위해서 튜플의 첫 번째 원소를 비용으로 설정
    edges.append((x[i + 1][0]- x[i][0], x[i][1], x[i + 1][1]))
    edges.append((y[i + 1][0]- y[i][0], y[i][1], y[i + 1][1]))
    edges.append((z[i + 1][0]- z[i][0], z[i][1], z[i + 1][1]))
    
# 간선을 비용순으로 정렬
edges.sort()

# 간선을 하나씩 확인하며
for edge in edges:
    cost, a, b = edge
    # 사이클이 발생하지 않는 경우에만 집합에 포함
    if find_parent(parent, a) != find_parent(parent, b):
        union_parent(parent, a, b)
        result += cost
        
print(result)

Leave a comment