今天真是醉了,来看一下除了穷举之外的求最大公约数算法.
来源:互联网 发布:嵌入式软件有哪些 编辑:程序博客网 时间:2024/03/29 19:18
简直了,想到这里感觉白活了。
下面来看一下吧。
1.欧几里得(Euclid)算法
又称辗转相除法,依据定理gcd(a,b)=gcd(b,a%b)
实现过程演示: sample:gcd(15,10)=gcd(10,5)=gcd(5,0)=5
C语言实现:
1 int Euclid_GCD(int a, int b)2 {3 return b?Euclid_GCD(b, a%b):a;4 }
2.Stein 算法
一般实际应用中的整数很少会超过64位(当然现在已经允许128位了),对于这样的整数,计算两个数之间的模是很简单的。对于字长为32位的平台,计算两个不超过32位的整数的模,只需要一个指令周期,而计算64位以下的整数模,也不过几个周期而已。但是对于更大的素数,这样的计算过程就不得不由用户来设计,为了计算两个超过 64位的整数的模,用户也许不得不采用类似于多位数除法手算过程中的试商法,这个过程不但复杂,而且消耗了很多CPU时间。对于现代密码算法,要求计算 128位以上的素数的情况比比皆是,设计这样的程序迫切希望能够抛弃除法和取模。
依据定理:
gcd(a,a)=a,也就是一个数和其自身的公约数仍是其自身。
gcd(ka,kb)=k*gcd(a,b),也就是最大公约数运算和倍乘运算可以交换。特殊地,当k=2时,说明两个偶数的最大公约数必然能被2整除。
当k与b互为质数,gcd(ka,b)=gcd(a,b),也就是约掉两个数中只有其中一个含有的因子不影响最大公约数。特殊地,当k=2时,说明计算一个偶数和一个奇数的最大公约数时,可以先将偶数除以2。
C语言实现:
1 int Stein_GCD(int x, int y) 2 { 3 if (x == 0) return y; 4 if (y == 0) return x; 5 if (x % 2 == 0 && y % 2 == 0) 6 return 2 * Stein_GCD(x >> 1, y >> 1); 7 else if (x % 2 == 0) 8 return Stein_GCD(x >> 1, y); 9 else if (y % 2 == 0)10 return Stein_GCD(x, y >> 1);11 else12 return Stein_GCD(min(x, y), fabs(x - y));13 }
阅读全文
0 0
- 今天真是醉了,来看一下除了穷举之外的求最大公约数算法.
- 求最大公约数的算法
- 求最大公约数的算法
- 求最大公约数的算法
- 求最大公约数的算法
- 求最大公约数的算法
- 求最大公约数的算法
- 求最大公约数的算法
- 今天可真是够郁闷的了!!!
- 今天离开了一家公司,发现除了备份之外,还很需要清除信息工具
- 7.7.1 穷举法求两个正整数的最大公约数
- [算法]求最大公约数的算法
- 求最大公约数的快速算法
- 求最大公约数的快速算法
- 求最大公约数的欧几里德算法
- 求最大公约数的高效率算法
- 求GCD(最大公约数)的算法
- 求最大公约数,最小公倍数的算法
- 欢迎使用CSDN-markdown编辑器
- js开发:javascript操作数组
- hihoCoder #1015 KMP算法
- DLL 函数导出的规则和方法
- hdoj1049 Climbing Worm(小学数学题)
- 今天真是醉了,来看一下除了穷举之外的求最大公约数算法.
- Problem D: 平面上的点和线——Point类、Line类 (V)
- 优化Mysql数据库的8个方法
- svn中检出项目后无法部署
- Oracle日期数据查询
- JAVA基础觉得很好所以....
- 使用mysql profiling功能剖析单条查询
- 浅谈HTTP传输(一些基本概念)
- Android加载SVG实现交互式地图绘制