时间限制:1000 ms | 内存限制:65535 KB
难度:3
-
描述
-
求a的b次方对c取余的摸
-
输入
-
第一行输入一个整数n表示测试数据的组数(n<100)
每组测试只有一行,其中有三个正整数a,b,c(1=<a,b,c<=1000000000)
输出
- 输出a的b次方对c取余之后的结果 样例输入
-
3 2 3 5 3 100 10 11 12345 12345
样例输出
-
3 1 10481
-
算法:快速求幂:
求a^b%c(RSA公钥的加密方法)算法1:公式a*b%c=((a%c)*b)%c, 每一步都进行这种处理,这就解决了a^b可能太大存不下的问题,但这个算法的时间复杂度没有优势。#include<iostream> using namespace std; long long f(long long a,long long b,long long c) { long long t,r; if(b==1) return a%c; t=f(a,b/2,c); r=(t*t)%c; if(b%2==1) r=(r*a)%c; return r; } int main() { int n,m; long long a,b,c; cin>>n; while(n--) { cin>>a>>b>>c; cout<<f(a,b,c)<<endl; } }算法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就是结果,即二进制扫描从最高位一直扫描到最低位 我们是把b进行二进制展开(就是转换为二进制)比如这个b就变成了1010110111111, 这样每次只需要求s1 mod c,s2mod c,... sn^(k-1) mod c,然后把对应位置上的二进制是1的项乘起来。在每次乘完以后求除c的余 #include <stdio.h> #include <stdlib.h> #include <string.h> long long a,b,c,m,temp,x; int main() { int t,i; scanf("%d",&t); while(t--) { scanf("%lld%lld%lld",&a,&b,&c); temp=a%c;m=1;x=1; while(b!=0)//用二分快速求幂的方法将复杂度 下降到log(2)c { if(b%2) {m=(temp*x)%c;x=m;} temp=(temp*temp)%c; b/=2; } printf("%lld\n",m); } return 0; }
-
第一行输入一个整数n表示测试数据的组数(n<100)