<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Multifunctional CRS Encryption Scheme on Isogenies of Non-Supersingular Edwards Curves</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Anatoly Bessalov</string-name>
          <email>a.bessalov@kubg.edu.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Serhii Abramov</string-name>
          <email>s.abramov.asp@kubg.edu.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Volodymyr Sokolov</string-name>
          <email>v.sokolov@kubg.edu.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Pavlo Skladannyi</string-name>
          <email>p.skladannyi@kubg.edu.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Oleksii Zhyltsov</string-name>
          <email>o.zhyltsov@kubg.edu.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Borys Grinchenko Kyiv University</institution>
          ,
          <addr-line>18/2 Bulvarno-Kudriavska str., Kyiv, 04053</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>CEUR Workshop Proceedings</institution>
          ,
          <addr-line>CEUR-WS.org</addr-line>
        </aff>
      </contrib-group>
      <fpage>12</fpage>
      <lpage>25</lpage>
      <abstract>
        <p>A multifunctional cryptosystem RCNIE on isogenies of non-supersingular Edwards curves is proposed, which solves the problems of Diffie-Hellman secret sharing, digital signature, and public key encryption. The problems of choosing the parameters of non-supersingular Edwards curves forming pairs of quadratic twist with orders p + 1 ± t ≡ 0mod8 over a prime field Fp are considered. Encryption algorithms with mutual authentication of Alice and Bob based on the sharing of their secrets are given, while the length of the key and the size of the digital signature are minimally short and do not exceed the size of the field Fp element. An illustration is given of the operation of the cryptosystem model on 4 degrees of isogenies {3,5,7,37} over the field F863 for a pair of quadratic twist curves with orders 840 and 888. It is shown that for non-supersingular curves there are main and dual cryptosystems, each of which has also an isomorphic cryptosystem. This allows you to perform parallel computing and speed up algorithms. A comparative evaluation of the arithmetic and properties of CSIDH and RCNIE is given. It is noted that we have not found strong arguments for the slow implementation of the CRS scheme in comparison with CSIDH. Taking into account the peculiarities of each of them, both schemes are certainly promising. Curve in generalized Edwards form, complete curve, twisted curve, quadratic curve, non-supersingular curve, curve order, point order, isomorphism, isogeny.</p>
      </abstract>
      <kwd-group>
        <kwd>Edwards</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>transitions,
associated
arithmetic [5–7].</p>
      <p>The presentation [1] of the authors of the PQC</p>
      <sec id="sec-1-1">
        <title>CSIDH</title>
        <p>algorithm
names the</p>
        <p>CRS
scheme
(Couveignes-Rostovtsev-Stolbunov) known since
the beginning of the 21st century as the first
proposed
scheme
on
isogenies
of
nonsupersingular elliptic curves [2–4]. Its remarkable
properties are the commutativity of isogenic
the
flexibility,
and
simplicity
with the use of prime field</p>
        <p>Further, the appearance of the PQC SIDH
(2011) and CSIDH (2018) algorithms already
uses the technology of supersingular elliptic
curves over the fields   2 and   , respectively,
which is justified
by the
relatively
faster
implementation of the algorithms [8]. In [1] it is
noted that CRS encryption is unacceptably slow
and can take several minutes at a security level of
128 bits.</p>
        <p>Being engaged in recent years in the problems
of modeling and modifying CSIDH [9–12], we
became interested in the reasons for the above
fact, which is not substantiated by anything. One
of the goals of this paper is to try to compare the
difficulties that hinder the execution of the
CSIDH algorithm and our simple CRS-like model
on
non-cyclic
curves.</p>
        <p>The
main
difference from CRS [2] in this work is, as in
CSIDH, the use of pairs of quadratic twist curves,
which
uses
the
property
of
isogeny
bidirectionality. It should be noted that the set of
non-supersingular elliptic curves is wider than the
corresponding set of supersingular curves with a
rough estimate by a factor of √ , and, as a result,
contains many potentialities. We managed to find
some of them, which are discussed in this article.</p>
        <p>In [1] and early implementations of CSIDH,
supersingular curves in the Montgomery form [1]
were used, but soon in [13] using the (W:
Z)coordinates [14] of curves in the Edwards form, a
gain of 20% was obtained in comparison with [1]
in the speed of calculations. Further, generalizing
the formulas for calculating isogenies for
Edwards curves [15] to twisted Edwards curves in
[9], we illustrated the implementation of CSIDH
models on non-cyclic quadratic and twisted
Edwards curves [10–12]. The last curves were
first defined in the fundamental work [16], but
with unfortunate terminology, so we use the
classification of curves in the Edwards form [17,
18]. An analysis of the properties of supersingular
Edwards curves of all classes is given in [19, 20].
In this paper, we use non-supersingular Edwards
curves of two classes with the same constraint p ≡
3mod4. This allows one to express the equations
of curves of a quadratic twist pair using additively
inverse parameters.</p>
        <p>A well-known problem of the CSIDH
algorithm is the vulnerability to a side channel
attack based on the measurement of the
computation time of the chain of isogenies of each
degree, which is proportional to the secret
exponent е of the key. In a large number of
articles [21–24], the solution to this problem is
proposed by increasing the exponents by fictitious
ones to a known maximum (Constant time
CSIDH). It is clear that such redundancy
significantly reduces the speed of the algorithm.
In [11], we proposed randomization of the CSIDH
algorithm as a method to counter side-channel
attacks. It is also used in this article.</p>
        <p>The key encapsulation algorithms [12, 25–27]
are now preferred to the classic Diffie-Hellman
secret sharing scheme. Such an algorithm can also
be constructed using the CRS scheme.</p>
        <p>The order of an elliptic curve E over a prime
field   is defined as # =  + 1 −  , where t is
the trace of the Frobenius endomorphism
equation. For a quadratic twist curve   ,
respectively, this order #  =  + 1 +  is
symmetric concerning the mean value  + 1. For
a supersingular curve,  = 0 and the orders of
both curves  + 1 coincide, and the sets of
isogeny degrees are the same, but the signs of the
exponents of the degrees are inverse to each other,
as in CSIDH. In the case of non-supersingular
curves, the orders of quadratic twist pairs differ by
2t, then there are different degrees of isogenies on
curves of two classes related as quadratic twist
pairs with different orders. This is the main
specificity of non-supersingular curves. The
exponents of the isogeny degrees of these two
curves, as in CSIDH, have opposite signs. The
alternation of degrees of isogenies according to
the randomization method is random, and the
simplicity of transitions of the chain of isogenies
from one class of curves to another is achieved by
the fact that their parameters are additively
inverse: ( ,  ) ↔ (− , − ).</p>
        <p>In Section 1, a brief review of the properties of
quadratic and twisted Edwards curves [18] is
given, and methods and options for choosing the
parameters of a pair of these curves as a pair of
quadratic twists are proposed. Section 2 presents
algorithms of the Randomized Commutative
Nonsupersingular Isogeny Encryption (RCNIE)
scheme on quadratic and twisted
nonsupersingular Edwards curves (NECs). In
Section 3 an illustration of the computation of
isogenic chains in the corresponding algorithms
on a model with degrees {3,5,7,37} over the field
 863 is given. Here we discuss the existence of
isomorphic and dual cryptosystems on isogenies
of non-supersingular curves, expanding their
possibilities and applications. Finally, Section 4
compares the CSIDH algorithm and the RCNIE
scheme.</p>
      </sec>
    </sec>
    <sec id="sec-2">
      <title>2. Choice of Non-Supersingular Non</title>
    </sec>
    <sec id="sec-3">
      <title>Cyclic Edwards Curves Parameters</title>
      <p>
        The elliptic curve   , the equation in the
generalized Edwards form [17] with two
parameters a and d is written as
  , :  2 +   2 = 1 +   2 2,
 ,  ∈   ∗,   ≠  ,   ≠ 1.
(
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        )
      </p>
      <p>For the first time, such a curve was proposed
in the fundamental work [12] with the coefficient
a at  2 and the term “twisted Edwards curve.” For
the correct division of curves in the Edwards form
into non-intersecting classes, we use our
classification [17].</p>
      <p>
        If the quadratic character is  ( ) = −1, then
curve (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) is isomorphic to the complete Edwards
curve [16] with one parameter d,  ( ) = −1,
 = 1. Another case χ(ad) = 1 generates 2 classes
of non-cyclic curves: quadratic and twisted
Edwards curves. In particular, if   ( ) =  ( ) =
1, curve (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) is isomorphic to the quadratic
Edwards curve [17] with one parameter
  :  2 +  2 = 1 +   2 2,
 ( ) = 1,   ≠ 1.
(
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        )
      </p>
      <p>
        The twisted Edwards curve is defined in [17]
as a special case of the curve (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) with   ( ) =
 ( ) = −1. The introduction of the second
parameter  into equation (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) in [16] is necessary
only for these conditions.
      </p>
      <p>
        In [16], curve (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) together with the complete
curve are called Edwards curves. At the same
time, their properties and structure are cardinally
different [17, 18]. The controversial terminology
in [16] sometimes leads to misunderstandings and
errors in scientific articles, which is discussed in
[10]. In this paper, in particular, the following
theorem is proved.
      </p>
      <p>
        Theorem 2 [10]. For a curve   ,
(
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) in
generalized
      </p>
      <p>Edwards
form
  2 2 defined over a prime field   , there exists
 2 +   2 = 1 +
a unique quadratic twist curve   ⃗ ,
parameters  ̄ =  ,    ̄ =  ,   ∈   ∗.</p>
      <p>2</p>
      <p>
        Its proof is given in [10]. From it, in particular,
it follows that in the class of complete Edwards
curves, the quadratic twist curve    =   −1 lies
inside this class, while for the quadratic curve (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ),
quadratic twist gives a twisted curve    , =
   , , ( ) = −1. Each of the 3 classes contains
equal sets  −3 curves (d ≠ 0, ±1). Then the
replacement of the class of complete Edwards
curves by 2 classes of non-cyclic Edwards curves
doubles the space of pairs of quadratic twist
curves in the CSIDH algorithm.
      </p>
      <p>We define a quadratic and twisted Edwards
curve as a pair of quadratic twists with parameters
 (</p>
      <p>) = 1,  ̅ =  ,    ̅=  ,   ( ) = −1. Since
we accept the condition  ≡ 3 mod4, we can
accept  = −1,   ̅ = − = −1,   ̅= − , where
 = 1,</p>
      <p>
        are the parameters of the quadratic
curve (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ), respectively  ̅,  ̅are twisted curves. In
other words, the transition from
quadratic to
twisted curve and vice versa can be defined as
  =  1,
      </p>
      <p>
        ↔  −1,− . Then the twisted Edwards
curve equation (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) can be written as
 −1,− :  2 −  2 = 1 −   2 2,
 ∈   ∗,   ≠ 1,   ( ) = 1.
      </p>
      <p>
        The orders of quadratic (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) and twisted (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
Edwards curves are comparable to 0mod8, then
p  − 1mod8 [17]. Note that equation (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ) has a
fixed parameter
      </p>
      <p>
        = −1 after which all curves (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
are determined by one parameter (− ). Quadratic
residues of parameters ( = 1 and  ) of curve (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        )
become quadratic non-residues (
= −1, −  ) of
curve (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ). This property of additively inverse
parameters simplifies the illustration
operation of cryptosystems on the isogenies of
these curves. Equation (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) can be written in a
different form with the change d→ad and fixing a.
      </p>
      <p>
        It is interesting that in equation (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ), as well as
for the
parameter, one can single out the parameter d
 = ( −2 +  −2 −  −2 −2).
      </p>
      <p>This
means that using the
multiplicative
inversion of the coordinates of eight known points
of the curve P = (±x, ±y), Q = (±y, ±x) it is easy to
calculate the unique parameter d of the curve. This
is useful in cryptanalysis and should be taken into
account in algorithms to protect against such an</p>
      <p>By analogy with CSIDH, it is easy to form the
general parameters of a CRS-like cryptosystem on
isogenies of Non-supersingular Edwards Curves

(NECs). Let  0 = ∏</p>
      <p>=1   and N = 8 0 be the
order of a quadratic supersingular Edwards curve
over a field with modulus  0 = 
− 1. By setting
the values of the Frobenius trace  = ±8 , 
1,2,3, …</p>
      <p>
        we determine the sum  0 ± 8
equal to the prime number p. Then over the field

 there exists a quadratic NEC (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) of order #Ed =
8 0 and a twisted curve (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ) of order # −1,−
=
= 
= 
± 16
      </p>
      <p>= 8 1.</p>
      <p>Example 1. For a set of isogeny degrees
{  } = {3,5,7},  0 = 105, N = 840,  0 = 839,
for m = 3 we get a prime number p = 839 + 24 =
863. Then the orders of the curves of the quadratic
twist pair are #Ed =840 = 8 ∙ 3 ∙ 5 ∙ 7 and # −1,−
=</p>
      <p>+ 48 = 888 = 8 ∙ 3 ∙ 37,  1 = 111.</p>
      <p>Another method is possible, not always lead to
success. You can set two sets of degrees that differ
by one or more elements and calculate the
arithmetic mean of the products of the elements of
these sets u. If 8u−1= p is a prime number, it is
the desired modulus of p.</p>
      <p>Example 2. Let {  } = {5,7},  0 = 35, 8 0 =
380. In the second set {  }
 = {3,11},  1 =
33, 8 1 = 364. The mean of orders is 8u = 372,
but 372 – 1 = 371 is not a prime number.</p>
      <p>However, 16 − 1 = 743 is a prime number, and

with</p>
      <p>
        attack.
(
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
over the field  743 one can obtain, with double
redundancy, pairs of quadratic twist curves with
orders 760 and 728, containing isogeny kernels of
degrees {5,7} and {3,11}.
      </p>
      <p>Example 3.</p>
      <p>Let
{  } = {3,5,7,11},  0 =
1155, 8 0 = 9240. In the second set {  }
{3,5,7,13} we change element 11 of the first sets
by 13 and calculate  1 = 1365, 8 1 = 10920.
The average value of the orders of the two curves
 =
is 10080 =  + 1, since p = 10079 is a prime
number. Here the trace of the Frobenius equation
is  = 840, and the orders of the curves of the
quadratic twist pair are #Ed = 9240, # −1,− =
10920. Note that for p = 9239, it is possible to
build CSIDH, for p = 10079 only CRS, and the
number 10919 is not prime.</p>
      <p>An approach can be considered rational if the
sets of isogenies degrees of pair curves of
quadratic twist intersect as much as possible. For
example, in
the
second
curve
with
the
preservation of all lower degrees, the highest
degree</p>
      <p>1 can be replaced by a higher one
 max2, then the orders of the two curves differ by
a factor  max1/ max2. This difference cannot
exceed
boundaries.</p>
      <p>the
allowable</p>
      <sec id="sec-3-1">
        <title>For real</title>
        <p>limits
2√ differ for large p with an estimate of 1 − 2
.</p>
        <p>Meanwhile, for CSIDH, for example, the ratio of
the maximum orders  max1/ max2 differs little
from 0.99, which is less than 1 on about 1%. In
the alternative case, for non-intersecting two sets
of degrees, there are also no such sets that would
give
practically
equal products. It is
more
pragmatic to specify the first non-supersingular
curve with minimal isogeny degrees, followed by
an appropriate factorization of the order of the
twist curve (Example 1). It will only add new
degrees of isogenies with a smaller number of
them. The number of isogenies degrees in both
curves can be made approximately equal by
selecting from a table of primes with alternating
numbers for the first and second curves (as in</p>
      </sec>
      <sec id="sec-3-2">
        <title>Example</title>
        <p>2). The
problem
of choosing the
parameters of quadratic twist pairs of
nonsupersingular curves still needs to be studied.</p>
      </sec>
      <sec id="sec-3-3">
        <title>Quadratic and twisted</title>
        <p>NEC as a pair of
quadratic
twists
have
different
orders
and
different structures. Except for the two points (0,
±1), all their points are different. Both curves are
non-cyclic
concerning
points of even
order
(contain 3 points of the 2nd order each, two of

√
which are singular points  1,2 = (±√
[17]). Quadratic NEC, in addition, contains 2
singular points of the 4th order ± 1 = (∞, ± 1
).</p>
        <p>The presence of 3 points of the 2nd order limits
the number 8 to the minimum even cofactor of the
order</p>
        <p>Ord E = 8n (n is odd) of twisted
and
quadratic Edwards curves [17]. The maximum
order of the points of these curves is Ord E/2.
Points of even orders mustn't be involved in the
calculation of isogenies of odd degrees (the first
multiplication by four of a random point P gives a
random point R of odd order n or a divisor of n).</p>
        <p>The choice of two classes of non-cyclic NECs
for cryptosystems on isogenies is justified by their
advantages over complete NECs:</p>
        <p>1. The number of all quadratic and twisted
Edwards curves ( − 3) is twice the number  −3
of
all
complete
corresponding proportion is also valid for the
number of isogenic NECs and, therefore, the
2
the
security of the cryptosystem.</p>
        <p>2. The transition to the curve of quadratic twist
 
↔  −1,−</p>
        <p>does not require the laborious
inversion of the parameter  ↔  −1, which is
necessary for a complete NEC.</p>
        <p>Along with isogenic curves (with different
Jinvariants), there are isomorphic curves
with
equal J-invariants, which are defined [16, 28]
 ( ,  ) =
 ( −  ) ≠ 0.</p>
        <p>16( 2 +  2 + 14 )</p>
        <p>
          3
 ( −  )4
,
(
          <xref ref-type="bibr" rid="ref12 ref4 ref8">4</xref>
          )
This
        </p>
        <p>parameter, in particular, recognizes
isomorphic curves with different values of the
parameter d. As a result of calculations of chains
of isogenies, one usually
makes the change
 →  ( ). This parameter is also used in the
ElGamal encryption scheme. For quadratic and
twisted</p>
      </sec>
      <sec id="sec-3-4">
        <title>Edwards</title>
        <p>curves  ( ) =  ( −1), i.e.
inverting the parameters d gives an isomorphic
curve.</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>3. Algorithms of the RCNIE Scheme</title>
      <p>on Quadratic and Twisted Edwards</p>
    </sec>
    <sec id="sec-5">
      <title>Curves</title>
      <p>
        Instead of supersingular ones, we start here
with non-supersingular Edwards curves (NEСs)
of two classes (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) and (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ), connected as pairs of
quadratic twist with orders  + 1 ±  ≡ 0mod8.
the
application
of this
function for
nonsupersingular curves is the cardinal complication
of the multiplicative inversion operation for some
of the isogeny degrees that are different in the sets
{  }and {  } . For our tasks, this does not violate
the efficiency of encryption algorithms, while at
the
same
time
complicating
the
tasks
of
cryptanalysis. We will return to this problem in
      </p>
      <sec id="sec-5-1">
        <title>Section 3.</title>
        <p>In all algorithms, the CGA encryption function
( ) encrypts the secret key 
using the mapping  = [ 1 1,  2 2, . . ,  
= ( 1,  2, . . ,   )
  ] and
the starting curve  0 into an isogenic curve  ′ =
 0 ∗ , whose parameter is taken either as the
corresponding public key or as a new short secret
key. The parameters of the quadratic NEC #  (0)=

 + 1 −  = 2</p>
        <p>∏ =1   ,  
field modulus  ≡  − 1 
≥ 3 and the prime</p>
      </sec>
      <sec id="sec-5-2">
        <title>8 are given. The</title>
        <p>
          (
          <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
          ) [11] is given below.
# −(01),− =  + 1 +  and the factorization of
degrees {  } of twisted NEC are determined. A
randomized algorithm for calculating
        </p>
      </sec>
      <sec id="sec-5-3">
        <title>Alice’s</title>
        <p>public
key  А</p>
        <p>
          using the secret key   =
( 1,  2, . . .   ) on the isogenies of curves (
          <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
          ) and
        </p>
        <sec id="sec-5-3-1">
          <title>Randomized algorithm 1: Evaluating encryption function on quadratic and twisted NEC</title>
          <p>Input:   ∈   ,  ( ) = 1 and a secret key   =
( 1,  2, . . .   ).
  , where   , :</p>
          <p>2 +  2 = 1 +  А,В 2 2, 
Output:   such that [ 1 1,  2 2, . . .     ] ∗   =
Let  0 = { |  &gt; 0} ,  1 = { |  &lt; 0},  0 =
1.
2.</p>
          <p>∏ ∈ 0   ,</p>
          <p>,  1 = ∏ ∈ 1   ,
While some   ≠ 0 do
3. Sample a random  ∈   ,
A nonzero value of  by a factor of √ expands
the set of curves and offers interesting new
applications.</p>
          <p>In contrast to the CRS scheme [2] on curves in
the Weierstrass form with two parameters, which
does not use pairs of quadratic twist, we build
speeding up our algorithms are the rejection of the
very laborious calculation of the isogenic function
 ( ) of point  and randomization of algorithms
[10].</p>
          <p>In this paper, we propose a PQC RCNIE
scheme. It is distinguished from the known ones
by the existence of four parallel cryptosystems
(with the addition of dual and two isomorphic
ones) and, most importantly, multi-functionality.</p>
          <p>The PQC CSIDH [1] algorithm is based on the
CGA (class group action) encryption function
over a prime field   . The CGA function defines
an isogenic mapping  of a supersingular elliptic
curve E of order #E = p + 1 into a curve  ′ =  ∗

degrees of isogenies and   are isogeny exponents
(number of isogenic transitions). The sign of the

exponent is   &gt; 0 for the original curve, and
 &lt; 0 for the quadratic twist curve. The mapping
 is commutative and is equally valid for all
elliptic curves over the field   regardless of their
order.
exponent</p>
          <p>{  }
exponents 


[    ]*[</p>
          <p>A pair of non-supersingular quadratic twist
curves have different orders and different sets of
degrees {  } and {  } , which
intersect. For the intersection {  } ∩} {  } , as in
CSIDH, each l-isogeny has both signs and then
may partially
−  ] = 1. This means that chains of
lisogenies of different signs are built in reverse
order and cancel each other out. Therefore, the
keys in CSIDH give for each degree 

an
of only
one
sign.</p>
          <p>
            For
nonsupersingular curves, the union of sets {  }∪
is constructed, and the signs
are determined by whether the
isogeny belongs to one of the curves of the
quadratic twist pair. One should strive for the
equiprobable use of both curves, which is true for
sets {  }and {  } of equal power. Since isogeny
curves of corresponding degrees exist in the
classes of curves (
            <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
            ) and (
            <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
            ) connected as pairs of
quadratic twist, to construct commutative chains
of isogenies we use the CGA encryption function
 = [ 1 1,  2 2, . . ,  
  ] [1]. A specific feature of
 −  ,
14. Return  А.
          </p>
          <p>5. Else a  −1,   1
   2 2 If  (( 2 − 1)/(  2 − 1) = 1,
  :  2 −  2 = 1 −
6. Compute  -coordinate of the point  =
   2 2,
( ,  ) ∈   ,
7. Compute 
13. Skip  in   and set   ← (  /  ) If   =
This algorithm has important differences from
the original algorithm 2 [1], which are discussed
in [11]. In addition to modifications related to the
randomization method of the CSIDH algorithm,
here we refuse the redundant isogenic function
φ(R) of a random point R, which radically speeds
up the algorithm.</p>
          <p>
            The idea of randomization is that for any
random value of the variable x (xy ≠ 0, ), the
point P = (x,y) with a known parameter d always
belongs to one of the two curves (
            <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
            ) or (
            <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
            ). This
makes it possible to double the speed ap selection
of a random point P and complicate the side
channel attack. This is also facilitated by a random
choice of isogeny degrees (as they are exhausted).
Also, if P is not well chosen, moving to a new
degree often fixes the problem faster than varying x.
          </p>
          <p>At the beginning of Algorithm 1, two subsets
  ,  = 0,1, with degree numbers   , are formed,
together with two factors  0 and  1 of the number
 =  0 1: the index λ=0 (  &gt; 0) corresponds to
the choice of a quadratic NEC, and λ = 1 is twisted
NEC (  &lt; 0). Since the order of the curve is
#  = 8 0, then in line 7 of the algorithm for the
curve   the point</p>
          <p>= 4 1 of odd order  0 is
calculated, and the curve  −1,− the point 
=
4 0
minimizes
of
odd
the
order  1is</p>
          <p>calculated.
cost
of
the
next</p>
          <p>This
scalar
multiplication, which determines the point Q of
the isogeny kernel of the degree  
Further, in line 10 of the algorithm, by doubling
points of the kernel &lt;Q&gt; are calculated.
the points,  = (  − 1)/2 x-coordinates of the
(line 9).</p>
          <p>In line 7 of Algorithm 1, double doubling the
random point P immediately allows you to get rid
of points of an even order (including special
points of the 2nd and 4th order) and then the
calculation of scalar multiplications in subgroups
of points of an odd order of the curve. Their task
of x-coordinates   of the kernel
is to find (  −1)</p>
          <p>
            2
according to the formula [15]
points &lt;Q&gt; of prime order   . As a result,

 ′ =    8,  = ∏ =1   , s = (
 − 1)/2
(
            <xref ref-type="bibr" rid="ref13">5</xref>
            )
the parameter  ′ of the   -isogenic quadratic NEC
is calculated. Twisted NEC parameters (
            <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
            ) а′ =
−1,  ′→ −  ′ . We emphasize that the concept of
RCNIE is the construction of chains of isogenic
curves as
          </p>
          <p>Abelian</p>
          <p>groups, and not isogenic
functions  ( ) of a random point  . The
laborintensive calculations of the latter in [1] are
redundant.
3.1.</p>
        </sec>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>Diffie-Hellman Non-Interactive Secret Sharing Algorithm</title>
      <p>compute ∏
modulus 
the elliptic curve   .</p>
      <p>1. Choice of parameters. For odd primes   ,

 =1   , choose an appropriate field</p>
      <p>= 2</p>
      <p>∏ =1   − 1,  
≥ 3 and start
2. Calculation of public keys. Alice and Bob
use secret keys in the form of vectors   ,В =
( 1,  2, . . ,   ) construct isogenic maps   ,В =
[ 1 1,  2 2, . . ,     ] and calculate the isogenic
curves   ,В =   ,В ∗  0 as their public keys.
These curves are determined by their parameters
up to isomorphism.</p>
      <p>3. Key exchange. Here the protocol is similar
to item 2 with the replacement  0 →   for Alice
and  0 →   for Bob. Knowing Bob's public key,
Alice calculates  
=   ∗  
=   ∗   ∗  0.</p>
      <p>Similar actions Bob gives the result  
  =   ∗   ∗  0 which coincides with the first
one</p>
      <p>due to the commutativity of the group
operation.</p>
      <p>The</p>
      <p>J-invariant
of
the</p>
      <p>curve
=   ∗
  В ( ВА) is taken as a shared secret [1].
3.2.</p>
    </sec>
    <sec id="sec-7">
      <title>Digital Signature Algorithm</title>
      <p>The tasks of the digital signature are to
authenticate the sender of message M and to verify
the integrity of the transmitted message by the
recipient. Alice usually uses her private key for
this, and Bob uses her public key for verification.
In the previous secret-sharing problem, both of
these functions are performed: Alice encrypts
Bob's public key</p>
      <p>with the secret key   , and
computes the shared secret  В . Bob uses Alice's
public key  А and his secret  В, to calculate the
same curve   В. The fact that   В =  В
that the task of authenticating Alice by Bob has
means
been completed.</p>
      <p>When checking, it remains for Bob to make
sure that together with the equality   В =  В for
the sent M and received M' messages, the hash
codes are identical: H(M') = H(M). Both qualities
can be easily combined by concatenation into one
with the secret of the first of them preserved. The
most secure from a quantum computer for this is
hashing the combined equalities above. These
symmetric procedures are included in the digital
signature algorithm below.</p>
      <sec id="sec-7-1">
        <title>Precomputation</title>
        <p>Based on each other’s public keys   and  В
and their private keys   ,В = ( 1,  2, . . ,   ) Alice
and Bob perform calculations of the previous
Diffie-Hellman secret-sharing scheme and find
curves   А=</p>
        <p>=  . Next, Alice respectively
forms, and Bob verifies the digital signature DS:
A. DS formation
1. Calculation h = H(M) &lt; p.
2. Calculation J( ) = J(  А).
3. Calculation DS =H[J(κ)|| h] &lt; p.
4. Dispatch (M, DS)→ B.</p>
        <p>B. DS Verification
1. Calculation h' = H(M').
2. Calculation J( ) = J(  ).
3. Calculation 
4. Checking</p>
        <p>′ = H[J( )|| ℎ′] &lt;  .
′ =</p>
        <p>. In the case of  ′ ≠  ,
the signature is incorrect.</p>
        <p>It is a non-standard signature in the sense that
in the asymmetric scheme, it is generated by
Alice's private key and verified by Bob with her
public key, and Bob's corresponding keys are not
required. This is because a non-interactive
DiffieHellman
scheme
was
adopted
as a
basis,
equalizing the conditions for Alice and Bob: both
parties use their secret keys and the public keys of
the other party to form a shared secret. Such
conditions mean that they work in the scheme of
mutual trust of symmetric cryptography, although
on the technology of asymmetric cryptography.
Perhaps this is permissible since Bob can write to
Alice with the same rights, and she is obliged to
believe him. To do this, they created a shared
secret. The main thing is that a multifunctional
cryptosystem works on this basis.</p>
        <p>This original signature does not exceed the size
of a prime field module and is rather concise. The
signature algorithm is closer to RSA and half as
long as the ElGamal signature. Here, only hashing
and field operations are performed on both sides,
with element inversions as the most complex
operations. The algorithm can equally rely on
CSIDH. The level of quantum security of finding
the secret curve   in the Diffie-Hellman secret
sharing problem is estimated as 4√ . Calculations
in the procedures for generating and verifying a
digital signature have the maximal security level
for a quantum computer.
3.3.</p>
      </sec>
    </sec>
    <sec id="sec-8">
      <title>ElGamal Encryption Algorithm</title>
      <p>This asymmetric cryptography algorithm was
proposed in [2] without using quadratic twist
curves. Interestingly, it is also based on a
secretsharing scheme. In the staging part, of the public
keys, only the recipient's key is—the  В curve.
Alice, on the other hand, calculates her public key
  , in the encryption process and transmits each
session to Bob along with the ciphertext. In
essence, this means that she uses a one-time secret
  ,, which forms her one-time public key   ,. The
combination of Alice's one-time keys   ,,   ,, and
Bob’s long-term keys   ,,   , allows both parties
to compute the one-time shared secret   , =   В,
It is used by Alice to encrypt the short message</p>
      <p>, and by Bob to decrypt it. We describe the
algorithm from [2] with the notations adopted
here.
• 
• 
•
•</p>
      <p>isogenies.</p>
      <p>General system parameters:
 is a prime field.
•   ,  = 1,2, . . ,  are isogeny degrees.
is integers
m
≤  
≤m-exponents of
m is the boundary value of the exponent.
•  0 is starting the elliptic curve and its equation.
•  В is Bob’s public key.</p>
      <p>is plain text.
and</p>
      <p>generating
3. Calculation of the shared Diffie-Hellman secret
4. Calculation of the J-invariant of the curve   ,:
5. Calculation of the ciphertext S = (M⋅   ) mod p.
6. Sending (  ,,  ) to Bob.</p>
      <p>B. Decryption
1. Based on the secret key  В, = ( 1,  2, . . ,   )
function form 

= [ 1 1 ,  2 2, . . ,  
  ].</p>
      <p>2. Calculation of the shared Diffie-Hellman secret
  В, =   , ∗ </p>
      <p>=   .
  =  (</p>
      <p>,).
3. Calculation of the J-invariant of the curve   ,:
4. Calculation of the plain text M = (S/  ) mod p.</p>
      <p>The use of one-time secrets in such a scheme
makes it more secure than a non-interactive</p>
      <sec id="sec-8-1">
        <title>Diffie-Hellman key exchange.</title>
        <p>All three of the above algorithms solve the
main problems of asymmetric cryptography. The
last two algorithms include the first one, which
becomes the base one. As a result, we can state
that in this paper we propose a multifunctional
cryptosystem PQC RCNIE on isogenies of
nonsupersingular
evaluations of its properties are discussed in the
next section.</p>
      </sec>
    </sec>
    <sec id="sec-9">
      <title>4. Modeling RCNIE</title>
      <p>As in the works [2, 29, 30], whose authors
concluded
the
theoretical
provisions
of</p>
      <p>a
cryptosystem on isogenies of elliptic curves,
resistant to quantum attacks, with examples of
encryption
of
messages on
odd</p>
      <p>degrees of
isogenies, in our works [9–12] we use modeling
as a way to easily illustrate the properties of
algorithms. Understanding these properties opens
the way for something new and better. Sometimes
the new is the well-forgotten old, as is the case
with curves in the Edwards form.</p>
      <p>
        Based on the data of Example 1 of Section 1,
we obtain four degrees of isogenies {  } =
{3,5,7,37}, the first three of which are factors of
order 840 of the quadratic curve (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ), and degrees
3 and 37 divide order 888 of the twisted curve (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
over the field  863 and the trace of the Frobenius
equation t =24. For the first curve (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) the signs of
the isogeny exponents   &gt; 0, and for curve (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
      </p>
      <p>&lt; 0. Here degree 3 is bidirectional (allows
both signs), and degrees 5 and 7 (  &gt; 0) and 37
(  &lt; 0) are unidirectional. Below we discuss
these
features
of
the</p>
      <p>isogenies of
nonsupersingular curves.</p>
      <p>
        With a relatively small field modulus p = 863,
it is not difficult to find estimated √ of the
parameters d of all curves (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) with order 840.
      </p>
      <p>Since they are squares, then a full enumeration
modulo p of all с = 2,3, … ,431, and</p>
      <p>
        =  2 gives
the set of all 62 values of the parameters d of the
NEC (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) and (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ), shown in Table 1. All the curves
together, respectively, are 124. In the class of
complete Edwards curves, there would be 62 of
them. Here the number of parameters is even, so
for each curve, there is an isomorphic curve with
the parameter 
(
        <xref ref-type="bibr" rid="ref12 ref4 ref8">4</xref>
        ).
      </p>
      <p>For
example,
↔  −1 and the same J-invariant</p>
      <p>
        169−1 = 623,  (169) =
 (623) = 826. Then there are 31 non-isomorphic
curves (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ), and the same number of curves (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ).
      </p>
      <p>Isogenies of all degrees have a prime period
An array of values of 62 quadratic and twisted NEC parameters d at p = 863, #E = 840, #Et=888 (t = 24)
169
19
71
300
688
400
779
616
113
843
729
111
618
311
339
161
308
444
858
623
818
253
302
673
706
210
116
192
725
281
436
705
486
589
181
309
503
318
75
43
32
852
684
665
573
231
551
840
472
728
307</p>
      <p>
        The number above arrow 1 means one step of
the 3-isogeny chain of the quadratic NEC curve
(
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) with the exponent   &gt; 0. Under the value of
the parameter  ( ) we write the degree of isogeny
in brackets. For the twisted curve (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ) with   &lt;
0, there is also a 3-isogeny chain of period π = 31
having a reverse order of alternation of isogenic
curves (the last chain and (
        <xref ref-type="bibr" rid="ref14">6</xref>
        ) are read in reverse
or opposite order). The number above the arrow
(–1) means one step of the isogeny of the curve
(
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ) with negative parameters. Here the remarkable
property of the twofold finding of the
multiplicative inversion of an element of an
isogenic chain arises. On the one hand, it's true:
 (0) ∗  (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) ∗  (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) ∗ … .∗  ( ) =  (0) ⇒
[ (0)] ∗ [ (0)]−1 = 1,
[ (0)]−1 =  (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) ∗  (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) ∗ … .∗  ( ).
      </p>
      <p>
        (
        <xref ref-type="bibr" rid="ref15">7</xref>
        )
      </p>
      <p>
        On the other hand, for bidirectional l-isogenies
of a pair of quadratic twist curves with exponents
±1, we have
 (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        ) ∗  (−1) = 1 ⇒ [ (
        <xref ref-type="bibr" rid="ref1 ref5 ref9">1</xref>
        )]−1 =  (−1)
      </p>
      <p>
        (
        <xref ref-type="bibr" rid="ref16">8</xref>
        )
      </p>
      <p>
        In other words, to multiplicatively reverse one
step of an isogenic chain, it is required in the
general case to find a chain of period π (see (
        <xref ref-type="bibr" rid="ref14">6</xref>
        )
and (
        <xref ref-type="bibr" rid="ref15">7</xref>
        )). The same problem for bidirectional
isogenies is solved in one step instead of π steps
(
        <xref ref-type="bibr" rid="ref16">8</xref>
        ). The opposite signs of the exponents of such
isogenies cancel each other out: [ +1]*[ −1] = 1.
      </p>
      <p>
        The above case takes place in CSIDH, which
distinguishes it favorably from CRS. But in
CSIDH, the isogeny of each degree in the keys is
used as a unidirectional exponent, which is
understandable, since different signs of the
exponent only neutralize each other [31]. Property
(
        <xref ref-type="bibr" rid="ref16">8</xref>
        ) is only useful for protection against
sidechannel attacks [11], and in the CSIKE problem
[12, 32]. Together with non-supersingular curves,
the inversion of a unidirectional isogeny element
according to (
        <xref ref-type="bibr" rid="ref15">7</xref>
        ) requires knowledge of the
isogeny period and computation time, which is
unrealizable for real cryptosystems. In these
cryptosystems, handling tasks should be avoided.
      </p>
      <p>
        Characteristically, for the starting curve
  (0) =  169, the sequence  ( ) (
        <xref ref-type="bibr" rid="ref14">6</xref>
        ) does not
contain the element 169−1 = 623 with the same
J-invariant. It follows that this is also true for all
elements of this sequence of period 31, any of
which can be taken as the starting element with
the corresponding cyclic shift (as in a cyclic code
of length 31). All J-invariants of parameters (
        <xref ref-type="bibr" rid="ref14">6</xref>
        )
are different. Isogenic curves of other degrees 5,
7, and 11 contain the same parameters and the
same period as (
        <xref ref-type="bibr" rid="ref14">6</xref>
        ), alternating in a different order.
      </p>
      <p>
        Next, we will see that the same operating
parameters contain all the calculations in the
secret sharing scheme. If you invert the starting
curve 169−1 → 623, you do not need to build
new isogenic chains, it is enough to invert the
results. In this case, the other half of the
parameters of Table 1, not included in (
        <xref ref-type="bibr" rid="ref14">6</xref>
        ), will be
working. Thus, there are two isomorphic
cryptosystems with different mutually inverse
parameters d and coinciding sets of J-invariants.
      </p>
      <p>If the starting curve is given and does not change
all parameters  ( ) of isogenic chains are unique
and there is no need to pass to the J-invariant of
the resulting curve. An isomorphic cryptosystem
can solve other problems in parallel, which
doubles the performance of such a system.</p>
      <p>Further, we will see that in addition to the
isomorphic cryptosystem, there is also a dual
cryptosystem, which also has its isomorphic one.</p>
      <p>Overall, there is the potential to quadruple the
performance of the RCNIE scheme.</p>
    </sec>
    <sec id="sec-10">
      <title>4.1. Implementation of the Diffie</title>
    </sec>
    <sec id="sec-11">
      <title>Hellman Secret Sharing Algorithm</title>
      <p>
        In our model with isogenies of degrees
{3,5,7,37}, to equalize the probabilities of
choosing the curves of a pair of quadratic torsion,
we will take all the degrees to be unidirectional,
then in the secret keys the degrees {5,7} will be
assigned to the quadratic curve (  &gt; 0), and the
degree {3,37} to twisted (  &lt; 0). Let’s take the
secret keys of Alice   = (
        <xref ref-type="bibr" rid="ref1 ref13 ref5 ref9">−2,5, 1, −4</xref>
        ) and Bob
  = (
        <xref ref-type="bibr" rid="ref11 ref11 ref3 ref3 ref7 ref7">−1, 3, 3, −5</xref>
        ). Let's compute each of their
public keys in 12 randomly chosen isogeny steps.
      </p>
      <p>
        Alice’s public key with a random choice of
curves and degrees is defined as
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
Bob’s symmetric calculation with   = (
        <xref ref-type="bibr" rid="ref11 ref11 ref3 ref3 ref7 ref7">−1, 3, 3, −5</xref>
        ):
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref13">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref15">7</xref>
        )
(
        <xref ref-type="bibr" rid="ref15">7</xref>
        )
give the same result due to the commutativity of
the isogenies   =   = 161 which
determines the quadratic curve   of the shared
secret. As noted above, this value is unique (for a
given starting curve) and here it is not required to
pass to the J-invariant in the shared secret  =
161.
      </p>
    </sec>
    <sec id="sec-12">
      <title>4.2. Implementation of the Digital</title>
    </sec>
    <sec id="sec-13">
      <title>Signature Algorithm</title>
      <p>The previous problem is solved at the
precomputation stage and also performs mutual
authentication of Alice and Bob. They calculated
the shared secret key  =161.</p>
      <p>In conclusion of this section, we note that the
example of the implementation of the encryption
algorithm here is redundant since with other
conditions it is given by the CRS co-authors in [2].
We will only comment on this example.</p>
      <p>In [2], an example was constructed for 6
unidirectional isogenies of degrees
{3,5,7,11,13,17} with their product 255255 over a
prime field with modulus p = 2038074743 and the
order #E = 2038078635 of the curve in the
Weierstrass form with two parameters A and B.
There is a typo in the last number in the article
since it is not divisible by all degrees except 3 and
5. It is not clear why the model required a
redundancy of approximately 104 times for the
field modulus and curve order. In the
corresponding number of times, the calculations
become more complicated, and both parameters A
and B of the curve increase (these are ten-digit
decimal numbers instead of the necessary
sixdigit ones). In addition, two Weierstrass curve
parameters can now be easily replaced with one
Edwards curve parameter d and double the
computational speed. Perhaps, here and in other
algorithms [29, 30, 33], the reasons for the
slowness of the CRS scheme are associated with
excessive redundancy of parameters, which can
be eliminated. But this is not a reason to consider
the CRS scheme unacceptable. We will return in
the next section to a discussion of this issue.</p>
    </sec>
    <sec id="sec-14">
      <title>5. Comparative Evaluations of Properties CSIDH and RCNIE</title>
      <p>
        In [11], we proposed a randomized CSIDH
model with bidirectional isogenies of degrees
{3.5.7} on curves (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) and (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ) above the field  839.
These parameters are close to the parameters of
non-supersingular curves (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) and (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        ) over the field
 863 with different orders 840 and 888 and a trace
of Frobenius equation t = 24. These two models
are most convenient and correctly compared.
      </p>
      <p>It can already be argued that with the transition
from CSIDH to a non-supersingular curve, one or
more new degrees of isogenies always appear (in
our case,  = 37). Their insignificant disadvantage
is unidirectional isogenies, and for large
cryptosystems, it is practically the inability to
inverse the isogenic chains and such tasks should
be avoided in algorithms. Here we see mutual
advantages and disadvantages.</p>
      <p>If we now turn to problems related to speed,
then we have not found any new reasons
inhibiting the execution of the algorithm. Usually,
when choosing a random point P at the beginning
of each step in calculating the isogenic curve, the
point P may be unsuccessful with a certain degree
  . This means that the order of point P does not
contain a factor   . The probability of such an
event  −1the more, the lower the degree and
reaches the maximum value of 1/3. We do not
recommend taking too small degrees in
cryptosystems, they are the most problematic. In
the described case, randomization allows random
transitions to other degrees of isogenies. From our
experience, unsuccessful random points arise with
the same frequency, regardless of whether the
curve is supersingular or not. As noted in the
previous section, the slowness of the
implementation of the calculations of the
isogenies is most likely associated with the
exorbitant redundancy of the characterization of
the prime field   and the curve order in the
models used. The reason for this redundancy [2]
remains unclear.</p>
      <p>For NEC, there is a unique ability to build not
only a pair of quadratic twists with the orders
p + 1 ± t but also inside each class to find a pair of
curves with the same order as in a quadratic twist
curve. We will call the corresponding curves dual.
Their existence allows you to replace quadratic
curves with twisted ones and vice versa. For
example, the degree of isogenies of l = 37 in our
model belongs to the twisted curves  −1,− of
power 64 and order 888. Over the field  863, there
is a curve   of order 888 with a minimum
parameter d = 6. Calculate for curves  
parameters  ( ) chains of 37-isogenic curves on
the period  = 31:
678 1</p>
      <p>→
(37)
(37)
12
(37)
→</p>
      <p>1</p>
      <p>
        Here we see half of the dual curves  
parameters of order 888. As in (
        <xref ref-type="bibr" rid="ref14">6</xref>
        ), in this
sequence, no element d has the inverse  −1. The
second half of the parameters  ( )) is calculated
by the inversion (for isomorphic curves) of the
above. The corresponding twisted curves have
order 840. The existence of dual curves makes it
possible to build two cryptosystems over the same
field of  863: the main and dual, the signs of the
exponent isogenies of which change places. These
cryptosystems can work independently, and,
therefore, double the number of tasks to be solved.
If you add isomorphic to each of the two
cryptosystems mentioned, four parallel
cryptosystems are formed with different sets of
parameters d that allow parallel independent
calculations. It is still unclear whether there is a
simple (as for pairs of quadratic twist) relationship
between the parameters of the main and dual
curves. This question remains open. In any case,
the existence of dual cryptosystems, unique for
non-supersingular curves, promises a 4-fold
expansion of the capabilities of cryptosystems on
the isogenies of elliptic curves. This prospect
requires further research.
      </p>
      <p>
        The results of the implementation of the
Edwards-CSIDH model [13] in projective
coordinates ( :  ) claim that it is faster than the
Montgomery-CSIDH models in coordinates
( :  ) by 20%. Note that this model in [13] is
built on the complete Edwards curves   with the
order #  = p+1 = 4n (n- odd) and the inversion
of the parameter  ↔  −1 when the transition to
the curve of quadratic twist. [9–12] use the fastest
arithmetic of quadratic (
        <xref ref-type="bibr" rid="ref10 ref2 ref2 ref6 ref6">2</xref>
        ) and twisted curves (
        <xref ref-type="bibr" rid="ref11 ref3 ref7">3</xref>
        )
with additive inversion to the parameters of a pair
of quadratic twists. The main advantage of these
classes of Edwards curves over a prime field   is
the doubling of several curves in the algorithm
with the corresponding increase in safety.
      </p>
    </sec>
    <sec id="sec-15">
      <title>6. Conclusion</title>
      <p>A multifunctional cryptosystem RCNIE on
isogenies of non-supersingular Edwards curves is
proposed, which solves the problems of
DiffieHellman secret sharing, digital signature, and
public key encryption. It is built on two classes of
non-cyclic Edwards curves forming pairs of
quadratic twists.</p>
      <p>The basic RCNIE algorithm is a secret sharing
algorithm that also serves to mutually authenticate
users.</p>
      <p>A model for the execution of
cryptoalgorithms on isogenies of 4 degrees {3,5,7,37} is
constructed and an analysis of its properties is
given. Examples of calculations of curve
parameters in crypto-algorithms using the
randomization method are given.</p>
      <p>The existence of the main and dual
cryptosystems on non-supersingular curves is
illustrated, for each of which there are isomorphic
cryptosystems with inverted parameters. The
possibility of parallel computing algorithms in
these cryptosystems allows you to quadruple the
performance of a complex cryptosystem, or use
some of them for redundancy and updating.</p>
      <p>A comparative evaluation of the arithmetic of
cryptosystems on isogenies of supersingular and
non-supersingular elliptic curves is given. It is
noted that the authors found no reason to consider
the latter technology to be slower than that
adopted in CSIDH. Since the number of all
nonsupersingular curves is estimated to be √ times
greater than the number of supersingular ones, it
is reasonable to use a number of their advantages
mentioned above in future applications.</p>
      <p>We believe that CSIDH and CRS technologies
should not be opposed, but should be developed
as promising, taking into account the features and
advantages of each of them.</p>
      <p>Further studies, it is planned to study new
approaches to the formation of sets of isogeny
degrees in RCNIE, as well as the digital signatures
algorithm.</p>
    </sec>
    <sec id="sec-16">
      <title>7. References</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Calculation h = H(M). Let</surname>
          </string-name>
          h =
          <volume>852</volume>
          &lt; p.
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2. Calculation of J( ) =
          <fpage>583</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Calculation</surname>
            <given-names>DS</given-names>
          </string-name>
          = H(
          <volume>852</volume>
          ||583) =
          <fpage>796</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Sending</surname>
          </string-name>
          (M, DS)→B, DS=
          <fpage>796</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          1. Calculation ℎ′ =  ( ′). Let ℎ′ =
          <volume>852</volume>
          &lt; p.
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          2. Calculation of J( ) =
          <fpage>583</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          3.
          <string-name>
            <surname>Calculation</surname>
            <given-names>DS</given-names>
          </string-name>
          ′ = H(
          <volume>852</volume>
          ||583) =
          <fpage>796</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          4. Verification:  ′ =  =
          <fpage>796</fpage>
          . DS is correct.
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>W.</given-names>
            <surname>Castryck</surname>
          </string-name>
          , et al.,
          <source>CSIDH: An Efficient Post-Quantum Commutative Group Action, Advances in Cryptology (ASIACRYPT)</source>
          (
          <year>2018</year>
          )
          <fpage>395</fpage>
          -
          <lpage>427</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          - 03332-3_
          <fpage>15</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>A.</given-names>
            <surname>Rostovtsev</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Stolbunov</surname>
          </string-name>
          ,
          <article-title>Public-Key Cryptosystem based on Isogenies</article-title>
          ,
          <source>IACR Cryptology ePrint Archive</source>
          <year>2006</year>
          /145 (
          <year>2006</year>
          ). https://ia.cr/
          <year>2006</year>
          /145
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>P. W.</given-names>
            <surname>Shor</surname>
          </string-name>
          , Algorithms for Quantum Computation:
          <article-title>Discrete Logarithms and Factoring</article-title>
          ,
          <source>in: 35th Annual Symposium on Foundations of Computer Science</source>
          (
          <year>1994</year>
          )
          <fpage>124</fpage>
          -
          <lpage>134</lpage>
          . doi:
          <volume>10</volume>
          .1109/SFCS.
          <year>1994</year>
          .365700
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>J. M.</given-names>
            <surname>Couveignes</surname>
          </string-name>
          ,
          <article-title>Hard homogeneous spaces</article-title>
          ,
          <source>IACR Cryptology ePrint Archive</source>
          <year>2006</year>
          /291 (
          <year>2006</year>
          ). https://ia.cr/
          <year>2006</year>
          /291
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.</given-names>
            <surname>Sokolov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Skladannyi</surname>
          </string-name>
          ,
          <article-title>Modeling of 3- and 5-Isogenies of Supersingular Edwards Curves</article-title>
          ,
          <source>in: 2nd International Workshop on Modern Machine Learning Technologies and Data Science (MoMLeT&amp;DS)</source>
          , vol.
          <volume>2631</volume>
          (
          <year>2020</year>
          )
          <fpage>30</fpage>
          -
          <lpage>39</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , et al.,
          <source>Analysis of 2-Isogeny Properties of Generalized Form Edwards Curves, Cybersecurity Providing in Information and Telecommunication Systems (CPITS)</source>
          , vol.
          <volume>2746</volume>
          (
          <year>2020</year>
          )
          <fpage>1</fpage>
          -
          <lpage>13</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>J. H.</given-names>
            <surname>Silverman</surname>
          </string-name>
          , The Arithmetic of Elliptic Curves, Graduate Texts in Mathematics (
          <year>2009</year>
          ). doi:
          <volume>10</volume>
          .1007/978-0-
          <fpage>387</fpage>
          -09494-6
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>A. V.</given-names>
            <surname>Sutherland</surname>
          </string-name>
          , Identifying Supersingular Elliptic Curves,
          <source>LMS Journal of Computation and Mathematics</source>
          <volume>15</volume>
          (
          <year>2012</year>
          )
          <fpage>317</fpage>
          -
          <lpage>325</lpage>
          . doi:
          <volume>10</volume>
          .1112/s1461157012001106
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , et al.,
          <source>Computing of Odd Degree Isogenies on Supersingular Twisted Edwards Eurves, Cybersecurity Providing in Information and Telecommunication Systems</source>
          , vol.
          <volume>2923</volume>
          (
          <year>2021</year>
          )
          <fpage>1</fpage>
          -
          <lpage>11</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          ,
          <article-title>How to Construct CSIDH on Quadratic and Twisted Edwards Curves</article-title>
          .
          <source>Cybersecurity: Education, Science, Technique</source>
          <volume>3</volume>
          (
          <issue>15</issue>
          ) (
          <year>2022</year>
          )
          <fpage>148</fpage>
          -
          <lpage>163</lpage>
          . doi:
          <volume>10</volume>
          .28925/
          <fpage>2663</fpage>
          -
          <lpage>4023</lpage>
          .
          <year>2022</year>
          .
          <volume>15</volume>
          .148163
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Kovalchuk</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Abramov</surname>
          </string-name>
          ,
          <article-title>Randomization of CSIDH Algorithm on Quadratic and Twisted Edwards Curves</article-title>
          ,
          <source>Cybersecurity: Education, Science, Technique</source>
          <volume>1</volume>
          (
          <issue>17</issue>
          ) (
          <year>2022</year>
          )
          <fpage>128</fpage>
          -
          <lpage>144</lpage>
          . doi:
          <volume>10</volume>
          .28925/
          <fpage>2663</fpage>
          -
          <lpage>4023</lpage>
          .
          <year>2022</year>
          .
          <volume>17</volume>
          .128144
        </mixed-citation>
      </ref>
      <ref id="ref20">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , et al.,
          <string-name>
            <surname>Modeling</surname>
            <given-names>CSIKE</given-names>
          </string-name>
          <article-title>Algorithm on Non-Cyclic Edwards Curves</article-title>
          ,
          <source>in: Cybersecurity Providing in Information and Telecommunication Systems</source>
          , vol.
          <volume>3288</volume>
          (
          <year>2022</year>
          )
          <fpage>1</fpage>
          -
          <lpage>10</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref21">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>S.</given-names>
            <surname>Kim</surname>
          </string-name>
          , et al.,
          <source>Optimized Method for Computing Odd-Degree Isogenies on Edwards Curves, Advances in Cryptology (ASIACRYPT)</source>
          (
          <year>2019</year>
          )
          <fpage>273</fpage>
          -
          <lpage>292</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          -34621-8_
          <fpage>10</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref22">
        <mixed-citation>
          [14]
          <string-name>
            <given-names>R.</given-names>
            <surname>Farashahi</surname>
          </string-name>
          , S. Hosseini,
          <source>Differential Addition on Twisted Edwards Curves, Lecture Notes in Computer Science</source>
          (
          <year>2017</year>
          )
          <fpage>366</fpage>
          -
          <lpage>378</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>319</fpage>
          -59870- 3_
          <fpage>21</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref23">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>D.</given-names>
            <surname>Moody</surname>
          </string-name>
          , D. Shumow,
          <article-title>Analogues of Vélu's Formulas for Isogenies on Alternate Models of Elliptic Curves</article-title>
          ,
          <source>Mathematics of Computation</source>
          <volume>85</volume>
          (
          <issue>300</issue>
          ) (
          <year>2015</year>
          )
          <fpage>1929</fpage>
          -
          <lpage>1951</lpage>
          . doi:
          <volume>10</volume>
          .1090/mcom/3036
        </mixed-citation>
      </ref>
      <ref id="ref24">
        <mixed-citation>
          [16]
          <string-name>
            <given-names>D. J.</given-names>
            <surname>Bernstein</surname>
          </string-name>
          , et al.,
          <source>Twisted Edwards Curves, Lecture Notes in Computer Science</source>
          , vol.
          <volume>5023</volume>
          (
          <year>2008</year>
          )
          <fpage>389</fpage>
          -
          <lpage>405</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>540</fpage>
          -68164-9_
          <fpage>26</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref25">
        <mixed-citation>
          [17]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          ,
          <article-title>Elliptic Curves in Edwards Form and Cryptography (</article-title>
          <year>2017</year>
          )
          <article-title>(In Russian)</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref26">
        <mixed-citation>
          [18]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>O.</given-names>
            <surname>Tsygankova</surname>
          </string-name>
          ,
          <article-title>Number of Curves in the Generalized Edwards Form with Minimal Even Cofactor of the Curve Order</article-title>
          ,
          <source>Problems of Information Transmission</source>
          <volume>53</volume>
          (
          <issue>1</issue>
          ) (
          <year>2017</year>
          )
          <fpage>92</fpage>
          -
          <lpage>101</lpage>
          . doi:
          <volume>10</volume>
          .1134/S0032946017010082
        </mixed-citation>
      </ref>
      <ref id="ref27">
        <mixed-citation>
          [19]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , L. Kovalchuk,
          <article-title>Supersingular Twisted Edwards Curves Over Prime Fields. I. Supersingular Twisted Edwards Curves with j-Invariants Equal to Zero and 123, Cybernetics</article-title>
          and
          <source>Systems Analysis</source>
          <volume>55</volume>
          (
          <issue>3</issue>
          ) (
          <year>2019</year>
          )
          <fpage>347</fpage>
          -
          <lpage>353</lpage>
          . doi:
          <volume>10</volume>
          .1007/s10559-019- 00140-9
        </mixed-citation>
      </ref>
      <ref id="ref28">
        <mixed-citation>
          [20]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , L. Kovalchuk,
          <article-title>Supersingular Twisted Edwards Curves over Prime Fields.* II. Supersingular Twisted Edwards Curves with the j-Invariant Equal to 663, Cybernetics</article-title>
          and
          <source>Systems Analysis</source>
          <volume>55</volume>
          (
          <issue>5</issue>
          ) (
          <year>2019</year>
          )
          <fpage>731</fpage>
          -
          <lpage>741</lpage>
          . doi:
          <volume>10</volume>
          .1007/s10559-019- 00183-y
        </mixed-citation>
      </ref>
      <ref id="ref29">
        <mixed-citation>
          [21]
          <string-name>
            <given-names>H.</given-names>
            <surname>Onuki</surname>
          </string-name>
          , et al.,
          <source>A Faster Constant-Time Algorithm of CSIDH Keeping Two Points. Lecture Notes in Computer Science</source>
          (
          <year>2019</year>
          )
          <fpage>23</fpage>
          -
          <lpage>33</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          -26834-
          <issue>3</issue>
          _
          <fpage>2</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref30">
        <mixed-citation>
          [22]
          <string-name>
            <given-names>A.</given-names>
            <surname>Jalali</surname>
          </string-name>
          , et al.,
          <source>Towards Optimized and Constant-Time CSIDH on Embedded Devices, Lecture Notes in Computer Science</source>
          (
          <year>2019</year>
          )
          <fpage>215</fpage>
          -
          <lpage>231</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          - 16350-1_
          <fpage>12</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref31">
        <mixed-citation>
          [23]
          <string-name>
            <given-names>M.</given-names>
            <surname>Meyer</surname>
          </string-name>
          , S. Reith,
          <string-name>
            <given-names>A Faster</given-names>
            <surname>Way to the</surname>
          </string-name>
          <string-name>
            <surname>CSIDH</surname>
          </string-name>
          ,
          <source>Lecture Notes in Computer Science</source>
          , vol.
          <volume>11356</volume>
          (
          <year>2018</year>
          )
          <fpage>137</fpage>
          -
          <lpage>152</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>030</fpage>
          -05378-
          <issue>9</issue>
          _
          <fpage>8</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref32">
        <mixed-citation>
          [24]
          <string-name>
            <given-names>M.</given-names>
            <surname>Meyer</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Campos</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Reith</surname>
          </string-name>
          ,
          <source>On Lions and Elligators: An Efficient Constant-Time Implementation of CSIDH, Cryptology ePrint Archive</source>
          <year>2018</year>
          /1198 (
          <year>2018</year>
          ). https://ia.cr/
          <year>2018</year>
          /1198
        </mixed-citation>
      </ref>
      <ref id="ref33">
        <mixed-citation>
          [25]
          <string-name>
            <given-names>D.</given-names>
            <surname>Jao</surname>
          </string-name>
          , et al.,
          <source>Supersingular Isogeny Key Encapsulation (NIST Round 2)</source>
          (
          <year>2019</year>
          ).
          <source>doi: 10.13140/RG.2.2.26543.07847</source>
        </mixed-citation>
      </ref>
      <ref id="ref34">
        <mixed-citation>
          [26]
          <article-title>National Institute of Standards and Technology. Post-quantum cryptography standardization (</article-title>
          <year>2016</year>
          ). https://csrc.nist.gov/Projects/Post-QuantumCryptography/Post-QuantumCryptography-Standardization
        </mixed-citation>
      </ref>
      <ref id="ref35">
        <mixed-citation>
          [27]
          <string-name>
            <surname>M. Qi.</surname>
          </string-name>
          ,
          <article-title>An Efficient Post-Quantum KEM from CSIDH</article-title>
          ,
          <source>Journal of Mathematical Cryptology</source>
          <volume>16</volume>
          (
          <issue>1</issue>
          ) (
          <year>2022</year>
          )
          <fpage>103</fpage>
          -
          <lpage>113</lpage>
          . doi:
          <volume>10</volume>
          .1515/jmc-2022
          <source>-0007</source>
        </mixed-citation>
      </ref>
      <ref id="ref36">
        <mixed-citation>
          [28]
          <string-name>
            <given-names>L. C.</given-names>
            <surname>Washington</surname>
          </string-name>
          , Elliptic
          <string-name>
            <surname>Curves</surname>
          </string-name>
          (
          <year>2008</year>
          ).
          <source>doi: 10.1201/9781420071474</source>
        </mixed-citation>
      </ref>
      <ref id="ref37">
        <mixed-citation>
          [29]
          <string-name>
            <given-names>A.</given-names>
            <surname>Stolbunov</surname>
          </string-name>
          ,
          <source>Constructing Public-Key Cryptographic Schemes based on Class Group Action on a Set of Isogenous Elliptic Curves, Advances in Mathematics of Communications</source>
          <volume>4</volume>
          (
          <issue>2</issue>
          ) (
          <year>2010</year>
          )
          <fpage>215</fpage>
          -
          <lpage>235</lpage>
          . doi:
          <volume>10</volume>
          .3934/amc.
          <year>2010</year>
          .
          <volume>4</volume>
          .
          <fpage>215</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref38">
        <mixed-citation>
          [30]
          <string-name>
            <given-names>A.</given-names>
            <surname>Stolbunov</surname>
          </string-name>
          .
          <source>Cryptographic Schemes Based on Isogenies</source>
          ,
          <source>PhD Thesis</source>
          (
          <year>2011</year>
          ).
          <source>doi: 10.13140/RG.2.2.20826.44488</source>
        </mixed-citation>
      </ref>
      <ref id="ref39">
        <mixed-citation>
          [31]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , et al.,
          <article-title>Implementation of the CSIDH Algorithm Model on Supersingular Twisted and Quadratic Edwards Curves, Cybersecurity Providing in Information and Telecommunication Systems (CPITS-II)</article-title>
          , vol.
          <volume>3187</volume>
          (
          <year>2022</year>
          )
          <fpage>302</fpage>
          -
          <lpage>309</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref40">
        <mixed-citation>
          [32]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bessalov</surname>
          </string-name>
          , et al.,
          <article-title>CSIKE-ENC Combined Encryption Scheme with Optimized Degrees of Isogeny Distribution</article-title>
          ,
          <source>in: Cybersecurity Providing in Information and Telecommunication Systems (CPITS)</source>
          , vol.
          <volume>3421</volume>
          (
          <year>2023</year>
          )
          <fpage>36</fpage>
          -
          <lpage>45</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref41">
        <mixed-citation>
          [33]
          <string-name>
            <given-names>J.</given-names>
            <surname>Kieffer</surname>
          </string-name>
          ,
          <article-title>Accelerating the Couveignes Rostovtsev Stolbunov key exchange protocol</article-title>
          , Rapport de stage de Master 2
          <article-title>Mathématiques fondamentales</article-title>
          , Université Paris VIUniversité Paris VI (
          <year>2017</year>
          )
          <article-title>(In French)</article-title>
          .
          <source>doi: 10</source>
          .48550/arXiv.
          <year>1804</year>
          .10128
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>