193 lines
7.2 KiB
Markdown
193 lines
7.2 KiB
Markdown
deck:: Logseq/coding tip
|
|
|
|
- ## 1. 개념(Concept)
|
|
- 크고 복잡한 문제를 {{c1 작고 간단한 하위 문제(Sub-problem)}}들로 나눈 뒤, 각각 해결하여 합치는 알고리즘 설계 패러다임.
|
|
id:: 695bbd99-6b9b-4019-a174-9238a89758de
|
|
-
|
|
- ## 2. 동작 원리 (Algorithm Flow)
|
|
- ### 1단계 {{c1 분할(Divide)}} -> 2단계 {{c1 정복(Conquer)}} -> 3단계 {{c1 결합(Combine)}}
|
|
id:: 695bbe3e-8b85-4abf-a9ed-6f0c24c8f2dd
|
|
#+BEGIN_EXTRA
|
|
분할 : 원래 문제를 같은 유형의 더 작은 하위 문제로 쪼갠다.
|
|
정복 : 하위 문제를 재귀적으로 해결한다.(문제가 충분히 작다면 바로 해결)
|
|
결합 : 하위 문제의 해답들을 합쳐 원래 문제의 답을 만든다.
|
|
#+END_EXTRA
|
|
-
|
|
- ## 3. 시간복잡도
|
|
- ### 문제의 크기를 $N$ 이라고 할 때 {{c1 $O(\log N)$}} 이나 {{c1 $O(N \log N)$}} 으로 감소됨.
|
|
id:: 695bbeb8-ecdc-4d0b-a29b-43231188ba86
|
|
-
|
|
- ## 4. 예시
|
|
- ### 1) 행렬의 거듭제곱($N \times N$ 크기의 행렬 $A$를 $B$번 거듭제곱 하기)
|
|
- #### 점화식
|
|
- $n$이 짝수일 때: $A^n = (A^{n/2}) \times (A^{n/2})$
|
|
- $n$이 홀수일 때: $A^n = (A^{n/2}) \times (A^{n/2}) \times A$
|
|
- #### 행렬곱(multiply_matrix) 함수 코드(python)
|
|
id:: 695bc027-6274-498e-b798-0e94060eb147
|
|
$N \times N$ 크기의 행렬 m1, m2와 N을 입력받아서 두 행렬을 곱한 행렬을 반환하는 코드
|
|
#card
|
|
- ```python
|
|
def multiply_matrix(m1, m2, N):
|
|
result = [[0]* N for _ in range(N)]
|
|
for i in range(N):
|
|
for j in range(N):
|
|
for k in range(N):
|
|
result[i][j] += m1[i][k] * m2[k][j]
|
|
result[i][j] %= 1000 # 문제 조건에 따라 모듈러 연산 추가
|
|
return result
|
|
```
|
|
- #### 분할정복(power_matrix) 함수 코드(python)
|
|
거듭제곱 될 행렬 adj, 현재 곱해지는 지수 n, 행렬의 크기 N을 입력받아서 분할정복으로 adj^n 을 구하는 코드
|
|
- 1) n이 1일 경우(차수가 1일 경우) #card
|
|
id:: 695bc161-b88a-4f55-aaf7-7bc17ccd4dcc
|
|
- ```python
|
|
# Base case : 지수가 1이면 그대로 반환한다.
|
|
if n==1 :
|
|
return [[num % 1000 for num in row] for row in adj]
|
|
```
|
|
- 2) 분할 #card
|
|
id:: 695bc1d6-8a7c-475c-9c51-45a672597058
|
|
- ```python
|
|
# Divide: 지수를 절반으로 나누어 재귀 호출
|
|
temp = power_matrix(adj, n//2, N)
|
|
```
|
|
- 3) 정복 및 결합 #card
|
|
id:: 695bc232-7ee9-4c02-a3b1-1f65061737f9
|
|
- ```python
|
|
if n % 2 == 0:
|
|
# 짝수면: temp * temp
|
|
return multiply_matrix(temp, temp, N)
|
|
else:
|
|
# 홀수면: temp * temp * A
|
|
return multiply_matrix(multiply_matrix(temp, temp, N), adj, N)
|
|
```
|
|
- 전체코드
|
|
- ```python
|
|
# 행렬 곱셈 함수 (N*N 행렬 두 개를 곱함)
|
|
def multiply_matrix(m1, m2, N):
|
|
result = [ * N for _ in range(N)]
|
|
for i in range(N):
|
|
for j in range(N):
|
|
for k in range(N):
|
|
result[i][j] += m1[i][k] * m2[k][j]
|
|
result[i][j] %= 1000 # 문제 조건에 따라 모듈러 연산 추가
|
|
return result
|
|
|
|
# 분할 정복을 이용한 거듭제곱 함수
|
|
def power_matrix(adj, n, N):
|
|
# Base Case: 지수가 1이면 그대로 반환
|
|
if n == 1:
|
|
return [[elem % 1000 for elem in row] for row in adj]
|
|
|
|
# Divide: 지수를 절반으로 나누어 재귀 호출
|
|
temp = power_matrix(adj, n // 2, N)
|
|
|
|
# Conquer & Combine
|
|
if n % 2 == 0:
|
|
# 짝수면: temp * temp
|
|
return multiply_matrix(temp, temp, N)
|
|
else:
|
|
# 홀수면: temp * temp * A
|
|
return multiply_matrix(multiply_matrix(temp, temp, N), adj, N)
|
|
|
|
```
|
|
-
|
|
- ### 2) 피보나치 수열을 행렬로 표현해서 분할정복으로 구하기
|
|
- #### 피보나치 수열의 행렬화
|
|
- 피보나치 수열의 점화식 $F_n = F_{n-1} + F_{n-2}$를 행렬 곱셈 형태로 변환.
|
|
- 1) 기본 점화식의 행렬화 #card
|
|
id:: 69623201-075b-4955-9d5b-aedcbcf16170
|
|
- $$F_{n+1} = F_n + F_{n-1}$$
|
|
$$F_n = F_{n-1} + F_{n-2}$$
|
|
- $$
|
|
\begin{bmatrix} F_{n+1} \\ F_n \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} F_n \\ F_{n-1} \end{bmatrix}
|
|
$$
|
|
- 2) 재귀적으로 쭉 대입하면 최종적으로 다음과 같은 식이 나옴. #card
|
|
id:: 69622f8c-b595-40fc-83b5-7298ca660f3a
|
|
- $$
|
|
\begin{bmatrix} F_{n+1} \\ F_n \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}^n \begin{bmatrix} F_1 \\ F_0 \end{bmatrix}
|
|
$$
|
|
- $$
|
|
\begin{bmatrix} F_{n+1} \\ F_n \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}^n \begin{bmatrix} 1 \\ 0 \end{bmatrix}
|
|
$$
|
|
#+BEGIN_EXTRA
|
|
즉 $\begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}^n$ 값의 1행 0열의 값이 $F_n$ 값이 된다.
|
|
#+END_EXTRA
|
|
- #### 결론 : $\begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}^n$ 을 계산하면 n번째 피보나치 수열의 항을 구할 수 있다.
|
|
- #### 예시문제 : n번째 피보나치 수를 1000000007 로 나눈 나머지를 출력하라.(Python)
|
|
id:: 6962333f-ebe3-43a9-bca6-96e9442286e1
|
|
정의할 함수는 mul_mat(m1, m2) 와 power_mat(adj, n) 이다.
|
|
아래와 같이 코드가 주어질 때 두 함수를 정의하라.
|
|
```python
|
|
import sys
|
|
input = sys.stdin.readline
|
|
|
|
# 문제에서 요구하는 나누기 상수
|
|
MOD = 1_000_000_007
|
|
|
|
# 1. 행렬 곱셈 함수 (2x2 고정)
|
|
def mul_mat(m1, m2):
|
|
# 행렬 곱
|
|
|
|
# 2. 분할 정복을 이용한 거듭제곱
|
|
def power_mat(adj, n):
|
|
# 분할정복
|
|
|
|
# Main
|
|
n = int(input())
|
|
|
|
# 기본 행렬 [[1, 1], [1, 0]]
|
|
base_matrix = [[1, 1], [1, 0]]
|
|
|
|
# 행렬 제곱 계산
|
|
result_matrix = power_mat(base_matrix, n)
|
|
|
|
# 1행 0열의 값이 n번째 값
|
|
print(result_matrix[1][0])
|
|
|
|
```
|
|
#card
|
|
- ```python
|
|
import sys
|
|
input = sys.stdin.readline
|
|
|
|
# 문제에서 요구하는 나누기 상수
|
|
MOD = 1_000_000_007
|
|
|
|
# 1. 행렬 곱셈 함수 (2x2 고정)
|
|
def mul_mat(m1, m2):
|
|
res = [[0]*2 for _ in range(2)]
|
|
for i in range(2):
|
|
for j in range(2):
|
|
for k in range(2):
|
|
res[i][j] += m1[i][k] * m2[k][j]
|
|
res[i][j] %= MOD
|
|
return res
|
|
|
|
# 2. 분할 정복을 이용한 거듭제곱
|
|
def power_mat(adj, n):
|
|
if n == 1:
|
|
return [[x % MOD for x in row] for row in adj]
|
|
|
|
temp = power_mat(adj, n // 2)
|
|
|
|
# 짝수일 때: A^(n//2) * A^(n//2)
|
|
if n % 2 == 0:
|
|
return mul_mat(temp, temp)
|
|
# 홀수일 때: A^(n//2) * A^(n//2) * A
|
|
else:
|
|
return mul_mat(mul_mat(temp, temp), adj)
|
|
|
|
# Main
|
|
n = int(input())
|
|
|
|
# 기본 행렬 [[1, 1], [1, 0]]
|
|
base_matrix = [[1, 1], [1, 0]]
|
|
|
|
# 행렬 제곱 계산
|
|
result_matrix = power_mat(base_matrix, n)
|
|
|
|
# 1행 0열의 값이 n번째 값
|
|
print(result_matrix[1][0])
|
|
|
|
``` |