0

这是一个获取素数的代码，我已经尽可能地提高了它的效率，但问题是我不能将它转换为 BigInteger，只要不能保存那么多信息；这里的代码：

public class p3{
    static long perfectNumber;
    static long mersenne;

    public static void main(String[] args) {
        long p = 2;
        while (true) {
            if( p % 2 == 0&&p!=2){
                p++;
            }
            else{
                if (isPrime(p) == true) {
                    mersenne = (long) (Math.pow(2, p) - 1);
                    if (isPrime(mersenne) == true) {
                        perfectNumber = (long) Math.pow(2, (p - 1)) * mersenne;
                        System.out.println(perfectNumber);
                    }
                }
                p+=1;
            }
        }
    }
    private static boolean isPrime(long testPrime) {
        for (long i = 3; i < Math.sqrt(testPrime); i += 2) {
            if (testPrime % i == 0) {
                    return false;
            }
        }
        return true;
    }
}
4

1 回答 1

1

我尝试使用 BigInteger 但代码不起作用，因为我不能将 BigInteger 指数与 pow 一起使用

你不需要。指数不需要像梅森素数和完美数一样大。他们可以有自己的独立isPrime()测试。事实上，他们需要满足int，而不是。longBigInteger.pow()

以下是您要求的，但可能不是您想要的。由于时间限制，我怀疑您会在原始代码之外获得一个以上的完美数字，这就是@WJS 将您推向不同方向的原因。

import java.math.BigInteger; 

public class p3 {
    static BigInteger TWO = new BigInteger("2");
    static BigInteger THREE = new BigInteger("3");

    public static void main(String[] args) {
        int p = 2;

        while (true) {
            if (isPrime(p)) {
                BigInteger mersenne = TWO.pow(p).subtract(BigInteger.ONE);

                if (isPrime(mersenne)) {
                    System.out.println(TWO.pow(p - 1).multiply(mersenne));
                }
            }

            p += (p == 2) ? 1 : 2;
        }
    }

    private static boolean isPrime(BigInteger number) {
        if (number.mod(TWO).equals(BigInteger.ZERO)) {
            return number.equals(TWO);
        }

        for (BigInteger i = THREE; number.compareTo(i.multiply(i)) >= 0; i = i.add(TWO)) {
            if (number.mod(i).equals(BigInteger.ZERO)) {
                return false;
            }
        }

        return true;
    }

    private static boolean isPrime(int number) {
        if (number % 2 == 0) {
            return number == 2;
        }

        for (int i = 3; number >= i * i; i += 2) {
            if (number % i == 0) {
                return false;
            }
        }

        return true;
    }
}

输出

> java p3
6
28
496
8128
33550336
8589869056
137438691328
2305843008139952128
2658455991569831744654692615953842176

您的原始代码输出 0 代替上面的最终（37 位）数字。因此，当前的问题确实是long无法容纳足够的信息。但超出这一点，您根本无法使用上述算法进行足够快的计算。

如果我们对上面的代码做一些简单的事情，比如替换这一行：

 if (isPrime(mersenne)) {

和：

if (mersenne.isProbablePrime(10)) {

代码会在缓慢爬行之前吐出前 20 个完美数字。调整你认为合适的确定性论证。isProbablePrime()

于 2019-10-05T04:33:27.370 回答