python3 - <<'PY' p=996979 d=14641 B=list(range(509847,509896)) A1=list(range(221445,221491)) A2=list(range(884827,884872)) print(len(B),len(A1),len(A2),91*len(B)+len(A1)+len(A2)) C2=C3=S2=S3=0 for x in B: aa=(d*x)%p S2=(S2+x*x)%983; S3=(S3+x*x*x)%983 C2=(C2+aa*aa)%983; C3=(C3+aa*aa*aa)%983 C2=(C2*91)%983; C3=(C3*91)%983 S2=(S2*91)%983; S3=(S3*91)%983 for x in A1+A2: aa=(d*x)%p S2=(S2+x*x)%983; S3=(S3+x*x*x)%983 C2=(C2+aa*aa)%983; C3=(C3+aa*aa*aa)%983 print(C2,C3,S2,S3) PY