快速幂

快速幂

例题： 洛谷1226

题目描述

输入b，p，k的值，求b^p mod k的值。其中b，p，k*k为长整型数。

输入格式

三个整数b,p,k.

输出格式

输出“b^p mod k=s”
s为运算结果

输入输出样例

输入#1

2 10 9

输出#1

2^10 mod 9=7

//https://www.luogu.org/problem/P1226
#include<cmath>
#include<bits/stdc++.h>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<string>
using namespace std;
#define LL long long
#define MOD 10000000
inline LL read()
{
	LL x=0,w=1;
	char ch=0;
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')
			w=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		x=x*10+ch-'0';
		ch=getchar();
	}
	return w*x;
}
LL b,p,k,s;
LL quickpow(LL x,LL y,LL m)//二分思想  
{
	x=x%m;//防止第一个x*x就爆掉longlong 
	LL ans=1;//记录答案 
	while(y)// 注意如果y=0的话不会进入循环 
	{
		if(y&1)//判断y是否是单数 
		{
			ans=ans*x;
			ans=ans%m;
		}
		x=x*x;
		x=x%m;
		y=y>>1;
	}
	ans=ans%m;//防止计算1^0%1=0的情况 
	return ans;
}
int main()
{
	b=read();
	p=read();
	k=read();
	s=quickpow(b,p,k);
	printf("%lld^%lld mod %lld=%lld\n",b,p,k,s);
	return 0;
}