p=5: A~ = 1*X^1*Y^0  (mod 5); delta,eps,m = (1, 0, 0)
p=7: A~ = 1*X^0*Y^1  (mod 7); delta,eps,m = (0, 1, 0)
p=11: A~ = 1*X^1*Y^1  (mod 11); delta,eps,m = (1, 1, 0)
p=13: A~ = 8*X^0*Y^2 + 6*X^3*Y^0  (mod 13); delta,eps,m = (0, 0, 1)
# Serre's examples p = 5, 7, 11, 13 reproduced: True
# primes 5 <= p <= 1000: 166; A~ not squarefree for: []
