watch out

本文是博主的 csdn 上搬过来的,格式有点崩,看不下去的可以去 博主的 csdn上看(上面 格式会好很多,并且有些公式也用 $\LaTeX$  update 上去了,但是博主也在 cnblogs 上更新了一下...lateX 都用上了,应该不至于不能看)

最近有点颓废啊,写篇blog振作一下…(不过没图的数论blog真是不对我胃口)

emmm …首先介绍一下这是一篇关于数论中较为重要 (主要可能经常要用到) 的一个知识分支—逆元
相信你听到数论之后可能鼠标就想往网页上的 X 键上移了,但是还是劝你看一看, 总能有点收获的吧
(反正我觉得自己讲的相对网上其他 blog 还算是蛮清楚的了)

解释一下逆元:

首先同学们小学的时候都学过除法吧(这个问题有点奇怪,滑稽)?
那么你一定知道,“除以一个数等同于乘以这个数的倒数”对吧?
所以逆元是什么东西呢? 首先这里有个式子: $({a\over b}) \% p$,这个式子的答案怎么求?
没错,暴力求是一种方法,但是当 $b$ 非常大的时候呢 ? 这个时候就要用到逆元了

这里有一句口诀 “除以一个数再取模等同于乘以这个数的逆元再取模 ”(自创) 
什么意思? 设 $inv[b]$ 是 $b$ 的逆元, 那么 $({a\over b}) \%p = (a*inv[b]) \%p$  
这下就你能大概理解逆元是个什么东西(以及它大概是用来干什么的)了吧 ?

那么逆元怎么求 ? 首先这里得纠正一点, 我们不能直接说一个数的逆元是多少,
应该这么说: 一个数 $x$ 在模 $p$ 的条件下的逆元是多少.(这点概念还是蛮重要的)
其次,我们不难得知一个数的逆元有无穷多个,但是我们只需要求得一个数的最小正整数逆元就行了

另外,一个数 $x$ 在模 $p$ 的条件下不一定有逆元, $x$ 关于 $p$ 的逆元存在 当且仅当 $x$ 和 $p$ 互质
这里有一个推导: (设 $a$ 为 $x$ 的逆元, $b$ 为任意整数)

$x * a ≡ 1 (mod\ p)$ = 将p连入式子=> $x*a = 1 -b*p$ => $x*a+b*p=1$

那么我们就得到了一组新方程: $x*a+b*p=1$ 若$x$ 和 $p$不互质, 则 $x$和 $p$ 存在公约数 $d=gcd(x,p)>1$

提取出$d$, 得到: $d(x/d*a+p/d*b) = 1$ 移项得到: $(x/d * a+ p/d *b) = 1/d$

易知 $x$ , $p$ 能整除 $d$, 所以括号内定为整数, 又因为$d > 1$ ,等式右边必为真分数, 等式无解

证毕

那么说了半天,逆元到底怎么求呢 ? 别着急,往下看—-
求逆元有好几种方法, 但是在这里我们只讨论两种方法(一个偷懒用的 ,一个不容易翻车的),
不过在此之前我们先来讲讲费马小定理 (别害怕,这玩意儿知道怎么用就好,不会叫你推的)

这是费马小定理的原始状态: $a^p ≡ a (mod\ p)$ ,其中 $a$ , $p$ 互质 [感性理解一下 a 一直乘 a  会出现循环节,而循环节必然是 p-1 的因数,所以乘着乘着就变回来了]
然后, 原式 =移项 => $a^p - a ≡0 (mod\ p) $=提取 a => $a*(a^{p-1} - 1) ≡0 (mod\ p)$  
=两边同除以 a => $a^{p-1}-1 ≡0 (mod\ p)$ =再次移项 => $a^{p-1} ≡1 (mod\ p)$ 
于是 费马小定理的结论就产生了(注意,这个结论很重要),求逆元时经常要用(主要偷懒求逆元时必用)


好了,现在我们再来考虑逆元的求法

一、 蒙哥马利快速幂模(偷懒求逆元大法):

首先听到这么高大上(主要是长)的名字,相信你的内心一定是崩溃的,仿佛见到了高斯定理或者欧拉定理一样(别着急啊,扩欧等下才讲)
但其实蒙哥马利快速模这个算法并没有那么高大上,它的局限性很大,只有在 $p$ 是质数的情况下才可以使用
首先我们设 $inv(a)$ 是 $a$ 的逆元 那么根据定义, $inv(a) * a ≡1 (mod\ p)$ 
再根据 费马小定理 $a^{p-1} ≡1$ , 易得 $inv(a) * a ≡a^{p-1} (mod\ p)$ 
移项,得: $inv(a) ≡ a^{p-2}$ 
于是我们得到了快速幂模算法的一个前提条件: $inv(a) ≡ a^{p-2} (mod\ p)$

下面附上代码(自己即兴手打):
[这里是位运算处理的方法,当然也有递归求的方法,那个代码短一些,不过这个代码我打的顺手]

 //by Judge
inline ll quick_pow(ll x,rint p){
ll res=;
while(p){
if(p&) res=(res*x)%mod;
x=(x*x)%mod, p>>=;
}
return res;
}
inline inv(ll a){
inv_a=quick_pow(a,mod-);
return inv_a;
}

下面我们不着急讲第二种算法,先来讲讲这种快速模算法的应用:
求 $C(n,r) \% p$ 的值 , 其中 $p$ 是大质数(如$1e9+7$) ,那么最后的答案是: 原式 = $n! * (r! * (n-r)!)^{p-2} \%p$(也就相当于 ${n!\over r! *(n-r)!} \% p$ ,就是加了逆元)

推导过程:

$C(n,r) (mod\ p) = {n! \over r! * (n-r)!} (mod\ p)$=多乘一个 1 => $ {n! \over r! * (n-r)!} * 1 (mod\ p)$ ----1式

∵ p是质数,再根据费马小定理可得: $(r! * (n-r)!)^(p-1) ≡1 (mod\ p)$ ----2式

再将2式带入1式中得到: $ C(n,r) ≡ {n! \over r! * (n-r)!} * (r! *(n-r)!)^{p-1} (mod\ p) 指数相消得到 $C(n,r) ≡n! * (r! * (n-r)!)^{p-2} (mod\ p)$

其实换元一下推导会简洁得多...

令 tmp = r! * (n-r)! ,那么有:

$C(n,r)%p = ( n! / tmp ) \%p = ( n! / tmp ) * 1 \%p $

= $( n! * tmp^{-1} ) * tmp^{p-1} \%p$

= $( n! * tmp^{p-2} ) \%p$ = 上述的式子

关于费马小定理的清奇证明可以看这里

不过呢,博主好像有写过费马小定理的花式证明,在这里

代码实现:

 //by Judge
#define ll long long
const ll mod=; //1e9+7或者其他题目给定的大质数
inline ll quick_pow(ll x,rint p){
ll res=;
while(p){
if(p&) res=(res*x)%mod;
x=(x*x)%mod, p>>=;
}
return res;
}
inline C_mod(rint n,rint m){
ll a=,b=;
for(rint i=;i<=n;++i)
a=(a*i)%mod;
for(rint i=;i<=m;++i)
b=(b*i)%mod;
for(rint i=;i<=n-m;++i)
b=(b*i)%mod;
return (a*quick_pow(b,mod-))%mod;
}

然后您可以试试用这个算法 A 掉乘法逆元的模板题(放心没坑,只会T而已,要A的话得用线性推逆元):P3811 【模板】乘法逆元


然后是上文提到过的扩欧

二、扩展欧几里得算法(一般我们都讲扩欧) :

这个东西讲起来真是烦,一开始我还一直以为扩欧就是在求逆元…其实扩欧是求同余方程的解…
首先如果你不是萌新的话,那么你应该看到过ex_gcd这样(恶心)的字眼,
这玩意儿其实是扩欧…(ex 表示扩展, gcd 最大公约数是可以用欧拉求的嘛)

不废话了,进入正题, 首先扩欧公式如下: $ax+by = (a,b)$   [注意这里是等于]
[$(a,b)$ 代表 $gcd(a,b)$ ,$ a$ 和 $b$ 是已知数, 要求的是 $x$ 和 $y$ 的一组解]
那么很容易可以知道, 当 a 和 b 互质的情况下,原来的等式就变成了: $ax+by = 1$  
然后我们怎么用这玩意儿求逆元呢? 我们先考虑求逆元的式子: $a * inv[a] ≡1 (mod\ p)$ 
其中$inv[a]$ 同上表示 $a$ 的逆元, 现在我们令 $x = inv[a] , b = p$;
那么同余的式子就成了: $ax ≡1 (mod b)$ , 然后我们把 $b$ 加入式子中,就成了这样:
$ax+by = 1$ (这里我就不慢慢推了,看博客也是要动动脑子的嘛 )其中 $y$ 和 $x$ 一样,也是未知数
emmmm还是讲讲, $ax ≡1 (mod b)$ => $ax = by+1$ (然后就把同余去掉了,变成了等式)
移项 => $ax+by = 1$ (由于 $y$ 是不确定的,而且我们求的主要是 $x$ 所以 $y$ 的正负性没什么影响)

这样你应该就能看懂了, 扩欧公式中求解的 $x$ 就是 $a$ 的逆元(一切好像都非常的显然了)
ok, 我还是没有讲 扩欧中的 $x$ 和 $y$ 怎么解 (这个才是真正麻烦的地方,解释起来很麻烦) 对吧?
哎, 式子先列出来再说, 我再斟酌斟酌…

$ax+by =(a,b) =1 , bX+(a\%b)Y =(b,a\%b) =1$ (注意大小写,并且显然这里要满足 $a$  $b$ 互质的条件)
=> ………… => $x=Y , y=X-\lfloor {a\over  b}\rfloor*Y$ (其中 $X$ 和 $Y$ 已知)
好了,自己推导吧! (试图想混过去的样子) 咳咳… 首先 $ax+by = bX+(a\%b)Y$ 对吧? (两式合并)
然后我们把 $a\%b$ 转化一下,变成: $ax+by = bX+(a-b*\lfloor{ a\over b} \rfloor ) Y$ (其实就是用了模的定义式)
=> $ax+by = bX+aY-b* \lfloor{ a\over b} \rfloor Y$ =提取公因式 a,b => $a(x-Y) = b(X- \lfloor {a\over b}\rfloor Y-y) $

那么很显然, 我们令 $x = Y$ , 令 $y = X-\lfloor{ a\over b} \rfloor Y$ 就能轻松的构造出一组 x 和 y 的解了

然而现在又出现了一个问题: X 和 Y 怎么求? 其实我们同上的方法求出 X 和 Y 就行了:
$bX+(a\%b)Y=(a\%b)X’ + (b\%(a\%b))Y’$ 如果你还不能理解, 那么我们令 $p = b$, $q = a\%b$ ,
那么原式 => $pX+qY=qX’+(p%q)Y’$ 这不就和之前的式子差不多了吗?
现在你应该很容易就能看出来这是在用递归的方法求解原式的解

但是递归总是有有边界的,不能永远递归下去对吧(非常显然)
那么这个递归式子的边界是什么呢? emmmm 提示一下, 这个边界和 gcd 的边界有点关系
对, 边界就是 b 等于 0 的时候, 我们不能再递归下去了,那么原式就成了:
$1x+0y=1$ [a此时为 $(a,b)$ ,b为 0 ] 这应该很好解吧? 我们令 $x=1$,$y=0$ 就行了(y=0是为了…方便嘛)

然后递归求解同上,下面附上代码(也是即兴手打的)
[下面常规版的,偷懒版和这个其实本质上没有区别,只是简化了一下代码,那么这里就不附上了,有兴趣的同学可以自行百度]

//by Judge
#define ll long long
const int mod=; //同上
void ex_gcd(ll a,ll b,ll &x,ll &y){
if(!b){ x=,y=; return ; }
ex_gcd(b,a%b,x,y);
ll t=x; x=y,y=t-(a/b)*y;
}
//当然你也可以这么写,更能体现公式:
// ll X=x,Y=y;
// x=Y,y=X-(a/b)*Y;
inline ll inv(ll a){
ll inv_a,y;
ex_gcd(a,mod,inv_a,y);
return (inv_a%mod+mod)%mod;
}
 

另外这里补充一下线性推逆元的方法 ($inv[i]$ 表示 $i$ 的逆元 , $fac[i]$ 表示 $i$ 的阶乘 , $p$为模数)

三、线性推逆元

1.求 1! ~ n! 的逆元:(因为我觉的阶乘的逆元反倒比较好推所以就先讲)

公式: $inv[i] ≡inv[i+1]*(i+1) (mod\ p)$

推导过程:

$fac[i] * inv[i] ≡1 (mod\ p)$

$fac[i+1] * inv[i+1] ≡1 (mod\ p) => fac[i]* (i+1) * inv[i+1] ≡1(mod\ p)$

由 同余的除法原理可得 : $inv[i] ≡inv[i+1] * (i+1)$

推导完毕

inline void get_finv(){
fac[]=finv[]=;
for(int i=;i<=n;++i;++i)
fac[i]=fac[i-]*i%mod;
finv[n]=quick_pow(fac[n],mod-);
for(int i=n-;i;--i)
finv[i]=finv[i+]*(i+)%mod;
}
 

2.求 1 ~ n 的逆元(这个还稍微麻烦些)

公式:$inv[i] ≡inv[p%i] * (- p/i) (mod\ p)$

推导过程:

令 $s = p/i , t = p\%i$ , 则有:

$s*i + t = p$ (显然)

然后 $s*i + t ≡ 0 (mod\ p)$

移项得 $t ≡ -s*i (mod\ p)$

同除以 $t * i 得 t / (t*i) ≡ -s*i / (t*i) (mod\ p)$

将除法转化为乘法 => $inv[i] ≡ -s * inv[t] (mod\ p)$

于是将 $s=p/i , t=p\%i$ 带入就可以得到:

$$inv[i] ≡ inv[p\%i] * (-p/i) (mod\ p)$$

推导完毕

代码:

inline void get_inv(){
inv[]=inv[]=;
for(int i=;i<=n;++i)
inv[i]=inv[mod%i]*(mod-mod/i)%mod;
}

然后你还可以用线性逆元求阶乘逆元了

// 这里是一种用线性逆元解阶乘逆元的方法
inline void get_finv2(){
finv[]=finv[1]=;
for(int i=2;i<=n;++i)
finv[i]=finv[mod%i]*(mod-mod/i)%mod;
for(int i=;i<=n;++i)
finv[i]=finv[i]*finv[i-]%mod;
}


如果你看完这些不过瘾,给你几个网站用来生啃代码:
四法求逆元
逆元的几种求法

最后的最后,逆元的这些东西…随便学学就好,到时候用得到(也可能用不到)
于是Judge的课堂又这样水过去了一节…
欢迎下次光临! _(:з」∠)_

(求点赞,求推荐,求收藏! ღ( ´・ᴗ・` )比心 )

逆元知识普及(扫盲篇) —— from Judge的更多相关文章

  1. 逆元知识普及(进阶篇) ——from Judge

    关于一些逆元知识的拓展 刚艹完一道 提高- 的黄题(曹冲养猪) ,于是又来混一波讲解了 ——承接上文扫盲篇   四.Lucas定理(求大组合数取模)   题外话 这里Lucas定理的证明需要用到很多关 ...

  2. Httpd服务入门知识-http协议版本,工作机制及http服务器应用扫盲篇

    Httpd服务入门知识-http协议版本,工作机制及http服务器应用扫盲篇 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任. 一.Internet与中国 Internet最早来源于美 ...

  3. 代理服务器基本知识普及代理IP使用方法!

    本文并未从专业角度进行详细讲解,而是从应用的角度出发来普及一些代理服务器的基本知识.文章明显是搜集多方资料的拼凑,而且比较老了,但往往越老的东西越接近事物的本质,更容易窥探到原理,对于刚接触的人来说, ...

  4. 转摘 MySQL扫盲篇

    一下文章摘自:http://www.jellythink.com/archives/636 MySQL扫盲篇 2014-09-15 分类:MySQL / 数据库 阅读(1412) 评论(1)  为什么 ...

  5. 给 Android 初学者的 Gradle 知识普及

    给 Android 初学者的 Gradle 知识普及:http://gold.xitu.io/entry/5778f8bd165abd0054b443b0/promote?utm_source=bai ...

  6. 分布式协调服务Zookeeper扫盲篇

    分布式协调服务Zookeeper扫盲篇 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任. 身为运维工程师对kubernetes(k8s)可能比较熟,那么etcd(go语言实现)分布式协 ...

  7. 高级Linux运维工程师必备技能(扫盲篇)

    高级Linux运维工程师必备技能(扫盲篇) 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任. 在了解文件系统之前,我们要学习一下磁盘存储数据的方式,大家都知道文件从内存若要持久化存储的 ...

  8. MySQL数据库扫盲篇

    MySQL数据库扫盲篇 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任.   一.MySQL概述 1>.什么是MySQL MySQL是瑞典的MySQL AB公司开发的一个可用于各 ...

  9. 详细介绍android rom移植知识普及

    详细介绍android rom移植知识普及 最近接到很多兄弟们的求助,也回答过无数个和下面这个问题类似的问题: 如何编译android 原生代码得到一个rom,然后跑到某某手机上. 鉴于很多兄弟对这块 ...

随机推荐

  1. ubuntu安装SCrapy

    依次安装 sudo apt-get install build-essential; sudo apt-get install python-dev; sudo apt-get install lib ...

  2. C# ftp 图片上传多快好省

    前言 此篇讲到的是图片上传功能,每个网站必定会有这样类似的功能,上传文件.上传图片等等.那么接下来,看看我们EF+uploadfile+ftp如何玩转上传图片吧 效果预览 具体实现 一个简单数据库 只 ...

  3. 编写高质量代码:改善Java程序的151个建议(第二章:基本类型)

    编写高质量代码:改善Java程序的151个建议(第二章:基本类型) 目录 建议21:用偶判断,不用奇判断 建议22:用整数类型处理货币 建议23:不要让类型默默转换 建议24:边界还是边界 建议25: ...

  4. redis远程连接报错记录

    错误如下 redis可视化工具连接测试 telnet ip 6379 修改关键参数如下 #开通外网访问 # bind 127.0.0.1 #以后台方式运行 daemonize no #取消保护模式,保 ...

  5. Django的rest_framework的序列化组件之序列化多表字段的方法

    首先,因为我们安装了restframework,所以我们需要在django的settings中引入restframework INSTALLED_APPS = [ 'django.contrib.ad ...

  6. 关于界面绘制过程多次回调ondraw()方法产生的问题

    最近项目中,出现一个问题,要做成的效果是这样的,但是一进去就变成这样了, 后来发现,刚进去是正常的,一闪而过,就变成全部了. 界面绘制过程,ondraw() 会被多次回调. 就是说在第一次绘制的时候是 ...

  7. ssh登陆过程图示

  8. Opencv 直方图比较

    #include <iostream>#include <opencv2/opencv.hpp> using namespace std;using namespace cv; ...

  9. 一个文件查看你选择 Run as Android applications 都干了啥

    <?xml version="1.0" encoding="UTF-8"?> <project name="PushFastDemo ...

  10. k8s 部署kube-dns

    [root@k8s-master src]# kubectl create -f kube-dns.yaml service "kube-dns" createdserviceac ...