How to find coefficient polynomials in an ideal

savecancel

h.lift(I) should do the job. See documentation for details.

Since lift is not implemented for BooleanPolynomialRing, we can workaround by considering the ideal generated by $f,g$ and $x_i^2+x_i$ in $GF(2)[x_1,\dots,x_n]$. Here is a sample code:

B.<x,y> = BooleanPolynomialRing()
I = ideal(x+1,y+1)
pol = x*y + 1
#print( pol.lift(I) )      # this is not implemented

#workaround:
F = GF(2)[B.gens()]
I2 = ideal( [F(g) for g in I.gens()] + [g^2+g for g in F.gens()] )
coefs = [B(c) for c in F(pol).lift(I2)[:len(I.gens())]]
print( coefs )