辗转相除法(欧几里德算法)原理

原理:

假设有两个数x和y，存在一个最大公约数z=(x,y)，即x和y都有公因数z，
那么x一定能被z整除，y也一定能被z整除，所以x和y的线性组合mx±ny也一定能被z整除。（m和n可取任意整数）


对于辗转相除法来说，思路就是：若x>y，设x/y=n余c，则x能表示成x=ny+c的形式，将ny移到左边就是x-ny=c，由于一般形式的mx±ny能被z整除，所以等号左边的x-ny（作为mx±ny的一个特例）就能被z整除，即x除y的余数c也能被z整除。

算法实现:

辗转相除法

有两整数a和b：

① a%b得余数c

② 若c=0，则b即为两数的最大公约数

③ 若c≠0，则a=b，b=c，再回去执行①

#include<iostream>
#include<cstdio>
using namespace std;
int main(){
	int a,b,m,c;
	scanf("%d %d",&a,&b);
	m=a*b;
	while(b!=0){
		c=a%b; a=b; b=c;
	}
	printf("The largest common divisor:%d\n", a);  
    printf("The least common multiple:%d\n", m/a); 
} 

其他算法:

 相减法

有两整数a和b：

① 若a>b，则a=a-b

② 若a<b，则b=b-a

③ 若a=b，则a（或b）即为两数的最大公约数

④ 若a≠b，则再回去执行①

例如求27和15的最大公约数过程为：

27－15＝12( 15>12 ) 15－12＝3( 12>3 )

12－3＝9( 9>3 ) 9－3＝6( 6>3 )

6－3＝3( 3==3 )

因此，3即为最大公约数

[cpp]  view plain  copy

1. #include<stdio.h>  
2. void main ( )  /* 相减法求最大公约数 */  
3. {    
4.    int m, n, a, b, c;  
5.     printf("Input two integer numbers:\n");  
6.    scanf ("%d,%d", &a, &b); m=a; n=b;   
7.      /* a, b不相等，大数减小数，直到相等为止。*/    
8.    while ( a!=b)   
9.          if (a>b)  a=a-b;        
10.          else  b=b-a;   
11.    printf("The largest common divisor:%d\n", a);  
12.    printf("The least common multiple:%d\n", m*n/a);  
13. }  

⑶穷举法

有两整数a和b：

① i=1

② 若a，b能同时被i整除，则t＝i

③ i++

④ 若 i <= a(或b)，则再回去执行②

⑤ 若 i > a(或b)，则t即为最大公约数，结束

改进：

① i= a(或b)

② 若a，b能同时被i整除，则i即为最大公约数，

结束

③ i--，再回去执行②

有两整数a和b：

① i=1

② 若a，b能同时被i整除，则t＝i

③ i++

④ 若 i <= a(或b)，则再回去执行②

⑤ 若 i > a(或b)，则t即为最大公约数，结束

改进：

① i= a(或b)

② 若a，b能同时被i整除，则i即为最大公约数，

结束

③ i--，再回去执行②

[cpp]  view plain  copy

1. #include<stdio.h>  
2. void main ()  /* 穷举法求最大公约数 */  
3. {    
4.    int  m, n, a, b, i, t;   
5.    printf("Input two integer numbers:\n");  
6.    scanf ("%d,%d", &a, &b); m=a;  n=b;   
7.    for (i=1; i<= a; i++)    
8.        if ( a%i == 0 && b%i ==0 )    t=i;  
9.    printf("The largest common divisor:%d\n", t);  
10.    printf("The least common multiple:%d\n", m*n/t);  
11. }   
12. /*  改进后的 
13.    for (t= a; t>0; t-- )     
14.        if ( a%t == 0 && b%t ==0 )    break;  
15. */  

[cpp]  view plain  copy

1. //穷举法求最小公倍数  
2.      for (i= a; ; i++ )  
3.          if ( i % a == 0 && i % b ==0 )     break;  
4.      printf("The least common multiple:%d\n", i )  
5.    
6. //多个数的最大公约数和最小公倍数  
7.      for (i= a; i>0; i-- )  
8.          if (a%i==0&&b%i==0&&c%i==0)     break;  
9.      printf("The largest common divisor:%d\n", i);  
10.      for (i= a; ; i++ )  
11.          if (i%a==0&&i%b==0&&i% c==0)    break;  
12.      printf("The least common multiple:%d\n", i )  

其他算法转自http://blog.csdn.net/iwm_next/article/details/7450424