LCM与GCD算法
LCM(最小公倍数)和GCD最大公因数在做ACM题时经常会用到,求两个整数的LCM和GCD有两种方法。
1. 辗转相除法(欧几里得算法)
定理:对于任意的两个整数a、b (a
≥
\geq
≥b), 有(a,b) = (b, a%b) . ((a, b)表示a,b的最大公因数)
证明如下:
a = qb + r,其中q为整数,0
≤
\leq
≤a
<
<
<b .
设d = (a, b),则b = md,a = nd .
则a = qmd + r = nd,进一步推出r = (n-qm)d .
故d也是r的因数,即d
≤
\leq
≤(b, r) = (b, a%b) .
同理,设p = (a, a%b) = (a, r),则r = sp,b = tp,d
≤
\leq
≤p.
则a = qtp + sp = (qt+s)p .
故p也是a的因数, 即p
≤
\leq
≤(a, b) = d.
综上,d = p,原命题得证 .
所以要求两个数的最大公因数,只需根据递推式不断进行递推,并更新a = b, b = a%b, 直到a%b = 0为止,则此时的a即为 (a, b) . 求得 (a, b)以后,则[a, b] (最小公倍数)便可由 a*b/(a, b)求得 .
2. 素因子分解
定理:任意一个正整数都能分解成若干个素数的幂的乘积的形式 .
证明略 .
由此可知,a =
p
1
a
1
p
2
a
2
.
.
.
p
n
a
n
p^{a_1}_1p^{a_2}_2...p^{a_n}_n
p1a1p2a2...pnan,b =
p
1
b
1
p
2
b
2
.
.
.
p
n
b
n
p^{b_1}_1p^{b_2}_2...p^{b_n}_n
p1b1p2b2...pnbn . 其中
a
i
,
b
i
≥
0
a_i,b_i\geq0
ai,bi≥0 .
故 (a, b) =
p
1
m
i
n
(
a
1
,
b
1
)
p
2
m
i
n
(
a
2
,
b
2
)
.
.
.
p
n
m
i
n
(
a
n
,
b
n
)
p^{min(a_1,b_1)}_1p^{min(a_2,b_2)}_2...p^{min(a_n,b_n)}_n
p1min(a1,b1)p2min(a2,b2)...pnmin(an,bn) .
[a, b] =
p
1
m
a
x
(
a
1
,
b
1
)
p
2
m
a
x
(
a
2
,
b
2
)
.
.
.
p
n
m
a
x
(
a
n
,
b
n
)
p^{max(a_1,b_1)}_1p^{max(a_2,b_2)}_2...p^{max(a_n,b_n)}_n
p1max(a1,b1)p2max(a2,b2)...pnmax(an,bn) .