p = 7 F.<a> = GF(p^2) %timeit a*a %timeit a+a
K = GF(p).algebraic_closure() G = K.subfield(2)[0] b = G.gen() %timeit b*b %timeit b+b