EXGCD+逆元

求解一组解让ax+by=gcd(a,b);

> 已知当b=0 的时候gcd(a,0)=a;
> 所以这个时候当b=0的时候,存在一组解x=1,y==0
已知a=k*b+r;
k=a/b;
> 那么已知一组解 b*x1+(a-k*b)*y1=gcd(a,b);
那么a*y1+b(x1-y1*k)=gcd(a,b);

同余方程

假如我们要求解
a* x三等b%(p)
那么等同于 (ax-b)%p==0;不妨设为y倍
a
x+py=b;-_-!
就这样好像就可以了。
假如b==1
那么只需要求出一组解
a
x+p*y=1;
同余方程

其实就是让上式中的b=1;

#include<bits/stdc++.h>
using namespace std;
struct dd {
    int gcd,x1,y1;
    dd(int xx=0,int yy=0,int zz=0):gcd(xx),x1(yy),y1(zz) {}
};
dd exgcd(int x,int y) {
    //if(x<y)swap(x,y);
    if(y==0)return dd(x,1,0);//当b==0的时候,一组解是x==1;
    //因为0是所有数的倍数。
    else {
        dd ans=exgcd(y,x%y);
        //cout<<"x="<<y<<"y="<<x%y<<"x1="<<ans.x1<<"y1="<<ans.y1<<endl;
        int k=x/y;
        return dd(ans.gcd,ans.y1,ans.x1-k*ans.y1);
    }
}
int main() {
    int a,b;
    cin>>a>>b;
//a*x+b*p=1;
//假如说我们已知一组解
//那么a*x1+b*y1=1;
//那么已知a*(x-x1)=b(y1-y);
//那么a/b=(y1-y)/(x-x1);
//所以一组解就是x+上b,y-去a 反向
    int x=exgcd(a,b).x1;
    //cout<<"x="<<x<<"y="<<x<<endl;
    if(x<0) {
        x=-(abs(x)%b);
        while(x<0) {
            x+=b;
        }
    }
    cout<<x<<endl;
    return 0;
}

逆元

  • 其实就是上面的b==1的情况,非常好求。
  • 当p=质数的时候。 x的乘法逆元就是pow(x,p-2); 使用快速幂即可
  • 还有一个乘法逆元

拓展中国剩余定理

假如第一组解x1,y1已经知道。
那么
a1x1+b1=a2x2+b2;
那么
那么对于第一组来说,k*b+a是一组通解

他的通解是b1*k+a1

b1*k+a1=k1*b2+a2;
那么b1*x+(-b2)y=(a2-a1);
求出来一组解,y,那么可以知道他的一组通解是
y2+k
b1;
那么可以同时求出第一组解和第二组解。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值