D=pow(379153,-1,p) B2=0 B3=0 for X in Z: for x in X: a=D*x%p B2=(B2+a*a)%983 B3=(B3+a*a*a)%983 print('DIRECT',len(Z),len(next(iter(Z))),B2,B3)