Digital signature method and apparatus

Information

  • Patent Application
  • 20170250819
  • Publication Number
    20170250819
  • Date Filed
    February 24, 2017
    7 years ago
  • Date Published
    August 31, 2017
    7 years ago
Abstract
A method for signing and subsequently verifying a digital message, including the following steps: generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x]; generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y]; producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively; producing a secret isomorphism from the first finite field to the second finite field; producing and publishing a public key that depends on F(y); producing a private key that depends on the secret isomorphism; producing a message digest by applying a hash function to the digital message and the public key; producing a digital signature using the message digest and the private key; and performing a verification procedure utilizing the digital signature and the public key.
Description
FIELD OF THE INVENTION

This invention relates to the field of cryptography and, more particularly, to a public key digital signature technique and system.


BACKGROUND OF THE INVENTION

Public key digital signatures are important for secure exchange of information between plural parties, for example between computers or mobile devices, between a smart card and a terminal, etc.


An earlier digital signature and authentication method and apparatus was described in U.S. Pat. No. 7,308,097, assigned to the same assignee as the present Application. Reference can also be made to “NTRUSign: Digital Signatures Using the NTRU Lattice”, J. Hoffstein, N. Howgrave Graham, J. Pipher, J. Silverman, and W. Whyte, Topics In Cryptology-CT-RSA 2003, Lecture Notes in Computer Science, Vol. 2612, Springer, Berlin, 2003.


The signing technique in the '097 Patent uses a mixing system based on multiplication in a ring and reduction modulo an ideal q in that ring; while the verification technique uses special properties of products of elements whose validity depends on elementary probability theory. The security of the identification/digital signature scheme comes from the interaction of reduction modulo q and the difficulty of forming products with special properties. In an embodiment of the digital signature scheme of the '097 Patent, the security also relies on the experimentally observed fact that for most lattices, it is very difficult to find a vector whose length is only a little bit longer than the shortest vector, and it is also difficult to find a lattice vector that is quite close to a randomly chosen nonlattice vector.


An improvement over the technique of the '092 Patent, which had reduced complexity and computational requirements for key generation and signing, was disclosed in copending U.S. patent application Ser. No. 14/544,426, assigned to the same assignee as the present Application, published as U.S. Patent Application Publication No. US2015/0229478, incorporated herein by reference. In a form of that invention, sometimes referred to a “pqNTRUSign” (mark of Security Innovation, Inc.), a method is set forth for signing and subsequently verifying a digital message, comprising the following steps implemented using at least one processor-based subsystem: selecting parameters including an integer q and a relatively smaller integer p that is coprime with q; generating random polynomial f relating to p and random polynomial g relating to q; producing a public key that includes h, where h is equal to a product that can be derived using g and the inverse off mod q; producing a private key from which f and g can be derived; storing the private key and publishing the public key; producing a message digest by applying a hash function to the digital message; producing a digital signature using the message digest and the private key; and performing a verification procedure utilizing the digital signature and the public key to determine whether the signature is valid.


It is among the objectives hereof to devise a digital signature method and system that has advantages over existing digital signature techniques, including those of the type described hereinabove.


SUMMARY OF THE INVENTION

The present invention utilizes, inter alia, a secret isomorphism between two finite fields to create a signature scheme.


In accordance with an embodiment of the invention, a method is set forth for signing and subsequently verifying a digital message, comprising the following steps implemented using at least one processor-based subsystem: generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x]; generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y]; producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively; producing a secret isomorphism from the first finite field to the second finite field; producing and publishing a public key that depends on F(y); producing a private key that depends on said secret isomorphism; producing a message digest by applying a hash function to the digital message and the public key; producing a digital signature using the message digest and the private key; and performing a verification procedure utilizing the digital signature and the public key to determine whether the signature is valid.


In a form of the embodiment just set forth, the first and second finite fields are designated, respectively, as x-space and y-space, and a further step comprises generating a specified lattice in x-space and using said isomorphism to generate a corresponding lattice in y-space, and the step of producing said digital signature includes initially producing a signature in the x-space lattice and then producing the digital signature in the corresponding y-space lattice.


In accordance with another embodiment of the invention, a method is set forth signing and sending a digital message, including the following steps implemented using at least one processor-based subsystem: generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x]; generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y]; producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively; producing a secret isomorphism from the first finite field to the second finite field; producing and publishing a public key that depends on F(y); producing a private key that depends on said secret isomorphism; producing a message digest by applying a hash function to the digital message and the public key; producing a digital signature using the message digest and the private key; and transmitting the digital signature.


In accordance with a further form of the invention, a system is set forth for signing and subsequently verifying a digital message, including: at least one processor-based subsystem that is programmed with instructions that cause the at least one processor subsystem to implement the following steps: generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x]; generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y]; producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively; producing a secret isomorphism from the first finite field to the second finite field; producing and publishing a public key that depends on F(y); producing a private key that depends on said secret isomorphism; producing a message digest by applying a hash function to the digital message and the public key; producing a digital signature using the message digest and the private key; and performing a verification procedure utilizing the digital signature and the public key to determine whether the signature is valid.


As described herein (see e.g. Appendix I), for reasons including the non-linear nature of the homomorphic encryption map, the effectiveness of attacks on the signature scheme of the invention is lessened, which enhances security and allows for smaller signatures and improved operating characteristics compared to existing signature schemes.


Further features and advantages of the invention will become more readily apparent from the following detailed description when taken in conjunction with the accompanying drawings.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a block diagram of a system that can be used in practicing embodiments of the invention.



FIG. 2 shows an example of a first finite field (also sometimes designated as x-space) and a second finite field (also sometimes designated as y-space), that is useful in understanding operation of the invention.



FIG. 3 again illustrates the finite fields of the FIG. 2 example, and also illustrates an exemplary isomorphism between polynomials of the respective finite fields, and shows some of the mappings between polynomials of the respective finite fields for an example in which the isomorphism is defined by a particular function that maps from polynomials of the first finite field (x-space) to polynomials of the second finite field (y-space).



FIG. 4A shows, on the left, a table of polynomial values in x-space, evaluated at an arbitrary value of x and, on the right, a table of corresponding (that is, mapped in accordance with the exemplary isomorphism) polynomial values in y-space. The tables illustrate the scrambling of the values that occurs in y-space.



FIG. 4B shows tables similar to those of FIG. 4A, but for a different group of polynomials. Again, the tables illustrate the scrambling of values that occurs in y-space.



FIG. 5 is a flow diagram of a public key digital signature technique which, when taken with the subsidiary flow diagrams referred to therein, can be used in implementing embodiments of the invention.



FIG. 6 is a flow diagram, in accordance with an embodiment hereof, of a routine for key generation.



FIG. 7 is a flow diagram, in accordance with an embodiment hereof, of a routine for signing a digital message.



FIG. 8 is a flow diagram, in accordance with an embodiment hereof, of a routine for verification of a digital signature.



FIG. 9, which includes FIGS. 9A and 9B placed one below another, is a flow diagram, in accordance with an embodiment hereof, of a further, and more detailed, routine for signing and subsequently verifying a digital message.





DETAILED DESCRIPTION


FIG. 1 is a block diagram of an example of a system that can be used in practicing embodiments of the invention. Two processor-based subsystems 105 and 155 are shown as being in communication over a channel 50, which may be, for example, any wired or wireless communication channel such as a telephone or internet communication channel in, for example, a cloud-based system. In the example hereof, the channel can be considered a secure or an insecure channel. The subsystem 105 includes processor 110 and the subsystem 155 includes processor 160. The subsystems can typically comprise mobile devices, computers, or terminals. When programmed in the manner to be described, the processors 110 and 160 and their associated circuits can be used to implement an embodiment of the invention and to practice an embodiment of the method of the invention. The processors 110 and 160 may each be any suitable processor, for example an electronic digital processor or microprocessor. It will be understood that any general purpose or special purpose processor, or other machine or circuitry that can perform the functions described herein, electronically, optically, or by other means, can be utilized. The subsystem 105 will typically include memories 123, clock and timing circuitry 121, input/output functions 118 and display 125, which may all be of conventional types. Inputs can include a touchscreen/keyboard input as represented at 103. Communication is via transceiver 135, which may comprise any suitable device for communicating signals.


The subsystem 155 in this illustrative embodiment can have a similar configuration to that of subsystem 105. The processor 160 has associated input/output circuitry 164, memories 168, clock and timing circuitry 173, and a display 176. Inputs include a touchscreen/keyboard 155. Communication of subsystem 155 with the outside world is via transceiver 162.


Co-inventors hereof, Jeffrey Hoffstein and Joseph Silverman, are inventors of copending PCT International Application No. PCT/US2016/041598, published as PCT Publication No. WO/2017/008043, entitled “Homomorphic Encryption”, and said PCT Publication is incorporated herein by reference. The '598 Application discloses systems, methods, and computer-readable storage devices storing instructions for homomorphic encryption via finite ring isomorphisms that includes constructing an isomorphism and an inverse isomorphism and using these for encryption that is homomorphic, in that certain operations can be performed on encrypted data without the need for first decrypting the data.


Consider f=Z/qZ[x]/(f(x)), a finite field Fqn of order qn. f(x) and F(y) define two copies of fqn, and Z/qZ[x]/(f(x))≅Z/qZ[y]/(F(y)) is a finite field isomorphism under a secret mapping x→φ(y).



FIG. 2 illustrates an example of the first finite field Fq(x)/f(x) (represented in conjunction with the lefthand box of the Figure) and the second finite field fq(y)/F(y) (represented in conjunction with the righthand box of the Figure). [In mathematical language, the first finite field is the quotient of the ring of polynomials in x, with coefficients mod q, by the ideal generated by the irreducible polynomial f(x), and the second finite field is the quotient of the ring of polynomials in y, with coefficients mod q, by the ideal generated by the irreducible polynomial F(y).] In this example, n=5 and there are 3125 (equals 55) polynomials listed in each box in an indexed order, with increasing coefficient values. In this example, the polynomial f(x) is x5+x4+4x3+x2+4x+1, and the polynomial F(y) is y5+y4+3y3+2y2+2y+4.


In the illustrative example of FIG. 3, which again shows the finite fields, the secret isomorphism between the finite fields is selected to be represented by the polynomial 4(y)=4y4+4y2+4y+3; that is, x in the finite field Fq(x)/f(x) (which is designated as x-space) maps to 4y4+4y2+4y+3 in the finite field Fq(y)/F(y) (which is designated as y-space). The arrows shown between x-space and y-space in FIG. 3 illustrate some of the mappings for this isomorphism.


As an example, the following calculations show how the index #726 (polynomial 1+0x+4x2+0x3+1x4) of the x-space finite field, maps to the index #837 (polynomial 2+2y+3y2+1y3+1y4) of the y-space finite field. Using 3+4y+4y2+4y4 for x in the expression 1+4x2+x4 gives 1+4(3+4y+4y2+4y4)2+(3+4y+4y2+4y4)4. When this is expanded and coefficients are reduced mod 5, the result is the long polynomial:






A(y)=3+3y+y2+4y3+y4+4y5+4y6+y7+y9+4y10+2y11+4y12+4y13+4y14+y16


and this can be written as:






A(y)=(1+y)2(4+2y+2y2+3y3+y4+y5)(4+4y+2y2+3y3+2y4+4y5+2y6+2y7+2y8+y9)+2+2y+3y2+y3+y5 mod 5.


This expression [A(y)] must be modded by F(y), which, in this example, is y5+y4+3y3+2y2+2y+4, so:






A(y)=(4+4y+2y2+3y3+2y4+4y5+2y6+2y7+2y8+y9)(1+y)2(F(y))+2+2y+3y2+y3+y4 mod 5.


So, in other words, when A(y) is divided by F(y) mod 5, the remainder is 2+2y+3y2+y3+y4, which, in the example of FIG. 3, is the polynomial index #837 in y-space. Ergo, polynomial index #726 in x-space maps to polynomial index #837 in y-space.


The foregoing also demonstrates that 1+4x2+x4 maps to 1+4 (image of x)2+(image of x)4 which demonstrates the homomorphic mapping property; that is, additive and multiplicative structure is preserved under the mapping.


Also, x maps to image of x=3+4y+4y2+4y4, which means that f(x)=x5+x4+4x3+x2+4x+1 maps to F(y)=(3+4y+4y2+4y4)5+(3+4y+4y2+4y4)4+4(3+4y+4y2+4y4)3+(3+4y+4y2+4 y4)2+4(3+4y+4y2+4y4)+1.


Since f(x) is zero in x-space (since Fq[x] is modded by f(x)), it should map to zero in y-space. The following demonstrates that it does.





Expanding F(y)mod 5 gives 4+4y+y3+3y4+2y5+4y6+2y7+y8+4y9+y10+2y11+4y13+y16+4y20


and if this is factored mod 5, one gets





4(4+2y+2y2+3y3+y4+y5)(4+2y2+2y3+4y4+y5)(1+3y+2y3+y4+3y5+2y6+2y7+y8+y10).






F(y)=4+2y+2y2+3y3+y4+y5





So A(y)=4(4+2y2+2y3+4y4+y5)(1+3y+2y3+y4+3y5+2y6+2y7+y8+y10)F(y)mod 5.


Stated another way, when the polynomial that is mapped to [called A(y)] is divided by F(y), the remainder is zero, so zero maps to zero, as it should.


The tables of FIGS. 4A and 4B illustrate the property hereof that the isomnrphism doesn't “see” the sizes of the coefficients of the x polynomials, which is an indication that images in y-space are not revealing of relevant information about x-space. FIGS. 4A and 4B show two examples of the principle. Both examples involve the same finite fields and the same secret mapping isomorphism as in FIGS. 2 and 3.


In FIG. 4A, twenty five polynomials in x-space are defined by the exemplary general polynomial f(x)=a4x4+a3x3+2x2+2x+2, where a4 and a3 can each take on any of the five values 0, 1, 2, 3, or 4, so that there are 25 possible polynomials. In the FIG. 4A table on the left a4 varies from 0 to 4 on the horizontal axis and a3 varies from 0 to 4 on the vertical axis. The numerical value in each of the 25 boxes is the value of the polynomial for the particular box, evaluated at x=5. For example, the upper leftmost box of the table, where a4=0 and a3=0, accordingly is associated with the polynomial f(x)=2x2+2x+2 which, evaluated at x=5, is (2)(25)+(2)(5)+2=62. As another example, consider the lower rightmost box of the table, where a4=4 and a3=4, so the associated polynomial is f(x)=4x4+4x3+2x2+2x+2 which, evaluated at x=5, is 4(625)+4(125)+2(25)+2(5)+2=3062.


In the FIG. 4A table on the right, each of the 25 boxes is associated with the polynomial F(y) that is the y-space image of the x-space polynomial f(x) in the corresponding box of the table on the left. The number in each box is the value of the associated polynomial, evaluated at y=5.


In FIG. 4B, twenty five polynomials in x-space are defined by another exemplary polynomial f(x)=2x4+2x3+a2x2+a1x+2, where a2 and a1 can, again, take on any of the five values, 0, 1, 2, 3, or 4 so that there are 25 possible polynomials. In the FIG. 4B table on the left a2 varies from 0 to 4 on the horizontal axis and a1 varies from 0 to 4 on the vertical axis. As before, the numerical value in each of the 25 boxes is the value of the polynomial for the particular box, evaluated at x=5. In this case, for example, the upper leftmost box of the table, where a2=0 and a1=0, accordingly is associated with the polynomial f(x)=2x4+2x3+2 which, evaluated at x=5, is (2)(625)+(2)(125)+2=1502. As another example, consider the lower rightmost box of the table, where a2=4 and a1=4, so the associated polynomial is f(x)=2x4+2x3+4x2+4x+2 which, evaluated at x=5, is 2(625)+2(125)+4(25)+4(5)+2=1622.


In the FIG. 4B table on the right, each of the 25 boxes is again associated with the polynomial F(y) that is the y-space image of the x-space polynomial f(x) in the corresponding box of the table on the left. The number in each box is the value of the associated polynomial, evaluated at y=5.


Accordingly, FIGS. 4A and 4B illustrate the scrambling that occurs in y-space as a result of the defined process using the secret isomorphism. As seen in the examples of these Figures, although the values of the selected polynomials in x-space are generally ordered in their tables, the y-space polynomials have no discernable order in their respective tables.



FIG. 5 illustrates a basic procedure that can be utilized with a public key digital signature technique, and refers to routines illustrated by other referenced flow diagrams which describe features in accordance with an embodiment of the invention. Reference can also be made to Appendix I for further details of embodiments of the invention. The block 510 represents the generating of the public key and private key signals and data, and the publishing of the public key. The routine of an embodiment thereof is described in conjunction with the flow diagram of FIG. 6. In the present example, this operation can be performed, for example, at the processor-based subsystem 105 of FIG. 1. The public key information can be published; that is, made available to any member of the public or to any desired group to whom the private key holder desires to send the digital signatures. Typically, although not necessarily, the public key may be made available at a central public key library facility or website where a directory of public key holders and their public keys are maintained.


The block 550 represents a routine that can be employed (that is, in this example, by the user of processor-based subsystem 105 of FIG. 1) for signing the digital message. This routine, in accordance with an embodiment of the invention, is described in conjunction with the flow diagram of FIG. 7. In this example, the digital signature is then transmitted over the channel 50 (FIG. 1), although it will be understood that the digital signature can, if desired, be stored and subsequently retrieved and verified, depending on the application or task being undertaken.


The block 570 represents a routine that can be employed (that is, in this example, by the user of processor-based subsystem 155 of FIG. 1) for using, inter alia, the public key to implement a verification procedure to either accept or reject the digital signature. This routine, in accordance with an embodiment of the invention, is described in conjunction with the flow diagram of FIG. 8.



FIG. 6 is a flow diagram of a routine, represented generally by the block 510 of FIG. 5, in accordance with an embodiment of the invention, for implementing key generation. Reference can also be made to Appendix I. The block 605 represents the inputting of initial parameters as, for example, described in Appendix I. The block 610 represents generating of an irreducible monic polynomial f(x) of degree n in a ring Fq(x), and the block 620 represents generating of an irreducible monic polynomial F(y) of degree n in a ring Fq(y). Then, as represented by the block 630, first and second finite fields are produced as Fq(x)/f(x) and Fq(y)/F(y), respectively. Next, as represented by the block 640, an isomorphism is produced from the first finite field to the second finite field. This can be, for example, an isomorphism e.g. in the form of a matrix or a polynomial of the type illustrated in conjunction with FIG. 3 as the mapping polynomial φ, but of a magnitude commensurate with a desired degree of security. A public key that depends on F(y) can then be output (block 650) and, typically, published. A private key, which depends on the isomorphism, is stored (block 660) for subsequent use.



FIG. 7 is a flow diagram of a routine, represented by the block 550 of FIG. 5, in accordance with an embodiment of the invention, for implementing the signing of a digital message that is to be stored and/or transmitted. Reference can also be made to Appendix I. The block 750 represents the inputting of the digital message that is to be encoded. Then, as represented by the block 760, a message digest is produced by applying a hash function to the digital message and the public key. The encoded digital signature is then produced (block 770) using the message digest and the private key. The encoded digital signature can then be output (block 780). If desired, an optional rejection sampling technique, represented by loop 775, can be employed, as set forth in copending U.S. patent application Ser. No. 14/121,041 (published as U.S. Patent Application Publication No. US2015/0033025) to check a candidate signature to determine if it may leak information that could be used in a transcript of signatures to reveal information about the public key and, if so, modifying one or more parameters used in producing the candidate digital signature, the procedure being repeated until an acceptable signature is obtained.



FIG. 8 is a flow diagram of a routine, represented by the block 570 of FIG. 5, in accordance with an embodiment of the invention, for implementing verification to accept or reject the digital signature. The block 810 represents the receiving of the digital signature, and the block 820 represents the inputting of the public key. Then, the decision block 830 represents the determination of whether the collective components of these and other inputs satisfy predetermined size and lattice presence conditions. If not, the signature is rejected. If so, however determination is also made (decision block 840) as to whether a predetermined congruence condition is met. If so, the digital signature is considered verified.



FIG. 9 is a flow diagram of a routine in accordance with a further embodiment of the invention for encoding and decoding a digital message. Reference can be made to Appendix I for further details. The block 910 represents inputting parameters such as p, q, n, and other parameters as described in Appendix I. The finite fields Fq(x)/(f(x)) (x-space, as described above) and Fq(y)/(F(y)) (y-space, as described above) are produced, as represented by the block 920. Next, the block 930 represents the selection of secret a(x), b(x) in x-space with small coefficients, and the computation of a secret h(x)=(pa(x))−1b(x). (As used herein in this context, a small coefficient is an integer with absolute value less than c, where c is significantly less than q, an example being q=2048, c=4.) Then, as represented by the block 940, a secret isomorphism (such as φ in FIG. 3) is produced, mapping h(x) (in x-space) to H(y) (in y-space). A public key that includes H(y) can then be published (block 950). Then, as represented by the block 955, a lattice Lh is produced in x-space and a lattice LH is produced in y-space. Lh is the set of all pairs of polynomials (s(x), t(x)) of degree less than or equal to n−1 such that t(x) is equivalent to S(x)h(x) mod(f(x),q), and LH is the set of all pairs of polynomials (S(y), T(y)) of degree less than or equal to n−1 such that T(y) is equivalent to S(y)H(y) mod (F(y),q).


A message digest is obtained from a hash of the digital message and the public key (block 960). In this embodiment, the message digest is a 2n dimensional polynomial (d0, d1, . . . dn-1, e0, e1 . . . en-1), with coefficients chosen from the interval [−p/2, p/2]. The block 970 represents finding (s(x), t(x)) in Ln such that the 2n dimensional vector of coefficients from (s(x), t(x)) satisfies the correct bound on absolute values, and is congruent to a specified 2n dimensional vector determined by the message digest and the polynomials Ci(y). The coefficients of this linear combination are output as the digital signature. Then, as represented by the block 985, the digital signature is received and subject to a verification procedure. In this embodiment, a determination is made (decision block 995) as to whether the linear combination of Ci(y) is in LH and, also, whether the 2n dimensional vector of coefficients has absolute values in a predetermined range and is congruent to (d0, d1 . . . dn-1, e0, e1, . . . en-1)mod p. If both criteria are satisfied, the digital signature is verified. If not, it is rejected.


The invention has been described with reference to particular preferred embodiments, but variations within the spirit and scope of the invention will occur to those skilled in the art. For example, while a digital signature technique has been described, it will be understood that an authentication procedure of the challenge-response-verification type can alternatively be implemented, using the technique hereof and employing the challenge as the message to be signed. Also, it will be understood that coefficients of polynomials can alternatively be represented in other forms including, but not limited to, matrices.


APPENDIX I

There are a number of possible methods of using a secret isomorphism between two finite fields to create a signature scheme. In the following, a method we call “pqFF-Sign” is set forth.


The approach utilizes the existing transcript-secure signature scheme pq-NTRUSign and also utilizes a finite field isomorphism. The pq-NTRUSign signature scheme is employed in a finite field custom-characterqn, and then the signature is homomorphically mapped to a different copy of the field custom-characterqn. Verification is still possible due to the homomorphic property of the map, but various lattice attacks that may have been possible on pq-NTRUSign are blunted or eliminated due to the non-linear nature of the homomorphic encryption map.


n a dimension parameter.


q a prime (or prime power) greater than a constant times n.


p a small prime (often taken to be 2).


f(x)εcustom-characterq[x] a short irreducible monic polynomial of degree n.


F(y)εcustom-characterq[y] a random irreducible monic polynomial of degree n.


ψ(x)εF(y)εcustom-characterq[y]. The map xcustom-characterφ(y) induces an isomorphism custom-characterq[x]/(f(x))→custom-characterq[y]/(F(y)).


ψ(x)εcustom-characterq[x] the map ycustom-characterψ(x) induces the inverse isomorphism custom-characterq[y]/(F(y))→custom-characterq[x]/(f(x)).


a(x), b(x)εcustom-characterq[x] short irreducible monic polynomials of degree n.


h(x)εcustom-characterq[x]≡b(x)·(pa(x))−1 (mod q).


H(y)εcustom-characterq[y]≡h(φ(y)) (mod q, F(y)).


U is an n-by-n matrix, small entries, invertible mod q.


We will lift mod q polynomials to polynomials having integer coefficients in the range (−½q,½q]. Define rings and fields








f

=





[
x
]



(

f


(
x
)


)



,


F

=





[
y
]



(

F


(
y
)





,



f


,
q



=




q



[
x
]



(

f


(
x
)


)



,



F
,
q


=





q



[
y
]



(

F


(
y
)


)


.






And define lattices






L
h={(u,vcustom-characterf2:v≡h·u(mod q)},






L
H={(U,Vcustom-characterF2:V=H·U(mod q)}.


We use U to define polynomials c1(x), . . . , cn(x)εcustom-characterq[x] of degree less than n by







(





c
1



(
x
)








c
2



(
x
)













c
n



(
x
)





)





U

-
1




(



x





x
2











x
n




)





(


mod











q

,

f


(
x
)



)

.






For 1≦j≦n we let






C
j(y)=cj(φ(y))εcustom-characterF,q


be the corresponding polynomials in the y-field.


Security Assumption


As U ranges over matrices with small coefficients that are invertible modulo q, the coefficients of U−1 mod q are uniformly distributed.


The polynomials c1(x), c2(x), . . . , cn(x) form a basis for custom-characterf,q and C1(y), C2(y), . . . , Cn(y) form a basis for custom-characterF,q. Each Cj(y) is the image of the corresponding cj(x) under the isomorphism that sends xcustom-characterφ(y). This same isomorphism preserves the coefficients of linear combinations of the cj(x), that is,





Σαjcj(x)custom-characterΣαjCj(y).


A key property that the scheme is based on is the fact that as








U






(





c
1



(
x
)








c
2



(
x
)













c
n



(
x
)





)





(



x





x
2











x
n




)



(


mod











q

,

f


(
x
)



)



,




and the coefficients of U are small, that each xi, for 1≦i≦n is expressible as a linear combination of the cj(x) with small coefficients. From this it follows that any polynomial in x with small coefficients, of the form t(x)=Σi=1ntixi, with the ti small, can in turn be written as a polynomial in cj(x) with small coefficients.


We will find it convenient to use polynomials reduced modulo f(x), which will necessarily have degree less than or equal to n 1, that is, of the form r(x)=Σi=0n-1rixi, with the ri short. Recall that f(x)=xn+fn-1xn-1+ . . . +f1x±1, where fiε{1, 0, −1}. Consequently





±1≡f1x− . . . −xn(mod f(x)).


and thus







r


(
x
)







i
=
1

n




r
i




x
i






±


r
0



(



-

f
1



x

-

-

x
n


)



+




i
=
1


n
-
1





r
i



x
i







(


mod






(

f


(
x
)


)


,









where the ri′ are also short, is expressible as a short linear combination of the cj(x).


Importantly, as the product of any two short polynomials in x remains short, such a product will also be writeable as a short linear combination of the cj(x). The reverse does not hold: with high probability a random linear combination of cj(x) with small coefficients will not equal a polynomial in x with small coefficients.


It is this property, that a product of short linear combinations of the cj(x) that correspond to short polynomials in x can again be written as a short linear combination of the cj(x), that allows the signer to solve a congruential lattice problem in Lh (just as in pq-NTRUSign) and then map the corresponding solution, with the same coefficients, back to LH.


There are two main security concerns that determine parameters in pq-NTRUSign. One is the problem of recovering the private key from the public NTRU key, and the other is the problem of forgery. Of these, the one that has the biggest impact on parameter size is the public key to private key problem. This is because, to make rejection sampling efficient, the q needs to be chosen large compared to n. This makes the lattice problem somewhat easier and forces an increase in the size of n. The forgery problem requires smaller parameters to achieve the same security levels.


In this context there appear at first to be two NTRU-type problems: Recovering a(x), b(x) from h(x), and recovering the corresponding polynomials A(y), B(y) from H(y).


The h(x) is private, and only revealed if the underlying isomorphism is discovered, in which case the scheme is considered broken. So this lattice problem does not arise.


On the other hand, the H(y) is public, but the corresponding problem of recovering A(y), B(y) from H(y) is not a lattice reduction problem as A(y), B(y) are generic polynomials with coefficients mod q, and not short. And as they are not short, recovery of them would not lead to any advantage.


There is a lattice attack to recover the matrix U from the Cj(y), which would suffice to break the scheme, but the dimension of the lattice required to accomplish this is greater than n2.


For this reason, it appears that it will suffice to set parameters to avoid forgery attacks, which should allow for smaller signatures and better operating characteristics.


Private Information: f(x), a(x), b(x), h(x), c1(x), . . . , cn(x) and U.


Public Information: F(y), H(y) and C1(y), . . . , Cn(y).


Digital Document Hash: A pair of mod p vectors δ, ε obtained by applying a hash function to the document being signed:






δ=δ1, . . . ,δn)ε(−½p,½p]n,






ε=ε1, . . . ,εnε(−½p,½p]n,


Signature: A pair of vectors





δ=(δ1, . . . ,δn)ε(−½q,½q]n,





ε=(ε1, . . . ,εn)ε(−½q,½q]n,


Verification: A signature on the document hash (δ,ε) is valid if it satisfies the following three conditions:

    • (1) δiδi (mod p) and εiεi (mod p) for all 1≦i≦n.
    • (2) |δi|≦½q−B and |εi|≦½q−B for all 1≦i≦n.
    • (3) Let






S(y)=δ1C1(y)+ . . . +δnCn(y),





and






T(y)=ε1C1(y)+ . . . +εnCn(y).


Then (S,T)εLH, i.e.,






T(y)=S(y)H(y)(mod q,F(y)).


Here B is a fixed small integer used to enable rejection sampling.


Signatures are created as in pq-NTRUSign working in the ring custom-characterf=custom-character[x]/(f(x)), with one change. Rather than creating polynomials with small coefficients relative to the standard basis 1, x, . . . , xn-1, we instead create polynomials with small coefficients relative to the basis c1(x), . . . , cn(x).


Step 1: Choose δj at random mod q such that q/2<δj≦q/2 and δjδj (mod p) and set








s
0



(
x
)


=




j
=
1

n




δ
j





c
j



(
x
)


.







Step 2: Define t0(x) by






t
0(x)≡s0(x)h(x)(mod q)


and write








t
0



(
x
)


=




i
=
0


n
-
1





t
i




x
i

.







Then (s0(x),t0(x))εLh.


Step 3: Write








t
0



(
x
)


=




j
=
1

n




η
j




c
j



(
x
)








for some η1, . . . , ηn.


To accomplish this, as described previously, write








t
0



(
x
)


=





i
=
0


n
-
1





t
i



x
i








i
=
1

n




t
i




x
i







(

mod






f
(
x
)


)








as described previously, and set





1, . . . ,ηn)≡(t′1, . . . ,t′n)U(q),


and select representatives for the ηj in the interval (−q/2/q/2]. The ηj will appear to be randomly and uniformly distributed mod q.


Step 4: Construct (u(x), v(x))εLh such that








u


(
x
)


=





j
=
1

n




δ
j

(
u
)





c
j



(
x
)







and






υ


(
x
)




=




j
=
1

n




δ
j

(
υ
)





c
j



(
x
)






,




with δj(u), δj(v) small, δj(u)≡0 (mod p), and δj(v)j(mod p) for all j.


To construct the desired (u(x), v(x)), search for an appropriate r(x) which is short, and set






u(x)=pr(x)a(x) and v(x)=r(x)v(x).


Such an r(x) must satisfy









r


(
x
)




b


(
x
)



=




j
=
1

n




δ
j

(
υ
)





c
j



(
x
)





,




with the δj(v) small and δj(v)jj(mod p),


and also satisfy









pr


(
x
)




a


(
x
)



=




j
=
1

n




δ
j

(
u
)





c
j



(
x
)





,




with the δj(u) small and δj(u)≡(mod p).


As r(x), a(x) are short, r(x)a(x) is also short, and we may write









r


(
x
)




a


(
x
)



=





i
=
0


n
-
1





d
i



x
i






f



,




with the di small. Then the δj(u) of








pr






(
x
)



a


(
x
)



=




j
=
1

n




δ
j






(
u
)






c
j



(
x
)





,




are given by





1(u), . . . ,δn(u))=p(d0, . . . ,dn-1)U.


As all the di and entries of U are small there is no wraparound mod q and each δj(u)≡0 (mod p). Thus for whatever short r(x) we find, the (δj(u)=0 (mod p) condition will hold.


We turn now to finding r(x) short, such that









r


(
x
)




b


(
x
)



=




j
=
1

n




δ
j






(
v
)






c
j



(
x
)





,




with δj(v) short and δj(v)εj−ηj(mod p).


To accomplish this, write b(x)=Σin-1bixi, set





(b0,0,b0,1, . . . ,b0,n-1)=(b0,b1, . . . ,bn-1),


and define (bi,0, bi,1, . . . , bi,n-1) by






x
i
b(x)=bi,0+bi,1x+ . . . +bi,n-1xn-1.


Let B denote the matrix whose i, j entry is bi,j, and let





β=BU.


Note that the entries βi,j of β are small because the bi,j and the entries of U are small.


For any







r


(
x
)


=





i
=
0


n
-
1





r
i



x
i








i
=
1

n




r
i




x
i







(

mod






(

f


(
x
)


)










we have








r


(
x
)




b


(
x
)



=


(


r
1


,

r
2


,





,

r
n



)


β







(





c
1



(
x
)








c
2



(
x
)













c
n



(
x
)





)

.






To solve for r(x), first define





(r′1,r′2, . . . ,r′n)≡(δ1(v),δ1(v), . . . ,δ1(v)−1(mod p)


and lift each r′j to r′jε(−p/2,p/2].


Then define δj9(v) by





1(v), . . . ,δn(v))≡(r′1, . . . ,r′n)β.


This accomplishes the goal









r


(
x
)




b


(
x
)



=




j
=
1

n




δ
j






(
v
)






c
j



(
x
)





,




with δj(v)=δj(v) (mod p).


After accomplishing Step 4, we have found a short pair of vectors (u(x), v(x))εLh with the appropriate congruential properties.


Having done so, set






s(x)=s0(x)+u(x) and t(x)=t0(x)+v(x).


Then






s


(
x
)


=





j
=
1

n




(


δ
j

+

δ
j

(
u
)



)




c
j



(
x
)







and






t
(
x
)



=




j
=
1

n




(


η
j

+

δ
j

(
v
)



)





c
j



(
x
)


.








By our construction, δjj(u) and ηjj(v). satisfy the required congruences mod p.


If, in addition, for an appropriate choice of custom-character, |δjj(u)|<q/2−custom-character, and |ηjj(v)|<q/2−custom-character, then we accept the signature and release it. If not, we repeat the process.


An argument very similar to that in pq-NTRUSign shows that this guarantees an information free transcript.


Why does short times short=pretty short?:


We will investigate the size of the coefficients of the remainder when a polynomial b(x) is divided by some other polynomial f(x), and in particular, how the coefficients of the remainder depend on the magnitude of the roots of f(x).


Example

Take n=150 and choose f(x) randomly to have the form






f(x)=x150+(random trinary polynomial of degree 90).


There are roughly 2142 such f. We next choose random trinary polynoimals g1(x), g2(x) of degree 150 and compute g1(x)g2(x) mod f(x). A sequence of random trials of f, g1, g2 produced polynomials whose (maximum coefficient, minimum coefficient) were





(343,−484),(255,−235),(607,−486),(411,−441),(552,−560), . . .


The spread of the coefficient range is very dependent on the size of the largest complex root of f(x). These roots will in general be considerably smaller if there is a large gap between the leading coefficient of highest degree (150 in the Example) and the non-zero coefficient of highest degree below the leading coefficient (90 in the Example).


Fix integers m≧n>0. Fix a polynomial







f


(
x
)


=





i
=
1

n







(

x
-

θ
i


)








[
x
]


.






Let






b


(
x
)


=




i
=
0


m
-
1





b
i



x
i







be chosen with each bi satisfying some probability distribution. Different coefficients may have different distributions, but we assume that they are independent and have mean 0. (In practice, our b(x) will be a product of plaintexts, so it will be a product of t polynomials whose coefficients are independent and more-or-less uniform in some interval. This means that the coefficients of b(x) each satisfy some sort of t-fold hypergeometric distribution, but note that the middle coefficients will be much larger than the ones near the top and the bottom. That is why we allow the coefficients of our b to have different distributions.) The independence means that






E(bibj)=E(bi)E(bj)=0 if iδj,


while the numbers E(bi2) depend on the distributions satisfied by the various bi.


We perform division with remainder,






b(x)=f(x)q(x)+r(x) with 0≦deg r<n.


As usual, we view the polynomials as vectors,






b=(b0, . . . ,bm) and r=(b0, . . . ,bn)


We let V denote the vanderMonde matrix of the θi's,







V
=



(

θ
i
j

)



1

i

n







0

j
<
n



=

(



1



θ
1







θ
1

n
-
1






1



θ
2







θ
2

n
-
1






















1



θ
n







θ
n

n
-
1





)



,




and we set







θ

(
j
)


=


(




θ
1
j






θ
2
j











θ
n
j




)

.





Then we set








b


(
θ
)


=


(




b


(

θ
1

)







b


(

θ
2

)












b


(

θ
n

)





)

=




j
=
0


m
-
1





b
j



θ

(
j
)






,




and similarly for r(θ).


We take the relation b(x)=f(x)q(x)+r(x) and substitute x=θ1, . . . , θn. Since f(θi)=0, this gives






ri)−bi) for all 1≦i≦n.


With our earlier notation, this is simply the equality of vectors






r(θ)=b(θ).


Now we observe that since r has degree at most n−1, we can write r(θ) as







r


(
θ
)


=





j
=
0


n
-
1





r
j



θ

(
j
)




=

Vr
.






Hence





r=V
−1
b(θ).


We now compute the expected value of ∥r∥2 as b(x) varies.













E


(



r


2

)


=



E


(





V

-
1




b


(
θ
)





2

)








=



E


(
t





b


(
θ
)


t



V

-
1




V

-
1




b


(
θ
)



)







=



E


(




j
,

k
=
0



m
-
1





b
k
l



θ


(
k
)


t




V

-
1




V

-
1




b
j



θ

(
j
)




)








=






j
,

k
=
0



m
-
1






E


(


b
k



b
j


)


t



θ


(
k
)


t




V

-
1




V

-
1




θ

(
j
)










=






j
=
0


m
-
1






E


(

b
j
2

)


t



θ


(
j
)


t




V

-
1




V

-
1




θ

(
j
)










=






j
=
0


m
-
1





E


(

b
j
2

)









V

-
1




θ

(
j
)





2

.










(
1
)







This last formula explains what's going on. If we assume that f(x) is fixed and that deg b(x) is large compared to n=deg f(x), then we obtain the rough, but useful, estimate







E


(



r


2

)










max

0

j
<
m









(


E


(

b
j
2

)


·


max

1

i

n







θ
i



j



)

.






Which term dominates will depend on the relative size of E(bj2) and max |θi|j for 0≦j<m. In our scenario, we have b(x)=ai(x) . . . at(x) with deg ai≈n, so m≈nt. The coefficients of the ai are uniform and small, so most of the coefficients of b are roughly Ct. Then E(∥r∥2) is roughly Ct max |θi|nt. So in order for decryption to work, we need roughly






g>(C max|θi|n)t.


As expected, we get exponential growth in t. But this shows very clearly how the largest root of f(x) has a major influence on the required size of q. Definintion:


Let f(x)εC[x] be a monic polynomial and let θ1, . . . , θn be the roots of f. We let









(
f
)


=


max

1

i

n







θ
i



.






Experiments clearly reveal the effect of the size of the roots of f(x). We fixed an f(x) of degree 11, chose 100 polynomials g(x) of degree 32 with random coefficients in [−2,2] and computed the largest coefficients of g(x) modulo f(x). We used the polynomials






f
1(x)=x11−x10+x9+x6−x5+x2−x−1






f
2(x)=x11+x10+x5−x4+x3−x2−x−1






f
3(x)=x11−x10+x7+x6+x5−x3−x2−1.


Then


















ƒ

custom-character  (ƒ)

Avg |g mod ƒ|∞
St. Dev. |g mod ƒ|∞





















ƒ1
1.1835
43.420
16.226



ƒ2
1.3511
352.250
191.452



ƒ3
1.4307
1167.720
666.196










We now consider if there is an advantage in taking the non-zero coefficients of f(x) to be in the lower degree terms. So we take f(x) to have the form






f(x)=xn+{tilde over (f)}(x),


where {tilde over (f)}(x) is random trinary of small degree. Simple estimates make it clear that such polynomials tend to have smaller roots than polynomials whose non-zero monomials have higher degree. In order to compare with our experiments, we took polynomials f(x) of degree 11 with non-zero coefficients only on monomials of degree at most 4, more precisely, we took






f(x)=x11+a4x4+a3x3+a2x2+a1x−1


with the ai randomly chosen from {±1}. The polynomial






f
4(x)=x11−x4+x3−x2++x−1


has






custom-character(f4)=1.18225,


so custom-character(f4) is comparable to custom-character(f1) for the f1(x) above.


For f4 and 100 samples, we found





Avg |g mod f4|=28.450 and St.Dev. |g mod f4|=15.658.


These may be compared with the roughly similar values 43.4 and 16.2 for f1. A likely reason for the difference is due to secondary effects due to the other roots. Thus the magnitudes of the roots of f1 are





1.18,1.18,1.15,1.15,1.08,1.08,1.00,1.00,0.890,0.890,0.578,


while the magnitudes of the roots of f4 are





1.18,1.18,1.00,1.00,1.00,1.00,1.00,0.953,0.953,0.888,0.888.


So the second largest root of f1 is significantly larger than the second largest root of f4.


As the above formula makes clear, the size of the inverse of the vanderMonde matrix Vf also has an effect. We list the sup norm and the spectral radius of Vf−1 for our two example polynomials.
















ƒ1
ƒ4




















Spectral Radius of Vƒ−1
7.766
5.522



Sup Norm of Vƒ−1
0.666
0.263











We note that the remainder coefficients for division by f1 and f4 resemble one another much more closely than do the remainder coefficients for division by f2 or f4. This suggests that it is not so much the distribution of non-zero monomials that affects the remainder coefficients as it is the size of the roots of f. However, if one desires to find an f with comparatively small roots, it is definitely advantageous to select f with non-zero monomials only in the lower degree terms.

Claims
  • 1. A method for signing and subsequently verifying a digital message, comprising the following steps implemented using at least one processor-based subsystem: generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x];generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y];producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively;producing a secret isomorphism from the first finite field to the second finite field;producing and publishing a public key that depends on F(y);producing a private key that depends on said secret isomorphism;producing a message digest by applying a hash function to the digital message and the public key;producing a digital signature using the message digest and the private key; andperforming a verification procedure utilizing the digital signature and the public key to determine whether the signature is valid.
  • 2. The method as defined by claim 1, wherein said first and second finite fields are designated, respectively, as x-space and y-space, and further comprising generating a specified lattice in x-space and using said isomorphism to generate a corresponding lattice in y-space, and wherein said step of producing said digital signature includes initially producing a signature in the x-space lattice and then producing the digital signature in the corresponding y-space lattice.
  • 3. The method as defined by claim 1, wherein said secret isomorphism comprises a polynomial.
  • 4. The method as defined by claim 1, further comprising transmitting the digital signature, and wherein said step of performing a verification procedure includes receiving the transmitted digital signature and performing the verification procedure on the received digital signature.
  • 5. The method as defined by claim 2, further comprising transmitting the digital signature, and wherein said step of performing a verification procedure includes receiving the transmitted digital signature and performing the verification procedure on the received digital signature.
  • 6. The method as defined by claim 1, wherein said digital message comprises a challenge communication from a verifier entity, and wherein said digital signature is transmitted to said verifier entity.
  • 7. The method as defined by claim 2, wherein said digital message comprises a challenge communication from a verifier entity, and wherein said digital signature is transmitted to said verifier entity.
  • 8. A method for signing and sending a digital message, comprising the following steps implemented using at least one processor-based subsystem: generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x];generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y];producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively;producing a secret isomorphism from the first finite field to the second finite field;producing and publishing a public key that depends on F(y);producing a private key that depends on said secret isomorphism;producing a message digest by applying a hash function to the digital message and the public key;producing a digital signature using the message digest and the private key; andtransmitting the digital signature.
  • 9. The method as defined by claim 8, wherein said first and second finite fields are designated, respectively, as x-space and y-space, and further comprising generating a specified lattice in x-space and using said isomorphism to generate a corresponding lattice in y-space, and wherein said step of producing said digital signature includes initially producing a signature in the x-space lattice and then producing the digital signature in the corresponding y-space lattice.
  • 10. The method as defined by claim 8, wherein said secret isomorphism comprises a polynomial.
  • 11. The method as defined by claim 8, further comprising receiving the digital signature, and employing a further process-based subsystem for performing a verification procedure on the received digital signature.
  • 12. The method as defined by claim 9, further comprising receiving the digital signature, and employing a further process-based subsystem for performing a verification procedure on the received digital signature.
  • 13. A system for signing and subsequently verifying a digital message, comprising: at least one processor-based subsystem that is programmed with instructions that cause the at least one processor subsystem to implement the following steps:generating an irreducible monic polynomial f(x) of degree n in a ring Fq[x];generating an irreducible monic polynomial F(y) of degree n in a ring Fq[y];producing first and second finite fields as Fq[x]/(f(x)) and Fq[y]/(F(y)), respectively;producing a secret isomorphism from the first finite field to the second finite field;producing and publishing a public key that depends on F(y);producing a private key that depends on said secret isomorphism;producing a message digest by applying a hash function to the digital message and the public key;producing a digital signature using the message digest and the private key; andperforming a verification procedure utilizing the digital signature and the public key to determine whether the signature is valid.
  • 14. The system as defined by claim 13, wherein said encoded digital signature is transmitted to a further processor-based subsystem, and said further processor-based system is programmed with instructions that cause the further processor-based system to implement said step of performing said verification procedure utilizing the encoded digital signature and the public key to determine whether the signature is valid.
  • 15. The system as defined by claim 13, wherein said first and second finite fields are designated, respectively, as x-space and y-space, and further comprising said at least one processor-based subsystem being programmed with instructions to implement: generating a specified lattice in x-space and using said isomorphism to generate a corresponding lattice in y-space, and wherein said step of producing said digital signature includes initially producing a signature in the x-space lattice and then producing the digital signature in the corresponding y-space lattice.
  • 16. The system as defined by claim 14, wherein said first and second finite fields are designated, respectively, as x-space and y-space, and further comprising said at least one processor-based subsystem being programmed with instructions to implement: generating a specified lattice in x-space and using said isomorphism to generate a corresponding lattice in y-space, and wherein said step of producing said digital signature includes initially producing a signature in the x-space lattice and then producing the digital signature in the corresponding y-space lattice.
RELATED APPLICATION

This application claims priority from U.S. Provisional Patent Application No. 62/389,390 filed Feb. 25, 2016, and said Provisional Patent Application is incorporated herein by reference.

Provisional Applications (1)
Number Date Country
62389390 Feb 2016 US