快速幂输出最后三位数的问题
给定一个A,B,求A的B次方的值的最后三位
快速幂 解题思路:
单纯的靠暴力去循环,复杂度是O(B),如果B=100000000,那么就要算很久,现在来假设3^10
3^10=3*3*3*3*3*3*3*3*3*3
//尽量想办法把指数变小来,这里的指数为10
3^10=(3*3)*(3*3)*(3*3)*(3*3)*(3*3)
3^10=(3*3)^5
3^10=9^5
//此时指数由10缩减一半变成了5,而底数变成了原来的平方,求3^10原本需要执行10次循环操作,求9^5却只需要执行5次循环操作,但是3^10却等于9^5,我们用一次(底数做平方操作)的操作减少了原本一半的循环量,特别是在幂特别大的时候效果非常好
//现在我们的问题是如何把指数5变成原来的一半,5是一个奇数,5的一半是2.5,但是我们知道,指数不能为小数,因此我们不能这么简单粗暴的直接执行5/2,然而,这里还有另一种方法能表示9^5
9^5=(9^4)*(9^1)
//此时我们抽出了一个底数的一次方,这里即为9^1,这个9^1我们先单独移出来,剩下的9^4又能够在执行“缩指数”操作了,把指数缩小一半,底数执行平方操作
9^5=(81^2)*(9^1)
//把指数缩小一半,底数执行平方操作
9^5=(6561^1)*(9^1)
//此时,我们发现指数又变成了一个奇数1,按照上面对指数为奇数的操作方法,应该抽出了一个底数的一次方,这里即为6561^1,这个6561^1我们先单独移出来,但是此时指数却变成了0,也就意味着我们无法再进行“缩指数”操作了。
9^5=(6561^0)*(9^1)*(6561^1)=1*(9^1)*(6561^1)=(9^1)*(6561^1)=9*6561=59049
我们能够发现,最后的结果是9*6561,而9是怎么产生的?是不是当指数为奇数5时,此时底数为9。那6561又是怎么产生的呢?是不是当指数为奇数1时,此时的底数为6561。所以我们能发现一个规律:最后求出的幂结果实际上就是在变化过程中所有当指数为奇数时底数的乘积
#include <iostream>
using namespace std;
int ksm(long long a, long long b)
{
long long m = 1;
while(b > 0)
{
/*if(b & 1)//if(b % 2 ==0)
{
//这里是记录对底数进行平方的操作,然后可以将指数除以二
m = m * a % 1000;
b /= 2;
}
else
{
b -= 1;
b /= 2;
a = a * a % 1000;
}
}*/
if(b & 1)//if(b % 2 ==0)
{
//这里是记录对底数进行平方的操作,然后可以将指数除以二
m = m * a % 1000;
}
//这里的优化是因为 b-1 后 b/2 等于 b/2
//b /= 2;
//这里的优化是因为b往后移动一位就是除以2
b >>= 1;
a = a * a % 1000;
}
//m是指数为奇数的情况下,所有底的乘积,且每次相乘都取模1000保留后三位
return m;
}
int main()
{
int a, b;
cin >> a >> b;
cout << ksm(a, b);
return 0;
}