Tripling-oriented Doche–Icart–Kohel curve explained

The tripling-oriented Doche–Icart–Kohel curve is a form of an elliptic curve that has been used lately in cryptography; it is a particular type of Weierstrass curve. At certain conditions some operations, as adding, doubling or tripling points, are faster to compute using this form.The Tripling oriented Doche–Icart–Kohel curve, often called with the abbreviation 3DIK has been introduced by Christophe Doche, Thomas Icart, and David R. Kohel in ^[1]

Definition

Let

K be a field of characteristic different form 2 and 3.

An elliptic curve in tripling oriented Doche–Icart–Kohel form is defined by the equation:

T _{a} : y ^{2} = x ^{3} + 3 a(x +1 ) ^{2}

with

a\in K .

A general point P on

T _{a} has affine coordinates (x,y) . The "point at infinity" represents the neutral element for the group law and it is written in projective coordinates as O = (0:1:0). The negation of a point P = (x, y) with respect to this neutral element is -P = (x, -y).

The group law

Consider an elliptic curve in the Tripling-oriented Doche-Icart-Kohel form in affine coordinates:

T _{a:} y ^{2=} x ^{3} + 3 a(x +1 ) ^{2,} a ≠ 0 , \tfrac{ 9 }{ 4 }.

As in other forms of elliptic curves, it is possible to define some "operations" between points, such as adding points, or doubling (See also The group law). In the following sections formulas to add, negate and doubling points are given. The addition and doubling formulas are often used for other operations: given a point P on an elliptic curve it is possible to compute [n]P, where n is an integer, using addition and doubling; computing multiples of points is important in elliptic curve cryptography and in Lenstra elliptic curve factorization.

Addition

Given

P _{1=(x} _{1,y} _{1)} and P _{2=(x} _{2,y} _{2)} on T _{a} , the point P _{3=(x} _{3,y} _{3)=P} _{1+P} _{2} has coordinates:\begin{align} x _{3} & = (1)/((x _{2} - x ^{2} _{1 )}) \left \{ - x ^{3+ (x} _{2-3 a)x} ^{2+6 ax} ^{2 )x} _{1+(y} ^{2-2 y} ^{2 y} _{1+(-x} ^{2 ))} ^{2} \right \} \\ y _{3} & = (1)/((x _{2} - x ^{3} _{1 )}) \left \{ ( - y _{1+2y} _{2)x} ^{3+ ( -3 ay} ^{1-3 y} _{2x} _{2+3ay} _{2)x} ^{2 y} ^{1+6 ax} _{2y} _{1-6ay} _{2x} _{2)x} _{1} \right.\\ & \left. + (y ^{3-3 y} _{2 y} ^{2 )y} ^{1+ (y} _{2x} ^{3+3 ay} ^{2 x} ^{3 ))} ^{2} \right \} \end{align}

Doubling

Given a point

P _{1=(x} _{1,y} _{1)} on T _{a} , the point P _{3=(x} _{3,y} _{3)=2P} _{1} has coordinates:\begin{align} x _{3} & = (9)/(+ 4 y ^{2} _{1} x ^{4} _{1}) (9)/(y ^{3} _{1}) + \left ( (9)/(+ y ^{2 a} _{1} ^{2}) (9)/(y ^{2 a} _{1}) \right )x ^{2+} _{1} \left ( (18)/(y ^{2 a} _{1} ^{2}) -2 \right ) x _{1+ (9)/(y ^{2 a} _{1} ^{2-3a})} \\ y _{3} & =- (27)/(- 8 y ^{6} _{1}) (81)/(4 y ^{5} _{1}) + \left ( - (81)/(- 2 y ^{3 a} _{1} ^{2}) (81)/(4 y ^{3 a} _{1}) \right )x ^{4+} _{1} \left ( - (27)/(- y ^{3 a} _{1} ^{3}) (81)/(+ y ^{3 a} _{1} ^{2}) (9)/(2 y _{1}) \right )x ^{3+} _{1} \left ( - (81)/(- y ^{3 a} _{1} ^{3}) (81)/(2 y ^{3} _{1}) a ^{2+ (27)/(2 y _{1a})} \right)x ^{2} _{1} \\ & + \left ( - (81)/(+ y ^{3 a} _{1} ^{3}) (9)/(+ y ^{2} _{1 a}) (9)/(y _{1a}) \right )x _{1} + \left ( - (27)/(+ y ^{3 a} _{1} ^{3}) (9)/(y ^{2} _{1 a}) - y _{1} \right ) \end{align}

Negation

Given a point

P _{1=(x} _{1,y} _{1)} on T _{a} , its negation with respect to the neutral element ( 0 : 1 : 0 ) is - P _{1=(x} _{1,-y} _{1)} .

There are also other formulas given in ^[2] for Tripling-oriented Doche–Icart–Kohel curves for fast tripling operation and mixed-addition.

New Jacobian coordinates

For computing on these curves points are usually represented in new Jacobian coordinates (Jⁿ):

a point in the new Jacobian coordinates is of the form

P = (X : Y : Z : Z ^{2)} ; moreover:P = (X:Y:Z:Z ^{2} ) = ( λ ^{2} X: λ ^{3} Y: λ Z: λ ^{2} Z ^{2} ),

for any

λ \in K .

This means, for example, that the point

P = (X:Y:Z:Z ^{2)} and the point Q = ( 4 X: 8 Y: 2 Z: 4 Z ^{2)} (for λ=2 ) are actually the same.P = (x,y) on T _{a} is written in the new Jacobian coordinates as P = (X : Y : Z : Z ^{2)} , where x = X/Z ^{2} and y = Y/Z ^{3} ; in this way, the equation for T _{a} becomes:T _{a} : Y ^{2} = X ^{3} + 3 aZ ^{2(X+Z} ^{2)} ^{2.}

The term

Z ^{2} of a point on the curve makes the mixed addition (that is the addition between two points in different systems of coordinates) more efficient.

The neutral element in new Jacobian coordinates is

( 1 : 1 : 0 : 0 ) .

Algorithms and examples

Addition

The following algorithm represents the sum of two points

P _{1} and P _{2} on an elliptic curve in the Tripling-oriented Doche-Icart-Kohel form. The result is a point P _{3=(X} _{3,Y} _{3,Z} _{3,ZZ} _{3)} .It is assumed that Z _{2=1} and that a _{3=3a} .The cost of this implementation is 7M + 4S + 1*a3 + 10add + 3*2 + 1*4, where M indicates the multiplications, S the squarings, a3 indicates the multiplication by the constant a₃, add represents the number of additions required.A = X _{2ZZ} _{1} B = Y _{2ZZ} _{1Z} _{1} C = X _{1-A} D = 2 (Y _{1-B)} F = C ^{2} F _{4} = 4 F Z _{3} = (Z ^{2- ZZ} _{1- F} E = {Z ^{2} _{3 }} G = CF _{4} H = AF _{4} X _{3} = D ^{2- G -2 H - a} _{3 E} Y _{3} = D(H - X _{3)-2BG} ZZ _{3} = E

Example

Let

P _{1=(1,\sqrt{13})} and P _{2=(0,\sqrt{3})} affine points on the elliptic curve over R :y ^{2} = x ^{3} +3 (x +1 ) ^{2} .

Then:

A = X _{2ZZ} _{1} = 0 B = Y _{2ZZ} _{1Z} _{1} = \sqrt{ 3 } C = X _{1-A} = 1 D = 2 (Y _{1-B)} = 2 (\sqrt{ 13 } - \sqrt{ 3 }) F = C ^{2} = 1 F _{4} = 4 F = 4 Z _{3} = (Z ^{2- ZZ} _{1- F} = 2 E = {Z ^{2} _{3 }} = 4 G = CF _{4} = 4 H = AF _{4} = 0 X _{3} = D ^{2- G -2 H - a} _{3 E} = 48-8 \sqrt{ 39 } Y _{3} = D(H - X _{3)-2BG} = 296 \sqrt{ 3 } -144 \sqrt{ 13 } ZZ _{3} = E = 4

Notice that in this case

Z _{1=Z} _{2=1} .The resulting point is P _{3} = (X _{3,Y} _{3,Z} _{3,ZZ} _{3)} = ( 48-8 \sqrt{ 39 }, 296 \sqrt{ 3 } -144 \sqrt{ 13 }, 2 , 4 ) , that in affine coordinates is P _{3=(12-2\sqrt{39},37\sqrt{3}-18\sqrt{13})} .

Doubling

The following algorithm represents the doubling of a point

P _{1} on an elliptic curve in the Tripling-oriented Doche-Icart-Kohel form.It is assumed that a _{3=3a} , a _{2=2a} .The cost of this implementation is 2M + 7S + 1*a2 + 1*a3 + 12add + 2*2 + 1*3 + 1*8; here M represents the multiplications, S the squarings, a2 and a3 indicates the multiplications by the constants a₂ and a₃ respectively, and add indicates the additions.A = {X ^{2} _{1 }} B = a _{2ZZ} _{1(X} _{1+ZZ} _{1)} C = 3 (A + B) D = {Y ^{2} _{1 }} E = D ^{2} Z _{3} = (Y _{1+Z} ^{2- D - ZZ} ^{1} ZZ _{3} = Z ^{2} _{3} F = 2 ((X ^{2- A - E)} _{1+ D)} X _{3} = C ^{2- a} _{3 ZZ} _{3-2F} Y _{3} = C(F - X _{3)-8E}

Example

Let

P _{1=(0,\sqrt{3})} be a point on y ^{2} = x ^{3} + 3 (x +1 ) ^{2} .

Then:

A = {X ^{2} _{1 }} = 0 B = a _{2ZZ} _{1(X} _{1+ZZ} _{1)} = 2 C = 3 (A + B) = 6 D = {Y ^{2} _{1 }} = 3 E = D ^{2} = 9 Z _{3} = (Y _{1+Z} ^{2- D - ZZ} ^{1} = 2 \sqrt{ 3 } ZZ _{3} = Z ^{2} _{3} = 12 F = 2 ((X ^{2- A - E)} _{1+ D)} = 0 X _{3} = C ^{2- a} _{3 ZZ} _{3-2F} = 0 Y _{3} = C(F - X _{3)-8E} = -72

Notice that here the point is in affine coordinates, so

Z _{1=1} .The resulting point is P _{3=(0,-72,2\sqrt{3},12)} , that in affine coordinates is P _{3} = ( 0 , - \sqrt{ 3 }) .

Equivalence with Weierstrass form

Any elliptic curve is birationally equivalent to another written in the Weierstrass form.

The following twisted tripling-oriented Doche-Icart-Kohel curve:

T _{a,λ} : y ^{2} = x ^{3} + 3λ a(x +λ ) ^{2}

can be transformed into the Weierstrass form by the map:

(x,y)\mapsto (x -λ a, y).

In this way

T _{a,λ} becomes:y ^{2} = x ^{3} - 3 { λ } ^{2a(a-2)x} + { λ } ^{3a(2a} ^{2-6a+3)} .

Conversely, given an elliptic curve in the Weierstrass form:

E _{c,d} : y ^{2} = x ^{3} + cx ^{2} + d ,

it is possible to find the "corresponding" Tripling-oriented Doche–Icart–Kohel curve, using the following relation:

λ = (-3 d(a -2 ))/(a( 2 a ^{2-6a+3)})

where is a root of the polynomial

6912 a(a -2 ) ^{3-j(4a-9),}

where

j = (6912 c ^{3})/(4 c ^{3+27d} ^{2})

is the j-invariant of the elliptic curve

E _{c,d} .

Notice that in this case the map given is not only a birational equivalence, but an isomorphism between curves.

Internal link

For more informations about the running-time required in a specific case, see Table of costs of operations in elliptic curves

External links

• http://hyperelliptic.org/EFD/g1p/index.html

References

• Book: Christophe Doche . Thomas Icart . David R. Kohel . amp . 2006 . appeared at PKC 2006, part of LNCS (Lecture Series in Computer Science) volume number 3958 . Efficient Scalar Multiplication by Isogeny Decompositions . Springer Verlag . 285–352 .
• Book: Daniel J. Bernstein, Tanja Lange . 2007 . appeared in G.L. Mullen, D. Panario, I.E. Shparlinski (eds.), Finite Fields and Applications (Proceedings 8th International Conference, Fq8, Melbourne, Australia, July 9–13, 2007) . Analysis and optimization of elliptic-curve single-scalar multiplication . Mathematics Subject Classification .
• Book: D.J.Bernstein, P.Birkner, T.Lange, and C.Peters . 2007 . appeared in K. Srinathan, C. Pandu Rangan, M. Yung (Eds.), Proceedings of the 8th International Conference on Cryptology in India: Progress in Cryptology (Indocrypt 2007) 9–13 December 2007, Chennai, India . Optimizing Double-Base Elliptic-Curve Single-Scalar Multiplication . Springer .
• http://hyperelliptic.org/EFD/g1p/auto-3dik-standard.html

Notes and References

1. Christophe Doche, Thomas Icart, and David R. Kohel, Efficient Scalar Multiplication by Isogeny Decompositions
2. Christophe Doche, Thomas Icart, and David R. Kohel, Efficient Scalar Multiplication by Isogeny Decompositions, page 198-199