题目:牛牛和牛可乐的赌约
题目分析:这是一个快速幂+逆元的板子题,逆元就是除法,因为计算机不擅长做除法,所以我们需要把它转化成乘法逆元,求逆元可以用费马小定理,这个定理证明了一个公式,我们可以用这个公式求得逆元。
费马小定理:假如p是质数,且(a,p)=1,那么 a(p-1)≡1(mod p)。即:假如a是整数,p是质数,且a,p互质(即两者只有一个公约数1),那么a的(p-1)次方除以p的余数恒等于1。
当涉及取模运算的计算中,如果有除法,不能直接除以一个数,而应该变成乘以它的乘法逆元。
当我们除以一个数n时,也就是乘上1/n,若x是1/n关于模N的逆元,则x=1/n (mod N),即 x*n=1(mod N)。由于我们做题时N常常为1000000007,而1000000007是个素数,
所以它满足了费马小定理,而满足费马小定理说明解唯一,所以我们可以直接得出x*n=n ^ (N-1)。那么x=n^(N-2),即为1/n关于模N的乘法逆元
看懂了上面那么题目就是一个简单的套题了
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7;
const int maxn = 1e6 + 10;
ll fast(ll a, ll b)
{
ll ans = 1;
while (b)
{
if (b & 1)
ans = (ans * a) % mod;
a = (a * a) % mod;
b >>= 1;
}
return ans;
}
ll inv(ll b)
{
return fast(b, mod - 2) % mod;
}
int main()
{
ll t, n, m;
scanf("%lld", &t);
while (t--)
{
scanf("%lld%lld", &n, &m);
ll x = inv(n);
printf("%lld\n", (mod + 1 - fast(x, m)) % mod);
}
}
时间并不会因为你的迷茫和迟疑而停留,就在你看这篇文章的同时,不知道有多少人在冥思苦想,在为算法废寝忘食,不知道有多少人在狂热地拍着代码,不知道又有多少提交一遍又一遍地刷新着OJ的status页面…… 没有谁生来就是神牛,而千里之行,始于足下! |
---|