HDU 1576 A/B 扩展欧几里德 Java

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

import java.io.BufferedReader;
import java.io.BufferedWriter;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;
import java.io.StreamTokenizer;

/**
 * 题意:中文题,不解析!!!
 *
 * 扩展欧几里德原理求解模的逆元,当然可以当作模线性方程去求解。
 *
 * 解法:由 (A/B)%9973 = x 得 A = B * x
 *       所以 n = A % 9973 --> n = B*x % 9973 --> B*x ≡ n(mod 9973)
 *       接下来典型的求解模线性方程。
 *
 * @author TinyDolphin
 */
public class Main {

    private static int x;
    private static int y;
    private static int exgcd(int a, int b) {
        if (b == 0) {
            x = 1;
            y = 0;
            return a;
        }
        int d = exgcd(b, a % b);
        int temp = x;
        x = y;
        y = temp - a / b * y;
        return d;
    }

    public static void main(String[] args) throws IOException {
        StreamTokenizer in = new StreamTokenizer(new BufferedReader(new InputStreamReader(System.in)));
        PrintWriter out = new PrintWriter(new BufferedWriter(new OutputStreamWriter(System.out)));
        int T;
        int N;
        int B;
        int a;
        int b;
        int n;
        int d;
        while (in.nextToken() != StreamTokenizer.TT_EOF) {
            T = (int) in.nval;
            while (T-- != 0) {
                in.nextToken();
                N = (int) in.nval;
                in.nextToken();
                B = (int) in.nval;
                a = B;
                b = N;
                n = 9973;
                d = exgcd(a, n);   // 题意已给出,d = gcd(B,9973) = 1
                // 扩展欧几里德算法求解模线性方程中的 定理一
                x = x * (b / d) % n;
                // 扩展欧几里德算法求解模线性方程中的 定理二
                // 注意:得出的 x 可能是负数,需要求最小的正整数解,所以需要 ((x + n/d ) % n/d
                x = (x + n / d) % (n / d);
                out.println(x);
            }
        }
        out.flush();
    }
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值