def shoup_poly(n, p, K):
    f = [K.one] * (n + 1)
    for i in range(1, n + 1):
        f[i] = (f[1 - i]**2 + K.one) % p
    return f
dataset/ETHPy150Open sympy/sympy/sympy/polys/benchmarks/bench_galoispolys.py/shoup_poly/swapped operands of `-`