链接:https://ac.nowcoder.com/acm/contest/275/A
来源:牛客网
转自:出处
复习内容:逆元的相关知识
在开始之前我们先介绍3个定理:
1.乘法逆元(在维基百科中也叫倒数,当然是 mod p后的,其实就是倒数不是吗?):
如果ax≡1 (mod p),且gcd(a,p)=1(a与p互质),则称a关于模p的乘法逆元为x。
2.费马小定理(定义来自维基百科):
如果a不是p的倍数,这个定理也可以写成
3.
(定义来自维基百科):
已知整数a、b,扩展欧几里得算法可以在求得a、b的最大公约数的同时,能找到整数x、y(其中一个很可能是负数),使它们满足贝祖等式。
好了,在明白上面的定理后我们开始分析乘法逆元:ax≡1 (mod p) 这个等式用中文描述就是 a乘一个数x并模p等于1,即 a%p*x%p=res,res%p=1;看上去就是同余定理的一个简单等式- -。那么问题来了。
为什么可以用费马小定理来求逆元呢?
由费马小定理 ap-1≡1 , 变形得 a*ap-2≡1(mod p),答案已经很明显了:若a,p互质,因为a*ap-2≡1(mod p)且a*x≡1(mod p),则x=ap-2(mod p),用快速幂可快速求之。
为什么可以用扩展欧几里得求得逆元?
我们都知道模就是余数,比如12%5=12-5*2=2,18%4=18-4*4=2。(/是程序运算中的除)
那么ax≡1 (mod p)即ax-yp=1.把y写成+的形式就是ax+py=1,为方便理解下面我们把p写成b就是ax+by=1。就表示x是a的模b乘法逆元,y是b的模a乘法逆元。然后就可以用扩展欧几里得求了。
知道逆元怎么算之后,那么乘法逆元有什么用呢?
做题时如果结果过大一般都会让你模一个数,确保结果不是很大,而这个数一般是1e9+7,而且这个数又是个素数,加减乘与模运算的顺序交换不会影响结果,但是除法不行。有的题目要求结果mod一个大质数,如果原本的结果中有除法,比如除以a,那就可以乘以a的逆元替代。(除一个数等于乘它的倒数,虽然这里的逆元不完全是倒数,但可以这么理解,毕竟乘法逆元就是倒数的扩展)。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
void exgcd(ll a,ll b,ll& d,ll& x,ll& y)
{
if(!b) { d = a; x = 1; y = 0; }
else{ exgcd(b, a%b, d, y, x); y -= x*(a/b); }
}
ll inv(ll a, ll p)
{
ll d, x, y;
exgcd(a, p, d, x, y);
return d == 1 ? (x+p)%p : -1;
}
int main()
{
ll a,p;
while(1)
{
scanf("%lld %lld",&a,&p);
printf("%lld\n",inv(a,p));
}
}
题目描述
你在一栋楼房下面,楼房一共有n层,第i层每秒有pi的概率会扔下一个东西并砸到你
求第一秒内你被砸到的概率
输入描述:
第一行一个整数n 之后有n行,第i+1行有两个整数ai,bi,表示
输出描述:
设答案为,你只需要找到一个最小的非负整数T,使得 输出这个T就行了
示例1
输入
复制
2 1 2 1 2
输出
复制
750000006
说明
一共只有如下状态: 1. 第一层和第二层都扔了下来 2. 第一层扔了下来 3. 第二层扔了下来 4. 第一层和第二层都没有扔下来 以上四种都是等概率发生的 除了第四种情况外,都会被砸到 因此被砸到的概率是 3/4,这个值在模1e9+7意义下就是750000006
备注:
数据范围 0 ≤ n ≤ 105 1 ≤ ai ≤ bi ≤ 105
此代码为学长的,自己的代码总是不对,就一点一点改,到头来忘了取模了。小辣鸡“”“”“”“”“”
#include<bits/stdc++.h>
using namespace std;
int mod=1e9+7;
typedef long long ll;
inline ll pow_mod(ll n,ll k)//n^k%mod
{
ll res=1;
n=n%mod;
while(k>0)
{
if(k&1)
res=res*n%mod;
n=n*n%mod;
k>>=1;
}
return res;
}
int main()
{
ll n,p=1,q=1,a,b;
scanf("%lld",&n);
for(int i=1;i<=n;i++)
{
scanf("%lld%lld",&a,&b);
p=p*(b-a)%mod;
q=q*b%mod;
}
ll ans=((q-p+mod)*pow_mod(q,mod-2))%mod;
printf("%lld\n",ans);
return 0;
}
同上:不用看
#include<bits/stdc++.h>
using namespace std;
int mod=1e9+7;
typedef long long ll;
inline ll pow_mod(ll n,ll k)//n^k%mod
{
ll res=1;
n=n%mod;
while(k>0)
{
if(k&1)
res=res*n%mod;
n=n*n%mod;
k>>=1;
}
return res;
}
int main()
{
ll n,a=1,b=1,a1,b1;
scanf("%lld",&n);
for(int i=1;i<=n;i++)
{
scanf("%lld%lld",&a1,&b1);
a=a*(b1-a1)%mod;
b=b*b1%mod;
}
ll temp=b-a+mod;
printf("%lld\n",((temp*pow_mod(b,(mod-2)))%mod));
return 0;
}