快速幂
1、快速幂的作用
快速幂可以帮助我们在 O ( log k ) O(\log k) O(logk)的时间复杂度内计算出 a k m o d p a^k\quad mod \quad p akmodp的结果。
2、算法原理
首先快速幂我们要处理出
log
k
+
1
\log k+1
logk+1个数,
a
2
0
m
o
d
p
,
a
2
1
m
o
d
p
,
a
2
2
m
o
d
p
,
.
.
.
,
a
2
log
k
m
o
d
p
a^{2^0}mod\quad p,a^{2^1}mod\quad p,a^{2^2}mod\quad p,...,a^{2^{\log k}}mod\quad p
a20modp,a21modp,a22modp,...,a2logkmodp处理完这些数之后就可以计算
a
k
a^k
ak.我们可以通过二进制的形式去凑。先把k转化成二进制形式,然后在出现1的位置乘上相应的我们前面预处理的数即可。
举个例子吧:
假设要我们算一下
4
5
m
o
d
10
4^5\quad mod \quad 10
45mod10的结果。
首先我们处理出
4
2
0
m
o
d
10
=
4
4^{2^0}\quad mod \quad 10=4
420mod10=4,
4
2
1
m
o
d
10
=
6
4^{2^1}\quad mod \quad 10=6
421mod10=6,
4
2
2
m
o
d
10
=
6
4^{2^2}\quad mod \quad 10=6
422mod10=6.然后5的二进制表示是101,我们在是1的位置上乘上相应的数,及
4
2
0
+
2
2
4^{2^0+2^2}
420+22,展开就是
4
2
0
+
4
2
2
4^{2^0}+4^{2^2}
420+422查表得
(
4
∗
6
)
m
o
d
10
=
4
(4*6)\quad mod\quad 10=4
(4∗6)mod10=4.然后就可以啦~
例题
快速幂例题
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
long long qmi(int a, int k, int p)
{
long long res = 1;
while (k)
{
if (k & 1) //如果这一位是1的话我们就乘上相应的数
res = (long long)res * a % p;
k >>= 1; //去掉最后一位
a = (long long)a * a % p; //每次算出相应的a值,及上一个值的二次方
}
return res;
}
void solve()
{
int a, k, p;
cin >> a >> k >> p;
cout << qmi(a, k, p) << endl;
}
int main()
{
int T;
cin >> T;
while (T--)
{
solve();
}
return 0;
}
3、逆元
我们要求
a
b
≡
a
×
x
(
m
o
d
p
)
{a\over b}\equiv a\times x \quad (mod \quad p)
ba≡a×x(modp)我们两边同时乘上一个b得
a
≡
a
×
x
×
b
(
m
o
d
p
)
a \equiv a\times x \times b\quad(mod\quad p)
a≡a×x×b(modp)我们约去a得到:
b
×
x
≡
1
(
m
o
d
p
)
b\times x\equiv 1\quad(mod\quad p)
b×x≡1(modp),然后b和p互质我们可以想到费马小定理那么x就等于
b
p
−
2
b^{p-2}
bp−2,这样的话就转化成了一个快速幂问题。如果大家不清楚费马小定理的话可以看这篇文章。欧拉定理及费马小定理
例题
逆元例题
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
long long pmi(int a, int p)
{
long long res = 1;
int k = p - 2;
while (k)
{
if (k & 1)
res = (long long)a * res % p;
k >>= 1;
a = (long long)a * a % p;
}
return res % p;
}
void solve()
{
int n, p;
cin >> n >> p;
if (n % p)
cout << pmi(n, p) << endl;
else
cout << "impossible" << endl;
}
int main()
{
ios::sync_with_stdio(false);
int T;
cin >> T;
while (T--)
{
solve();
}
return 0;
}