快速幂,积取模总结
来源:互联网 发布:ubuntu 优麒麟区别 编辑:程序博客网 时间:2024/06/16 03:37
求a^b%c(这就是著名的RSA公钥的加密方法)
当a,b很大时,直接求解这个问题不太可能
当a,b很大时,直接求解这个问题不太可能
算法1:利用a*b%c=((a%c)*b)%c,这样每一步都进行这种处理,这就解决了a^b可能太大存不下的问题,但这个算法的时间复杂度依然没有得到优化
算法2:另一种算法利用了二分的思想,可以达到O(logn)。
可以把b按二进制展开为:b = p(n)*2^n + p(n-1)*2^(n-1) +…+ p(1)*2 + p(0)
其中p(i) (0<=i<=n)为 0 或 1
这样 a^b = a^ (p(n)*2^n + p(n-1)*2^(n-1) +...+ p(1)*2 + p(0))
= a^(p(n)*2^n) * a^(p(n-1)*2^(n-1)) *...* a^(p(1)*2) * a^p(0)
对于p(i)=0的情况, a^(p(i) * 2^(i-1) ) = a^0 = 1,不用处理
我们要考虑的仅仅是p(i)=1的情况
化简:a^(2^i) = a^(2^(i-1) * 2) = ( a^( p(i) * 2^(i-1) ) )^2(这里很重要!!)
利用这一点,我们可以递推地算出所有的a^(2^i)
当然由算法1的结论,我们加上取模运算a^(2^i)%c = ( (a^(2^(i-1))%c) * a^(2^(i-1))) %c
于是再把所有满足p(i)=1的a^(2^i)%c按照算法1乘起来再%c就是结果
即二进制扫描从最高位一直扫描到最低位。
代码:
算法2:另一种算法利用了二分的思想,可以达到O(logn)。
可以把b按二进制展开为:b = p(n)*2^n + p(n-1)*2^(n-1) +…+ p(1)*2 + p(0)
其中p(i) (0<=i<=n)为 0 或 1
这样 a^b = a^ (p(n)*2^n + p(n-1)*2^(n-1) +...+ p(1)*2 + p(0))
= a^(p(n)*2^n) * a^(p(n-1)*2^(n-1)) *...* a^(p(1)*2) * a^p(0)
对于p(i)=0的情况, a^(p(i) * 2^(i-1) ) = a^0 = 1,不用处理
我们要考虑的仅仅是p(i)=1的情况
化简:a^(2^i) = a^(2^(i-1) * 2) = ( a^( p(i) * 2^(i-1) ) )^2(这里很重要!!)
利用这一点,我们可以递推地算出所有的a^(2^i)
当然由算法1的结论,我们加上取模运算a^(2^i)%c = ( (a^(2^(i-1))%c) * a^(2^(i-1))) %c
于是再把所有满足p(i)=1的a^(2^i)%c按照算法1乘起来再%c就是结果
即二进制扫描从最高位一直扫描到最低位。
代码:
算法3:与算法二类似,直接利用位运算将b进行移位,该算法减少了递归资源的占用,比算法二要好一点
代码:
代码:
快速求积,快速求幂,大指数取模
传说中的O(lgn)时间的快速算术算法和超大整数的取模算法。
1.快速求积,a*b=a*2*b/2
int fast_mul(int a, int b){
int m = 0;
while(b){
if(b & 0x01){
//a*b = a+a(b-1)
m += a;
--b;
}else{
//a*b = a*2*b/2
a <<= 1;
b >>= 1;
}
}
return m;
}
2.快速求幂,a^e=a^(2*e/2)
int fast_exp(int a, int e){
int exp = 1;
while(e){
if(e & 0x01){
//a^e = a*a^(e-1)
exp *= a;
--e;
}else{
//a^e = a^(2*e/2)
a *= a;
e >>= 1;
}
}
return exp;
}
3.大整数取模,(a*b)%m = ((a%m)*(b%m))%m
int bigint_mod(int a, int n, int m){
int mod = 1;
while(n){
if(n & 0x01){
//(a*b)%m = ((a%m)*b)%m = (a*(b%m))%m = ((a%m)*(b%m))%m
mod = (mod*a)%m;
}
a = (a*a)%m;
n >>= 1;
}
return mod;
}
传说中的O(lgn)时间的快速算术算法和超大整数的取模算法。
1.快速求积,a*b=a*2*b/2
int fast_mul(int a, int b){
int m = 0;
while(b){
if(b & 0x01){
//a*b = a+a(b-1)
m += a;
--b;
}else{
//a*b = a*2*b/2
a <<= 1;
b >>= 1;
}
}
return m;
}
2.快速求幂,a^e=a^(2*e/2)
int fast_exp(int a, int e){
int exp = 1;
while(e){
if(e & 0x01){
//a^e = a*a^(e-1)
exp *= a;
--e;
}else{
//a^e = a^(2*e/2)
a *= a;
e >>= 1;
}
}
return exp;
}
3.大整数取模,(a*b)%m = ((a%m)*(b%m))%m
int bigint_mod(int a, int n, int m){
int mod = 1;
while(n){
if(n & 0x01){
//(a*b)%m = ((a%m)*b)%m = (a*(b%m))%m = ((a%m)*(b%m))%m
mod = (mod*a)%m;
}
a = (a*a)%m;
n >>= 1;
}
return mod;
}
- 快速幂,积取模总结
- 总结--快速幂
- 矩阵快速幂总结
- 快速幂总结
- 快速幂总结(一)
- 矩阵快速幂总结
- 矩阵快速幂 总结
- 矩阵快速幂总结
- 矩阵快速幂(总结)
- 快速幂总结,最大公约数(递归版本)
- 最简短总结-整数快速幂
- JQuery快速学习方法总结
- 总结的快速排序
- 快速排序总结
- Python快速教程总结
- 【排序总结--插入、快速】
- 快速排序总结
- 快速排序总结_legend
- 【重新上本科】再写main函数
- 凡要成为高手的人,首先要有骨气
- 重定向
- PS图片去除文字的六个方法(图文)
- 双数据源多线程图片下载(google、flickr)
- 快速幂,积取模总结
- 关于php学习笔记四
- Servlet的生命周期
- PHP的三种安装模式的比较
- 数据库主动式通知编程
- Java使用MulticastSocket实现多点广播
- 通过Eclipse的插件自动下载和更新功能
- Linux安装Tomcat6.0
- HDU 4261 Estimation (DP+优先队列优化)