2022-07-24 23:03:35 +08:00
|
|
|
//! This script uses the Sieve of Eratosthenes to calculate prime numbers.
|
2020-03-12 23:46:40 +08:00
|
|
|
|
2020-04-11 16:06:57 +08:00
|
|
|
let now = timestamp();
|
|
|
|
|
2020-05-31 22:47:39 +08:00
|
|
|
const MAX_NUMBER_TO_CHECK = 1_000_000; // 9592 primes <= 100000
|
2020-03-12 23:46:40 +08:00
|
|
|
|
|
|
|
let prime_mask = [];
|
2021-12-15 12:06:17 +08:00
|
|
|
prime_mask.pad(MAX_NUMBER_TO_CHECK + 1, true);
|
2020-03-12 23:46:40 +08:00
|
|
|
|
2020-04-22 19:28:08 +08:00
|
|
|
prime_mask[0] = false;
|
|
|
|
prime_mask[1] = false;
|
2020-03-12 23:46:40 +08:00
|
|
|
|
|
|
|
let total_primes_found = 0;
|
|
|
|
|
2021-12-15 12:06:17 +08:00
|
|
|
for p in 2..=MAX_NUMBER_TO_CHECK {
|
2020-05-03 22:12:10 +08:00
|
|
|
if !prime_mask[p] { continue; }
|
2020-03-12 23:46:40 +08:00
|
|
|
|
2020-05-31 15:57:53 +08:00
|
|
|
//print(p);
|
2020-03-12 23:46:40 +08:00
|
|
|
|
2020-05-03 22:12:10 +08:00
|
|
|
total_primes_found += 1;
|
|
|
|
|
2021-12-15 12:06:17 +08:00
|
|
|
for i in range(2 * p, MAX_NUMBER_TO_CHECK + 1, p) {
|
2020-05-03 22:12:10 +08:00
|
|
|
prime_mask[i] = false;
|
2020-03-12 23:46:40 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2021-04-04 23:22:45 +08:00
|
|
|
print(`Total ${total_primes_found} primes <= ${MAX_NUMBER_TO_CHECK}`);
|
|
|
|
print(`Run time = ${now.elapsed} seconds.`);
|
2020-05-13 11:56:48 +08:00
|
|
|
|
2020-05-31 22:47:39 +08:00
|
|
|
if total_primes_found != 78_498 {
|
|
|
|
print("The answer is WRONG! Should be 78,498!");
|
2021-04-09 22:48:47 +08:00
|
|
|
}
|