'Polynomial exact pseudo-quotient of ``f`` by ``g``.
Examples
>>> from sympy import Poly
>>> from sympy.abc import x
>>> Poly(x**2 - 1, x).pexquo(Poly(2*x - 2, x))
Poly(2*x + 2, x, domain=\'ZZ\')
>>> Poly(x**2 + 1, x).pexquo(Poly(2*x - 4, x))
Traceback (most recent call last):
ExactQuotientFailed: 2*x - 4 does not divide x**2 + 1'

def pexquo(f, g):

    (_, per, F, G) = f._unify(g)
    if hasattr(f.rep, 'pexquo'):
        try:
            result = F.pexquo(G)
        except ExactQuotientFailed as exc:
            raise exc.new(f.as_expr(), g.as_expr())
    else:
        raise OperationNotSupported(f, 'pexquo')
    return per(result)
