求最大公约数--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)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值