r= x^16 - x^8 + 1
q= (x^22 - 2*x^19 + x^16 - x^14 + 2*x^11 - x^8 + x^6 +x^3 + 1)  / 3
P= q_k16method66.numerator()
v= x
curve_side_security=128
field_side_security=146
gfield_side_security=144
best_results_file=d-k16method66-128-4.txt
best parametes :
5608.      18.36      80.60      146.6      603.1      998.8         128

curve_side_security=192
field_side_security=192
gfield_side_security=194
best_results_file=d-k16method66-192-4.txt
best parametes :
10090.     24.70      104.2      192.4      812.5      1661.         230

curve_side_security=256
field_side_security=256
gfield_side_security=256
best_results_file=d-k16method66-256-4.txt
best parametes :
18940.     32.00      136.9      256.5      1041.      2882.         431