乘法逆元--模板

最新推荐文章于 2024-07-24 19:35:30 发布

一只蓝色小鲨鱼

最新推荐文章于 2024-07-24 19:35:30 发布

阅读量214

点赞数 1

分类专栏：总结的语法文章标签：算法 c++ 图论

本文链接：https://blog.csdn.net/m0_64267361/article/details/136314804

版权

总结的语法专栏收录该内容

23 篇文章 7 订阅

订阅专栏

本文介绍了如何使用快速幂算法和扩展欧几里得方法计算大整数a在模p下的逆元，以及逆元在解决a/bmodp问题中的作用，特别提到当p为质数时，逆元为a^(p-2)，非质数情况下的处理方式。

摘要由CSDN通过智能技术生成

876. 快速幂求逆元

//b*x≡1(modp)；如果是逆元；
//如果 b和p互质，那么有b^(p-1)≡(1mod p);
//所以 b*b^(p-2)≡1(mod p);
//那么b的逆元就是b^(p-2);
//p为质数，可知a的p-2次方为a的逆元；

2.p当不是质数时,则需要用扩展欧几里得求逆元
3.逆元的作用
当 a,b,很大时 求 a/b mod p 的值
而 a/b mod p≠((a mod p)/(b mod p))mod p
因此可以借助逆元转化为乘法,再算,设b的逆元为b−1
则  a/b mod p= a*b^-1 mod p

#include <iostream>
#include <algorithm>
using namespace std;

typedef long long LL;

int qmi(int a,int k,int p)
{ 
    LL res=1;
    while(k)
    {
        if(k&1) res=(LL)res*a%p;
        k>>=1;
        a=(LL)a*a%p;
    }
    return res;
}

int main()
{
    int n;
    cin>>n;
    while(n--)
    {
        int a,p;
        scanf("%d%d",&a,&p);
        if(a%p==0)printf("impossible");//特判一下；a是p的倍数。余数为0；
        else 
        printf("%d",qmi(a,p-2,p));a不是p的倍数；则有费马小定理；
        puts("");
    }
    return 0;
}