>
 

MATH 392 -- Seminar in Computational Commutative Algebra 

Generalized Resultants for Elimination 

October 25, 2006 

 

When we have an ideal generated by several polynomials, 

the generalized resultant is computed by introducing  

new u-variables as follows. 

> FList := [x^4*y^2+z, x^2*z-y+1, x^2+4*y*z+2]; 1
 

[x^4*y^2+z, x^2*z-y+1, x^2+4*y*z+2] 

> gRes := resultant(FList[1], u[2]*FList[2]+u[3]*FList[3], x); 1
 

(16*y^4*z^2*u[3]^2-8*y^4*u[3]*z*u[2]+8*y^3*u[2]*z*u[3]+16*y^3*u[3]^2*z+u[2]^2*y^4-2*y^3*u[2]^2-4*y^3*u[2]*u[3]+u[2]^2*y^2+4*u[2]*y^2*u[3]+4*y^2*u[3]^2+u[2]^2*z^3+2*u[2]*z^2*u[3]+z*u[3]^2)^2
(16*y^4*z^2*u[3]^2-8*y^4*u[3]*z*u[2]+8*y^3*u[2]*z*u[3]+16*y^3*u[3]^2*z+u[2]^2*y^4-2*y^3*u[2]^2-4*y^3*u[2]*u[3]+u[2]^2*y^2+4*u[2]*y^2*u[3]+4*y^2*u[3]^2+u[2]^2*z^3+2*u[2]*z^2*u[3]+z*u[3]^2)^2
(16*y^4*z^2*u[3]^2-8*y^4*u[3]*z*u[2]+8*y^3*u[2]*z*u[3]+16*y^3*u[3]^2*z+u[2]^2*y^4-2*y^3*u[2]^2-4*y^3*u[2]*u[3]+u[2]^2*y^2+4*u[2]*y^2*u[3]+4*y^2*u[3]^2+u[2]^2*z^3+2*u[2]*z^2*u[3]+z*u[3]^2)^2
 

>
 

We can expand this polynomial, and collect like terms in the u-variables: 

> gRes := collect(expand(gRes), [u[2], u[3]]); 1
 

(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
(-4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3)*u[2]^4+(-8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8*y^4-16*y^4*z^4+8*y^2*z^3-16*y^8*z+4*y^4*z^2-8*y^3*z^2+16*y^3*z^4+48*y^7*z-24*y^...
 

>
 

We claim that each of the coefficients of the monomials u^alpha 

here is an element of the elimination ideal To verify that 

claim, we will compute a Grobner basis of  I = <FList>: 

 

> with(Groebner); -1
 

> B := Basis(FList, plex(x, y, z)); 1
 

[256*z^9+16*z^3-32*z^2+z+64*z^4+96*z^5+4+256*z^7, 16*y-128*z^8-20*z^3+14*z^2+25*z-24*z^4-96*z^6+64*z^7+48*z^5-12, 8*x^2+23*z-18*z^2-64*z^7-44*z^3-48*z^5+12-128*z^8-24*z^4-96*z^6]
[256*z^9+16*z^3-32*z^2+z+64*z^4+96*z^5+4+256*z^7, 16*y-128*z^8-20*z^3+14*z^2+25*z-24*z^4-96*z^6+64*z^7+48*z^5-12, 8*x^2+23*z-18*z^2-64*z^7-44*z^3-48*z^5+12-128*z^8-24*z^4-96*z^6]
[256*z^9+16*z^3-32*z^2+z+64*z^4+96*z^5+4+256*z^7, 16*y-128*z^8-20*z^3+14*z^2+25*z-24*z^4-96*z^6+64*z^7+48*z^5-12, 8*x^2+23*z-18*z^2-64*z^7-44*z^3-48*z^5+12-128*z^8-24*z^4-96*z^6]
[256*z^9+16*z^3-32*z^2+z+64*z^4+96*z^5+4+256*z^7, 16*y-128*z^8-20*z^3+14*z^2+25*z-24*z^4-96*z^6+64*z^7+48*z^5-12, 8*x^2+23*z-18*z^2-64*z^7-44*z^3-48*z^5+12-128*z^8-24*z^4-96*z^6]
 

Then "pick off" the coefficients of each monomial in the  u-variables. 

Note that each contains only  y,z 

> cList := [coeffs(gRes, [u[2], u[3]])]; 1
 

[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
[z^2+32*y^3*z^2+128*y^5*z+8*y^2*z+16*y^4+512*y^7*z^3+384*y^6*z^2+256*y^8*z^4+32*y^4*z^3, -4*y^5+y^4+y^8+6*y^6-4*y^7-4*y^3*z^3+z^6+2*y^4*z^3+2*y^2*z^3, -8*y^7+16*y^5*z-48*y^6*z+24*y^6+4*z^5+4*y^2*z^2+8...
 

>
 

Then we divide the coefficients by the Grobner basis: 

> for i to nops(cList) do NormalForm(cList[i], B, plex(x, y, z)) end do; 1
 

0 

0 

0 

0 

0 

This shows that every one of the coefficients of the generalized resultant is 

in the elimination ideal.    

 

CAUTION:  While these polynomials are in the elimination ideal, they do  

not necessarily generate the elimination ideal.  For instance, in this case,  

we can see that the coefficients generate an ideal strictly contained in the  

elimination ideal: 

> BB := Basis(cList, plex(y, z)); 1
 

[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
[-255*z^4+16*z^2+8*z^3+64*z^5+131072*z^18-2048*z^11+114688*z^16-2048*z^9-128*z^7+1568*z^6-3648*z^8+16384*z^13+57344*z^14+17920*z^12+65536*z^20+7680*z^10+32768*z^15, 16*y*z^3+64*y*z^2+465*z^4-48*z^2+88...
 

Recall, reduced Grobner bases are unique.  These polynomials are not the 

same as the Grobner basis for the elimination ideal obtained from B above. 

For instance, the univariate polynomial in z  in BB is a multiple of the  

univariate polynomial in z from B: 

> rem(BB[1], B[1], z); 1
 

0 

> factor(BB[1]); 1
 

z^2*(256*z^9+16*z^3-32*z^2+z+64*z^4+96*z^5+4+256*z^7)^2 

> factor(B[1]); 1
 

256*z^9+16*z^3-32*z^2+z+64*z^4+96*z^5+4+256*z^7 

>