Calculating Groebner base with constants

edit

I want to see what the Groebner base of an ideal is, but the polynomials generating it have some constants inside. Is there any way to do that with Sage?

var('a')
R.<x,y> = PolynomialRing(CC, 2, order='lex')
I = ideal(y^4*a, y^2*x - y*a - x, y - x*a)
B = I.groebner_basis()
B

This doesn't seem to work. a is the constant. Also, this example is just a dummy, I want to know in general how to have constant coefficient in polynomials.