求最大公约数--Stein算法

Stein算法

伪代码：

Stein算法(假设0<=b<a):
r←0
while b>0
   do if a偶,b偶 then a←a>>1   b←b>>1   r←r+1
        else if a偶,b奇 then a←a>>1
        else if a奇,b偶 then b←b>>1
        else if a奇,b奇 then a←(a-b)>>1
        if a<b then 交换a,b
return a<<r

1. unsigned MaxDivisor(unsigned a, unsigned b)   
2. {   
3.     unsigned c = 0;   
4.     while(1)  
5.     {   
6.     // 退出条件   
7.         if(a==0)   
8.             return b << c;  
9.         else if(b == 0)   
10.             return a << c;  
11.     // 为提高速度，采用位的与运算，避免用取模判断奇偶   
12.         if(!(a & 1) && !(b & 1)) //a,b 都是偶数   
13.         {   
14.             a >>= 1; b >>= 1; ++c;   
15.         }   
16.         else if(!(a & 1) && (b & 1)) //a偶 b奇   
17.         {   
18.             a >>= 1;   
19.         }   
20.         else if((a & 1) && !(b & 1)) //a奇 b偶   
21.         {  
22.              b >>= 1;   
23.         }   
24.         else if((a & 1) && (b & 1)) //a,b都是奇数   
25.         {   
26.             unsigned tmp = a>b?b:a; //取较小的一个   
27.             a = a>b?a-b:(b-a); //绝对差值  
28.             b = tmp;   
29.         }   
30.     }  
31. }   

1.a,b均为偶，明显有公约数2

2.a,b奇偶性不同，则偶数的2可以去掉，因为不可能和奇数有公约数

3.a,b均为奇，此时可能有奇数的公约数

a = k * x,b = k* y (k为奇)

a = a - b = k * (x - y)

不断减小其他因子，直到最后a = b，相减a = 0,则b为k，再左移（即乘上之前的2^c）