// CANONICAL LIFTING -- CHARACTERISTIC p = 2


a1 := a0^2;

x1 := a0;

y1 := (x0^2 + x0)*y0 + x0^3 + a0*x0^2 + a0;

a2 := a0^4;

x2 := a0*x0^2 + x0^6;

y2 := (a0^2 + a0^4) + (a0 + a0^3)*x0^2 + 
      a0^4*x0^4 + (1 + a0^2)*x0^5 + (1 + a0 + a0^2)*
        x0^6 + ((1 + a0^2)*x0^3 + a0*x0^4 + x0^6)*y0;

a3 := a0^8;

x3 := a0^8 + a0*x0^6 + (a0+a0^2+a0^4)*x0^8 + x0^10 + x0^120;

y3 := a0^4 + a0^6 + a0^8 + a0^5*x0^2 + 
  (a0^2 + a0^3 + a0^4 + a0^8 + a0^10)*x0^4 + a0^4*x0^5 + 
  (a0 + a0^2 + a0^3 + a0^7 + a0^9)*x0^6 + 
  (a0^2 + a0^3 + a0^4)*x0^7 + 
  (a0 + a0^2 + a0^3 + a0^4 + a0^6)*x0^8 + 
  (1 + a0 + a0^2 + a0^4 + a0^6 + a0^8)*x0^9 + 
  (1 + a0 + a0^2 + a0^4 + a0^6 + a0^8)*x0^10 + 
  (1 + a0^2 + a0^3)*x0^11 + (1 + a0^2 + a0^3 + a0^4 + a0^8)*
   x0^12 + (1 + a0^2 + a0^4)*x0^13 + 
  (1 + a0 + a0^2 + a0^4)*x0^14 + x0^16 + 
  (a0^4*x0^3 + a0^5*x0^4 + a0^2*x0^5 + (a0^2 + a0^3)*x0^6 + 
    (1 + a0^3 + a0^4 + a0^6 + a0^8)*x0^7 + 
    (a0 + a0^3 + a0^5)*x0^8 + (1 + a0)*x0^9 + 
    (1 + a0^3 + a0^4)*x0^10 + (a0 + a0^4)*x0^11 + 
    (a0 + a0^2)*x0^12 + x0^13 + x0^14)*y0;
