そこで、素数を生成するアトキンアルゴリズムのふるいを作成しました (プロジェクト オイラーの問題用)。という素数 (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;
}
}