multi_types_solver.py 文件源码

def shift_polynomial(coef, d):
    # given polynomial sum_k a_k x^k -> sum_k a_k (x+d)^k
    coef2 = np.zeros(5, dtype="float64")
    for i in range(5):
        for j in range(i, 5):
            coef2[i] = coef2[i] + coef[j] * (d**(j - i)) * sp.binom(j, i)
    return coef2