所以我做了一个阿特金算法的筛子来生成素数(针对欧拉计划问题)。它返回一个名为 的素数 (int[]) 数组primes
。问题是,每当我尝试从某个索引访问该数组时(我在我的main
方法中这样做),它都会抛出一个java.lang.ArrayIndexOutOfBounds
异常。请帮忙(如果您认为我的逻辑与仙女不符,并且有更简单的方法可以做到这一点,请告诉我)!
import java.lang.*;
import java.util.*;
public class SieveOfAtkin {
public static void main(String[] stuff) {
int limit = getInt("Limit?");
int[] p = getPrimes(limit);
for(int i = 0; i < p.length; i++) {
System.out.println(p[i]);
}
}
public static int[] getPrimes(int limit) {
boolean[] isPrime = new boolean[limit + 1];
double sqrt = Math.sqrt(limit);
int n = 0;
for(int i = 0; i < limit + 1; i++) {
isPrime[i] = false;
}
for(int x = 0; x < sqrt; x++) {
for(int y = 0; y < sqrt; y++) {
n = 4 * x * x + y * y;
if(n <= limit && (n % 12 == 1 || n % 12 == 5)) {
isPrime[n] = !isPrime[n];
}
n = 3 * x * x + y * y;
if(n <= limit && n % 12 == 7) {
isPrime[n] = !isPrime[n];
}
n = 3 * x * x - y * y;
if(n <= limit && n % 12 == 11 && x > y) {
isPrime[n] = !isPrime[n];
}
}
}
for(int i = 5; i < sqrt; i++) {
if(isPrime[i]) {
for(int j = i * i; j < limit; j = j + (i * i)) {
isPrime[j] = false;
}
}
}
int count = 0;
for(int i = 0; i < isPrime.length; i++) {
if(isPrime[i]) {
count++;
}
}
int[] primes = new int[count];
int found = 0;
if (limit > 2) {
primes[found++] = 2;
}
if (limit > 3) {
primes[found++] = 3;
}
for (int p = 5; p <= limit; p += 2) {
if (isPrime[p]) {
primes[found++] = p;
}
}
return primes;
}
public static int getInt(String prompt) {
System.out.print(prompt + " ");
int integer = input.nextInt();
input.nextLine();
return integer;
}
}