polynomialQuotientRemainder[b^4--2*b^3--4*b^2--4*b--1,b^2--1,b]