欧几里得算法求解x,y。和扩展欧几里得

/*
裴蜀定理:若ax+by = z,则 gcd(a,b)| z

如何求解 (以下讨论a>b)
显然当 b=0,gcd(a,b)=a。此时 x=1,y=0;
当a>b>0 时
设 ax1+ by1= gcd(a,b);
bx2+ (a mod b)y2= gcd(b,a mod b);
根据欧几里德原理有 gcd(a,b) = gcd(b,a mod b);
则:ax1+ by1= bx2+ (a mod b)y2;
即:ax1+ by1= bx2+ (a - [a / b] * b)y2 = ay2+ bx2- [a / b] * by2;(a mod b = a - [a / b]*b;[a / b]代表a整除b)
也就是ax1+ by1 = ay2 + b(x2- [a / b] *y2);
根据恒等定理得:x1=y2;y1=x2- [a / b] *y2;
这样我们就得到了求解 x1,y1 的方法:x1,y1 的值基于 x2,y2
由引理我们知道:ax+by = z,z为gcd(a,b)若干倍,所以我们先求解ax+by = gcd(a,b),再将求出的解乘以 z/gcd(a,b)就好了。
 */


public class 欧几里得算法 {

    static long x;
    static long y;
    public static void main(String[] args) {

        printxy(2,3,1);
        System.out.println(x);   //定义的x和y
        System.out.println(y);


    }

    //辗转相除法求最大公约数,就是欧几里得算法
   public int gcd(int a,int b){
        if(b==0) return a;
        else return gcd(b,a%b);
    }

    //扩展欧几里得算法的功能就更强大了,它可以用来求二元一次方程的通解.
    public static long exGcd(int a,int b)//网上有很多&x,&y 的方法,在Java中没有,只能定义全局变量
    {
        if(b==0)
        {
            x = 1;
            y = 0;
            return a;
        }
        long r = exGcd(b,a%b);
        long t;
        t = (int) x;
        x = y;            //这一步其实很好理解,就是a换成b,b换成a%b
        y = t-a/b*y;      //这一步是a%b =  a- (a/b)*b
        return r;
    }

    //求求二元一次方程的通解
    public static void printxy(int a,int b,int m){
        long l = exGcd(a, b);
        long beishu = m/l;    //获取倍数
        if(m%l != 0) System.out.println("无解");
        x = x * l;  //通解x
        y = y * l;  //通解y
        System.out.println(l);
    }

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值