题目:http://acm.hdu.edu.cn/showproblem.php?pid=2028
学习了求最大公约数的stein方法,一个不容易理解的地方是当a和b都为奇数时,stein(a,b) = stein(abs(a - b)/2, min(a,b)):
假设a、b的最大公约数为c,则a = mc,b = nc,从而abs(a - b) = abs(m-n)c,min(a,b) = min(m,n)c,所以abs(a - b)和min(a,b)的最大公约数即是a和b的最大公约数,另一方面,由于a和b都为奇数,abs(a-b)为偶数,min(a,b)为奇数,所以可以先将abs(a-b)除以2在和min(a,b)求最大公约数,因为偶数和奇数的最大公约数一定不会被2整除
#include <stdio.h>
#include <stdlib.h>
int stein(int a, int b)
{
if(a == 0) return b;
if(b == 0) return a;
if(a & 1){
if(b & 1) return stein(abs(a - b) >> 1, a > b ? b : a);
else return stein(a, b >> 1);
}
else{
if(b & 1) return stein(a >> 1, b);
else return stein(a >> 1, b >> 1) << 1;
}
}
int lowestCommonMultiple(int arr[], int n)
{
int lcm = arr[0], i = 1, gcd;
for(; i < n; ++i){
if(lcm > 1){
gcd = stein(lcm, arr[i]);
lcm = lcm / gcd * arr[i];
}
else lcm = arr[i];
}
return lcm;
}
int main()
{
int i, n, a[200];
while(scanf("%d", &n) == 1){
for(i = 0; i < n; ++i) scanf("%d", a + i);
printf("%d\n", lowestCommonMultiple(a, n));
}
return 0;
}