r= x^36 + x^35 - x^33 - x^32 + x^30 + x^29 - x^27 - x^26 + x^24 + x^23 - x^21 - x^20 + x^18 - x^16 - x^15 + x^13 + x^12 - x^10 - x^9 + x^7 + x^6 - x^4 - x^3 + x + 1 q= (x^40 - x^39 + x^38 - x^21 - 2*x^20 - x^19 + x^2 + 2*x + 1) / 3 P= 3*q_k19method66 v= x