[PUBLISHER] upload files #129

_posts/Algorithms/BOJ/2023-11-02-23-random-ps-1.md

@@ -0,0 +1,89 @@
+---
+share: true
+toc: true
+math: true
+categories:
+  - Algorithms
+  - BOJ
+tags:
+  - algorithms
+  - boj
+title: 랜덤 PS일지 (1)
+date: 2023-11-02
+github_title: 2023-11-02-23-random-ps-1
+---
+
+## 30323번
+
+- [BOJ 30323](https://www.acmicpc.net/problem/30323): Exponentiation
+
+주어진 $\alpha = x + x^{-1}$와 $\beta$에 대하여 $x^\beta + x^{-\beta} \pmod m$을 구하면 된다.
+
+먼 옛날 곱셈 공식의 변형
+
+$$
+x^2 + \frac{1}{x^2} = \left( x + \frac{1}{x} \right)^2 - 2
+$$
+
+을 공부하던 기억을 떠올리며...
+
+$\beta$가 짝수인 경우에는 쉽게 처리할 수 있다.
+
+$$
+x^\beta + x^{-\beta} = (x^{\beta/2} + x^{-\beta/2})^2 - 2.
+$$
+
+홀수인 경우에는[^1] 지수의 크기를 절반으로 줄이는 것이 핵심일테니,
+
+$$
+x^\beta + x^{-\beta} = \left( x^{\frac{\beta+1}{2}} + x^{-\frac{\beta+1}{2}}\right)\left( x^{\frac{\beta-1}{2}} + x^{-\frac{\beta-1}{2}}\right) - (x + x^{-1})
+$$
+
+를 생각하면 된다. 같은 지수에 대해서 여러 번 계산을 할 수도 있으므로, 결과를 저장해뒀다가 꺼내서 쓰도록 하자.
+
+## 12728, 12925번
+
+- [BOJ 12728](https://www.acmicpc.net/problem/12728): $n$제곱 계산
+- [BOJ 12925](https://www.acmicpc.net/problem/12925): Numbers
+
+두 문제가 같은 문제이다. $(3 + \sqrt{5})^n$의 정수부분의 마지막 $3$자리를 구하면 된다.
+
+중학교 시절 에이급 수학에서 $(3 + 2\sqrt{2})^5$의 정수부분을 구하라는 문제를 봤었는데 이 때 사용했던 아이디어가 켤레무리수를 생각하는 것이었다. 비슷한 아이디어를 2017학년도 서울대학교 공과대학 수시 일반 심층 면접에서도 $(2 + \sqrt{5})^n$이 나와 사용했었다. 그리고...
+
+> **정리.** $\alpha = 3 + \sqrt{5}$, $\beta = 3 - \sqrt{5}$ 일 때, $\alpha^n + \beta^n \in \mathbb{N}$.[^2]
+
+여기서 핵심은 $0 < \beta < 1$ 임을 이용하는 것이다. 따라서, $\alpha^n$의 정수부분은 $\alpha^n + \beta^n - 1$이 된다. 이제 $\alpha^n + \beta^n$만 구하면 된다. 근과 계수의 관계를 이용하면 수열 $s_n = \alpha^n + \beta^n$에 대한 귀납적 정의를 얻을 수 있다.
+
+$n$이 크기 때문에 계산은 행렬을 빠르게 거듭제곱하는 형태를 사용하면 된다.
+
+## 4215번
+
+- [BOJ 4215](https://www.acmicpc.net/problem/4215): To Add or to Multiply
+
+우선 주어진 범위 내에서 가능한 곱셈의 횟수는 최대 $29$번이다. ($2^{30} > 10^9$) 그러므로 모든 $m$에 대해 조사할 시간은 충분하다.
+
+곱셈을 $g$번 한다고 가정하고, 이 때 $0, \dots, g$ 번째 덧셈 횟수를 $h_i$라 하자. 그렇다면 입력이 $x$일 때, 우리가 생성한 프로그램의 최종 값은
+
+$$
+x \cdot m^g + a \sum_{i=0}^g m^{g-i} h_i
+$$
+
+가 된다. 입력 범위 $[p, q]$를 $[r, s]$로 보내야 하므로, 부등식을 세워보면
+
+$$
+ \tag{1}
+a\sum_{i=0}^g m^{g-i} h_i \in \left[ r - p \cdot m^g, s - q\cdot m^g \right]
+$$
+
+여야 함을 알 수 있다. 그러면 이제 프로그램의 길이를 최소화 해야하는데, 길이가 $\sum_{i=0}^g h_i + g$ 이므로 $h_i$의 합을 최소화하는 방향으로 생각해본다. 좌변의 식을 조금 풀어서 써보면
+
+$$
+a\sum_{i=0}^g m^{g-i} h_i =(am^g)h_0 + (am^{g-1})h_1 + \cdots + (a)h_g
+$$
+
+이므로 $\sum h_i$를 최소화 하려면 $h_0$부터 $h_g$까지 순서대로 잡아주면 된다. 이 때 $(1)$의 범위를 만족해야 하므로, 범위를 만족하는 최소의 $h_0$가 잡히면 바로 끝. 만약 불가능하다면 $r - p\cdot m^g$를 넘지 않는 최대의 $h_0$를 잡고 이번에는 구간 양 끝에서 $am^gh_0$를 뺀 뒤 같은 방법으로 $h_1$을 잡으면 된다.
+
+모든 가능한 프로그램의 후보를 얻었다면, 가장 짧은 것을 찾고 사전 순으로 제일 먼저 오는 것을 찾으면 된다. 사전 순 정렬의 경우 귀납적으로 생각하면 쉽게 구현할 수 있다. 앞에서부터 연산의 종류와 횟수를 비교하면 된다.
+
+[^1]: 원래 빠른 거듭제곱을 할 때는 $a^n = a \cdot (a^2)^{(n-1)/2}$ 으로 했던 것 같은데 이 경우에는 잘 안되므로...
+[^2]: 증명은 귀납법. 이항정리를 써도 좋고, 수열의 귀납적 정의를 사용해도 좋다.