// CANONICAL LIFTING -- CHARACTERISTIC p = 5


a1:= (a0^3*b0^2 + b0^4)/a0;

b1:= 4*a0^6*b0 + a0^3*b0^3 + b0^5;

x1:=  4/a0*x0^7 + 4*b0/a0*x0^4 + a0*x0^3 + 3*b0*x0^2 + 3*b0^2/a0*x0 + a0*b0;

P1:=  1/a0*x0^8 + 2*x0^6 + 3*b0/a0*x0^5 + 2*a0*x0^4 + 3*b0*x0^3 + a0^2*x0^2 + 
    3*a0*b0*x0 + 3*a0^3 + 3*b0^2;

a2:= (2*a0^36 + a0^33*b0^2 + a0^30*b0^4 + 3*a0^27*b0^6 + 2*a0^24*b0^8 + 
    a0^18*b0^12 + 4*a0^12*b0^16 + 3*a0^9*b0^18 + 4*a0^6*b0^20 + 4*a0^3*b0^22 + 
    4*b0^24)/a0^11;

b2:= a0^36*b0 + 4*a0^33*b0^3 + 3*a0^27*b0^7 + 4*a0^21*b0^11 + 4*a0^15*b0^15 + 
    a0^12*b0^17 + 3*a0^6*b0^21 + b0^25;

x2 :=  2/a0^10*x0^45 + 2/a0^6*x0^37 + (a0^6 + a0^3*b0^2 + b0^4)/a0^11*x0^35 + 
    2*b0/a0^6*x0^34 + 3/a0^4*x0^33 + (4*a0^3 + 4*b0^2)/a0^6*x0^31 + 
    4*b0^5/a0^10*x0^30 + 4*b0^2/a0^5*x0^29 + (a0^3 + 3*b0^2)/a0^4*x0^27 + 
    2*b0/a0^2*x0^26 + (4*a0^3*b0 + 3*b0^3)/a0^4*x0^24 + (4*a0^6 + 3*a0^3*b0^2 + 
    b0^4)/a0^5*x0^23 + (2*a0^3*b0^3 + 4*b0^5)/a0^6*x0^22 + (a0^3*b0^2 + 
    b0^4)/a0^4*x0^21 + (2*a0^12*b0 + 4*a0^9*b0^3 + 2*a0^6*b0^5 + a0^3*b0^7 + 
    b0^9)/a0^11*x0^20 + (a0^9 + 3*a0^6*b0^2 + 4*b0^6)/a0^6*x0^19 + (2*a0^6*b0 + 
    a0^3*b0^3 + 4*b0^5)/a0^4*x0^18 + (2*a0^6 + a0^3*b0^2 + b0^4)/a0^2*x0^17 + 
    (3*a0^6*b0^3 + 4*a0^3*b0^5 + 3*b0^7)/a0^6*x0^16 + (4*a0^9 + 4*a0^6*b0^2 + 
    3*b0^6)/a0^4*x0^15 + (4*a0^6*b0 + 2*b0^5)/a0^2*x0^14 + (3*a0^9 + 4*a0^6*b0^2
    + 3*a0^3*b0^4 + 2*b0^6)/a0^3*x0^13 + (3*a0^9*b0 + 2*a0^6*b0^3 + a0^3*b0^5 + 
    2*b0^7)/a0^4*x0^12 + (2*a0^9*b0^2 + 3*a0^6*b0^4 + 4*a0^3*b0^6 + 
    2*b0^8)/a0^5*x0^11 + (2*a0^6*b0^3 + 4*a0^3*b0^5 + 3*b0^7)/a0^3*x0^10 + 
    (3*a0^9*b0^2 + 2*a0^6*b0^4 + 4*a0^3*b0^6 + 4*b0^8)/a0^4*x0^9 + (2*a0^9*b0^3 
    + 2*a0^6*b0^5 + 4*a0^3*b0^7 + b0^9)/a0^5*x0^8 + (3*a0^9*b0^4 + 3*a0^6*b0^6 +
    a0^3*b0^8 + 2*b0^10)/a0^6*x0^7 + (3*a0^9*b0^3 + 3*a0^3*b0^7 + 
    2*b0^9)/a0^4*x0^6 + (4*a0^21 + 3*a0^15*b0^4 + a0^12*b0^6 + 2*a0^9*b0^8 + 
    a0^6*b0^10 + 2*a0^3*b0^12 + 2*b0^14)/a0^11*x0^5 + (a0^9*b0^5 + 3*a0^6*b0^7 +
    4*a0^3*b0^9 + 2*b0^11)/a0^6*x0^4 + (4*a0^9*b0^4 + 4*a0^6*b0^6 + a0^3*b0^8 + 
    3*b0^10)/a0^4*x0^3 + (2*a0^9*b0^5 + 2*a0^6*b0^7 + 4*b0^11)/a0^5*x0^2 + 
    (2*a0^9*b0^6 + 4*b0^12)/a0^6*x0 + (4*a0^21*b0 + a0^18*b0^3 + 4*a0^15*b0^5 + 
    a0^12*b0^7 + a0^9*b0^9 + 3*a0^6*b0^11 + 3*b0^15)/a0^10;

P2 :=  4/a0^10*x0^56 + 3/a0^9*x0^54 + 3*b0/a0^10*x0^53 + 4/a0^8*x0^52 + 
    3*b0/a0^9*x0^51 + 4*b0^2/a0^10*x0^50 + 3/a0^6*x0^48 + (4*a0^3*b0^2 + 
    4*b0^4)/a0^11*x0^46 + 4*b0/a0^6*x0^45 + (3*a0^6 + 3*a0^3*b0^2 + 
    3*b0^4)/a0^10*x0^44 + (3*a0^3*b0^3 + 3*b0^5)/a0^11*x0^43 + (2*a0^6 + 
    4*a0^3*b0^2 + 4*b0^4)/a0^9*x0^42 + (3*a0^6*b0 + 3*a0^3*b0^3 + 
    4*b0^5)/a0^10*x0^41 + (a0^9 + 3*a0^6*b0^2 + 4*a0^3*b0^4 + 
    4*b0^6)/a0^11*x0^40 + (2*a0^6*b0 + 2*b0^5)/a0^9*x0^39 + (2*a0^9 + a0^6*b0^2 
    + 2*b0^6)/a0^10*x0^38 + (2*a0^6*b0 + 4*a0^3*b0^3 + b0^5)/a0^8*x0^37 + 
    (a0^6*b0^2 + 4*a0^3*b0^4 + 2*b0^6)/a0^9*x0^36 + (4*a0^6*b0^3 + 
    b0^7)/a0^10*x0^35 + 3*b0^4/a0^5*x0^34 + (4*a0^3*b0^3 + 3*b0^5)/a0^6*x0^33 + 
    (2*a0^6 + 2*b0^4)/a0^4*x0^32 + (a0^12*b0 + a0^9*b0^3 + 3*a0^6*b0^5 + 
    2*a0^3*b0^7 + 2*b0^9)/a0^11*x0^31 + (4*a0^9 + 4*b0^6)/a0^6*x0^30 + 
    (2*a0^12*b0 + 2*a0^9*b0^3 + 3*a0^6*b0^5 + 4*a0^3*b0^7 + 4*b0^9)/a0^10*x0^29 
    + (2*a0^12*b0^2 + a0^9*b0^4 + 4*a0^6*b0^6 + 4*a0^3*b0^8 + 
    4*b0^10)/a0^11*x0^28 + (4*a0^12*b0 + a0^9*b0^3 + 4*a0^6*b0^5 + 2*a0^3*b0^7 +
    2*b0^9)/a0^9*x0^27 + (4*a0^15 + 2*a0^12*b0^2 + 4*a0^3*b0^8 + 
    4*b0^10)/a0^10*x0^26 + (4*a0^15*b0 + a0^12*b0^3 + 2*a0^9*b0^5 + 2*a0^6*b0^7 
    + 2*a0^3*b0^9 + 2*b0^11)/a0^11*x0^25 + (4*a0^6 + 4*a0^3*b0^2 + 2*b0^4)*x0^24
    + (3*a0^9*b0 + 3*a0^6*b0^3 + 2*b0^7)/a0^4*x0^23 + (3*a0^9 + 3*a0^6*b0^2 + 
    4*a0^3*b0^4 + 2*b0^6)/a0^2*x0^22 + (2*a0^9*b0 + b0^7)/a0^3*x0^21 + 
    (3*a0^9*b0^2 + 4*a0^3*b0^6 + 2*b0^8)/a0^4*x0^20 + (4*a0^7*b0 + 
    a0^4*b0^3)*x0^19 + (3*a0^6*b0^2 + 3*a0^3*b0^4 + 2*b0^6)*x0^18 + (3*a0^9*b0 +
    4*a0^6*b0^3 + 3*a0^3*b0^5 + 3*b0^7)/a0*x0^17 + (3*a0^12*b0^2 + 
    2*b0^10)/a0^5*x0^16 + (a0^12*b0 + 2*a0^9*b0^3 + 4*a0^6*b0^5 + 4*a0^3*b0^7 + 
    4*b0^9)/a0^3*x0^15 + (4*a0^15 + a0^9*b0^4 + 2*a0^3*b0^8 + 
    4*b0^10)/a0^4*x0^14 + (2*a0^15*b0 + 3*a0^12*b0^3 + 3*a0^9*b0^5 + a0^6*b0^7 +
    a0^3*b0^9 + 4*b0^11)/a0^5*x0^13 + (4*a0^15 + 4*a0^9*b0^4 + 2*a0^6*b0^6 + 
    2*a0^3*b0^8 + 2*b0^10)/a0^3*x0^12 + (2*a0^18*b0^3 + a0^15*b0^5 + a0^12*b0^7 
    + a0^9*b0^9 + 4*a0^6*b0^11 + 4*b0^15)/a0^10*x0^11 + (a0^18 + 4*a0^15*b0^2 + 
    3*a0^9*b0^6 + 2*a0^6*b0^8 + 4*a0^3*b0^10 + 2*b0^12)/a0^5*x0^10 + (a0^21*b0 +
    4*a0^18*b0^3 + 4*a0^12*b0^7 + 4*a0^9*b0^9 + 3*b0^15)/a0^9*x0^9 + (4*a0^24 + 
    a0^21*b0^2 + 2*a0^18*b0^4 + a0^15*b0^6 + 4*a0^12*b0^8 + a0^9*b0^10 + 
    3*b0^16)/a0^10*x0^8 + (a0^21*b0 + a0^18*b0^3 + a0^15*b0^5 + 3*a0^12*b0^7 + 
    2*a0^9*b0^9 + 4*b0^15)/a0^8*x0^7 + (4*a0^24 + 2*a0^21*b0^2 + 3*a0^18*b0^4 + 
    3*a0^15*b0^6 + 4*a0^12*b0^8 + 3*a0^9*b0^10 + 3*b0^16)/a0^9*x0^6 + 
    (2*a0^24*b0 + 2*a0^21*b0^3 + a0^18*b0^5 + 4*a0^15*b0^7 + 2*a0^9*b0^11 + 
    4*b0^17)/a0^10*x0^5 + (2*a0^16 + a0^13*b0^2 + 4*a0^4*b0^8 + 4*a0*b0^10)*x0^4
    + (a0^15*b0 + 3*a0^12*b0^3 + 4*a0^6*b0^7 + 2*a0^3*b0^9 + b0^11)*x0^3 + 
    (4*a0^17 + 4*a0^14*b0^2 + 4*a0^11*b0^4 + a0^8*b0^6 + 4*a0^5*b0^8 + 
    2*a0^2*b0^10)*x0^2 + (a0^16*b0 + 4*a0^13*b0^3 + 4*a0^10*b0^5 + 4*a0^7*b0^7 +
    3*a0^4*b0^9 + a0*b0^11)*x0 + 4*a0^18 + a0^12*b0^4 + 3*a0^9*b0^6 + a0^6*b0^8 
    + 3*b0^12;

a3 := (2*a0^225 + 2*a0^219*b0^4 + a0^216*b0^6 + 3*a0^213*b0^8 +
        a0^210*b0^10 + a0^207*b0^12 + a0^204*b0^14 + 3*a0^198*b0^18 +
        4*a0^195*b0^20 + 4*a0^192*b0^22 + 4*a0^189*b0^24 +
        3*a0^186*b0^26 + 4*a0^183*b0^28 + a0^180*b0^30 +
        3*a0^177*b0^32 + 3*a0^174*b0^34 + 3*a0^171*b0^36 +
        4*a0^168*b0^38 + 2*a0^165*b0^40 + 2*a0^162*b0^42 +
        2*a0^159*b0^44 + 3*a0^156*b0^46 + 3*a0^153*b0^48 +
        a0^150*b0^50 + 2*a0^147*b0^52 + 3*a0^141*b0^56 + a0^138*b0^58
        + 2*a0^132*b0^62 + 4*a0^129*b0^64 + 3*a0^126*b0^66 +
        4*a0^123*b0^68 + 2*a0^120*b0^70 + 3*a0^117*b0^72 +
        a0^114*b0^74 + 3*a0^111*b0^76 + 4*a0^108*b0^78 + a0^102*b0^82
        + 4*a0^99*b0^84 + a0^96*b0^86 + 4*a0^93*b0^88 + a0^90*b0^90 +
        2*a0^87*b0^92 + 2*a0^84*b0^94 + 4*a0^81*b0^96 + 4*a0^78*b0^98
        + a0^75*b0^100 + 2*a0^72*b0^102 + a0^63*b0^108 + a0^60*b0^110
        + a0^57*b0^112 + a0^54*b0^114 + 2*a0^51*b0^116 +
        2*a0^48*b0^118 + 3*a0^45*b0^120 + a0^42*b0^122 + a0^39*b0^124
        + 2*a0^30*b0^130 + a0^15*b0^140 + 3*b0^150)/a0^100;

b3 := (4*a0^261*b0 + a0^258*b0^3 + 4*a0^255*b0^5 + 2*a0^252*b0^7 +
        4*a0^249*b0^9 + 2*a0^246*b0^11 + 3*a0^243*b0^13 +
        4*a0^240*b0^15 + 2*a0^237*b0^17 + 3*a0^228*b0^23 +
        2*a0^225*b0^25 + a0^222*b0^27 + 2*a0^216*b0^31 +
        4*a0^213*b0^33 + 3*a0^207*b0^37 + a0^204*b0^39 +
        2*a0^201*b0^41 + 4*a0^198*b0^43 + a0^192*b0^47 +
        3*a0^189*b0^49 + 3*a0^186*b0^51 + 4*a0^183*b0^53 +
        a0^180*b0^55 + 4*a0^177*b0^57 + 2*a0^174*b0^59 + a0^168*b0^63
        + 2*a0^165*b0^65 + 3*a0^162*b0^67 + 2*a0^159*b0^69 +
        3*a0^156*b0^71 + a0^153*b0^73 + 2*a0^147*b0^77 +
        3*a0^144*b0^79 + 3*a0^141*b0^81 + 2*a0^135*b0^85 +
        3*a0^132*b0^87 + a0^129*b0^89 + 3*a0^126*b0^91 +
        2*a0^123*b0^93 + a0^117*b0^97 + 3*a0^111*b0^101 +
        3*a0^105*b0^105 + 2*a0^99*b0^109 + 4*a0^96*b0^111 +
        4*a0^93*b0^113 + 3*a0^90*b0^115 + 3*a0^87*b0^117 +
        2*a0^75*b0^125 + a0^66*b0^131 + a0^63*b0^133 + 3*a0^57*b0^137
        + a0^54*b0^139 + 4*a0^51*b0^141 + 4*a0^48*b0^143 +
        a0^45*b0^145 + 4*a0^42*b0^147 + 2*a0^39*b0^149 +
        2*a0^36*b0^151 + a0^33*b0^153 + 2*a0^30*b0^155 +
        2*a0^27*b0^157 + a0^24*b0^159 + 2*a0^21*b0^161 +
        3*b0^175)/a0^75;
