# Find expansion of polynomial in an ideal

I have a polynomial `p`

and some other polynomials `p_1,...,p_k`

which are elements of a multivariate polynomial ring. Say something like `P = PolynomialRing(QQ,'a,b,c,d,e,f,g')`

. I know that `p`

belongs to the ideal generated by `p_1,..,p_k`

because when I ask for a groebner basis of `I`

I can see explicitly `p`

there. How do I find the expression of `p`

as a linear combination of the `p_i`

's with coefficients in `P`

?