def inv(z):
<EOL>z = qdiv(z)<EOL>z2 = z * z % PRIME  <EOL>z9 = pow2(z2, <NUM_LIT:2>) * z % PRIME  <EOL>z11 = z9 * z2 % PRIME  <EOL>z2_5_0 = (z11*z11) % PRIME * z9 % PRIME  <EOL>z2_10_0 = pow2(z2_5_0, <NUM_LIT:5>) * z2_5_0 % PRIME  <EOL>z2_20_0 = pow2(z2_10_0, <NUM_LIT:10>) * z2_10_0 % PRIME  <EOL>z2_40_0 = pow2(z2_20_0, <NUM_LIT:20>) * z2_20_0 % PRIME<EOL>z2_50_0 = pow2(z2_40_0, <NUM_LIT:10>) * z2_10_0 % PRIME<EOL>z2_100_0 = pow2(z2_50_0, <NUM_LIT:50>) * z2_50_0 % PRIME<EOL>z2_200_0 = pow2(z2_100_0, <NUM_LIT:100>) * z2_100_0 % PRIME<EOL>z2_250_0 = pow2(z2_200_0, <NUM_LIT:50>) * z2_50_0 % PRIME  <EOL>return pow2(z2_250_0, <NUM_LIT:5>) * z11 % PRIME<EOL>
$= z^{-1} \mod q$, for z != 0