CMakeLists.txt

pojcs101_M03247_palindromic_prime_init1.py

@@ -0,0 +1,29 @@
+import math
+
+MAX_N = math.ceil(1e9 ** 0.5) + 5
+
+primes: list[int] = [2, 3, 5, 7, 11]
+
+
+def append_prime():
+    for next_number in range(primes[-1]+2, MAX_N + 10000, 2):
+        for p in primes:
+            if next_number % p == 0:
+                break
+            if p * p > next_number:
+                primes.append(next_number)
+                print(next_number)
+                return
+
+
+def write_primes(file_name):
+    with open(file_name, 'w') as f:
+        f.write(", ".join(map(str, primes)))
+
+
+if __name__ == '__main__':
+    while primes[-1] < MAX_N:
+        append_prime()
+
+    print(len(primes))
+    write_primes('pojcs101_M03247_palindromic_prime_init1.out')