POJ 1061 青蛙的约会 扩展欧几里德 Java

典型的利用扩展欧几里德算法求解模线性方程!!!
【请点击蓝色字体,查看算法详情】

import java.io.BufferedReader;
import java.io.BufferedWriter;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;
import java.math.BigInteger;
import java.util.Scanner;

/**
 * 题意:中文题,不解析!
 *
 * 分析:典型的模线性方程求解。
 *
 * 解法:①、已知n,m,L,x,y 设步数为 s,由题意得出方程:(ms-ns) % L = y-x
 *           化简得:(m-n)s ≡ (y-x)(mod L)
 *       ②、根据扩展欧几里德算法求解模线性方程。
 *
 * @author TinyDolphin
 */
public class Main {

    private static BigInteger x;
    private static BigInteger y;

    // 扩展欧几里德算法
    public static BigInteger exGcd(BigInteger a, BigInteger b) {
        if (b.compareTo(BigInteger.valueOf(0)) == 0) {
            x = BigInteger.valueOf(1);
            y = BigInteger.valueOf(0);
            return a;
        } else {
            BigInteger d = exGcd(b, a.mod(b));
            BigInteger temp = x;
            x = y;
            y = temp.subtract(a.divide(b).multiply(y));
            return d;
        }
    }

    public static void main(String[] args) {
        Scanner in = new Scanner(new BufferedReader(new InputStreamReader(System.in)));
        PrintWriter out = new PrintWriter(new BufferedWriter(new OutputStreamWriter(System.out)));
        BigInteger X;
        BigInteger Y;
        BigInteger M;
        BigInteger N;
        BigInteger L;
        BigInteger gcd;
        while (in.hasNext()) {
            X = in.nextBigInteger();
            Y = in.nextBigInteger();
            M = in.nextBigInteger();
            N = in.nextBigInteger();
            L = in.nextBigInteger();
            BigInteger b = Y.subtract(X);
            BigInteger a = M.subtract(N);
            BigInteger n = L;
            gcd = exGcd(a, n);
            if (b.mod(gcd).compareTo(BigInteger.valueOf(0)) == 0) {
                // 扩展欧几里德算法求解模线性方程中的 定理一
                x = x.multiply(b.divide(gcd)).mod(n);
                // 扩展欧几里德算法求解模线性方程中的 定理二
                // 注意:得出的 x 可能是负数,需要求最小的正整数解,所以需要 ((x + n/d ) % n/d
                x = x.add(n.divide(gcd)).mod(n.divide(gcd));
                out.println(x);
            } else {
                out.println("Impossible");
            }
        }
        out.flush();
    }
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值