유클리드 호제법

섬네일
**유클리드 호제법**(-互除法, Euclidean algorithm) 또는 **유클리드 알고리즘**은 2개의 자연수 또는 정식(整式)의 최대공약수를 구하는 알고리즘의 하나이다. 호제법이란 말은 두 수가 서로(互) 상대방 수를 나누어(除)서 결국 원하는 수를 얻는 알고리즘을 나타낸다. 2개의 자연수(또는 정식) a, b에 대해서 a를 b로 나눈 나머지를 r이라 하면(단, a>b), a와 b의 최대공약수는 b와 r의 최대공약수와 같다. 이 성질에 따라, b를 r로 나눈 나머지 r'를 구하고, 다시 r을 r'로 나눈 나머지를 구하는 과정을 반복하여 나머지가 0이 되었을 때 나누는 수가 a와 b의 최대공약수이다. 이는 명시적으로 기술된 가장 오래된 알고리즘으로서도 알려져 있으며, 기원전 300년경에 쓰인 《원론》 제7권, 명제 1부터 3까지에 해당한다. 최소공배수, 최대공약수를 구하기 위해서는 다음과 같은 표를 쓴다.

## 예시
1071과 1029의 최대공약수를 구하면,
* 1071은 1029로 나누어 떨어지지 않기 때문에, 1071을 1029로 나눈 나머지를 구한다. ≫ 42
* 1029는 42로 나누어 떨어지지 않기 때문에, 1029를 42로 나눈 나머지를 구한다. ≫ 21
* 42는 21로 나누어 떨어진다.
따라서, 최대공약수는 21이다.

78696과 19332의 최대공약수를 구하면,
     **78696** ＝ **19332**×4 ＋ 1368
     19332 ＝ 1368×14 ＋ 180
      1368 ＝ 180×7 ＋ 108
       180 ＝ 108×1 ＋ 72
       108 ＝ 72×1 ＋ 36
        72 ＝ 36×2 ＋ 0

따라서, 최대공약수는 36이다.

## 정수에서 반영됨으로써의 유클리드 호제법
본문에 있는 내용을 다시 한 번 상기해 보면, 다음 정리를 이용하여 유클리드 호제법이 성립하게 된다.

### 정리
a, b \in \mathbb{Z}이고, a를 b로 나눈 나머지가 r이라고 하자. (여기서 a\ge b이고, r은 0 \le r 인 정수.)

a와 b의 최대공약수를 \left(a, b\right)라고 하면, 다음이 성립한다.
:\left(a, b\right) = \left(b, r\right).

상기한 1071과 1029에 대한 예시를 위와 같은 표현으로 적어보면,
:\left(1071, 1029\right) = \left(1029, 42\right) = \left(42, 21\right) = \left(21, 0\right)=21
처럼 쓸 수 있다.

### 증명
a, b \in \mathbb{Z}이고, a\ge b라고 하자.

그러면, a=bq+r을 만족하는 유일한 정수 q, r이 존재한다. 이때, 0 \le r 이다.

\left(a, b\right) = d,
a=d\alpha,
b=d\beta라고 하자. 즉, \alpha와 \beta는 서로소이다.
:a=bq+r.
:\Rightarrow d\alpha =d\beta q+r.
:\Rightarrow d|r.(즉, r은 d의 배수)

이제, r=d\rho.라고 하자.

여기서 β와 ρ가 서로소라면, b(= dβ)와 r(=dρ)의 최대 공약수도 d가 된다.

만약 \left(\beta, \rho\right) = d' > 1(서로소가 아닌 수, 즉 다른 공약수를 가지는 수)라면, \beta = d'\beta', \rho=d'\rho'으로 두었을 때,
:a=bq+r.
:\Rightarrow d\alpha = d\beta q+d\rho = d d'\beta' q + d d' \rho'= d d'\left(\beta' q + \rho'\right).
:\Rightarrow \alpha = d'\left(\beta' q + \rho'\right).

이 되므로, d'|\alpha이다. (즉, α는 d'의 배수 )

즉, d'|\alpha, d'|\beta가 되어 \alpha와 \beta는 서로소라는 것에 모순이다.

이는 \left(\beta, \rho\right) = d' > 1 라는 가정에서 나타나는 모순이므로 \left(\beta, \rho\right) = 1(즉,  β와 ρ는 서로소)이다.

따라서 곧 \left(b, r\right) = d라는 것이다.

쉽게 풀어주면,

A를 B로 나눈 몫을 Q, 나머지를 R이라고 할 때, A=BQ+R로 나타낼 수 있다.A와 B의 최대 공약수를 G라고 할 때,R=A-BQ이므로 R도 G 로 나눌 수 있게 된다. 즉, G는 R과 B의 공약수가 된다. 따라서 R=GR_1,B=GB_1라고 쓸 수 있다.이때 R_1, B_1은 서로소, 즉 공약수를 갖지 않는다. R_1, B_1가 공약수 C를 가지면 R_1=CR_2, B_1=CB_2가 되어,R=GR_1=GCR_2,B=GB_1=GCB_2라고 쓸 수 있다. 그러면,A=BQ+R=GCB_2Q+GCR=GC(B_2Q+R_2 )가 된다. 이때 A, B의 최대 공약수가 CG가 되어 G가 최대공약수라는 것과 모순된다. 따라서 A와 B의 최대 공약수 G는 R과 B의 최대 공약수이다.

## 다항식에 대한 유클리드 호제법
임의의 두 다항식에 대해서도 유클리드 호제법과 같은 원리를 사용하여, 두 다항식의 최대차수공약다항식을 구할 수 있다. 즉, 다음과 같은 정리가 성립한다.

### 정리
f, g \in \mathbb{C}[x]이고, f를 g로 나눈 나머지가 r(x)이라고 하자. (여기서 \operatorname{deg}f\ge \operatorname{deg}g이고, r은 0 \le \operatorname{deg}r 인 다항식, \operatorname{deg}f는 f의 차수를 뜻한다.)

f와 g의 최대차수공약다항식을 \left(f, g\right)라고 하면, 다음이 성립한다.
:\left(f, g\right) = \left(g, r\right).

#### 다항식의 경우의 예시
f(x)=x^3+2x^2+x, g(x)=x^2-1을 예로 들어 보자.

:x^3+2x^2+x=(x^2-1)(x+2)+(2x+2).
:x^2-1=(2x+2)\left(\frac{x-1}{2}\right).

위 식으로부터, \left(x^3+2x^2+x, x^2-1\right) = \left(x^2-1, 2x+2\right)=\left(2x+2, 0\right)=x+1이다.

(\mathbb{C}[x]에서, 상수는 가역원이므로, 공약다항식에 대해서 말할 때에는 곱해주어서 계수를 조절해줄 수 있다.)

### 증명
f\in\mathbb{C}[x],g \in \mathbb{C}[x]\setminus\{0\}이고, \operatorname{deg}f\ge \operatorname{deg}g라고 하자.

그러면, \mathbb{C}[x]는 유클리드 정역이므로, f=gq+r을 만족하는 유일한 다항식 q, r이 존재한다. 이때, r=0이거나 \operatorname{deg}r 이다.

\left(f, g\right) = d, f=df', g=dg'라고 하자. 즉, \operatorname{deg}(f',g')=0이다.
:f=gq+r.
:\Rightarrow df' =dg'q+r.
:\Rightarrow d|r.

이제, r=dr'.라고 하자.

만약 \left(g', r'\right) = d'이 상수가 아니라면, g' = d'g*, r'=d'r*으로 두었을 때,
:f=gq+r.
:\Rightarrow df' = dg'q+dr' = d d'g* q + d d' r*= d d'\left(g* q + r*\right).
:\Rightarrow f' = d'\left(g* q + r*\right).
이 되므로, d'|f'이다.

즉, d'|f', d'|g'이 되어 \operatorname{deg}(f',g')=0이라는 것에 모순이다.

그러므로 \operatorname{deg}\left(g', r'\right) = 0이고, 그것은 곧 \left(g, r\right) = d라는 것이다.

## 알고리즘
# 입력으로 두 수 m,n(m>n)이 들어온다.
# n이 0이라면, m을 출력하고 알고리즘을 종료한다.
# m이 n으로 나누어 떨어지면, n을 출력하고 알고리즘을 종료한다.
# 그렇지 않으면, m을 n으로 나눈 나머지를 새롭게 m에 대입하고, m과 n을 바꾸고 3번으로 돌아온다.

위 과정은 “n, m에 대해서 나머지 연산을 실시할 수 있다”라는 조건에만 의존하므로, 정수환뿐 아니라, 임의의 유클리드 정역에 대해도 똑같은 과정을 거치면 공약인자가 구해진다.
### 소스 코드
* C

```c
int gcd(int a, int b)
{
	return b ? gcd(b, a%b) : a;
}
```


```c
int gcd(int a, int b)
{
    int c;
	while(b)
	{
		c = a % b;
		a = b;
		b = c;
	}
    return a;
}
```

* 파이썬

```python
def gcd(m,n):
	if m < n:
		m, n = n, m
	if n == 0:
		return m
    if m % n == 0:
		return n
	else:
		return gcd(n, m%n)
```


```python
def gcd(m,n):
    while n != 0:
       t = m%n
       (m,n) = (n,t)
    return abs(m)
```


```python
def gcd(m,n):
    while n! = 0:
	    if m < n:
		    m, n = n, m
	    if n == 0:
		    return m
	    if m % n == 0:
		    return n
```

```python3
def gcd(m,n):
    if n == 0:
        return m
    mod = m % n
    if mod != 0:
        m, n = n, mod
        return gcd(m, n)
    else:
        return n
```

* 자바

```java
public static int gcd(int p, int q)
 {
	if (q == 0) return p;
	return gcd(q, p%q);
 }
```


* F#

```fsharp
let rec gcd m n =
  match n with
  | 0 -> m
  | _ -> gcd n (m % n)
```


* 프롤로그

```prolog
gcd(A, 0, D) :- A >= 0, D is A, !.
gcd(A, 0, D) :- A < 0, D is -1 * A, !.

gcd(A, B, D) :- A < 0, B > 0, !, gcd(-1 * A, B, D).
gcd(A, B, D) :- A >= 0, B < 0, !, gcd(A, -1 * B, D).
gcd(A, B, D) :- A < 0, B < 0, !, gcd(-1 * A, -1 * B, D).

gcd(A, B, D) :- A >= 0, B > 0, A >= B, R is A mod B, !, gcd(B, R, D).
gcd(A, B, D) :- A >= 0, B > 0, B > A, !, gcd(B, A, D).
```


## 호제법의 확장
### 정수의 경우
정수 m, n 의 최대공약수(Greatest Common Divisor)를 gcd(m,n)와 같이 나타낼 때, 확장된 유클리드 호제법을 이용하여, am + bn = gcd(m,n)의 해가 되는 정수 a, b 짝을 찾아낼 수 있다.(a, b 중 한개는 보통 음수가 된다.) 이 식은 베주의 정의라고 한다. 위에서 든 예의 경우,

 1071 = 1 * 1029 + 42
 1029 = 24 * 42 + 21
 42 = 2 * 21
특히, m, n이 서로소(gcd(m,n) = 1)인 경우 유용한데, 그럼 위의 식은 am + bn = 1이 되고, 여기서 a는 모듈로 연산의 곱의 역원(modular multiplicative inverse)이 되기 때문이다.

### 다항식의 경우
f, g \in \mathbb{C}[x]의 최대차수공약다항식을 \left(f, g\right)라고 하면, fp + gq = \left(f, g\right)인 p, q \in \mathbb{C}[x]를 찾아낼 수 있다. 위에서 든 예의 경우,

f(x)=x^3+2x^2+x, g(x)=x^2-1

:x^3+2x^2+x=(x^2-1)(x+2)+(2x+2).
:x^2-1=(2x+2)\left(\frac{x-1}{2}\right).
이므로

:2x+2=x^3+2x^2+x-(x^2-1)(x+2)
가 된다.

특히, \left(f, g\right)=1인 경우 fp + gq = 1인 p, q \in \mathbb{C}[x]를 찾아낼 수 있다.

## 계산량
유클리드 호제법을 이용하여 나머지를 취하는 조작을, 최악의 경우라도 작은 수를 십진법으로 표시한 자리수의 5배를 반복하기 전에 최대공약수에 이른다(라메의 정리).

최대공약수를 구할 때, 소인수 분해를 하는 것이 더 쉽다고 생각하는 사람이 있을지도 모르지만, 위 정리는 소인수 분해를 이용하는 것보다도 유클리드 호제법을 실시하여 구하는 것이 훨씬 빠르다는 것을 말한다.

실제로 계산량 이론에서 최대공약수를 요구하는 것은 “쉬운 문제”로서 알려져 있고, 소인수 분해는 “어려운 문제”로 여겨지고 있다. 입력된 2개의 정수의 비트수를 n이라고 하면, 유클리드 알고리즘은 Θ(n) 회의 나눗셈으로 최대공약수를 구할 수 있다.

## 연분수
위에 든 예에서 나온 1071과 1029의 최대공약수를 구하는 과정은 다음과 같이 나타낼 수 있다.

:\begin{matrix}1071 &=& 1029 \times 1 + 42\\
1029 &=& 42 \times 24 + 21\\
42 &=& 21 \times 2\end{matrix}

즉,

:\begin{matrix}1071 / 1029 &=& 1 + 42 / 1029\\
1029 / 42 &=& 24 + 21 / 42\\
42 / 21 &=& 2\end{matrix}

따라서,

:\frac{1071}{1029} = 1 + \frac{1}{24 + \frac{1}{2}}

이와 같이, n와 m(n>m)의 최대공약수를 구할 때의 유클리드 호제법의 나눗셈의 모양은 유리수 n/m 의 연분수 전개와 같다.

## 같이 보기

* 최소공배수

분류:수론 알고리즘
분류:에우클레이데스