CPLIB_FPS_POLYNOMIAL_INTERPOLATION = 1
proc polynomialInterpolation[T: BarrettModint or MontgomeryModint]( xs, ys: seq[T]): seq[T]