Improved version of poly_mulmod.

Polynomial division is replaced by some subtractions, exploiting the fact that x^n = -1 (mod (x^n + 1)).
Additionally this allows to get rid of a temporal buffer which held the intermediate result.