乘法逆元模板
来源:互联网 发布:中国网络图书馆 编辑:程序博客网 时间:2024/05/17 11:05
我们知道,在模意义下,加减乘运算是封闭的;
例如在模17意义下,8+11=((8%17)+(11%17))%17=2,19+21=((19%17)+(21%17))=6;
减法同理,但是在C++编译器中取模的意义是与数学上不同的,负数mod正数应该是非负数,但是C++却可能模出负数,所以我们要多加一些处理,对于x,可以(x%mod+mod)%mod保证其正确性,比如在模17意义下,21-19=((21%17)-(19%17))=2,35-21=((35%17)-(21%17))=14.
乘法也是这样,先模再乘,乘完再模,其意义不会变化。
但是除法不是。
于是我们需要用到除法逆元。
关于逆元的定义。
定义 对a∈Zm,存在b∈Zm,使得a+b ≡ 0 (mod m),则b是a的加法逆元, 记b= - a。定义 对a∈Zm,存在b∈Zm,使得a×b ≡1 (mod m),则b为a的乘法逆元 ,记b=1/a;
加法逆元很少用啊,那么怎么求乘法逆元呢?
方法一:我们要找到a×b ≡1 (mod m)的b,所以很容易想到用扩展欧几里得求出a*b+m*y=1的一组解,b就是a在模m意义下的乘法逆元。当然前提是a与m互质,因为只有这样gcd(a,m)才是1
ll exgcd(ll a,ll b,ll &x,ll &y){if(b==0){x=1;y=0;return a;}int r=exgcd(b,a%b,x,y);ll t=y;y=x-(a/b)*y;x=t;return r;}ll inv(ll a){ll x,y;ll r=exgcd(b,mod,x,y);if(r==1) return (x%mod+mod)%mod;return -1;}
方法二:根据费马小定理,a^(p-1)≡ 1(mod p)
所以稍加变形,a*a^(p-2)≡ 1(mod p)
即a^(p-2)是a在模p意义下的乘法逆元。
快速幂一下就好了,复杂度O(log p)
方法三:O(n)线性预处理
inv[1]=1;for(int i=2;i<=n;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;可以再推一遍求出阶乘的逆元。
求出逆元之后,要求(a/b)%mod, 设b`为b在模mod意义下的逆元,直接(a%mod*b`%mod)%mod即可。
1 0
- 乘法逆元模板
- 乘法逆元 模板
- 乘法逆元模板
- [模板]乘法逆元
- 【模板】乘法逆元
- [模板]乘法逆元
- ACM 乘法逆元 模板
- [P3811][模板]乘法逆元
- 洛谷 P3811 【模板】乘法逆元
- 洛谷 P3811 【模板】乘法逆元
- 51nod 1256 乘法逆元 (模板)
- HDU 3037:Saving Beans Lucas模板题+乘法逆元
- Loj-110乘法逆元(线性模板题)
- codeforces 300C 乘法逆元 (乘法逆元模为素数的模板)
- 乘法逆元
- 乘法逆元
- 关于乘法逆元
- 乘法逆元
- angularJS依赖注入(二)服务对比
- 二维map找出值中值得最大值,借助二维数组
- CSS3 Filter的十种特效
- 中文字体的英文名称 css 使用
- Android:仿微博Listview与多线程功能(JSON、数据下载、图片加载、Listview、异步任务)
- 乘法逆元模板
- Codeforces Round #287 (Div. 2)-D. The Maths Lecture(数位dp)
- container_of函数原理分析
- Android三种播放视频的方式
- UIWebView、WKWebView使用详解及性能分析
- LeetCode 9. Palindrome Number
- Android:控件WebView显示网页
- ObjectAid UML Explorer for Eclipse
- 【Caffe安装】Caffe安装(有GPU)--wanglei