# Patent

## Title
Pairing computation device, pairing computation method and recording medium recording pairing computation program

## Abstract

There is provided a pairing computation device, a pairing computation method and a recording medium recording a pairing computation program all enabling a pairing computation at high speed. In the case where the Ate pairing e(Q,P) is defined as in equation (1) and k is any one of even, a multiple of 3, a multiple of 4, and a multiple of 6, arithmetic operations of a rational function which is required to derive the Miller function fs,Q(P) are performed on the proper subfield specified by a twist curve using a quadratic and cubic non-residue v which becomes 1 when an exponentiation of this fs,Q(P) to the power of (qk−1)/r is performed.

## Background
FIELD OF THE INVENTION
The present invention relates to a pairing computation device, a pairing computation method, and a non-transitory recording medium recording a paring computation program all enabling a pairing computation at high speed.
DESCRIPTION OF THE RELATED ART
Recently, since there have been available high speed telecommunication lines at low prices, on a network such as the Internet, it has been possible to provide various services such as distribution of music or video, Internet banking, and electronic application to administrative agencies.
Further, corporate business styles have changed, such that it has been made possible, using so-called mobile terminals such as a notebook PC or a mobile phone, to access the company's server system from an outside location to input or get various kinds of information.
Particularly in Japan, services in which mobile phones are used as connecting terminals to a network are satisfactory, so that advent of so-called ubiquitous society in which users can access the required network from anywhere to get various kinds of information or to use services has become more likely.
In accessing a network in this way to get various kinds of information or to use services, it has been made possible to get information or to use services by the fact that authentication processing is performed through an authentication server equipped with required authentication means and a user who has accessed the network is authenticated as one of the certain users registered in advance.
Particularly in recent years, since, using digital signature technology, it has been made possible to safely deal with even highly-confidential information on a network by enabling to guarantee that the information dealt with is not altered by a third party or is not leaked to a third party, more active use of a network has become popular.
However, in this digital signature, since an individual user is identified, a history of the user is to be accumulated as information based on an authentication processing in an authentication server. This history information is a kind of private information, and hence, there has been a problem of protection of private information.
Accordingly, there has been proposed to use digital group signature which is extended from the digital signature. In the case where the digital group signature is used, a user sends anonymously to an authentication server a signature data which only certifies that the user belongs to a certain group, then the authentication server authenticates that the user belongs to the prescribed group based on the received signature data without identifying the individual user. Therefore, while it is possible to prevent accumulating history information on each individual user's usage in the authentication server, it is also possible to block unfair use by users who do not belong to the group.
Here, in the anonymous authentication in the digital group signature, a pairing computation has been used. In the pairing computation, arithmetic operations using a function of 2 inputs and 1 output are performed. It takes advantage of the fact that, for example, letting P be a rational point over a prime field Fq, and Q be a rational point over a finite extension field Fqk of degree k, in the case of pairing where P and Q are inputted and an element z in the finite extension field F*qk is outputted, when a times of P and b times of Q are inputted, an ab-th power of z is calculated. The property like this is called bilinearity. In addition, here, “k” is called an embedding degree, and “F*qk” is meant to be correctly displayed as in the following representation, but due to display restrictions, it is denoted as “F*qk”.

F*qk 

In general, as rational points P and Q, points on an elliptic curve are used respectively. A pairing computation on the elliptic curve like this consists of a step of performing arithmetic operations using Miller's algorithm and a step of performing arithmetic operations at the final exponentiation.
For example, in a digital group signature whose group consists of 10000 members, there has been known a method which, when the digital group signature performs verification, performs a pairing computation for members who invalidated their right in order to make it possible to flexibly correspond to validation and invalidation of access right of each member. In this case, supposing 100 members who invalidated their right, there requires 100 times of the pairing computation. Since one pairing computation by a currently general electronic computer requires approximately 0.1 second, 100 times of the pairing computation comes to require 10 seconds, and hence, use of the pairing computation in a large scale digital group signature system has not been practical.
Accordingly, various kinds of development have been made to improve pairing computation speed. For example, there has been proposed an acceleration technique for pairing computation defined on an elliptic curve over a finite field (see Patent Document 1 for example).
Patent Document 1: JP-A-2005-316267.
SUMMARY OF THE INVENTION
However, currently proposed acceleration techniques for pairing computation have not yet come to the speed suited for the use in a large scale digital group signature, and hence, there has been a problem that the number of members of the digital group signature must be limited.
The inventors, in view of the present situation, have made a study to speed up pairing computation and have made the present invention.
According to a first aspect of the present invention, there is provided a pairing computation device, wherein an elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq (q: a power of a prime number greater than 3), letting an embedding degree be 2h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q2h/(F*q2h)r,

letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    2
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E14
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3 using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q2h−1)/r.

In particular, the pairing computation device is characterized in that it includes:
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1, m←ypv−3/2 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−3/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′,T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−3/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q2h−1)/r.
According to a second aspect of the present invention, there is provided a pairing computation device, wherein an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 3h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map,

e:G2×G1→F*q3h/(F*q3h)r 

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    3
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E15
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1 using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q3h−1)/r.

In particular, the pairing computation device is characterized in that it includes:
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′,T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q3h−1)/r.
According to a third aspect of the present invention, there is provided a pairing computation device, wherein an elliptic curve equation is given as y2=x3+ax, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 4h-th degree (h: a natural number), qh−1 be dividable by 4, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q])

as a non-degenerate bilinear map

e:G2×G1→F*q4h/(F*q4h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    4
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E16
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q4h−1)/r.

In particular, the pairing computation device is characterized in that it includes:
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1/2, m←ypv−3/4 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−3/4 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′,T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−3/4 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q4h−1)/r.
According to a fourth aspect of the present invention, there is provided a pairing computation device, wherein an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 6h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q6h/(F*q6h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    6
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E17
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1 using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q6h−1)/r.

In particular, the pairing computation device is characterized in that it includes:
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′,T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q6h−1)/r.
According to a fifth aspect of the present invention, there is provided a pairing computation method, wherein an elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq (q: a power of a prime number greater than 3), letting an embedding degree be 2h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q2h/(F*q2h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    2
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E18
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q2h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3, and by comprising:

a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1, m←ypv−3/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q2h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a sixth aspect of the present invention, there is provided a pairing computation method, wherein an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 3h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q3h/(F*q3h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as,





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    3
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E19
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q3h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by comprising:

a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q3h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a seventh aspect of the present invention, there is provided a pairing computation method, wherein an elliptic curve equation is given as y2=x3+ax, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 4h-th degree (h: a natural number), qh−1 be dividable by 4, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q4h/(F*q4h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    4
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E20
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q4h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x, and by comprising:

a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/2, m←ypv−3/4 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/4 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/4 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q4h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a eighth aspect of the present invention, there is provided a pairing computation method, wherein an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 6h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q6h/(F*q6h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    6
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E21
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q6h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by comprising:

a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)vv−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q6h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a ninth aspect of the present invention, there is provided a non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+ax+b, aεFq bεFq (q: a power of a prime number greater than 3), letting an embedding degree be 2h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q])

as a non-degenerate bilinear map

e:G2×G1→F*q2h/(F*q2h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    2
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E22
        ]
      
    
  




the pairing computation program characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q2h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1, m←ypv−3/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q2h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a tenth aspect of the present invention, there is provided a non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 3h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q3h/(F*q3h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    3
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E23
        ]
      
    
  




the pairing computation program characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q3h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by causing the electronic computer to perform:

a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m) by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q3h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a eleventh aspect of the present invention, there is provided a non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+ax, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 4h-th degree (h: a natural number), qh−1 be dividable by 4, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q4h/(F*q4h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    4
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E24
        ]
      
    
  




the pairing computation program characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q4h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/2, m←ypv−3/4 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/4 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/4 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q4h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
According to a twelfth aspect of the present invention, there is provided a non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting an embedding degree be 6h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, and φq be the Frobenius endomorphism, the Ate pairing e is defined, using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q]),

as a non-degenerate bilinear map

e:G2×G1→F*q6h/(F*q6h)r,

and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as





  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                F
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    6
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E25
        ]
      
    
  




the pairing computation program characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q6h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of the tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of the straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q6h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.
The present invention, in performing pairing computation, reduces a burden of arithmetic operations by performing arithmetic operations of a rational function which is required to derive fs,Q(P) on the lower degree proper subfield using a twist curve, thus enabling to speedup the pairing computation.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is an explanatory view of conventional Miller's algorithm;
FIG. 2 is an explanatory view of Miller's algorithm according to the present invention;
FIG. 3 is a flowchart of a pairing computation program according to the present invention; and
FIG. 4 is a block diagram of a pairing computation device according to an embodiment of the present invention.



EXPLANATION OF SYMBOLS
10 electronic computer
11 CPU
12 storage unit
13 memory unit
14 bus
15 input output control part
20 telecommunication lines
30 client device
DESCRIPTION OF THE PREFERRED EMBODIMENTS
A paring computation device, a pairing computation method, and a pairing computation program according to the present invention make possible high speed computation by speeding up an arithmetic operation using Miller's algorithm in a paring computation which is constituted of a step of performing an arithmetic operation using Miller's algorithm and a step of performing the final exponentiation operation.
The Miller's algorithm is known as an algorithm shown in FIG. 1. Supposing an elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq (q: a power of a prime number greater than 3), letting an embedding degree be 2h-th degree (h: a natural number), E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, φq be the Frobenius endomorphism, the Miller's algorithm inputs s=t−1 using a trace t of the Frobenius endomorphism φq, P which satisfies PεE(Fq), and Q which satisfies QεE(Fq2h), and outputs fs,Q(P).
In addition, an embedding degree is not limited to even degree. In the case where the embedding degree is a multiple of 3, it is desirable that an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), an embedding degree is 3h-th degree (h: a natural number), and an additive group consisting of rational points on a pairing-enabled elliptic curve defined over Fq3h is used.
Further, particularly in the case where the embedding degree is a multiple of 4, it is desirable that an elliptic curve equation is given as y2=x3+ax, aεFq (q: a power of a prime number greater than 3), an embedding degree is 4h-th degree (h: a natural number), 4 divides qh−1, and additive group consisting of rational points on a pairing-enabled elliptic curve defined over Fq4h is used.
Still further, particularly in the case where the embedding degree is a multiple of 6, it is desirable that an elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), an embedding degree is 6h-th degree (h: a natural number), and an additive group consisting of rational points on a pairing-enabled elliptic curve defined over Fq6h is used.
In the first step of the Miller's algorithm shown in FIG. 1, substitution operations f←1, T←Q are performed and the following arithmetic operations are iterated based on a condition of iteration in the second step. In the second step, to be more specific, the number of bits obtained by binary representation of s is set as an initial value of i, an iteration count is controlled by decrementing i by 1 for each iteration of arithmetic operations. s[i] in the fifth step denotes a value of the i-th bit from the least bit of s in binary representation, and s[i] is any one of “0” and “1”.
In the third step of the Miller's algorithm shown in FIG. 1, an arithmetic operation of a value of lT,T(xp,yp) substituted by coordinates of rational point P in the left hand side of lT,T(x,y)=0 which is the tangent line passing through rational point T is performed. In the fourth step, an elliptic curve doubling is performed as T←2T.
In the fifth step of the Miller's algorithm shown in FIG. 1, it is determined whether or not a bit value of s[i] in the case of binary representation of s is 1, and in the case of 1, in the sixth step, an arithmetic operation of a value of lT,Q(xp,yp) substituted by coordinates of rational point P in the left hand side of lT,Q(x,y)=0 which is the straight line passing through rational points T and Q is performed. In the seventh step, an elliptic curve addition is performed as T←T+Q.
In the tenth step of the Miller's algorithm shown in FIG. 1, arithmetic operations are iterated based on the condition of iteration in the second step, and a result of the arithmetic operations is outputted. In pairing computation, finally, an exponentiation operation to the power of (q2h−1)/r is performed with respect to fs,Q(P) which is a result of arithmetic operations outputted from the Miller's algorithm.
In addition, in the case of the embedding degree being a multiple of 3, the final exponentiation is to the power of (q3h−1)/r, in the case of the embedding degree being a multiple of 4, the final exponentiation is to the power of (q4h−1)/r, and in the case of the embedding degree being a multiple of 6, the final exponentiation is to the power of (q6h−1)/r.
In what follows, arithmetic operation methods of lT,T(xp,yp) in the third step and lT,Q(xp,yp) in the fifth step are explained. Here, let the embedding degree be even degree, that is, 2h-th degree (h: a natural number) and the elliptic curve be y2=x3+ax+b, aεFq, bεFq (q: a power of a prime number greater than 3).
Letting T=(xT, yT), and Q=(xQ, yQ), λT,T which is the gradient of lT,T(xp,yp) and λT,Q which is the gradient of lT,Q(xp,yp) are given by the following equation.




  
    
      
        
          λ
          
            T
            ,
            T
          
        
        =
        
          
            
              3
              ⁢
              
                x
                T
                2
              
            
            +
            a
          
          
            2
            ⁢
            
              y
              T
            
          
        
      
    
    
      
        [
        E26
        ]
      
    
  
  
    
      
        
          
            λ
            
              T
              ,
              Q
            
          
          =
          
            
              
                y
                Q
              
              -
              
                y
                T
              
            
            
              
                x
                Q
              
              -
              
                x
                T
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          (
          
            T
            ≠
            Q
          
          )
        
      
    
    
      
        [
        E27
        ]
      
    
  




A value of lT,Q(xp,yp) is calculated as follows using this gradient λT,Q.

λT,Q(xp,yp)=(xp−xT)λT,Q−(yp−yT)  [E28]

Although, usually, the calculation is performed on E(Fq2h) using [E26] and [E27], in the case where the embedding degree is an even degree, that is, 2h-th degree (h: a natural number) and the elliptic curve is given as y2=x3+ax+b, aεFq, bεFq (q: a power of a prime number greater than 3), there exists a twist curve of E(Fqh), that is E′(Fqh):y2=x3+av−2x+bv−3, vεFqh, and there exists a homomorphism from E′(Fqh) to E(Fq2h), that is, (x,y)→(xv,yv3/2). Here, a new parameter v is a quadratic non-residue element which becomes 1 by an exponentiation of itself to the power of (q2h−1)/r.
Therefore, it is possible to perform the elliptic curve addition on E(Fq2h) [r] substituting the elliptic curve addition on E′(Fqh) [r] which has a 1⁄2 extension degree, thus enabling to reduce a burden of arithmetic operations and to perform high speed computation.
In particular, since there exist rational points T′,Q′ on E′(Fqh) where relations T=(xT, yT)=(xT′v, yT′v3/2), Q=(xQ, yQ)=(xQ′v, yQ′v3/2) hold with respect to rational points T,QεE(Fq2h), first, λT,T, the gradient of lT,T(x,y) and λT,Q, the gradient of lT,Q(xp,yp) are transformed as follows.




  
    
      
        
          λ
          
            T
            ,
            T
          
        
        =
        
          
            
              
                3
                ⁢
                
                  x
                  T
                  2
                
              
              +
              a
            
            
              2
              ⁢
              
                y
                T
              
            
          
          =
          
            
              
                
                  (
                  
                    
                      3
                      ⁢
                      
                        x
                        T
                        2
                      
                    
                    +
                    
                      av
                      
                        -
                        2
                      
                    
                  
                  )
                
                ⁢
                
                  v
                  2
                
              
              
                2
                ⁢
                
                  y
                  T
                  ′
                
                ⁢
                
                  v
                  
                    3
                    2
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    T
                    ′
                  
                
              
              ⁢
              
                v
                
                  1
                  2
                
              
            
          
        
      
    
    
      
        [
        E29
        ]
      
    
  
  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            
              
                y
                Q
              
              -
              
                y
                T
              
            
            
              
                x
                Q
              
              -
              
                x
                T
              
            
          
          =
          
            
              
                
                  (
                  
                    
                      y
                      Q
                      ′
                    
                    -
                    
                      y
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    3
                    2
                  
                
              
              
                
                  (
                  
                    
                      x
                      Q
                      ′
                    
                    -
                    
                      x
                      T
                      ′
                    
                  
                  )
                
                ⁢
                v
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    Q
                    ′
                  
                
              
              ⁢
              
                
                  v
                  
                    1
                    2
                  
                
                ⁢
                
                  

                
                (
                
                  T
                  ≠
                  Q
                
                )
              
            
          
        
      
    
    
      
        [
        E30
        ]
      
    
  




The gradient λT,Q becomes as follows from [E29] and [E30], irrespective of T=Q, T≠Q.




  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            λ
            
              
                T
                ′
              
              ,
              
                Q
                ′
              
            
          
          ⁢
          
            v
            
              1
              2
            
          
        
      
    
    
      
        [
        
          E
          ⁢
          
            
          
          ⁢
          31
        
        ]
      
    
  




Therefore, from [E31], the equation [E28] becomes as follows.




  
    
      
        
          
            
              
                
                  l
                  
                    T
                    ,
                    Q
                  
                
                ⁡
                
                  (
                  
                    
                      x
                      P
                    
                    ,
                    
                      y
                      P
                    
                  
                  )
                
              
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        x
                        T
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      T
                      ,
                      Q
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      y
                      T
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        
                          x
                          
                            T
                            ′
                          
                        
                        ⁢
                        v
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        ′
                      
                    
                  
                  ⁢
                  
                    v
                    
                      1
                      2
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      
                        y
                        
                          T
                          ′
                        
                      
                      ⁢
                      
                        v
                        
                          3
                          2
                        
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    
                      
                        (
                        
                          
                            
                              x
                              P
                            
                            ⁢
                            
                              v
                              
                                -
                                1
                              
                            
                          
                          -
                          
                            x
                            
                              T
                              ′
                            
                          
                        
                        )
                      
                      ⁢
                      
                        λ
                        
                          
                            T
                            ′
                          
                          ,
                          
                            Q
                            ′
                          
                        
                      
                    
                    -
                    
                      (
                      
                        
                          
                            y
                            P
                          
                          ⁢
                          
                            v
                            
                              -
                              
                                3
                                2
                              
                            
                          
                        
                        -
                        
                          y
                          T
                          ′
                        
                      
                      )
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    3
                    2
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    l
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        ′
                      
                    
                  
                  (
                  
                    
                      
                        x
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          1
                        
                      
                    
                    ,
                    
                      
                        y
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            3
                            2
                          
                        
                      
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    3
                    2
                  
                
              
            
          
        
      
    
    
      
        [
        
          E
          ⁢
          
            
          
          ⁢
          32
        
        ]
      
    
  




Since the final part of v3/2 in the equation [E32] becomes 1 at the final exponentiation to the power of (q2h−1)/r in the pairing computation, in calculating the lT,Q(xp,yp) in the Miller's algorithm, it is enough to calculate lT′,Q′(xpv−1,ypv3/2) on the proper subfield instead of calculating lT,Q(xp,yp), thus enabling to reduce a burden of arithmetic operations significantly.
That is, the elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq (q: a power of a prime number greater than 3), letting the embedding degree be 2h-th degree (h: a natural number), E be an additive group constituted of rational points on the pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, φq be the Frobenius endomorphism, and using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q])

the Ate pairing e is defined as a non-degenerate bilinear map, e: G2×G1→F*q2h/(F*q2h)r, and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using Miller function fs,Q(•), in computing the Ate pairing e(Q,P) as





  
    
      
        
          e
          ⁡
          
            (
            
              Q
              ,
              P
            
            )
          
        
        =
        
          
            
              F
              
                s
                ,
                Q
              
            
            ⁡
            
              (
              P
              )
            
          
          
            
              
                q
                
                  2
                  ⁢
                  h
                
              
              -
              1
            
            r
          
        
      
    
    
      
        [
        E33
        ]
      
    
  




the arithmetic operations of the rational function lT,T(xp,yp) and lT,Q(xp,yp) which are required to derive fs,Q(P) are performed on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3, using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q2h−1)/r, as arithmetic operations lT′,T′(xpv−1,ypv−3/2) and lT′,Q′(xpv−1,ypv−3/2), and hence, it is possible to reduce a burden of arithmetic operations and to perform high speed computation.
In the case where the embedding degree is a multiple of 3, that is, 3h-th degree (h: a natural number) and the elliptic curve is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), there exists a twist curve of E(Fqh), that is E′(Fqh):y2=x3+av−1, vεFqh, and there exists a homomorphism from E′(Fqh) to E(Fq3h), that is, (x,y)→(xv1/3,yv1/2). Here, a new parameter v is a quadratic and cubic non-residue element which becomes 1 by an exponentiation of itself to the power of (q3h−1)/r.
Therefore, it is possible to perform the elliptic curve addition on E(Fq3h) [r] substituting the elliptic curve addition on E′(Fqh) [r] which has a 1⁄3 extension degree, thus enabling to reduce a burden of arithmetic operations and to perform high speed computation.
In particular, since there exist rational points T′,Q′ on E′(Fqh) where relations T=(xT, yT)=(xT′v1/3, yT′v1/2), Q=(xQ, yQ)=(xQ′v1/3, yQ′v1/2) hold with respect to rational points T,QεE(Fq3h), first, λT,T, the gradient of lT,T(x,y) and λT,Q, the gradient of lT,Q(xp,yp) are transformed as follows.




  
    
      
        
          λ
          
            T
            ,
            T
          
        
        =
        
          
            
              3
              ⁢
              
                x
                T
                2
              
            
            
              2
              ⁢
              
                y
                T
              
            
          
          =
          
            
              
                3
                ⁢
                
                  
                    x
                    ′
                  
                  T
                  2
                
                ⁢
                
                  v
                  
                    2
                    3
                  
                
              
              
                2
                ⁢
                
                  y
                  T
                  ′
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    T
                    ′
                  
                
              
              ⁢
              
                v
                
                  1
                  6
                
              
            
          
        
      
    
    
      
        [
        E34
        ]
      
    
  
  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            
              
                y
                Q
              
              -
              
                y
                T
              
            
            
              
                x
                Q
              
              -
              
                x
                T
              
            
          
          =
          
            
              
                
                  (
                  
                    
                      y
                      Q
                      ′
                    
                    -
                    
                      y
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
              
                
                  (
                  
                    
                      x
                      Q
                      ′
                    
                    -
                    
                      x
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    3
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    Q
                    ′
                  
                
              
              ⁢
              
                
                  v
                  
                    1
                    6
                  
                
                ⁢
                
                  

                
                (
                
                  T
                  ≠
                  Q
                
                )
              
            
          
        
      
    
    
      
        [
        E35
        ]
      
    
  




The gradient λT,Q becomes as follows from [E34] and [E35], irrespective of T=Q, T≠Q.




  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            λ
            
              
                T
                ′
              
              ,
              
                Q
                ′
              
            
          
          ⁢
          
            v
            
              1
              6
            
          
        
      
    
    
      
        [
        
          E
          ⁢
          
            
          
          ⁢
          36
        
        ]
      
    
  




Therefore, from [E36], the equation [E28] becomes as follows.




  
    
      
        
          
            
              
                
                  l
                  
                    T
                    ,
                    Q
                  
                
                ⁡
                
                  (
                  
                    
                      x
                      P
                    
                    ,
                    
                      y
                      P
                    
                  
                  )
                
              
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        x
                        T
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      T
                      ,
                      Q
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      y
                      T
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        
                          x
                          
                            T
                            ′
                          
                        
                        ⁢
                        
                          v
                          
                            1
                            3
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        ′
                      
                    
                  
                  ⁢
                  
                    v
                    
                      1
                      6
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      
                        y
                        
                          T
                          ′
                        
                      
                      ⁢
                      
                        v
                        
                          1
                          2
                        
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    
                      
                        (
                        
                          
                            
                              x
                              P
                            
                            ⁢
                            
                              v
                              
                                -
                                
                                  1
                                  3
                                
                              
                            
                          
                          -
                          
                            x
                            
                              T
                              ′
                            
                          
                        
                        )
                      
                      ⁢
                      
                        λ
                        
                          
                            T
                            ′
                          
                          ,
                          
                            Q
                            ′
                          
                        
                      
                    
                    -
                    
                      (
                      
                        
                          
                            y
                            P
                          
                          ⁢
                          
                            v
                            
                              -
                              
                                1
                                2
                              
                            
                          
                        
                        -
                        
                          y
                          T
                          ′
                        
                      
                      )
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    l
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        ′
                      
                    
                  
                  (
                  
                    
                      
                        x
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            1
                            3
                          
                        
                      
                    
                    ,
                    
                      
                        y
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            1
                            2
                          
                        
                      
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
          
        
      
    
    
      
        [
        
          E
          ⁢
          
            
          
          ⁢
          37
        
        ]
      
    
  




Since the final part of v1/2 in the equation [E37] becomes 1 at the final exponentiation to the power of (q3h−1)/r in the pairing computation, in calculating the lT,Q(xp,yp) in the Miller's algorithm, it is enough to calculate lT′,Q′(xpv−1/3,ypv−1/2) on the proper subfield instead of calculating lT,Q(xp,yp), thus enabling to reduce a burden of arithmetic operations significantly.
That is, the elliptic curve equation is given as y2=x3+a, aεFq, (q: a power of a prime number greater than 3), letting the embedding degree be 3h-th degree (h: a natural number), E be an additive group constituted of rational points on the pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, φq be the Frobenius endomorphism, and using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q])

the Ate pairing e is defined as a non-degenerate bilinear map, e: G2×G1→F*q3h/(F*q3h)r, and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), in computing the Ate pairing e(Q,P) as





  
    
      
        
          e
          ⁡
          
            (
            
              Q
              ,
              P
            
            )
          
        
        =
        
          
            
              F
              
                s
                ,
                Q
              
            
            ⁡
            
              (
              P
              )
            
          
          
            
              
                q
                
                  3
                  ⁢
                  h
                
              
              -
              1
            
            r
          
        
      
    
    
      
        [
        E38
        ]
      
    
  




the arithmetic operations of the rational function lT,T(xp,yp) and lT,Q(xp,yp) which are required to derive fs,Q(P) are performed on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q3h−1)/r, as arithmetic operations lT′,T′(xpv−1/3, ypv−1/2) and lT′,Q′(xpv−1/3,ypv−1/2) and hence, it is possible to reduce a burden of arithmetic operations and to perform high speed computation.
In the case where the embedding degree is a multiple of 4, that is, 4h-th degree (h: a natural number) and the elliptic curve is given as y2=x3+ax, aεFq, 4 divides qh−1 (q: a power of a prime number greater than 3), there exists a twist curve of E(Fqh), that is E′(Fqh):y2=x3+av−1x, vεFqh, and there exists a homomorphism from E′(Fqh) to E(Fq4h), that is, (x,y)→(xv1/2,yv3/4). Here, a new parameter v is a quadratic non-residue element which becomes 1 by an exponentiation of itself to the power of (q4h−1)/r.
Therefore, it is possible to perform the elliptic curve addition on E(Fq4h) [r] substituting the elliptic curve addition on E′(Fqh) [r] which has a 1⁄4 extension degree, thus enabling to reduce a burden of arithmetic operations and to perform high speed computation.
In particular, since there exist rational points T′,Q′ on E′(Fqh) where relations T=(xT, yT)=(xT′v1/2, yT′v3/4), Q=(xQ, yQ)=(xQ′v1/2, yQ′v3/4) hold with respect to rational points T,QεE(Fq4h), first, λT,T, the gradient of lT,T(x,y) and λT,Q, the gradient of lT,Q(xp,yp) are transformed as follows.




  
    
      
        
          λ
          
            T
            ,
            T
          
        
        =
        
          
            
              
                3
                ⁢
                
                  x
                  T
                  2
                
              
              +
              a
            
            
              2
              ⁢
              
                y
                T
              
            
          
          =
          
            
              
                
                  (
                  
                    
                      3
                      ⁢
                      
                        
                          x
                          ′
                        
                        T
                        2
                      
                    
                    +
                    
                      av
                      
                        -
                        1
                      
                    
                  
                  )
                
                ⁢
                v
              
              
                2
                ⁢
                
                  y
                  T
                  ′
                
                ⁢
                
                  v
                  
                    3
                    4
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    T
                    ′
                  
                
              
              ⁢
              
                v
                
                  1
                  4
                
              
            
          
        
      
    
    
      
        [
        E39
        ]
      
    
  
  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            
              
                y
                Q
              
              -
              
                y
                T
              
            
            
              
                x
                Q
              
              -
              
                x
                T
              
            
          
          =
          
            
              
                
                  (
                  
                    
                      y
                      Q
                      ′
                    
                    -
                    
                      y
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    3
                    4
                  
                
              
              
                
                  (
                  
                    
                      x
                      Q
                      ′
                    
                    -
                    
                      x
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    Q
                    ′
                  
                
              
              ⁢
              
                
                  v
                  
                    1
                    4
                  
                
                ⁢
                
                  

                
                (
                
                  T
                  ≠
                  Q
                
                )
              
            
          
        
      
    
    
      
        [
        E40
        ]
      
    
  




The gradient λT,Q becomes as follows from [E39] and [E40], irrespective of T=Q, T≠Q.




  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            λ
            
              
                T
                ′
              
              ,
              
                Q
                ′
              
            
          
          ⁢
          
            v
            
              1
              4
            
          
        
      
    
    
      
        [
        
          E
          ⁢
          
            
          
          ⁢
          41
        
        ]
      
    
  




Therefore, from [E41], the equation [E28] becomes as follows.




  
    
      
        
          
            
              
                
                  1
                  
                    T
                    ,
                    Q
                  
                
                ⁢
                
                  (
                  
                    
                      x
                      P
                    
                    ,
                    
                      y
                      P
                    
                  
                  )
                
              
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        x
                        T
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      T
                      ,
                      Q
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      y
                      T
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        
                          x
                          
                            T
                            ′
                          
                        
                        ⁢
                        
                          v
                          
                            1
                            2
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        ′
                      
                    
                  
                  ⁢
                  
                    v
                    
                      1
                      4
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      
                        y
                        
                          T
                          ′
                        
                      
                      ⁢
                      
                        v
                        
                          3
                          4
                        
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    
                      
                        (
                        
                          
                            
                              x
                              P
                            
                            ⁢
                            
                              v
                              
                                -
                                
                                  1
                                  2
                                
                              
                            
                          
                          -
                          
                            x
                            
                              T
                              ′
                            
                          
                        
                        )
                      
                      ⁢
                      
                        λ
                        
                          
                            T
                            ′
                          
                          ,
                          
                            Q
                            ′
                          
                        
                      
                    
                    -
                    
                      (
                      
                        
                          
                            y
                            P
                          
                          ⁢
                          
                            v
                            
                              -
                              
                                3
                                4
                              
                            
                          
                        
                        -
                        
                          y
                          T
                          ′
                        
                      
                      )
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    3
                    4
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  1
                  
                    
                      T
                      ′
                    
                    ,
                    
                      Q
                      ′
                    
                  
                
                ⁢
                
                  (
                  
                    
                      
                        x
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            1
                            2
                          
                        
                      
                    
                    ,
                    
                      
                        y
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            3
                            4
                          
                        
                      
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    3
                    4
                  
                
              
            
          
        
      
    
    
      
        [
        
          E
          ⁢
          
            
          
          ⁢
          42
        
        ]
      
    
  




Since the final part of v3/4 in the equation [E42] becomes 1 at the final exponentiation to the power of (q4h−1)/r in the pairing computation, in calculating the lT,Q(xp,yp) in the Miller's algorithm, it is enough to calculate lT′,Q′(xpv−1/2,ypv−3/4), on the proper subfield instead of calculating lT,Q(xp,yp), thus enabling to reduce a burden of arithmetic operations significantly.
That is, the elliptic curve equation is given as y2=x3+ax, aεFq, (q: a power of a prime number greater than 3), letting the embedding degree be 4h-th degree, qh−1 be dividable by 4 (h: a natural number), E be an additive group constituted of rational points on the pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, φq be the Frobenius endomorphism, and using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q])

the Ate pairing e is defined as a non-degenerate bilinear map, e: G2×G1→F*q4h/(F*q4h)r, and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using the Miller function fs,Q(•), in computing the Ate pairing e(Q,P) as





  
    
      
        
          e
          ⁡
          
            (
            
              Q
              ,
              P
            
            )
          
        
        =
        
          
            
              F
              
                s
                ,
                Q
              
            
            ⁡
            
              (
              P
              )
            
          
          
            
              
                q
                
                  4
                  ⁢
                  h
                
              
              -
              1
            
            r
          
        
      
    
    
      
        [
        E43
        ]
      
    
  




the arithmetic operations of the rational function lT,T(xp,yp) and lT,Q(xp,yp) which are required to derive fs,Q(P) are performed on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x, using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q4h−1)/r, as arithmetic operations lT′,T′(xpv−1/2,ypv−3/4) and lT′,Q′(xpv−1/2,ypv−3/4), and hence, it is possible to reduce a burden of arithmetic operations and to perform high speed computation.
In the case where the embedding degree is a multiple of 6, that is, 6h-th degree (h: a natural number) and the elliptic curve is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), there exists a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, vεFqh, and there exists a homomorphism from E′(Fqh) to E(Fq6h), that is, (x,y)→(xv1/3,yv1/2). Here, a new parameter v is a quadratic and cubic non-residue element which becomes 1 by an exponentiation of itself to the power of (q6h−1)/r.
Therefore, it is possible to perform the elliptic curve addition on E(Fq6h) [r] substituting the elliptic curve addition on E′(Fqh) [r] which has a 1⁄6 extension degree thus enabling to reduce a burden of arithmetic operations and to perform high speed computation.
In particular, since there exist rational points T′,Q′ on E′(Fqh) where relations T=(xT, yT)=(xT′v1/3, yT′v1/2), Q=(xQ, yQ)=(xQ′v1/3, yQ′v1/2) hold with respect to rational points T,QεE(Fq6h), first, λT,T, the gradient of lT,T(x,y) and λT,Q, the gradient of lT,Q(xp,yp) are transformed as follows.




  
    
      
        
          λ
          
            T
            ,
            T
          
        
        =
        
          
            
              3
              ⁢
              
                x
                T
                2
              
            
            
              2
              ⁢
              
                y
                T
              
            
          
          =
          
            
              
                3
                ⁢
                
                  
                    x
                    ′
                  
                  T
                  2
                
                ⁢
                
                  v
                  
                    2
                    3
                  
                
              
              
                2
                ⁢
                
                  y
                  T
                  ′
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    T
                    ′
                  
                
              
              ⁢
              
                v
                
                  1
                  6
                
              
            
          
        
      
    
    
      
        [
        E44
        ]
      
    
  
  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            
              
                y
                Q
              
              -
              
                y
                T
              
            
            
              
                x
                Q
              
              -
              
                x
                T
              
            
          
          =
          
            
              
                
                  (
                  
                    
                      y
                      Q
                      ′
                    
                    -
                    
                      y
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
              
                
                  (
                  
                    
                      x
                      Q
                      ′
                    
                    -
                    
                      x
                      T
                      ′
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    3
                  
                
              
            
            =
            
              
                λ
                
                  
                    T
                    ′
                  
                  ,
                  
                    Q
                    ′
                  
                
              
              ⁢
              
                
                  v
                  
                    1
                    6
                  
                
                ⁢
                
                  

                
                (
                
                  T
                  ≠
                  Q
                
                )
              
            
          
        
      
    
    
      
        [
        E45
        ]
      
    
  




The gradient λT,Q becomes as follows from [E44] and [E45], irrespective of T=Q, T≠Q.




  
    
      
        
          λ
          
            T
            ,
            Q
          
        
        =
        
          
            λ
            
              
                T
                ′
              
              ,
              
                Q
                ′
              
            
          
          ⁢
          
            
          
          ⁢
          
            v
            
              1
              6
            
          
        
      
    
    
      
        [
        E46
        ]
      
    
  




Therefore, from [E46], the equation [E28] becomes as follows.




  
    
      
        
          
            
              
                
                  l
                  
                    T
                    ,
                    Q
                  
                
                ⁡
                
                  (
                  
                    
                      x
                      P
                    
                    ,
                    
                      y
                      P
                    
                  
                  )
                
              
              =
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        x
                        T
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      T
                      ,
                      Q
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      p
                    
                    -
                    
                      y
                      T
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
                
                  
                    (
                    
                      
                        x
                        P
                      
                      -
                      
                        
                          x
                          
                            T
                            
                              ′
                              ⁢
                              
                                
                              
                            
                          
                        
                        ⁢
                        
                          v
                          
                            1
                            3
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    λ
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        
                          ′
                          ⁢
                          
                            
                          
                        
                      
                    
                  
                  ⁢
                  
                    v
                    
                      1
                      6
                    
                  
                
                -
                
                  (
                  
                    
                      y
                      P
                    
                    -
                    
                      
                        y
                        
                          T
                          
                            ′
                            ⁢
                            
                              
                            
                          
                        
                      
                      ⁢
                      
                        v
                        
                          1
                          2
                        
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
                
                  (
                  
                    
                      
                        (
                        
                          
                            
                              x
                              P
                            
                            ⁢
                            
                              v
                              
                                -
                                
                                  1
                                  3
                                
                              
                            
                          
                          -
                          
                            x
                            
                              T
                              ′
                            
                          
                        
                        )
                      
                      ⁢
                      
                        λ
                        
                          
                            T
                            ′
                          
                          ,
                          
                            Q
                            ′
                          
                        
                      
                    
                    -
                    
                      (
                      
                        
                          
                            y
                            P
                          
                          ⁢
                          
                            v
                            
                              -
                              
                                1
                                2
                              
                            
                          
                        
                        -
                        
                          y
                          T
                          ′
                        
                      
                      )
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    l
                    
                      
                        T
                        ′
                      
                      ,
                      
                        Q
                        ′
                      
                    
                  
                  (
                  
                    
                      
                        x
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            1
                            3
                          
                        
                      
                    
                    ,
                    
                      
                        y
                        P
                      
                      ⁢
                      
                        v
                        
                          -
                          
                            1
                            2
                          
                        
                      
                    
                  
                  )
                
                ⁢
                
                  v
                  
                    1
                    2
                  
                
              
            
          
        
      
    
    
      
        [
        E47
        ]
      
    
  




Since the final part of v1/2 in the equation [E47] becomes 1 at the final exponentiation to the power of (q6h−1)/r in the pairing computation, in calculating the lT,Q(xp,yp) in the Miller's algorithm, it is enough to calculate lT′,Q′(xpv−1/3, ypv−1/2) on the proper subfield instead of calculating lT,Q(xp,yp), thus enabling to reduce a burden of arithmetic operations significantly.
That is, the elliptic curve equation is given as y2=x3+a, aεFq (q: a power of a prime number greater than 3), letting the embedding degree be 6h-th degree (h: a natural number), E be an additive group constituted of rational points on the pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, φq be the Frobenius endomorphism, and using

G1=E[r]∩Ker(φq−[1])

G2=E[r]∩Ker(φq−[q])

the Ate pairing e is defined as a non-degenerate bilinear map, e: G2×G1→F*q6h/(F*q6h)r, and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using Miller function fs,Q(•), in computing the Ate pairing e(Q,P) as





  
    
      
        
          e
          ⁡
          
            (
            
              Q
              ,
              P
            
            )
          
        
        =
        
          
            
              f
              
                s
                ,
                Q
              
            
            ⁡
            
              (
              P
              )
            
          
          
            
              
                q
                
                  6
                  ⁢
                  h
                
              
              -
              1
            
            r
          
        
      
    
    
      
        [
        E48
        ]
      
    
  




the arithmetic operations of the rational function lT,T(xp,yp) and lT,Q(xp,yp) which are required to derive fs,Q(P) are performed on the proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to the power of (q6h−1)/r, as arithmetic operations lT′,T′(xpv−1/3,ypv−1/2) and lT′,Q′(xpv−1/3,ypv−1/2), and hence, it is possible to reduce a burden of arithmetic operations and to perform high speed computation.
In what follows, a pairing computation program in the case of extension degree of 12 which is the least common multiple of even degree, a degree of a multiple of 3, a degree of a multiple of 4, and a degree of a multiple of 6 is explained. FIG. 2 is Miller's algorithm corresponding to FIG. 1. Here, q, which has been a power of a prime number greater than 3, is assumed to be a prime number P.
In a pairing computation device constituted of an electronic computer, such as an authentication server or the like, a pairing computation is performed by a program based on a flowchart shown in FIG. 3, thus realizing a large scale digital group signature.
Here, the electronic computer which performs a pairing computation is the electronic computer which makes executable the pairing computation program integrating a program based on the Miller's algorithm shown in FIG. 2. As shown in FIG. 4, the electronic computer 10 includes: CPU 11 which performs arithmetic operations; memory unit 12 such as a hard disk or the like which stores various programs such as a pairing computation program and data used in the pairing computation program; and memory unit 13 constituted of RAM or the like which as well as expands the pairing computation program and makes it executable, temporally stores data generated in the course of execution of the pairing computation program. In FIG. 4, numeral 14 is a bus. In addition, the electronic computer 10 further includes a non-transitory recording medium drive not shown in the drawing, and hence, it is possible to allow memory unit 12 to store the non-transitory recording medium which records the pairing computation program via the non-transitory recording medium drive.
In the digital group signature technique which uses the pairing computation program according to the present invention, the electronic computer 10 makes it possible to connect to telecommunication lines 20 such as the Internet and to receive signature data of the digital group signature transmitted from client device 30 which is connected to the telecommunication lines 20. In FIG. 4, numeral 15 is an input output control part of the electronic computer 10. The electronic computer 10 temporally stores signature data of the digital group signature transmitted from the client device 30 in the memory unit 13.
The CPU 11 of the electronic computer 10, when signature data of the digital group signature is transmitted from the client device 30, temporally stores the transmitted signature data in the memory unit 13 and starts a pairing computation program (step S1).
By means of the started pairing computation program, the CPU 11 of the electronic computer 10 functions as an input unit, receives inputs of s=t−1 based on a trace t of the Frobenius endomorphism φp, P which satisfies PεE(Fp), and Q′ which satisfies Q′εE′(Fp2) (step S2). The signature data temporally stored in the memory unit 13 is any one of P and Q′ and is generally P. Q′ and the trace t of the Frobenius endomorphism φp or s=t−1 are stored in the predetermined address of the memory unit 13, and hence, they are read out from the predetermined address and inputted.
Next, the CPU 11 of the electronic computer 10 functions as a substitution unit and performs substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp) (step S3). Here, a value of v−1/3 and a value of v−1/2 are know and are store in the predetermined address in the memory unit 12 or the memory unit 13 and hence calculation of xpv−1/3 and ypv−1/2 are performed by reading them out from the predetermined address and inputting them.
In addition, these substitution operations, since the embedding degree is 12-th degree and a degree of a multiple of 6, are n←xpv−1/3, m←ypv−1/2, however in the case of the embedding degree being an even number degree, substitution operations are n←xpv−1, m←ypv−3/2, in the case of the embedding degree being a multiple of 3, substitution operations are n←xpv−1/3, m←ypv−1/2, and in the case of the embedding degree being a multiple of 4, substitution operations are n←xpv−1/2, m←ypv−3/4.
Next, the CPU 11 of the electronic computer 10 functions as the substitution unit and performs substitution operations f←1, T′←Q′ (step S4). This is a so-called initial value setting.
Next, the CPU 11 of the electronic computer 10 calculates a value of f, by setting the number of bits of s in binary representation as an initial value of i, and by decrementing i by 1 for each iteration of the following arithmetic operations (step S5). Here, as described above, s[i] denotes a value of the i-th bit counted from the least bit of s in binary representation.
In step S5, to be more specific, the CPU 11 of the electronic computer 10 functions as a first arithmetic operation unit, calculates lT,T(xp,yp)v−1/2 as lT′,T′(n,m) and performs a substitution operation of f←f2·lT′,T′(n,m) instead of calculating a value of lT,T(xp,yp) substituted by coordinates (xp,yp) of the rational point P in the left hand side of the tangent line lT,T(x,y)=0 passing through the rational point T, then next, the CPU 11 of the electronic computer 10 functions as a second arithmetic operation unit and performs a substitution operation of T′←2T′.
Further, in the case where a value of s[i] which is the i-th bit of s in binary representation is 1, the CPU 11 of the electronic computer 10 functions as the third arithmetic operation unit, calculates lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m) and performs a substitution operation of f←f·lT′,Q′(n,m) instead of calculating a value lT,Q(xp,yp) substituted by coordinates (xp,yp) of the rational point P in the left hand side of the straight line lT,Q(x,y)=0 passing through the rational points T and Q, then next, the CPU 11 of the electronic computer 10 functions as the fourth arithmetic operation unit and performs a substitution operation of T′←T′+Q′.
After calculation in step S5, the CPU 11 of the electronic computer 10 temporally stores a calculation result f in the memory unit 13 and finishes the calculation of the Miller function (step S6).
Next, the CPU 11 of the electronic computer 10 performs an exponentiation to the power of (q12−1)/r using the calculation result of f, and temporally stores the result as a result of pairing computation in the memory unit 13 (step S7).
The CPU 11 of the electronic computer 10 iterates the above described pairing computation the number of times corresponding to the number of invalidated users and finishes the pairing computation program.
In the case where the electronic computer 10 is an authentication server, after that, the CPU 11 of the electronic computer 10 performs authentication using the result of the pairing computation.
In this way, in a pairing computation device constituted of an electronic computer which performs a pairing computation, in the case where an embedding degree is 12, by performing pairing computation as an elliptic curve addition on E′(Fp2) [r] instead of an elliptic curve addition on E(Fp12) [r], it is possible to reduce the required time based on Miller's algorithm by approximately 30%, thus enabling to perform a pairing computation at higher speed.
Therefore, it is possible to utilize the pairing computation for the digital group signature with more members and to utilize the digital group signature in a broader range.
In the embodiment according to the present invention, the explanation about pairing computation used in digital group signature has been made, however, a pairing computation device, a pairing computation method and a non-transitory recording medium recording a pairing computation program are not limited to the case where they are used for digital group signature and are applicable to appropriate pairing computation as needed.

## Claims

1. A pairing computation device, wherein an elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq, q is a prime number greater than 3, letting an embedding degree be 2h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q2h/(F*q2h)r,


letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    2
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E1
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on a proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3 using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q2h−1)/r,
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1, m←ypv−3/2 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−3/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′,T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of a predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−3/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q2h−1)/r.

2. A pairing computation device, wherein an elliptic curve equation is given as y2=x3+a, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 3h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q])


as a non-degenerate bilinear map,


e:G2×G1→F*q3h/(F*q3h)r 


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    3
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E2
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on a proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1 using v εFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q3h−1 )/r,
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′, T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of a predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q3h−1)/r.

3. A pairing computation device, wherein an elliptic curve equation is given as y2=x3+ax, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 4h-th degree, h is a natural number, qh−1 be dividable by 4, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q])


as a non-degenerate bilinear map


e:G2×G1→F*q4h/(F*q4h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    4
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E3
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on a proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x using v εFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q4h−1)/r,
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1/2, m←ypv−3/4 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−3/4 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′, T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of a predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−3/4 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q4h−1)/r.

4. A pairing computation device, wherein an elliptic curve equation is given as y2=x3+a, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 6h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q6h/(F*q6h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P) is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    6
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E4
        ]
      
    
  




the pairing computation device characterized by comprising: in order to perform an arithmetic operation of a rational function which is required to derive fs,Q(P) on a proper subfield by a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1 using v εFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q6h−1)/r,
an input unit which is configured to receive inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh);
a substitution unit which is configured to perform substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp) and to perform substitution operations f←1, T′←Q′;
a first arithmetic operation unit which is configured to calculate lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and to perform a substitution operation f←f2·lT′,T′(n,m);
a second arithmetic operation unit which is configured to perform a substitution operation T′←2T′;
a third arithmetic operation unit which is configured, when a value of a predetermined bit of s in binary representation is 1, to calculate lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and to perform a substitution operation f←f·lT′,Q′(n,m);
a fourth arithmetic operation unit which is configured, when a value of the predetermined bit of s in binary representation is 1, to perform a substitution operation T′←T′+Q′; and
a fifth arithmetic operation unit which is configured to perform an exponentiation operation to the power of (q6h−1)/r.

5. A pairing computation method, wherein an elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq, q is a power of a prime number greater than 3, letting an embedding degree be 2h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q2h/(F*q2h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    2
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E5
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q2h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3, and by comprising:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1, m←ypv−3/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q2h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

6. A pairing computation method, wherein an elliptic curve equation is given as y2=x3+a, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 3h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q3h/(F*q3h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as,




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    3
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E6
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q3h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by comprising:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q3h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

7. A pairing computation method, wherein an elliptic curve equation is given as y2=x3+ax, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 4h-th degree, h is a natural number, qh−1 be dividable by 4, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q4h/(F*q4h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    4
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E7
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q4h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x, and by comprising:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/2, m←ypv−3/4 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/4 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/4 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q4h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

8. A pairing computation method, wherein an elliptic curve equation is given as y2=x3+a, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 6h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q6h/(Fq6h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), the Ate pairing e(Q,P), by using an electronic computer, is computed as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    6
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        [
        E8
        ]
      
    
  




the pairing computation method characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q6h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by comprising:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q6h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

9. A non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+ax+b, aεFq, bεFq, q is a power of a prime number greater than 3, letting an embedding degree be 2h-th degree, h is a natural number, E bean additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq2h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q2h/(F*q2h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    2
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        
      
    
  



the pairing computation program characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q2h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−2x+bv−3, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1, m←ypv−3/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q2h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

10. A non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+a, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 3h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq3h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q3h/(F*q3h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as




  
    
      e
      ⁡
      
        (
        
          Q
          ,
          P
        
        )
      
    
    =
    
      
        
          f
          
            s
            ,
            Q
          
        
        ⁡
        
          (
          P
          )
        
      
      
        
          
            q
            
              3
              ⁢
              h
            
          
          -
          1
        
        r
      
    
  
  ,




the pairing computation program characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q3h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m←ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m) by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q3h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

11. A non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+ax, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 4h-th degree, h is a natural number, qh−1 be dividable by 4, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq4h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q45/(F*q4h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as




  
    
      
        
          
            e
            ⁡
            
              (
              
                Q
                ,
                P
              
              )
            
          
          =
          
            
              
                f
                
                  s
                  ,
                  Q
                
              
              ⁡
              
                (
                P
                )
              
            
            
              
                
                  q
                  
                    4
                    ⁢
                    h
                  
                
                -
                1
              
              r
            
          
        
        ,
      
    
    
      
        
      
    
  



the pairing computation program characterized by using vεFqh which is a quadratic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q4h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1x, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/2, m←ypv−3/4 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−3/4 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−3/4 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q4h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.

12. A non-transitory recording medium recording a pairing computation program, wherein an elliptic curve equation is given as y2=x3+a, aεFq, q is a power of a prime number greater than 3, letting an embedding degree be 6h-th degree, h is a natural number, E be an additive group constituted of rational points on a pairing-enabled elliptic curve defined over Fq6h, E[r] be a set of rational points of prime order r, and φq be a Frobenius endomorphism, an Ate pairing e is defined, using


G1=E[r]∩Ker(φq−[1])



G2=E[r]∩Ker(φq−[q]),


as a non-degenerate bilinear map


e:G2×G1→F*q6h/(F*q6h)r,


and letting PεG1, QεG2, and letting s=t−1 using a trace t of the Frobenius endomorphism φq, and using a Miller function fs,Q(•), an electronic computer is caused to compute the Ate pairing e(Q,P) as




  
    
      e
      ⁡
      
        (
        
          Q
          ,
          P
        
        )
      
    
    =
    
      
        
          f
          
            s
            ,
            Q
          
        
        ⁡
        
          (
          P
          )
        
      
      
        
          
            q
            
              6
              ⁢
              h
            
          
          -
          1
        
        r
      
    
  
  ,



the pairing computation program characterized by using vεFqh which is a quadratic and cubic non-residue element and becomes 1 by an exponentiation of fs,Q(P) to a power of (q6h−1)/r and a twist curve of E(Fqh), that is, E′(Fqh):y2=x3+av−1, and by causing the electronic computer to perform:
a step of receiving inputs of said s=t−1, P which satisfies PεE(Fq), and Q′ which satisfies Q′εE′(Fqh), by causing the electronic computer to function as an input unit;
a step of performing substitution operations n←xpv−1/3, m→ypv−1/2 with respect to coordinates of the rational point P (xp,yp), by causing the electronic computer to function as a substitution unit;
a step of performing substitution operations f←1, T′←Q′, by causing the electronic computer to function as a substitution unit;
a step of calculating lT,T(xp,yp)v−1/2 as lT′,T′(n,m), instead of calculating lT,T(xp,yp) which is a value substituted by coordinates (xp,yp) of the rational point P in the left hand side of an equation lT,T(x,y)=0 of a tangent line passing through the rational point T, and performing a substitution operation f←f2·lT′,T′(n,m), by causing the electronic computer to function as a first arithmetic operation unit;
a step of performing a substitution operation T′←2T′, by causing the electronic computer to function as a second arithmetic operation unit;
a step of, when a value of a predetermined bit of s in binary representation is 1, calculating lT,Q(xp,yp)v−1/2 as lT′,Q′(n,m), instead of calculating lT,Q(xp,yp) which is a value substituted by coordinates of the rational point P in the left hand side of an equation lT,Q(x,y)=0 of a straight line passing through the rational points T and Q, and performing a substitution operation f←f·lT′,Q′(n,m), by causing the electronic computer to function as a third arithmetic operation unit;
a step of, when a value of the predetermined bit of s in binary representation is 1, performing a substitution operation T′←T′+Q′, by causing the electronic computer to function as a fourth arithmetic operation unit; and
a step of performing an exponentiation operation to the power of (q6h−1)/r, by causing the electronic computer to function as a fifth arithmetic operation unit.