乘法逆元(模板)

乘法逆元定义:

b存在乘法逆元的充要条件是b与模数m互质 原因:b * x ≡ 1 (mod m) 如果b和m不互质,则 b * x肯定是m的倍数,b * x%m=0

 所以b%m==0 ,b不存在乘法逆元

1.当n为质数时,可以用快速幂求逆元:

  a / b(整除) ≡ a * x (mod m)

  两边同乘b可得 a ≡ a * b * x (mod m) -> 1 ≡ b * x (mod m)

  同 b * x ≡ 1 (mod m)

  由费马小定理可知,当m为质数时:  b ^ (m - 1) ≡ 1 (mod m)

拆一个b出来可得 b * b ^ (m - 2) ≡ 1 (mod m)

     所以当n为质数时,b的乘法逆元 x = b ^ (m - 2)%m

求a的乘法逆元:

import java.util.*;

public class Main{
        static long quick_pow(long a,long b,long c){
                long res=1;
                while(b>0){
                        if((b&1)==1) res=res*a%c;
                        a=a*a%c;
                        b>>=1;
                }
                return res;
        }
        public static void main(String[] args) {
                Scanner scan=new Scanner(System.in);
                int t=scan.nextInt();
                while(t-->0){
                        long a=scan.nextLong();
                        long p=scan.nextLong();
                        if(a%p!=0) System.out.println(quick_pow(a,p-2,p));
                        else System.out.println("impossible");
                }
        }
}

原文地址:https://www.cnblogs.com/qdu-lkc/p/12262827.html

时间: 2024-10-07 05:04:31

乘法逆元(模板)的相关文章

乘法逆元模板

1.扩展欧几里得求逆元 typedef long long ll; //ax + by = gcd(a,b) //传入固定值a,b.放回 d=gcd(a,b), x , y void extendgcd(ll a,ll b,ll &d,ll &x,ll &y) { if(b==0){d=a;x=1;y=0;return;} extendgcd(b,a%b,d,y,x); y-=x*(a/b); } //Ax=1(mod M),gcd(A,M)==1 //输入:10^18>=A

求乘法逆元模板(扩展欧几里得)

void exgcb(LL a,LL b,LL &d,LL &x,LL &y){ if(!b){d=a;x=1;y=0;return;} exgcb(b,a%b,d,y,x);y-=x*(a/b); } LL ny(LL a,LL b){ ///求a关于b的逆元(要求a,b互质) LL d,x,y; exgcb(a,b,d,x,y); return d==1?(x+b)%b:-1; }

【codevs 1200】【NOIP 2012】同余方程 拓展欧几里德求乘法逆元模板题

模板,,, #include<cstdio> using namespace std; void exgcd(long long a,long long b,long long &x,long long &y){ if (b==0) {x=1; y=0;} else {exgcd(b,a%b,x,y); int t=y; y=x-a/b*y; x=t;} } int main(){ long long a,b,x,y; scanf("%lld %lld\n"

P3811 【模板】乘法逆元

P3811 [模板]乘法逆元 题目背景 这是一道模板题 题目描述 给定n,p求1~n中所有整数在模p意义下的乘法逆元. 输入输出格式 输入格式: 一行n,p 输出格式: n行,第i行表示i在模p意义下的逆元. 输入输出样例 输入样例#1: 10 13 输出样例#1: 1 7 9 10 8 11 2 5 3 4 说明 1≤n≤3×10?6??,n<p<20000528 输入保证 p 为质数. 我们有三种办法求逆元 由欧拉定理可知 当gcd(a,n)==1 时 我们有 Aφ(n-1)≡ 1(mod

洛谷 P3811 【模板】乘法逆元 如题

P3811 [模板]乘法逆元 时空限制1s / 256MB 题目背景 这是一道模板题 题目描述 给定n,p求1~n中所有整数在模p意义下的乘法逆元. 输入输出格式 输入格式: 一行n,p 输出格式: n行,第i行表示i在模p意义下的逆元. 输入输出样例 输入样例#1: 10 13 输出样例#1: 1 7 9 10 8 11 2 5 3 4 说明 1 \leq n \leq 3 \times 10 ^ 6, n < p < 200005281≤n≤3×106,n<p<20000528

题解 P3811 【模板】乘法逆元

题意求\(i\)在模\(p\)意义下的逆元\(\frac{1}{i}\)即\(inv(i)\).题目数据范围很明显规定了要求一个线性求逆元的算法. 令\(p=ai+b\),则有: \[ai+b\equiv 0(\mod p)\] 移项得: \[ai\equiv -b(\mod p)\] 系数化简得: \[i\equiv -\frac{b}{a}(\mod p)\] 取倒数得: \[\frac{1}{i} \equiv -\frac{a}{b}(\mod p)\] 即 \[inv(i) \equi

[洛谷P3811]【模板】乘法逆元

题目大意:给你n和质数p,求1~n在模p意义下的乘法逆元(n<p). 解题思路:由于$n<p<20000528$,所以扩展欧几里得是会超时的.这儿就要用到线性推逆元大法辣!→不懂戳这里← 注意乘法可能会超过int,所以计算时先转化为long long即可. C++ Code: #include<cstdio> int n,p,inv[20000529]; int main(){ scanf("%d%d",&n,&p); puts("

【模板】求1~n的整数的乘法逆元

洛谷3811 先用n!p-2求出n!的乘法逆元 因为有(i-1)!-1=i!-1*i (mod p),于是我们可以O(n)求出i!-1 再用i!-1*(i-1)!=i-1 (mod p)即是答案 1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cstdio> 5 #include<cmath> 6 #include<algorithm> 7

[模板]乘法逆元

本博客所有代码基于题目 luogu_P3811 逆元: 一般用于求 (a/b) mod p  定义: 若 a*x ≡ 1 (mod p) ,且 a 与 p 互质,那么我们就能定义: x 为 a 的逆元,记为 a^-1 ,所以我们也可以称 x 为 a 的倒数(mod p意义下). 所以对于 (a/b) mod p ,我们就可以求出 b 在 mod p 意义下的逆元,然后乘上 a ,再 mod p ,就是这个乘法逆元的值了. 求法: First:费马小定理 定理内容:如果 a , p 互质,那么 a

Luogu P3811 [模板]乘法逆元 题解报告

题目传送门 [题目大意] 给定$n$,求$1-n$在膜$p$意义下的乘法逆元. [思路分析] 好的原本我只会求单个数的逆元,然后被告知了这道题之后发现自己不会做(我果然还是太弱了),于是就学了一下递推求逆元. 设$p=k*i+r$,则可得$k*i+r\equiv0(mod\ p)$,然后乘上$i^{-1},r^{-1}$即可得到$k*r^{-1}+i^{-1}\equiv0(mod\ p)$ 由于$k=\lfloor \frac{p}{i}\rfloor,r=p\ mod\ i$,所以$i^{-