牛客小白月赛12 - B - 快速幂+快速乘
程序员文章站
2022-03-24 15:45:56
...
题目描述
找到了心仪的小姐姐月月后,华华很高兴的和她聊着天。然而月月的作业很多,不能继续陪华华聊天了。华华为了尽快和月月继续聊天,就提出帮她做一部分作业。
月月的其中一项作业是:给定正整数A、B、P,求ABmodP的值。华华觉得这实在是毫无意义,所以决定写一个程序来做。但是华华并不会写程序,所以这个任务就交给你了。
因为月月的作业很多,所以有T组询问。
输入描述:
第一行一个正整数T表示测试数据组数。 接下来T行,每行三个正整数A、B、P,含义如上文。
输出描述:
输出T行,每行一个非负整数表示答案。
示例1
输入
2 2 5 10 57284938291657 827493857294857 384729583748273
输出
2 18924650048745
备注:
1≤T≤103
思路:
快速幂+快速乘
快速乘:主要解决大数相乘取模,如两个1e18的数相乘,模数也是1e18级别的,就算取模,还是会爆long long。这时候就要用到快速乘。
说明一下,快速乘并不是很快(肯定没直接乘快O(1)),它是把a*b换成b个a相加,用来解决取模爆long long问题
代码如下:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>
using namespace std;
#define ll long long
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
typedef pair<int,int>P;
const int INF=0x3f3f3f3f;
const int N=1000010,mod=32767;
ll ksc(ll a,ll b,ll p){//a*b%p
ll res=0;
while(b){
if(b&1)res=(res+a)%p;
b>>=1;
a=(a+a)%p;
}
return res;
}
ll ksm(ll a,ll b,ll p){
ll res=1;
while(b){
if(b&1)res=ksc(res,a,p)%p;
b>>=1;
a=ksc(a,a,p)%p;
}
return res;
}
int main(){
ll a,b,p;
int t;
scanf("%d",&t);
while(t--){
scanf("%lld%lld%lld",&a,&b,&p);
ll ans=ksm(a,b,p);
printf("%lld\n",ans);
}
}
非常规操作:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>
using namespace std;
#define ll long long
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
typedef pair<int,int>P;
const int INF=0x3f3f3f3f;
const int N=1000010,mod=32767;
ll ksm(ll aa,ll b,ll p){
__int128 res=1,a=aa%p;
while(b){
if(b&1)res=((res%p)*(a%p))%p;
b>>=1;
a=((a%p)*(a%p))%p;
}
return res;
}
int main(){
int t;
ll a,b,p;
scanf("%d",&t);
while(t--){
scanf("%lld%lld%lld",&a,&b,&p);
printf("%lld\n",ksm(a,b,p));
}
}
上一篇: 基于分治法的---快速排序
下一篇: 说的就是你