python3 - <<'PY' N=7**6 M=lambda x:x-7*(x//7) D=lambda x:[(x//7**i)-7*((x//7**i)//7) for i in range(5,-1,-1)] ax=[0]*N qx=[0]*N for x in range(N): a,b,c,d,e,f=D(x) ax[x]=sum([M(-b)*16807,M(a-b)*2401,M(-d)*343,c*49,M(-f)*7,e]) qx[x]=M(sum([(2*a-b)**2,(c-d)**2,(e-f)**2,b*sum([c,d,e,f])])) ids=[] for x in range(N): y=x