python3 - <<'PY' p=996979 d=14641 def sums(r): s2=s3=0 for x in r: a=(d*x)%p s2=(s2+a*a)%983 s3=(s3+a*a*a)%983 return s2,s3 B=range(509847,509896) A1=range(221445,221491) A2=range(884827,884872) sb=sums(B); sa1=sums(A1); sa2=sums(A2) print(sb,sa1,sa2) print((91*sb[0]+sa1[0]+sa2[0])%983,(91*sb[1]+sa1[1]+sa2[1])%983) PY