Java Fermat因子分解算法BigInteger不起作用

Question

我正在使用BigInteger实现Fermat分解算法，因此可以分解。 但是目前，该代码无法正常工作。 它由于某种原因挂起。 有人可以将我定向到问题所在，还是让我知道我的算法是否错误？ BigInteger使生活变得困难，因此我不得不寻找平方根方法。

import java.math.BigInteger;
import java.util.Scanner;

public class Fermat
{
    /** Fermat factor **/
    public void FermatFactor(BigInteger N)
    {
        BigInteger a = sqrt(N);
        BigInteger b2 = a.multiply(a).subtract(N);

        while (!isSquare(b2)) {
            a = a.add(a);
            b2 = a.multiply(a).subtract(N);
        }

        BigInteger r1 = a.subtract(sqrt(b2));
        BigInteger r2 = N.divide(r1);
        display(r1, r2);
    }

    /** function to display roots **/
    public void display(BigInteger r1, BigInteger r2) {
        System.out.println("\nRoots = "+ r1 +" , "+ r2);    
    }

    /** function to check if N is a perfect square or not **/
    public boolean isSquare(BigInteger N) {
        BigInteger ONE = new BigInteger("1");
        BigInteger sqr = sqrt(N);

        if (sqr.multiply(sqr) == N  || (sqr.add(ONE)).multiply(sqr.add(ONE)) == N)
            return true;
        return false;
    }


    public static BigInteger sqrt(BigInteger x)
            throws IllegalArgumentException {
        if (x.compareTo(BigInteger.ZERO) < 0) {
            throw new IllegalArgumentException("Negative argument.");
        }
        // square roots of 0 and 1 are trivial and
        // y == 0 will cause a divide-by-zero exception
        if (x == BigInteger.ZERO || x == BigInteger.ONE) {
            return x;
        } // end if
        BigInteger two = BigInteger.valueOf(2L);
        BigInteger y;
        // starting with y = x / 2 avoids magnitude issues with x squared
        for (y = x.divide(two);
                y.compareTo(x.divide(y)) > 0;
                y = ((x.divide(y)).add(y)).divide(two));
        if (x.compareTo(y.multiply(y)) == 0) {
            return y;
        } else {
            return y.add(BigInteger.ONE);
        }
    } // end bigIntSqRootCeil


    /** main method **/
    public static void main(String[] args) 
    {
        Scanner scan = new Scanner(System.in);
        System.out.println("Fermat Factorization Test\n");
        System.out.println("Enter odd number");
        BigInteger N = scan.nextBigInteger();
        Fermat ff = new Fermat();
        ff.FermatFactor(N);
        scan.close();
    }
}

我知道我有很多错误，但是可以提供任何帮助。 谢谢。

Answer 1

您的“ for”循环：

for (y = x.divide(two);
    y.compareTo(x.divide(y)) > 0;
    y = ((x.divide(y)).add(y)).divide(two));

不终止。 也许您可以跟踪变量“ y”的值，以猜测何时必须停止。

编辑：那是错误的（请参阅评论）。 问题出在行中

a = a.add(a)

内部程序FermatFactor。 应该是

a = a.add(ONE)

在我的机器上，我还遇到了使用'A == B'测试相等性的麻烦。 方法“ A.equals（B）”修复了该问题。