乘法逆元笔记（等以后遇到例题再加，现在就6道^_^）

已于 2022-05-20 10:43:19 修改
阅读量2k 收藏 2
点赞数 2
分类专栏： acm 文章标签： c++
于 2022-04-03 08:54:33 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_64839851/article/details/123932297
版权

拓展欧几里得求逆元

求从1-n中mod p的逆元

inline ll exgcd(ll a,ll b,ll &x,ll &y){
  return !b?(x=1,y=0,a):(exgcd(b,a%b,y,x),y-= (a/b)*x);
}

int main(){
    ll n,p;
    cin>>n>>p;
    for(ll i=1;i<=n;i++){
        ll j,k;
        //ij≡1(mod p)  ij+p(-k)=1
        ll d=exgcd(i,p,j,k);
        k=-k;
        j=(j%p+p)%p;
        cout<<j<<endl;
    }
}

快速幂求逆元

ij≡1（mod p）由费马小定理知ij≡i^(p-1)（mod p)，即j≡i^(p-2)（mod p)

然后就开始快乐的快速幂把

！！！！！快速幂只限于p是素数的情况！！！！！

ll qpow(ll x,ll n,ll p){//x^n mod p
    x%=p;
    ll ans=1;
    while(n){
        if(n&1)ans=qmul(x,ans,p);
        n/=2;
        x=qmul(x,x,p);
    }
    return ans;
}

int main(){
    ll n,p;
    cin>>n>>p;
    for(ll i=1;i<=n;i++){
        ll j;
        j=qpow(i,p-2,p);
        cout<<j<<endl;
    }
}

上述两种方法求线性逆元都较慢，下给出线性逆元求法

这里记对于每个i的乘法逆元为i^{-1}（事实上也就是这么写的）,所以我们只需要求i^{-1}.

首先我们有一个,1^{-1}=1 ( modp )，然后设 p=k*i+r,(1<r<i<p)

则k*i+r≡0（mod p）

然后乘上i^{-1},r^{-1}就可以得到:

k*r^{-1}+i^{-1}≡0(mod p)

i^{-1}≡−k∗r^{-1}(modp),即i^{-1}≡−⌊p/i​⌋∗(p%i)^{-1}(modp)

于是，我们就可以从前面推出当前的逆元了。

1.递推公式为inv[i]=-(p/i)*inv[p%i]

(这里注意到前面一项为负数，进行一下体调整)

2.inv[i]=(p-p/i)*inv[p%i]

(这里加上了一个p*inv[p%i]，但由于是mod p不会影响)

3.inv[i]=(p-p/i)*inv[p%i]%p   //最终递推公式

根据题目的i的规模，可以采用记忆化递推进行优化

    ll p,n;//multiplicative inverse
    cin>>n>>p;
    inv[1]=1;inv[0]=0;
    for(ll i=2;i<=n;i++){
        inv[i]=inv[p%i]*(p-p/i)%p;
    }

ps:当 inv[i]=0时，表示不存在乘法逆元

下面是一些例题和题解

 P1082 [NOIP2012 提高组] 同余方程

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define LOCAL

inline ll exgcd(ll a,ll b,ll &x,ll &y){
  return !b?(x=1,y=0,a):(exgcd(b,a%b,y,x),y-= (a/b)*x);
}

int main(){
    ll a,b,x,k;
    cin>>a>>b;
    //ax≡1(mod b)  ax+bk=1
    //ax≡a^b-1(mod b) x≡a^b-2(mod b)
    ll d=exgcd(a,b,x,k);
    x=(x%b+b)%b;
    cout<<x;
}

P3811 【模板】乘法逆元（线性求逆元）

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll inv[3000010];
int main(){
    ll p,n;
    scanf("%lld%lld",&n,&p);
    inv[1]=1;inv[0]=0;
    for(ll i=2;i<=n;i++) inv[i]=inv[p%i]*(p-p/i)%p;
    for(ll i=1;i<=n;i++) printf("%lld\n",inv[i]);
}

第一次用cin,cout被卡了 。。。。。。 

P2054 [AHOI2005]洗牌 

根据题目意思有下图（李姐一下，蹩脚的画图技术）

 1.首先令n=2k，原位置L，洗牌后的位置L'；

2.1当L>k时 L'=(L-k)*2-1=2L-n-1

2.2当1=<L<=k时 L'=2L

3.则L'≡2L（mod n+1）

 1.设初始位置为P，末位置为L（已知）也就是P*2^{M}≡L（mod n+1）

2.P≡L*(2^{M})^{-1}（mod n+1），下面只用求2^{M}的乘法逆元即可

3.记2^{M}≡r（mod n+1）(这里使用快速幂求r)

4.然后就是最简单的乘法逆元了

 又是被文件重复定向坑的一天

用longlong洛谷只有84，最后一组数据会爆longlong 

这里用了__int128，注意输入输出要用快读快写

上述ll代表__int128，因为之前是用longlong 错了，懒得改了

 P4071 [SDOI2016]排列计数

题目：求有多少种 1 到 n的排列 a，满足序列恰好有 m 个位置，使得 a_i =i。答案对 10^9 + 7取模。

注意1e9+7是十位数中第一个素数（所以可以快速幂了）

考虑恰有m个位置a_i !=i，那么对于j+1个位置不同分3类，这里从n个数中选了m个位置不同，则共有C_{n}^{m}种取法，记m个位置均不满足a_i =i的取法有A_{m} 种

1.前j个元素均满足题设，将j+1与前j个数中的任意一个交换得到j+1个位置不同的序列，也就是j*A_{j}种

 

2. 前j个元素中仅j-1个满足题设,j*A_{j-1}种

 

 3. 前j个元素中小于j-1个满足题设，那么新引入的j+1，与前面不满足题设的位置交换，肯定至少还剩1个位置不满足，也就是这样的排列是不行的

综上A_{j+1}=j*(A_{j}+A_{j-1})，其中A_{1}=0，A_{2}=1

 超时代码！！！！！！洛谷60分

#include<bits/stdc++.h>
using namespace std;
#define ll long long

inline ll read(){
    char ch=getchar();ll s=0,w=1;
    while(ch<48||ch>57)
        if(ch=='-')w=-1;ch=getchar();
    while(ch>=48&&ch<=57)
        s=(s<<1)+(s<<3)+ch-48;ch=getchar();
    return s*w;
}

inline ll qmul(ll a,ll b,ll mod){//a*b mod mod
    ll ans=0;
    while(b){
        if (b&1)ans=(ans+a)%mod;
        b>>=1;
        a=a<<1%mod;
    }
    return ans;
}

inline void write(ll x){
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}

const ll p=1e9+7;
ll inv[1000005],A[1000005],step[1000005];

int main(){
    ll T;
    T=read();
    inv[1]=1,inv[0]=1,A[0]=1,A[2]=1,step[1]=1;

    //A[i]=(i-1)(A[i-1]+A[i-2]) i>=3;
    for(ll i=3;i<=1e6+5;i++)
        A[i]=qmul(i-1,A[i-1]+A[i-2],p);

    //求组合数的模
    //注意到组合数中有除法，由于除法的模时候没有意义的，所以这时候就要求逆元
    //所以只需求1到m的乘法逆元，由于是线性的，用线性求逆元的方法
    for(ll i=2;i<=1e6+5;i++)
        inv[i]=inv[p%i]*(p-p/i)%p;
    
    //这里求了inv[i]!的模
    for(ll i=2;i<=1e6+5;i++)
        inv[i]=inv[i]*inv[i-1]%p;
        
    //i!的模
    for(ll i=2;i<=1e6+5;i++)
        step[i]=step[i-1]*i%p;

    while(T--){
        ll n,m,c_mod=1;
        n=read(),m=read();
        c_mod=qmul(inv[n-m],inv[m],p);
        c_mod=qmul(c_mod,step[n],p);
        write(qmul(c_mod,A[n-m],p));putchar('\n');
    }
}

然后我把qmul删了 ，因为这题好像不会爆long long ，然后就快乐ac了哈哈哈

#include<bits/stdc++.h>
using namespace std;
#define ll long long
 
inline ll read(){
    char ch=getchar();ll s=0,w=1;
    while(ch<48||ch>57)
        if(ch=='-')w=-1;ch=getchar();
    while(ch>=48&&ch<=57)
        s=(s<<1)+(s<<3)+ch-48;ch=getchar();
    return s*w;
}
 
inline void write(ll x){
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
 
const ll p=1e9+7;
ll inv[1000005],A[1000005],step[1000005];
int main(){
    ll T;
    T=read();
    inv[1]=1,inv[0]=1,A[0]=1,A[2]=1,step[1]=1;
 
    //A[i]=(i-1)(A[i-1]+A[i-2]) i>=3;
    for(ll i=3;i<=1e6+5;i++)
        A[i]=(i-1)*(A[i-1]+A[i-2])%p;
 
    //求组合数的模
    //注意到组合数中有除法，由于除法的模时候没有意义的，所以这时候就要求逆元
    //所以只需求1到m的乘法逆元，由于是线性的，用线性求逆元的方法
    for(ll i=2;i<=1e6+5;i++)
        inv[i]=inv[p%i]*(p-p/i)%p;
    
    //这里求了inv[i]!的模
    for(ll i=2;i<=1e6+5;i++)
        inv[i]=inv[i]*inv[i-1]%p;
        
    //i!的模
    for(ll i=2;i<=1e6+5;i++)
        step[i]=step[i-1]*i%p;
 
    while(T--){
        ll n,m,c_mod=1;
        n=read(),m=read();
        c_mod=inv[n-m]*inv[m]%p;
        c_mod=c_mod*step[n]%p;
        write(c_mod*A[n-m]%p);putchar('\n');
    }
}

P4942 小凯的数字

首先一个数mod 9的余数等于它各个位上的和mod 9的余数，也就是这题求(l+r)(l-r+1)/2 mod 9（等差数列求和）注意这里由有2,运用逆元，可以将其变为乘5.

#include<bits/stdc++.h>
using namespace std;
#define ll long long

inline ll read(){
    char ch=getchar();ll s=0,w=1;
    while(ch<48||ch>57){
        if(ch=='-')w=-1;ch=getchar();
    }
    while(ch>=48&&ch<=57){
        s=(s<<1)+(s<<3)+ch-48;ch=getchar();
    }
    return s*w;
}

inline void write(ll x){
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}

int main(){
    ll t;
    t=read();
    while(t--){
        ll l,r,ans;
        l=read(),r=read();
        if(l>r)swap(l,r);
        ans=(l+r)%9*(r-l+1)%9*5%9;
        write(ans);putchar('\n');
    }
}

P2613 【模板】有理数取余

由于19260817是质数，所以不会angry  哈哈哈哈哈哈 

然后求乘法逆元直接上小费马就可以了

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll ans,a,b,p=19260817;
string s1,s2;
ll qpow(ll x,ll n){
    ll res=1;
    while(n){
        if(n&1)res=res*x%p;
        n>>=1;
        x=x*x%p;
    }
    return res;
}
int main(){
    cin>>s1>>s2;
    for(int i=0;i<s1.size();i++){
        a=a*10+s1[i]-'0';
        a%=p;
    }
    for(int i=0;i<s2.size();i++){
        b=b*10+s2[i]-'0';
        b%=p;
    }
    cout<<a*qpow(b,p-2)%p;
    return 0;
}

P5431 【模板】乘法逆元 2 题解

这道题是真的毒瘤，卡常凡死了；

首先将所求通分 \sum \frac{k^{i}\prod_{j!=i}^{} a_{j}}{\prod a_{i}} ，这样求就只用求一次逆元了

然后预处理a[i]的前缀积后缀积，O（1）时间求出分子中的连乘

代码如下：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=5*1e6;
ll n,k1=1,p,k,ans,s1[N+10],s2[N+10],a[N+10];
inline ll qpow(ll x,ll n){
    ll res=1;
    while(n){
        if(n&1)res=res*x%p;
        n>>=1;
        x=x*x%p;
    }
    return res;
}

inline int read(){
    char ch=getchar();int s=0,w=1;
    while(ch<48||ch>57)
        if(ch=='-')w=-1;ch=getchar();
    while(ch>=48&&ch<=57)
        s=(s<<1)+(s<<3)+ch-48;ch=getchar();
    return s*w;
}

void inti(){
    n=read(); p=read(); k=read();
    s1[0]=s2[n+1]=1;
    for(int i=1;i<=n;i++) a[i]=read();
    for(int i=1;i<=n;i++) s1[i]=s1[i-1]*a[i]%p;
    for(int i=n;i>=1;i--) s2[i]=s2[i+1]*a[i]%p;
}

int main(){
    inti();
    ll inv=qpow(s1[n],p-2);
    for(int i=1;i<=n;i++){
        k1=k1*k%p;
        ans=ans+s1[i-1]*s2[i+1]%p*k1%p;
        ans%=p;
    }
    cout<<ans*inv%p;
    return 0;
}

Cot287
关注
  • 2
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
ctf.show与学习笔记
Q221022的博客
01-29 3123
ctf.show 1.密码学签到 看到密文之后不难发现，flag就是倒着写的，直接将密文倒叙输出就是答案 2. 看到题目，全是由符号组成的，将题目内容全部复制到控制台，回车即可得到答案 3. 第三题与第二题看起来是一样的，所以先将其复制到控制台，发现报错了，有中文的引号，将所有的引号全部换成英文之后还是没有得到flag，········（不会了，没有思路了） 4. 直接用工具计算出私钥d即可 5. 根据rsa密码的解密规则解出m 6. 没有判断出...
读书笔记：从《Hotspot实战》一书中引出的想法^_^.zip
最新发布
07-25
读书笔记：从《Hotspot实战》一书中引出的想法^_^
乘法逆元
weixin_30765319的博客
08-21 102
重新整理了一下笔记 引入 一个数的逆元就是数论中，该数在模另一个数的意义下的倒数。 这么说或许太抽想了，我们可以看这个： 求组合数 $\Large{\left( ^a_b \right)}$ 取模 $10^9+7$ 的值。 这道题看上去简单，但是我们发现： 组合数是有乘有除的，而我们一般都是先乘完所有的再一个个除，这时就会出事了： 就比如 $\Large{\left( ...
CTF学习笔记一——RSA加密
qq_45198339的博客
01-20 6380
本文对ctf中的rsa题型进行了归纳总结
离散数学_笔记
weixin_42892844的博客
07-28 1785
内容大纲 ╭1> 命题与命题公式 |2> 命题逻辑的推理理论 |3> 谓词逻辑 |4> 集合 |5> 关系与函数 离散数学 < |6> 代数系统的一般概念 |7> 格与布尔代数 |8> 图 ╰9> 图的应用 第一章:命题与命题公式 ╭1.命题与命题连接词...
密码学-古典密码学习笔记
baola的博客
08-23 1742
本文作为学习笔记来使用，对正在上密码学课程的学生来说例题比较有价值，此外在图片和文字上大量参考了各种密码学书籍，并添加了自己的注解和理解；
2021涅普冬令营Misc笔记与题解
break_cat的博客
02-10 3283
前言： 在我进行CTF签到活动的不长时间里，也没有专门学过，平常也就应付应付签到题，正好最近有个挺好的入门课，结果有事绊着，到现在录播就看了三个，这两天趁没事了赶紧补一补，先前零碎而间断地写练习场的题，写到今天有点写不下去了，于是就先小结一下，有些题不会，求指点。 笔记大多是从幻灯片摘下来的，自己稍微又加了点儿。 题目来自练习场https://camp.hackingfor.fun/ 所有图片均为题原图。 题目前面有问号的为存在疑问或没能做出来的题，求指点。 Misc Misc简介 Misc为英文mi.
笔记
xgx984826498的博客
02-21 157
对sort 函数制定排序规则 例题 给出n个区间，判断时候有哪个区间是被包含在哪个区间的 ，如果有多个输出其中的一组。先输出被包含的区间（他的次叙数），再输出包含他的区间的次序输，如果没有输出的就输出-1 -1 思路 先按左边界从小到大排序，左边界相同的就再按右边界从大到小排序，如果再相同就按序号从小到大排序，可以给sort函数制定一种排序规则，这题我觉得我得背下来，既有结构体还有对sort...
基础数论知识复习笔记
qq_49593247的博客
03-13 3655
目录 欧拉函数 快速幂求逆元 拓展欧几里得算法 欧拉函数 在忘掉了很多学过的东西之后,我开始了漫长的复习之旅; 首先是欧拉函数,要搞清楚欧拉函数,要理解欧拉函数的含义.欧拉函数表示的是在1到n之中和n互质的数的个数.其通式如下:. 即phi[ n ] = n * ( 1 - ( 1 / p1 ) )...i * ( 1 - ( 1 / p1 ) );phi为欧拉函数,p为该数字i的每一个质因数. 欧拉函数性质,欧拉函数后面乘的( 1 - ( 1 / p1 ) )只和欧拉函数的质因.
离散数学 笔记整理
NorthStar的博客
06-22 3914
离散数学 内容大纲 ╭1> 命题与命题公式 |2> 命题逻辑的推理理论 |3> 谓词逻辑 |4> 集合 |5> 关系与函数 离散数学 < |6> 代数系统的一般概念 |7> 格与布尔代数 |8> 图 ╰9> 图的应用 第一章:命题与命题公式 ╭1.命题与命题连接词:
StudyWorkspace:这是一个iOS学习的Demo集, 会把经常用的一些小笔记记录在此 ^_^
04-08
"经常用的一些小笔记记录在此 ^_^" 暗示这个Demo集不仅包含了代码实例，还可能包含了作者在学习过程中积累的技巧、心得或提示，这对于初学者来说是非常宝贵的资源。 【描述分析】 描述部分 "#这是一个iOS学习的Demo...
Android应用开发笔记13^Android移动应用界面的模板化设计.pdf
11-18
Android应用开发笔记13^Android移动应用界面的模板化设计.pdf
要点、考点与例题——第6章 一阶电路-教程与笔记习题
05-19
要点、考点与例题——第6章 一阶电路
YUNBIJI.zip_spring ajax_云笔记_云笔记 项目_云笔记项目
09-20
云笔记项目主要用于客户在线学习时记录笔记，在该功能基础上又扩展了分享，收藏，活动等功能，该项目采用MVC设计模式，Spring+MyBatis，Ajax,jQuery框架开发，在此技术上实现项目的登录模块，笔记本模块，笔记模块，...
【C++】C++中的find方法介绍
qq_35902025的博客
07-25 395
find 函数是 std::string 类中非常有用的成员函数之一，它提供了灵活的方式来查找子字符串或字符。使用 npos 可以检查查找操作是否成功。
《Cross-Modal Dynamic Transfer Learning for Multimodal Emotion Recognition》
buyaotutou的博客
07-25 1134
例如：第一章 Python 机器学习入门之pandas的使用文章链接Multimodal Emotion Recognition is an important research area for developing human-centric applications, especially in the context of video platforms. Most existing models have attempted to develop sophisticated fusion tech
背单词工具(C++)
zeyeqianli的博客
07-24 3343
【代码】背单词工具(C++)
C++ //练习 15.28 定义一个存放Quote对象的vector，将Bulk_quote对象传入其中。计算vector中所有元素总的net_price。
Made in Code
07-25 613
C++ Primer（第5版） 练习 15.28
"动力节点老杜Spring6配套笔记PDF文件分享
动力节点老杜Spring6配套的笔记PDF文件包括了对Spring6的基本概念、核心功能、使用方法、最佳实践等方面的详细介绍和讲解。这些笔记对于想要学习和掌握Spring6框架的开发者来说是非常有帮助的。通过阅读这些笔记，...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值