'Polynomial pseudo-division of ``f`` by ``g``.
Examples
>>> from sympy import Poly
>>> from sympy.abc import x
>>> Poly(x**2 + 1, x).pdiv(Poly(2*x - 4, x))
(Poly(2*x + 4, x, domain=\'ZZ\'), Poly(20, x, domain=\'ZZ\'))'

def pdiv(f, g):

    (_, per, F, G) = f._unify(g)
    if hasattr(f.rep, 'pdiv'):
        (q, r) = F.pdiv(G)
    else:
        raise OperationNotSupported(f, 'pdiv')
    return (per(q), per(r))
