扩展欧几里得求乘法逆元

好久没发博客了,总结下今晚学的扩展欧几里得吧

虽然现在也没学会

ax + by = c

扩展欧几里得知道a b c 求x y

求乘法逆元

a*x = 1(mod m)

知道a和m 求x

先复习一下欧几里得求最大公约数吧

inr gcd(int x,int y)
{
    if(y==0) return x;
    else return gcd(y,x%y);
}

求 ax + by = c。

int exgcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int  r = exgcd(b,a%b,x,y);
    int t = x;x = y;
    y = t - a/b*y;
    return r;
}

int main()
{
    int a,b,c,x,y;
    while(cin>>a>>b>>c){
        int d = exgcd(a,b,x,y);
        if(c%d)
            printf("-1\n");
        int k = c/d;
        x*=k;
        y*=k;
        printf("%d %d\n",x,y);
    }
    return 0;
}

 求乘法逆元

int exgcd(int a,int b,int &x,int &y)
{
    if(a==0&&b==0) return -1;
    if(b==0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int  r = exgcd(b,a%b,y,x);
    y  -= a/b*x;
    return r;
}
int cal(int a,int m)
{
    int x,y;
    int gcd=exgcd(a,m,x,y);
    if(gcd==1) return (x%m+m)%m;

    return -1;
}
int main()
{
    int a,m;
    while(cin>>a>>m){
    int n = cal(a,m);
    cout<<n<<endl;
    }
    return 0;
}

 

 

 

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值