【位运算|矩阵快速幂思想|倍增】

题目1

求 a 的 b 次方对 p 取模的值。

输入格式
三个整数 a,b,p ,在同一行用空格隔开。

输出格式
输出一个整数，表示a^b mod p的值。

数据范围
0≤a,b,p≤109
数据保证 p≠0
输入样例：
3 2 7
输出样例：
2

一、思想

底数倍增
把指数缩小，底数放大，从而进行优化
例如
3^10=3333333333
3^10=(33)(33)(33)(33)(33)
3^10=(3*3)5
3¹⁰⁼⁹5
9^5=（94）（9^1）
9^5=（94）（9^1）
9^5=（65611）*(9^1)
规律：结果等于指数为1时，放大的底数相乘
时间复杂度为log2（b）

二、代码

*/ 
//用更快的位运算 
int main()
{
	long long res =1;
	long long base;
	long long x;
	cin>>base>>x;
	
	if(x==0)
	{
		cout<<base<<endl;
	}
	
	while(x>0)
	{
		if(x&1) //为奇数 分离出一个base*到结果上 位与运算奇数为1 
		{
			res=res*base%1000;
		}
		x = x>>1 ; //相当于 x-1 后再x/2 
		base=(base*base)%1000;  //指数缩小两倍 底数扩大 
	} 
	cout<<res<<endl;
	return 0;
}

题目2

求 a 乘 b 对 p 取模的值。

输入格式
第一行输入整数a，第二行输入整数b，第三行输入整数p。

输出格式
输出一个整数，表示a*b mod p的值。

数据范围
1≤a,b,p≤1018
输入样例：
3
4
5
输出样例：
2

一、思想

乘数倍增
例如
a为乘数 b为被乘数
b=10时：
a * 10
=(2a)*5
=(2a)+(2a)*4
=(2a)+(4a)*2
=(2a)+(8a)*1
=(2a)+(8a)+(16a)*0
=10a

总结
当b为奇数时，b=b-1 然后 结果加上当前被放大的乘数
当b为偶数时，b=b/2 然后 乘数放大2倍

时间复杂度为log2（b）

二、代码

*/ 
//用更快的位运算 
#include<iostream>
using namespace std;
int main()
{
	long long a;
	long long b;
	long long p;
	cin>>a>>b>>p;
	long long res;

	while(b)
	{
		if(b & 1)
		{
			res = (res+a) %p;
		}
		b=b>>1;
		a=a*2%p;
	}
	cout<<res<<endl;
 }