PolynomialQuotientRemainder[x^3 - 3 x - 1, x^2 + 2, x]