def _dmp_zz_gcd_interpolate(h, x, v, K):
    f = []

    while not dmp_zero_p(h, v):
        g = dmp_ground_trunc(h, x, v, K)
        f.insert(0, g)

        h = dmp_sub(h, g, v, K)
        h = dmp_quo_ground(h, x, v, K)

    if K.is_negative(dmp_ground_LC(f, v + 1, K)):
        return dmp_neg(f, v + 1, K)
    else:
        return f
Interpolate polynomial GCD from integer GCD.
dataset/ETHPy150Open sympy/sympy/sympy/polys/euclidtools.py/_dmp_zz_gcd_interpolate