# Patent

## Title
Cryptographic method and apparatus for non-linearly merging a data block and a key

## Abstract

A method and apparatus for cryptographically converting a digital input data block into a digital output data block. The apparatus has an input for supplying the input data block and a further input for supplying a code conversion digital key K1. Cryptographic processing merges a selected part M1 of the digital input data block with the key K1 to produce a data block B1 that is non-linearly dependent on M1 and K1. The merging is performed in one sequentially inseparable step. The digital output block is derived from a selected part of the data block B1.

## Background
This application is a continuation of application Ser. No. 08/859,591, filed May 20, 1997.


FIELD OF THE INVENTION
The invention relates to a method for converting a digital input block into a digital output block; said conversion comprising the step of merging a selected part M1 of the digital input block with a first key K1 and producing a data block B1 which non-linearly depends on the selected part M1 and the first key K1, and where a selected part of the digital output block is derived from the data block B1.
The invention further relates to an apparatus for cryptographically converting a digital input block into a digital output block; the apparatus comprising first input means for supplying the digital input block; second input means for supplying a first key K1; cryptographic processing means for converting the digital input block into the digital output block; such conversion comprising merging a selected part M1 of the digital input block with the first key K1 and producing a data block B1 which non-linearly depends on the selected part M1 and the first key K1, and where a selected part of the digital output block is derived from the data block B1; and output means for outputting the digital output block.
BACKGROUND OF THE INVENTION
The Data Encryption Standard (DES) of the National Bureau of Standard [FIPS publication 46, 1977 Jan. 15] describes a widely used algorithm for converting a digital input block into a digital output block. Such an algorithm is generally referred to as a block cipher. The DES algorithm is used for encrypting (enciphering) and decrypting (deciphering) binary coded information. Encrypting converts intelligible data, referred to as plaintext, into an unintelligible form, referred to as ciphertext. Decrypting the ciphertext converts the data back to its original form. In the so-called electronic code book mode, DES is used to encrypt blocks of 64 bits of plaintext into corresponding blocks of 64 bits of ciphertext. In this mode, the encryption uses keys which are derived from a 64 bit key, of which 56 bits may be freely selected. FIG. 1 shows the overall structure of DES during encrypting. In the encrypting computation, the input (64 bit plaintext) is first permuted using a 64 bit fixed permutation IP. The result is split into 32 left bits L0 and 32 right bits R0. The right bits are transformed using a cipher function f(R0,K1), where K1 is a sub-key. The result f(R0,K1) is added (bit-wise modulo 2) to the left bits, followed by interchanging the two resulting 32 bit blocks L0□f(R0,K1) and R0. This procedure is continued iteratively for a total of 16 rounds. At the end of the last round the inverse permutation of the initial permutation IP is applied.
In the calculation of f(Ri,Ki+1) the 32 right bits Ri are first expanded to 48 bits in the box E, as illustrated in FIG. 2. According to a given table this expansion is performed by taking some input bits twice as an output bit and others only once. Then, the expanded 48 bits are added (bit-wise modulo 2) to the 48 key bits Ki. The resulting 48 bits are split into 8 groups of 6 bits each. Each of these groups is processed by an S box (Si), which reduces the 6 bits to 4 bits in a non-linear operation. The eight Si boxes are given in the form of a table. The total output is 32 bits, which is permuted in the box P. P is also given in the form of a table.
FIG. 3 illustrates the key schedule calculation. The key consists of 64 bits, of which only 56 are used in the algorithm. Those 56 bits should be chosen randomly. Eight complementing error detecting bits are used to make the parity of each byte of the key odd. The selection of the 56 bits is performed in box PC1, together with a permutation. The result is split into two 28 bit words C0 and D0. To obtain the 48 key bits for each round, first the words C0 and D0 are left shifted once or twice. A selection and a permutation PC2 are then applied to the result. The output of PC2 is the 48 bit sub-key K1 which is used in f(R0,K1). The process of shifting, selecting and permutating is repeated to generate a sub-key for each round. A table specifies how many shifts must be performed to obtain the next 48 bits of the sub-key for the following round.
The same algorithm and key can be used for decrypting a ciphertext. The initial permutation for the decrypting cancels the inverse permutation of the encrypting. Each round consists of a, so-called, Feistel cipher. It is well-known that for Feistel-ciphers the inverse operation consists of using the same rounds as used for encrypting but applying the sub-keys in inverse order. As such, the first decrypting round must be supplied with the same sub-key as used for the sixteenth encrypting round, the second decrypting round must be supplied with the same sub-key as used for the fifteenth encrypting round, etc. It is also well-known how the DES algorithm can be used in other encryption modes, such as the cipher feedback mode. In this mode, the DES algorithm is used to generate a stream of statistically random binary bits, which are combined with the plaintext, using, for instance, an exclusive-or logic operation.
The DES algorithm, in essence, comprises an initial permutation, followed by sixteen key-dependent computations on part of the data and terminated with an inverse permutation. Each key dependent computation comprises adding (modulo 2) key-dependent bits to the data part, followed by a non-linear operation on sub-blocks of the data part, and terminated by a permutation (linear operation) of the data part.
In general, DES is considered to be a good encryption/decryption tool. It is, however, an open question whether or not DES has remained secure over the past years, particularly in view of the recent very powerful differential cryptanalytic attacks.
SUMMARY OF THE INVENTION
It is an object of the invention to provide a cryptographic method and apparatus of the kind set forth which is more robust against cryptanalytic attacks.
To achieve this object, the cryptographic method according to the invention is characterised in that the step of merging the data and the key is performed by executing a non-linear function g for non-linearly merging said selected part M1 of the data and said first key K1 in one, sequentially inseparable step. In the DES system, as shown in FIG. 2, in a first processing step the R data is bit-wise added to the key, followed by a second processing step of non-linearly processing the result (S-boxes) . According to the invention, an algorithm is used which non-linearly merges data with a key in one step (i.e. one, sequentially inseparable step). As such, adding the key bits to the data is an integrated part of the nonlinear operation, making the system more immune against modern attacks, such as differential cryptanalysis.
In an embodiment of the method according to the invention, in each round both parts of the digital input block are processed, giving a better encryption result than for conventional Feistel ciphers, such as DES, where during each round only half of the digital input block is being processed. To ensure that the same system can be used for both encryption and decryption, one part of the data is processed using an operation g, whereas the other half is processed using the inverse operation g−1. Using this scheme, decrypting is performed by using the same system but supplying the keys in reverse order to the rounds (during decryption the first non-linear step is supplied with the key which, during encryption, was supplied to the last non-linear step, etc). Compared to a conventional implementation of a Feistel cipher with twice as many rounds, the system according to the invention is faster.
The measure of splitting a relatively large data block and key, of for instance 64 bits, into smaller sub-blocks and sub-keys simplifies realtime nonlinear processing.
In an embodiment of the method according to the invention, a constant is used to enhance the quality of the encryption. Advantageously, the constant is predetermined per system, forming, for instance, a customer-specific constant. Alternatively, the constant is generated using a pseudorandom generator.
The invention provides a way for nonlinearly merging the data sub-block and the sub-key in one step. Additionally, different inputs all result in different outputs. This increases the immunity of the system against cryptanalytic attacks, compared to DES where the non-linear operation reduces the 6bit input sub-block to a 4-bit output sub-block, implying that the same output is produced for four different inputs.
In an embodiment of the method according to the invention a constant is used to enhance the quality of the encryption. Advantageously, the constant is predetermined per system, forming, for instance, a customerspecific constant. Alternatively, the constant is generated using a pseudo-random generator.
In an embodiment of the method according to the invention individual sub-blocks corresponding to different parts of the digital input block are swapped to improve the quality of the encryption.
Preferably, the sub-block mi comprises eight data bits. This further improves the quality of the non-linear operation compared to DES, where the non-linear operation converts six to four bits.
Another embodiment has the advantage of reducing the multiplication in GF(28) to operations in GF(24), making it possible to achieve a simpler or more cost-effective implementation.
The multiplication in GF(28) may be reduced to operations in GF(24).
An embodiment of the method according to the invention is characterised in that β is a root of an irreducible polynomial h(x)=x4+x3+x2+x+1 over GF(2). This is a preferred choice for β, allowing the use of the so-called shifted polynomial base.
An embodiment of the method according to the invention is characterised in that calculating the inverse b in an element of GF(28) comprises performing a series of calculations in GF(24). By reducing the inverse operation in CF(28) to operations in CF(24) a simpler or more-cost effective implementation can be achieved.
An embodiment of the method according to the invention is characterised in that calculating the inverse of said element b comprises calculating (a02+a0a1+a12β)−1((a0+a1)+a1D). This is an effective way of reducing the inverse operation in GF(28) to operations in GF(24).
An embodiment of the method according to the invention is characterised in that said first key K1 comprises 64 data bits and wherein each of said sub-keys k1 comprises eight data bits. By using a large key the quality of the encryption is increased.
To achieve the object of the invention, the apparatus according to the invention is characterised in that said cryptographic processing means is arranged to perform said merging by executing a non-linear function g for non-linearly merging said selected part M1 and said first key K1 in one, sequentially inseparable step.
These and other aspects of the invention will be apparent from and elucidated with reference to the embodiments shown in the drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows the processing steps for the DES system,
FIG. 2 illustrates details of merging the data with the key and the non-linear operation in DES,
FIG. 3 illustrates details of the key calculation in DES,
FIG. 4 shows a block diagram of the cryptographic apparatus of the invention,
FIG. 5 illustrates separate processing of two parts of the digital input block,
FIG. 6 illustrates processing of a part of the digital input block in the form of sub-blocks,
FIG. 7 illustrates processing of two parts in the form of subblocks, and
FIG. 8 shows an overall encryption system.



DETAILED DESCRIPTION OF THE INVENTION
FIG. 4 shows a block diagram of the cryptographic apparatus 400 according to the invention. For the purpose of explaining the invention, the system is described in the electronic code book mode. Persons skilled in the art will be able to use the system in other modes as well. The apparatus 400 comprises first input means 420 for providing a digital input block M. The digital input block N may be any suitable size. Preferably, M is sufficiently large, for instance 128 bits, to obtain a reasonably secure encryption result. The apparatus 400 further comprises cryptographic processing means 420 for converting the digital input block into a digital output block. Advantageously, the digital output block has substangtially equal length as the digital input block. The apparatus 400 comprises output means 430 for outputting the digital output block. Basically, the cryptographic processing means 420 converts the digital input block M into the digital output block by merging a selected part M1 of the digital input block N with a first key K1, producing a data block B1 which non-linearly depends on M1 and K1. The merging is performed in one, sequentially inseparable step. The digital output block is derived from B1 and the remaining part of M, which is not part of M1. To obtain the first key K1, the cryptographic apparatus 400 comprises second input block 440. As will be described in more details below, a second part M2 of the digital input block may be nonlinearly merged with a second key K2, preferably, using an operation inverse to the operation for merging M1 and K1, producing a data block B2. In this case, the digital output block also depends on B2. To obtain the second key K2, the cryptographic apparatus 400 comprises third input block 450.
It will be appreciated that the cryptographic apparatus 400 may be implemented using a conventional computer, such as a PC, or using a dedicated encryption/decryption device. The digital input block may be obtained in various ways, such as via a communication network, from a data storage medium, such as a harddisk or floppy disk, or directly being entered by a user. Similarly, the digital output block may be output in various ways, such as via a communication network, stored on a data storage medium, or displayed to a user. Preferably, secure means are used to this end. The cryptographic processing means 420 may be a conventional processor, such as for instance used in personal computers, but may also be a dedicated cryptographic processor. The cryptographic apparatus 400 may, in part or in whole, be implemented on a smart-card.
Details of the cryptographic conversion process will now be described for encrypting blocks of 128 bits of plaintext into corresponding blocks of 128 bits of ciphertext. Persons skilled in the art will be able to use the system for other block sizes as well. Data sizes shown in the Figures are given for reasons of clarity and should be treated as examples only. The description focuses on the non-linear processing of the data and the merging of the key with the data as performed in one round. As such the invention can be applied in a system as shown in FIG. 1, comprising multiple rounds and also including a linear operation on the data block in each round.
As shown in FIG. 5, the message block M of 128 bits is divided into a first part M1 and a second part M2 (a left and a right block). Preferably, both parts are of equal size, 64 bits. It will be appreciated that M1 and M2 may also be derived from M using a more complicated selection process. M1 is processed using a non-linear function g. In principle, it is not required to process M2 during the same round. Advantageously, M2 is processed in the same round using the inverse function g−1. Each of the functions g and g−1 non-linearly merges, M1 or, respectively, M2 with a key K1 or, respectively K2. Preferably, the data parts and the keys have the same size. Since it is difficult to implement a good non-linear operation on a large data block and non-linearly processing a large data block is time consuming, the data parts M1 and M2 are split into sub-blocks. FIG. 6 illustrates this for M1. FIG. 7 illustrates the splitting of M1 and M2. Using 64-bit data parts M1 and M2, advantageously, the parts are each split into eight 8-bit elements, where M1=(m0, m1, . . . , m7) and M2=(m8, m9, . . . , m15). The two keys K1 and K2 may be derived from a larger key, for instance, by splitting a 128 bit key into two 64-bit keys K1 and K2. The two keys K1 and K2 may be split further. Using 64-bit keys, advantageously, each key is split into 8-bit sub-keys, giving a total of sixteen 8-bit sub-keys kj, j=0 . . . 15. Each of the sub-keys kj is associated with the corresponding sub-block mj. Each sub-block is processed separately. Preferably, the sub-blocks are processed in parallel. If preferred, the sub-blocks relating to one round may also be serially processed. The first group of sub-blocks, forming M1, are each processed by a cipher function f. The second group of sub-blocks are each processed by the inverse function f−1.
For the cryptographic operations, an n-bit sub-block or sub-key is considered to represent an element of GF(2n) (Galois Field). All operations are, therefore, in GF(2n).
In its basic form, the cipher function f has two inputs mj and kj and one output tj as also illustrated in FIGS. 6 and 7, where tj=f(mj, kj), for j=0 to 7. In the basic form, the cipher function f involves one operation h(bj, kj) with an output of substantially equal size as bj. The function h has a data sub-block bj and a sub-key kj as input, where bj=mj for the basic form of the cipher function f. The function f (in this embodiment the same as the function h) is defined as follows for j=0 . . . 7:




  
    
      
        
          
            h
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          
            k
            j
          
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                k
                j
              
              )
            
            
              -
              2
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              2
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        
          k
          j
        
      
    
  




Similarly, in its basic form the inverse cipher function f−1 has two inputs mj and kj and one output tj as also illustrated in FIGS. 6 and 7, where tj=f−1(mj, kj), for j=8 to 15. The inverse cipher function f−1 involves also one operation, h−1(bj.kj) with an output of substantially equal size as bj. The function h−1 is the inverse of h. As before, bj=mj in the basic form of the cipher function f−1. The function f−1 (in this embodiment the same as the functions h−1) is defined as follows for j=8 . . . 15:




  
    
      
        
          
            
              h
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              
                b
                j
              
              ·
              
                k
                j
                2
              
            
          
          ≠
          1
        
      
    
  
  
    
      
        
        ⁢
        
          
            k
            j
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              
                1
                2
              
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            
              b
              j
            
            ·
            
              k
              j
              2
            
          
        
        =
        1
      
    
  




In a further embodiment, the outputs tj of the cipher functions f(tj=f(mj, kj), for j=0 to 7) and the outputs of the inverse cipher function f−1(tj=f−1(mj, kj), for j=8 to 15) are swapped in the following manner: tj<−>t15−j for j=0 to 7 and the swapped subblocks are concatenated to form a single result. This is illustrated in FIG. 7.
In a further embodiment, a constant is added (bit-wise module 2) to each data sub-block mj before executing the function h. Preferably, eight independent constants pj (j=0 . . . 7) are used, each being added to the corresponding data sub-block mj. The same function h is used as before, now operating on bj=mj⊕pj. The cipher function f is now defined as follows:



  
    
      b
      j
    
    =
    
      
        m
        j
      
      ⊕
      
        p
        j
      
    
  




  
    
      
        
          
            
              h
              ⁡
              
                (
                
                  
                    b
                    j
                  
                  ,
                  
                    k
                    j
                  
                
                )
              
            
            =
            
            ⁢
            
              
                (
                
                  
                    b
                    j
                  
                  ·
                  
                    k
                    j
                  
                
                )
              
              
                -
                1
              
            
          
          ,
        
      
      
        
          
            
              if
              ⁢
              
                
              
              ⁢
              
                b
                j
              
            
            ≠
            0
          
          ,
          
            
              k
              j
            
            ≠
            0
          
          ,
          
            
              and
              ⁢
              
                
              
              ⁢
              
                b
                j
              
            
            ≠
            
              k
              j
            
          
        
      
    
    
      
        
          
          ⁢
          
            
              
                (
                
                  k
                  j
                
                )
              
              
                -
                2
              
            
            ,
          
        
      
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          =
          0
        
      
    
    
      
        
          
          ⁢
          
            
              
                (
                
                  b
                  j
                
                )
              
              
                -
                2
              
            
            ,
          
        
      
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              k
              j
            
          
          =
          0
        
      
    
    
      
        
          
          ⁢
          
            0
            ,
          
        
      
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          =
          
            k
            j
          
        
      
    
  



Similarly, for the inverse cipher function f−1 also a constant is added (bit-wise module 2) to each data sub-block mj. To allow the inverse function f−1 to be used to decrypt text encrypted using the cipher function f, the constant is added after the function h. Preferably, the same eight independent constants pj (j=0 . . . 7) are used as used for the cipher function f. Now, the constants pj are being added to the 15−j-th stream (j=0 . . . 7). As a consequence, the inverse cipher function f−1 involves the following two operations (j=8 . . . 15):





  
    
      
        
          
            
              h
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              
                b
                j
              
              ·
              
                k
                j
                2
              
            
          
          ≠
          1
        
      
    
  
  
    
      
        
        ⁢
        
          
            k
            j
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              
                1
                2
              
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            
              b
              j
            
            ·
            
              k
              j
              2
            
          
        
        =
        1
      
    
  




tj=h−1(bj.kj)⊕p15−j  2.

Finally, tj and t15−j are swapped (j=0 . . . 7).
In a further embodiment, a further constant is added (bit-wise module 2) to each data sub-block mj after executing the function h. Preferably, eight independent constants dj (j=0 . . . 7) are used, each being added to the corresponding data sub-block mj. The same function h is used as before. The cipher function f is now defined as follows:

bj=mj⊕pj  1.





  
    
      
        
          
            h
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          
            k
            j
          
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                k
                j
              
              )
            
            
              -
              2
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              2
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        
          k
          j
        
      
    
  




tj=h(bj, kj)⊕dj  3.

Similarly, for the inverse cipher function f−1 also a constant is added (bit-wise module 2) to each data sub-block mj. To allow the inverse function f−1 to be used to decrypt text encrypted using the cipher function f, the constant is added before executing the function h. Preferably, the same eight independent constants dj(j=0 . . . 7) are used as used for the cipher function f. Now, the constants dj are being added to the 15−j-th stream (j=0 . . . 7). The same function h−1 is used as before, now operating on bj=mj⊕d15−j. As a consequence, the inverse cipher function f−1 involves the following three operations (j=8 . . . 15):

bj=mj⊕d15−j.  1.





  
    
      
        
          
            
              h
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              
                b
                j
              
              ·
              
                k
                j
                2
              
            
          
          ≠
          1
        
      
    
  
  
    
      
        
        ⁢
        
          
            k
            j
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              
                1
                2
              
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            
              b
              j
            
            ·
            
              k
              j
              2
            
          
        
        =
        1
      
    
  




tj=h−1(bj.kj)⊕p15−j  3.

Finally, tj and t15−j are swapped (j=0 . . . 7).
It will be appreciated that it is also possible to use the constants dj without using constants pj.
In a further embodiment, the cipher function f raises the outcome of the function h to a power of two. The same function h is used as before. The cipher function f is now defined as follows:

bj=mj⊕pj  1.





  
    
      
        
          
            h
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          
            k
            j
          
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                k
                j
              
              )
            
            
              -
              2
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              2
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        
          k
          j
        
      
    
  




sj=h(bj, kj)2fi  3.

tj=sj⊕dj  4.

Similarly, the inverse cipher function f−1 also raises a data sub-block to a power of 2. To allow the inverse function f−1 to be used to decrypt text encrypted using the cipher function f, the additional operation is performed before executing the function h. The same function h−1 is used as before, now operating on bj=mj⊕d15−j. As a consequence, the inverse cipher function f−1 involves the following four operations (j=8 . . . 15):

qj=mj⊕d15−j  1.

bj=qj2f(j−7)  2.





  
    
      
        
          
            
              h
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  b
                  j
                
                ,
                
                  k
                  j
                
              
              )
            
          
          =
          
          ⁢
          
            
              (
              
                
                  b
                  j
                
                ·
                
                  k
                  j
                
              
              )
            
            
              -
              1
            
          
        
        ,
      
    
    
      
        
          
            if
            ⁢
            
              
            
            ⁢
            
              b
              j
            
          
          ≠
          0
        
        ,
        
          
            k
            j
          
          ≠
          0
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              
                b
                j
              
              ·
              
                k
                j
                2
              
            
          
          ≠
          1
        
      
    
  
  
    
      
        
        ⁢
        
          
            k
            j
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            b
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                b
                j
              
              )
            
            
              -
              
                1
                2
              
            
          
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            k
            j
          
        
        =
        0
      
    
  
  
    
      
        
        ⁢
        
          0
          ,
        
      
    
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            
              b
              j
            
            ·
            
              k
              j
              2
            
          
        
        =
        1
      
    
  




tj=h−1(bj.kj)⊕p15−j  3.

Finally, tj and t15−j are swapped (j=0 . . . 7).
It will be appreciated that it is also possible to use the operation of raising to a power of 2 without using one or both of the constants dj and pj.
For decrypting the same algorithm is used as for encrypting, but the sub-keys are swapped: instead of kj, k15−j is used, j=0 . . . 15.
The Multiplication in GF(28)
In principle, for the invention any multiplication in GF(28) may be used. An example of a VLSI implementation of multiplications in GF(2m) is given in [P. A. Scott, “A fast VLSI multiplier for GF(2m)”, IEEE Journal on selected areas in communications, Vol. SAC-4, No. 1, January 1986, pages 62–66]. Advantageously, the following mechanism is used to reduce the multiplication in GF(28) to a series of multiplications and additions in GF(24) . As is known in the art, in finite fields with a characteristic of 2 (e.g. GF(2n) ) and the Galois field represented in binary arithmetic, the subtraction operation (i.e. the inverse of addition) is the same as the addition operation. For convenience, the “+” symbol is used herein for this addition/subtraction operation, although a “−” symbol may be equivalently substituted for ease of understanding, as required.
Let in GF(24), β be the non-trivial root of β5=1 (non-trivial means β≠1 or, equally, β is the root of the irreducible polynomial h(x)=x4+x3+x2+x+1 over GF(2), since: x5+1=(x+1)(x4+x3+x2+x+1) ). The normal base β, β2, β4, β8 is taken as the base in GF(24). Since according to the polynomial β8=β3, this is the same as the so-called shifted polynomial base: β, β2, β3, β4.
Let D be an element of GF(28), defined as a root of the irreducible polynomial k(x)=x2+x+β over GF(24). Every element of GF(28) can be represented as a0+a1.D, with a0 and a1 being elements of GF(24). In binary terminology, the number b of GF(28) can be represented using eight bits, arranged as a vector (a0, a1), with a0, a1 having four bits, representing numbers of GF(24). As such, the base in GF(28) is: β, β2, β3, β4, Dβ, Dβ2, Dβ3, Dβ4. Two elements b and c of GF(28), represented as b=a0+a1.D and c=a2+a3.D, with aiεGF(24), can be multiplied as follows:

b.c=(a0+a1.D).(a2+a3.D)=a0a2+(a1a2+a0a3).d+a1a3.D2.

Using the fact that D is a root of k(x), which implies that: D2=D+β, this gives the multiplication result:

b.c=(a0a2+a1a3β)+(a1a2+a0a3+a1a3).D.

This has reduced the multiplication of two elements of GF(28) to a series of multiplications and additions in GF(24).

The Inverse in GF(28)

In principle any known method may be used to calculate the inverse of an element in GF(28). Advantageously, if the previous method has been used to reduce the multiplication in GF(28) to a multiplication in GF(24), then the following method is used to reduce the inverse operation in GF(28) to an inverse operation in GF(24).
The inverse b−1 of an element b in GF(28), where b is represented as b=a0+a1.D, with a1εGF(24), is given by:



  
    
      
        b
        
          -
          1
        
      
      =
      
        
          
            (
            
              
                a
                0
                2
              
              +
              
                
                  a
                  0
                
                ⁢
                
                  a
                  1
                
              
              +
              
                
                  a
                  1
                  2
                
                ⁢
                β
              
            
            )
          
          
            -
            1
          
        
        ·
        
          (
          
            
              a
              0
            
            +
            
              a
              1
            
            +
            
              
                a
                1
              
              ⁢
              D
            
          
          )
        
      
    
    ,
    
      
    
    ⁢
    
      since:
    
  




  
    
      
        
          
            b
            
              -
              1
            
          
          ⁣
          
            b
            =
            
              
                
                  (
                  
                    
                      a
                      0
                      2
                    
                    +
                    
                      
                        a
                        0
                      
                      ⁢
                      
                        a
                        1
                      
                    
                    +
                    
                      
                        a
                        1
                        2
                      
                      ⁢
                      β
                    
                  
                  )
                
                
                  -
                  1
                
              
              ·
              
                (
                
                  
                    a
                    0
                  
                  +
                  
                    a
                    1
                  
                  +
                  
                    
                      a
                      1
                    
                    ⁢
                    D
                  
                
                )
              
              ·
              
                (
                
                  
                    a
                    0
                  
                  +
                  
                    
                      a
                      1
                    
                    ⁢
                    D
                  
                
                )
              
            
          
        
      
    
    
      
        
          
            =
            
              
                
                  (
                  
                    
                      a
                      0
                      2
                    
                    +
                    
                      
                        a
                        0
                      
                      ⁢
                      
                        a
                        1
                      
                    
                    +
                    
                      
                        a
                        1
                        2
                      
                      ⁢
                      β
                    
                  
                  )
                
                
                  -
                  1
                
              
              ·
              
                (
                
                  
                    a
                    0
                    2
                  
                  +
                  
                    
                      a
                      0
                    
                    ⁢
                    
                      a
                      1
                    
                  
                  +
                  
                    
                      a
                      1
                      2
                    
                    ⁢
                    D
                  
                  +
                  
                    
                      a
                      1
                      2
                    
                    ⁢
                    
                      D2
                      2
                    
                  
                
                )
              
            
          
          ,
        
      
    
  



and since D2+D=β, this gives: b−1.b=1.

In this way the inverse operation in GF(28) is reduced to an inverse operation in GF(24) and a series of multiplications and additions in GF(24).

Multiplication in GF(24)

In principle, any multiplication in GF(24) may be used. Advantageously, as described before, the shifted polynomial base β, β2, β3, β4 is taken as the base in GF(24), where β is the root of the irreducible polynomial h(x)=x4+x3+x2+x+1 over GF(2), and β5=1 in GF(24). Since β is a root of h, this implies:

    β4+β3+β2+β=1. Assuming that the base elements are named e1, e2, e3 and e4, with ei=βi, the base elements are multiplied in the following way, using the definition of β:

e1.e1=β.β=β2=e2

e1.e2=β.β2=β3=e3

e1.e3=β.β3=β4=e4

e1.e4=β.β4=β5=1=e1+e2+e3+e4

e2.e2=β2.β2=β4=e4

e2.e3=β2.β3=β5=1=e1+e2+e3+e4

e2.e4=β2.β4=β6=β=e1

e3.e3=β3.β3=β6=β=e1

e3.e4=β3.β4=β7=β2=e2

e4.e4=β4.β4=β8=β3=e3

This in principle defines the multiplication in GF(24). In binary terms the multiplication can be seen as follows. With respect to the base, each element b in GF(24) can be represented as b=b0e1+b1e2+b2e3+b3e4, with biεGF(2). As such, the element b can be represented by a 4-dimensional vector with binary components (b1, b1, b2, b3). On a micro-processor this can be represented using a nibble. In binary terms, the multiplication of two elements b and c in GF(24) can be seen as follows, assuming the two elements are represented by b=(b0, b1, b2, b3) and c=(c0, c1, c2, c3). Multiplying the two elements in the normal way gives:







  
    
      
        
          b
          ·
          c
        
        =
        
        ⁢
        
          
            
              (
              
                
                  b
                  0
                
                ⁢
                
                  c
                  0
                
              
              )
            
            ⁢
            
              β
              2
            
          
          +
          
            
              (
              
                
                  
                    b
                    
                      0
                      ⁢
                      
                        
                      
                    
                  
                  ⁢
                  
                    c
                    1
                  
                
                +
                
                  
                    b
                    1
                  
                  ⁢
                  
                    c
                    0
                  
                
              
              )
            
            ⁢
            
              β
              3
            
          
          +
          
            
              (
              
                
                  
                    b
                    0
                  
                  ⁢
                  
                    c
                    2
                  
                
                +
                
                  
                    b
                    1
                  
                  ⁢
                  
                    c
                    1
                  
                
                +
                
                  
                    b
                    2
                  
                  ⁢
                  
                    c
                    0
                  
                
              
              )
            
            ⁢
            
              β
              4
            
          
          +
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                
                  
                    b
                    0
                  
                  ⁢
                  
                    c
                    3
                  
                
                +
                
                  
                    b
                    1
                  
                  ⁢
                  
                    c
                    2
                  
                
                +
                
                  
                    b
                    2
                  
                  ⁢
                  
                    c
                    1
                  
                
                +
                
                  
                    b
                    3
                  
                  ⁢
                  
                    c
                    0
                  
                
              
              )
            
            ⁢
            
              β
              5
            
          
          +
          
            
              (
              
                
                  
                    b
                    1
                  
                  ⁢
                  
                    c
                    3
                  
                
                +
                
                  
                    b
                    2
                  
                  ⁢
                  
                    c
                    2
                  
                
                +
                
                  
                    b
                    3
                  
                  ⁢
                  
                    c
                    1
                  
                
              
              )
            
            ⁢
            
              β
              6
            
          
          +
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                
                  
                    b
                    2
                  
                  ⁢
                  
                    c
                    3
                  
                
                +
                
                  
                    b
                    3
                  
                  ⁢
                  
                    c
                    2
                  
                
              
              )
            
            ⁢
            
              β
              7
            
          
          +
          
            
              (
              
                
                  b
                  3
                
                ⁢
                
                  c
                  3
                
              
              )
            
            ⁢
            
              β
              8
            
          
        
      
    
  




Using the definition of β to replace β5 by β4+β3+β2+β, β6 by β, β7 by β2, and β8 by β3, gives the following four components:

b.c=(b1c3+b2c2+b3c1+b0c3+b1c2+b2c1+b3c0)β+(b0c0+b2c3+b3c2+b0c3+

b1c2+b2c1+b3c0)β2+(b0c1+b1c0+b3c3+b0c3+b1c2+b2c1+b3c0)β3+(b0c2+b1c1+

b2c0+b0c3+b1c2+b2c1+b3c0)β4

The result of the multiplication, in binary terms, is, therefore, given by:





  
    b
    .
    c
  
  =
  
    (
    
      
        
          
            b
            1
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            
              0
              ⁢
              
                
              
            
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            1
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            0
          
        
      
      ,
      
        
          
            b
            0
          
          ⁢
          
            c
            0
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            
              0
              ⁢
              
                
              
            
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            1
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            0
          
        
      
      ,
      
        
          
            b
            0
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            1
          
          ⁢
          
            c
            0
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            0
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            1
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            0
          
        
      
      ,
      
        
          
            b
            0
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            1
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            0
          
        
        +
        
          
            b
            0
          
          ⁢
          
            c
            3
          
        
        +
        
          
            b
            1
          
          ⁢
          
            c
            2
          
        
        +
        
          
            b
            2
          
          ⁢
          
            c
            1
          
        
        +
        
          
            b
            3
          
          ⁢
          
            c
            0
          
        
      
    
    )
  




Inverse Operation in GF(24)

Using the normal base β, β2, β4, β8, each element x of GF(24) can be written as b=a.β+b.β2+c.β4+d.β8, with a, b, c, d εGF(2). As such, each element can be represented by a 4-dimensional vector (a, b, c, d).
In order to obtain the inverse of b (b−1):


    
    
        calculate the following intermediate results: ab, a{overscore (b)}, āb, bc, {overscore (b)}c, b{overscore (c)}, cd, {overscore (c)}d, c{overscore (d)}, da, {overscore (d)}a, dā, where ab is the binary AND of a and b (a AND b) and ā is the binary complement of a (NOT a).
        calculate the first bit of b−1 by using cd, {overscore (c)}d, c{overscore (d)}, āb, b{overscore (c)}, and {overscore (d)}a as follows:
        (cd) OR (ā AND c{overscore (d)}) OR ({overscore (c)}d AND āb) OR (b{overscore (c)} AND {overscore (d)}a)
        calculate the second bit of b−1 by using da, {overscore (d)}a, dā, {overscore (b)}c, c{overscore (d)}, āb as follows:
        (da) OR ({overscore (b)} AND dā) OR ({overscore (d)}a AND {overscore (b)}c) OR (c{overscore (d)} AND āb)
        calculate the third bit of b−1 by using ab, āb, a{overscore (b)}, {overscore (c)}d, dā, {overscore (b)}c as follows:
        (ab) OR ({overscore (c)} AND a{overscore (b)}) OR (āb AND {overscore (c)}d) OR (dā AND {overscore (b)}c)
        calculate the fourth bit of b−1 by using bc, {overscore (b)}c, b{overscore (c)}, {overscore (d)}a, a{overscore (b)}, {overscore (c)}d as follows:
        (bc) OR ({overscore (d)} AND b{overscore (c)}) OR ({overscore (b)}c AND {overscore (d)}a) OR (a{overscore (b)} AND {overscore (c)}d)
    
    


Besides being used in a DES-like system as shown in FIG. 1, a dedicated system can be built around the non-linear algorithm of the invention. Such a system is shown in FIG. 8. In this system, the blocks are processed using the non-linear operation NL of the invention and a linear operation LIN. The first step is the non-linear operation. This is followed by an iteration of the linear operation followed by the non-linear operation. It is expected that a sufficiently safe system is achieved by performing six non-linear operations (i.e. using five rounds), provided that the linear operation mixes the data bits thoroughly. Preferably, 15 rounds are used. Each of the linear operations is the same. Also, each of the non-linear operations is the same, but each non-linear operation uses a different key of 128 bits. Advantageously, keys are derived from one global key of, for instance, 256 bits, using a key schedule calculation. The same key is used for encryption as well as decryption. In most cases the key is provided using a smart-card. For the linear operation, advantageously, instead of a permutation a more complex matrix is used. As described before, in addition to the key, each non-linear operation may, optionally, use a constant C of 128 bits, which is split in the constants pj and dj. The constant may be the same for each operation. Advantageously, each non-linear operations is provided with a separate constant. The constants may be predetermined per system (e.g. a customer-specific constant). Alternatively, the constant is generated using a pseudo-random generator.

## Claims

1. A program segment stored on a computer readable medium for cryptographically converting a digital input data block M into a digital output data block; said program segment comprising:
a program portion for merging a selected part M1 of said digital input data block M with a first digital key B1 to produce a data block B1 which non-linearly depends on said selected part M1 and said first digital key K1; and
a program portion for deriving said digital output block from said data block B1 end the remaining part of the digital input data block M, wherein said merging step is performed by executing a non-linear function g for non-linearly merging said selected part M1 and said first key K1 in a single step.

2. A program segment as claimed in claim 1, comprising:
a program portion for splitting said digital input block into said selected part M1 and a second part M2 before executing said program portion for merging;
a program portion for executing a non-linear function g−1 to merge said second block M2 with a second key K2 in one step, producing a data block B2 as output; said non-linear function g−1 being the inverse of said non-linear function g; and
a program portion for forming combined data from data in said data black B1 and in said data block B2; said digital output block being derived from said combined data.

3. A program segment as claimed in claim 1, wherein said program portion for merging comprises:
a program portion for splitting said selected part M1 in a first plurality n of sub-blacks m0, . . . , mn−1 of equal length;
a program portion for splitting said first key K1 in said first plurality n of sub-keys k0, . . . , kn−1 having equal length, the sub-key K1 corresponding to the sub-block m1, for i=0 to n−1;
a program portion for separately processing each of said sub-blocks mi by executing for each of said sub-blocks mi a same non-linear function h for non-linearly merging a sub-block bi derived from said sub-block mi with said corresponding sub-key ki in one, sequentially inseparable step and producing said first plurality of output sub-blocks h(bi, ki); and
a program portion for combining sub-blocks ti derived from said first plurality of said output sub-blocks h(bi, ki) to form said data block B1.

4. A program segment as claimed in claim 2, wherein said program portion for executing said non-linear function g−1 comprises:
a program portion for splitting said second part M2 in said first plurality n of sub-blocks mn, . . . , m2n−1, having equal length;
a program portion for splitting said key K2 in said first plurality n of sub-keys kn, . . . , k2n −1, having equal length, the sub-key ki corresponding to the sub-block mi, for i=n to 2n−1;
a program portion for executing for each of said sub-blocks mi a same non-linear function h−1 for non-linearly merging a sub-block bi derived from said sub-block mi with said corresponding sub-key ki and producing said first plurality of an output sub-block h−1(bi, ki); said function h−1 being the inverse of said function h; and
a program portion for combining sub-blocks ti derived from said first plurality of output sub-blocks h−1(bi, ki) to form said data block B2.

5. A program segment as claimed in claim 3, wherein said sub-block bi is derived from said sub-block mi by bit-wise adding a constant pi to said sub-block mi, said constant pi having equal length as said sub-block mi.

6. A program segment as claimed in claim 3, characterised in that said function h(bi,ki) is defined by:


h(bi,ki)=(bi,ki)−1, if bi≠0, ki≠0, and bi≠ki



h(bi,ki)=(ki)−2, if bi=0



h(bi,ki)=(bi)−2, if ki=0



h(bi,ki)=0, if bi=ki,


where the multiplication and inverse operations are predetermined Galois Field multiplication and inverse operations.

7. A program segment as claimed in claim 6, wherein deriving said sub-blocks ti from said output sub-blocks h(bi, ki) comprises bit-wise adding a constant di to said output sub-block h(bi,ki), said constant di having equal length as said sub-block mi.

8. A program segment as claimed in claim 7, wherein deriving said sub-blocks ti from said output sub-blocks h(bi, ki) further comprises raising h(bi,ki)⊕di to a power 21, using said predetermined Galois Field multiplication.

9. A program segment as claimed in claim 6, wherein deriving said sub-blocks ti from said output sub-blocks h(bi, ki) comprises raising said output sub-block h(bi,ki) to a power 2i, using said predetermined Galois Field (GF) multiplication.

10. A program segment as claimed in claim 4, wherein said combined data is formed by:
swapping the sub-blocks ti and t2n−1−i, for i =0 to n−1.

11. A program segment as claimed in claim 6, wherein said sub-block mi comprises eight data bits, and wherein said multiplying of two elements b and c of GF(28) comprises executing a series of multiplications and additions in GF(24).

12. A program segment as claimed in claim 11, wherein said multiplying of said two elements b and c comprises:
representing b as a0+a1.D and c as a2+a3.D, where a0, a1, a2 and a3 are elements of GF(24), and where D is an element of GF(28) defined as a root of an irreducible polynomial k(x)=x2+x+β over GF(24) where β is an element of GF(24); and
calculating (a0a2+a1a3β)+(a1a2+a0a3+a1a3).D.

13. A program segment as claimed in claim 12, wherein β is a root of an irreducible polynomial h(x)=x4+x3+x2+x+1 over GF(2).

14. A program segment as claimed in claim 6, wherein said sub-block mi comprises eight data bits, and wherein calculating the inverse of an element b of GF(28) comprises performing a series of calculations in GF(24).

15. A program segment as claimed in claim 14, wherein calculating the inverse of said element b comprises:
representing b as a0+a1.D, where a0 and a1 are elements of GF(24), and where D is an element of GF(28) defined as a root of an irreducible polynomial k(x)=x2+x+β over GF(24), where β is an element of GF(24); and
calculating (a02+a0a1+a12β−1 ((a0+a1)+a1D).

16. A processor for cryptographically converting a digital input block into a digital output block, said processor comprising:
a first input means for obtaining said digital input block;
a second input means for obtaining a first key K1; and
a cryptographic processing means for converting the digital input block into the digital output block by executing a non-linear funtion g for non-linearly merging said selected part M1 and said first key K1 in one step and producing a data block B1 which non-linearly depends on said selected part M1 and said first key K1, where a selected part of said digital output block is derived from said data block B1.

17. A processor as claimed in claim 16, wherein said processor comprises a third input means for obtaining a second key K2, and wherein said processor is arranged to:
split said digital input block into said selected part M1 and a second part M2 before performing said merging;
perform a non-linear function g−1 to merge said second block M2 with said second key K2 in one step, producing a data block B2 as output; said non-linear function g−1 being the inverse of said non-linear function g; and
combine data from data in said data block B1 and in said data block B2; said digital output block being derived from said combined data.

18. A processor as claimed in claim 16, wherein said merging comprises:
splitting said selected part M1 in a first plurality n of sub-blocks m0, . . . , mn−1 of equal length;
splitting said first key K1 in said first plurality a of sub-keys k0, . . . , kn−1, having equal length, the sub-key ki corresponding to the sub-block mi, for i=0 to n−1; and
separately processing each of said sub-blocks mi by executing for each of said sub-blocks mi a same non-linear function h for non-linearly merging a sub-block bi derived from said sub-block mi with said corresponding sub-key ki in one, sequentially inseparable step and producing said first plurality of output sub-blocks h(bi, ki); and
combining sub-blocks ti derived from said first plurality of said output sub-blocks h(bi, ki) to form said data block B1.

19. A processor as claimed in claim 18, wherein said function


h(bi,ki)=(bi,ki)−1, if bi≠0, ki≠0, and bi≠ki



h(bi,ki)=(ki)−2, if bi=0



h(bi,ki)=(bi)−2, if ki=0



h(bi,ki)=0, if bi=ki,


where the multiplication and inverse operations are predetermined Galois Field multiplication and inverse operations.

20. A processor as claimed in claim 19, wherein said sub-block mi comprises eight data bits, and wherein said multiplying of two elements b and c of GF(28) comprises:
representing b as a0+a1.D and c as a2+a3.D, where a0, a1, a2 and a3 are elements of GF(24), and where D is an element of GF(28) defined as a root of an irreducible polynomial k(x)=x2+x+β over GF(24), where β is an element of GF(24); and
calculating (a0a2+a1a3β)+(a1a2+a0a3+a1a3).D;

and wherein calculating the inverse of an element b of GF(28) comprises calculating (a02+a0a1+a12β)−1(a0+a1D).