def METHOD_NAME(alpha_s, p):
    lambda_s = np.zeros((1, len(alpha_s)), dtype="int64")
    for i in range(len(alpha_s)):
        cur_alpha = alpha_s[i]
        den = PI([cur_alpha - o for o in alpha_s if cur_alpha != o], p)
        num = PI([0 - o for o in alpha_s if cur_alpha != o], p)
        lambda_s[0][i] = divmod(num, den, p)
    return lambda_s.astype("int64")