최소 공통 조상 알고리즘

HeeSeong·2021년 3월 22일
0

파이썬 코딩테스트

목록 보기
12/12
post-thumbnail

Ⅰ. LCA 알고리즘


  • 두 노드의 공통된 조상 중에서 가장 가까운 조상을 찾는 문제입니다.

LCA 알고리즘 동작순서


  1. 모든 노드에 대한 깊이를 계산합니다.

  2. 최소 공통 조상을 찾을 두 노드를 확인합니다.

  3. 먼저 두 노드의 깊이가 동일하도록 거슬러 올라갑니다.

  4. 부모가 같아질 때까지 반복적으로 두 노드의 부모 방향으로 거슬러 올라갑니다.

  5. 모든 LCA(a, b) 연산에 대하여 3~4번의 과정을 반복합니다.


LCA 알고리즘 시간 복잡도


  • 매 쿼리마다 부모 방향으로 거슬로 올라가기 위해 최악의 경우 O(N)O(N)의 시간 복잡도가 요구됩니다. (한쪽에 치우쳐서 일렬로 된 노드인 경우)

  • 따라서 모든 쿼리(M)를 처리할 때의 시간 복잡도는 O(NM)O(NM)입니다.


Ⅱ. 개선된 LCA 알고리즘



  • 각 노드가 거슬러 올라가는 속도를 빠르게 만드는 방법입니다.

  • 15칸을 올라가면 8칸 🠖 4칸 🠖 2칸 🠖 1칸

  • 2의 제곱 형태로 거슬러 올라가도록 하면 O(logN)O(logN)의 시간 복잡도를 보장할 수 있습니다.

  • 메모리를 더 사용하여 각 노드에 대하여 2i2^i번째 부모에 대한 정보를 기록합니다.

  • 즉, 공간 복잡도를 늘리고 시간 복잡도를 줄이는 trade-off


개선된 LCA 알고리즘 시간 복잡도


  • 다이나믹 프로그래밍(DP)를 이용해 시간 복잡도를 개선할 수 있습니다.

  • 매 쿼리(M)마다 부모를 거슬러 올라가기 위해 O(logN)O(logN)의 복잡도가 필요합니다.

  • 모든 쿼리를 처리할 때의 시간 복잡도는 O(MlogN)O(MlogN)입니다.


Ⅲ. 예제


① LCA


❔ 문제 설명


N(2 ≤ N ≤ 50,000)개의 정점으로 이루어진 트리가 주어진다. 트리의 각 정점은 1번부터 N번까지 번호가 매겨져 있으며, 루트는 1번이다.

첫째 줄에 노드의 개수 N이 주어지고, 다음 N-1개 줄에는 트리 상에서 연결된 두 정점이 주어진다. 그 다음 줄에는 가장 가까운 공통 조상을 알고싶은 쌍의 개수 M이 주어지고, 다음 M개 줄에는 정점 쌍이 주어진다.

M개의 줄에 차례대로 입력받은 두 정점의 가장 가까운 공통 조상을 출력한다.


💡 문제 해결 아이디어


  • 문제 이름 그대로 최소 공통 조상 찾기 알고리즘을 이용해 푸는 문제입니다.

✔️ 풀이

Python

import sys
sys.setrecursionlimit(int(1e5)) # 런타임 오류를 피하기 위한 재귀 깊이 제한 설정

n = int(input())

parent = [0] * (n + 1) # 부모 노드 정보
d = [0] * (n + 1) # 각 노드까지의 깊이
c = [0] * (n + 1) # 각 노드의 깊이가 계산되었는지 여부
graph = [[] for _ in range(n + 1)] # 그래프(graph) 정보

for _ in range(n - 1):
    a, b = map(int, input().split())
    graph[a].append(b)
    graph[b].append(a)

# 루트 노드부터 시작하여 깊이(depth)를 구하는 함수
def dfs(x, depth):
    c[x] = True
    d[x] = depth
    for y in graph[x]:
        if c[y]: # 이미 깊이를 구했다면 넘기기
            continue
        parent[y] = x
        dfs(y, depth + 1)

# A와 B의 최소 공통 조상을 찾는 함수
def lca(a, b):
    # 먼저 깊이(depth)가 동일하도록
    while d[a] != d[b]:
        if d[a] > d[b]:
            a = parent[a]
        else:
            b = parent[b]
    # 노드가 같아지도록
    while a != b:
        a = parent[a]
        b = parent[b]
    return a

dfs(1, 0) # 루트 노드는 1번 노드

m = int(input())

for i in range(m):
    a, b = map(int, input().split())
    print(lca(a, b))

② LCA 2


❔ 문제 설명


N(2 ≤ N ≤ 100,000)개의 정점으로 이루어진 트리가 주어진다. 트리의 각 정점은 1번부터 N번까지 번호가 매겨져 있으며, 루트는 1번이다.

두 노드의 쌍 M(1 ≤ M ≤ 100,000)개가 주어졌을 때, 두 노드의 가장 가까운 공통 조상이 몇 번인지 출력한다.

첫째 줄에 노드의 개수 N이 주어지고, 다음 N-1개 줄에는 트리 상에서 연결된 두 정점이 주어진다. 그 다음 줄에는 가장 가까운 공통 조상을 알고싶은 쌍의 개수 M이 주어지고, 다음 M개 줄에는 정점 쌍이 주어진다.

M개의 줄에 차례대로 입력받은 두 정점의 가장 가까운 공통 조상을 출력한다.


💡 문제 해결 아이디어


  • 일반 LCA 알고리즘 로직으로는 이 문제는 시간 초과 판정을 받습니다.

  • 개선된 LCA 알고리즘을 사용해서 풀어야 합니다.


✔️ 풀이

Python

import sys
input = sys.stdin.readline # 시간 초과를 피하기 위한 빠른 입력 함수
sys.setrecursionlimit(int(1e5)) # 런타임 오류를 피하기 위한 재귀 깊이 제한 설정
LOG = 21 # 2^20 = 1,000,000

n = int(input())
parent = [[0] * LOG for _ in range(n + 1)] # 부모 노드 정보
d = [0] * (n + 1) # 각 노드까지의 깊이
c = [0] * (n + 1) # 각 노드의 깊이가 계산되었는지 여부
graph = [[] for _ in range(n + 1)] # 그래프(graph) 정보

for _ in range(n - 1):
    a, b = map(int, input().split())
    graph[a].append(b)
    graph[b].append(a)

# 루트 노드부터 시작하여 깊이(depth)를 구하는 함수
def dfs(x, depth):
    c[x] = True
    d[x] = depth
    for y in graph[x]:
        if c[y]: # 이미 깊이를 구했다면 넘기기
            continue
        parent[y][0] = x
        dfs(y, depth + 1)

# 전체 부모 관계를 설정하는 함수
def set_parent():
    dfs(1, 0) # 루트 노드는 1번 노드
    for i in range(1, LOG):
        for j in range(1, n + 1):
            parent[j][i] = parent[parent[j][i - 1]][i - 1]

# A와 B의 최소 공통 조상을 찾는 함수
def lca(a, b):
    # b가 더 깊도록 설정
    if d[a] > d[b]:
        a, b = b, a
    # 먼저 깊이(depth)가 동일하도록
    for i in range(LOG - 1, -1, -1):
        if d[b] - d[a] >= (1 << i):
            b = parent[b][i]
    # 부모가 같아지도록
    if a == b:
        return a;
    for i in range(LOG - 1, -1, -1):
        # 조상을 향해 거슬러 올라가기
        if parent[a][i] != parent[b][i]:
            a = parent[a][i]
            b = parent[b][i]
    # 이후에 부모가 찾고자 하는 조상
    return parent[a][0]

set_parent()

m = int(input())

for i in range(m):
    a, b = map(int, input().split())
    print(lca(a, b))
profile
끊임없이 성장하고 싶은 개발자

0개의 댓글