Revision history [back]

Perform these two steps at once:

    F = E.division_polynomial(N)
    ...
        z = F(x)

By passing x to division_polynomial. Make sure that x is in the field R, and move the line after where x is defined:

z = E.division_polynomial(N, R(x))

With this change, your program runs within 2 seconds on my computer.