python3 - <<'PY' def mod(a,b): return a-b*(a//b) def T(p,b): return mod(sum((73,29*p,31*b,17*p*b,11*p*p,13*b*b)),997) for p in range(5): row=[] for b in range(7): t=T(p,b) k=mod(t,4) u=mod(sum((t,p,2*b)),4) v=mod(sum((2*t,p,b)),4) w=mod(sum((t,3*p,b)),4) row.append((b,t,k,u,v,w)) print(p,row) PY