GROS BUG PRGA !!!!!!!!!!!!!!!!!!!!!!!!! A CORRIGER PARTOUT !!!!!!!!!

[Cipher_code.git] / OneRoundIoT / OneRound / one_round_par2.cpp

diff --git a/OneRoundIoT/OneRound/one_round_par2.cpp b/OneRoundIoT/OneRound/one_round_par2.cpp
index e42a6d5251273cf19c30a20967c9dea7aa64a47f..4a069b4653a8766b76f8fe8804d3814866f33c38 100644 (file)
--- a/OneRoundIoT/OneRound/one_round_par2.cpp
+++ b/OneRoundIoT/OneRound/one_round_par2.cpp
@@ -12,9 +12,9 @@
 #include <omp.h>
 
 /*#include <cryptopp/hex.h>
 #include <omp.h>
 
 /*#include <cryptopp/hex.h>

-#include <cryptopp/sha.h>
-#include <cryptopp/osrng.h>
-#include <cryptopp/secblock.h>
+  #include <cryptopp/sha.h>
+  #include <cryptopp/osrng.h>
+  #include <cryptopp/secblock.h>

 */
 
 
 */
 
 


@@ -136,11 +136,11 @@ void prga(uchar *sc, int ldata, uchar *r) {
 void rotate(uchar *RM1, uchar *RM2, int size, int n)
 {
   int i;
 void rotate(uchar *RM1, uchar *RM2, int size, int n)
 {
   int i;

-   for (i = 0; i< size-n; i++)
-     RM2[i+n] = RM1[i];
+  for (i = 0; i< size-n; i++)
+    RM2[i+n] = RM1[i];

 
 

-   for (i = 0; i< n; i++)
-     RM2[i] = RM1[size-n-1+i];
+  for (i = 0; i< n; i++)
+    RM2[i] = RM1[size-n-1+i];

 }
 
 
 }
 
 


@@ -148,16 +148,6 @@ template<int h2>
 void encrypt_ctr(uchar* seq_in, uchar *seq_out, int len,uchar* RM1,int *Pbox, int *PboxRM, uchar *Sbox1, uchar *Sbox2, int enc, int num) {
 
 
 void encrypt_ctr(uchar* seq_in, uchar *seq_out, int len,uchar* RM1,int *Pbox, int *PboxRM, uchar *Sbox1, uchar *Sbox2, int enc, int num) {
 
 

-  uchar X[h2];
-  uchar fX[h2*num];
-
-
-  
-
-
-  for(int a=0;a<h2;a++) {
-     X[a]=Sbox1[a&0xFF];           //Warning according to the size of h2, we can be outsize of Sbox1[a]
-   }

 
 
   int loc_len=len/num;
 
 
   int loc_len=len/num;


@@ -166,125 +156,112 @@ void encrypt_ctr(uchar* seq_in, uchar *seq_out, int len,uchar* RM1,int *Pbox, in
   for(int p=0;p<num;p++) {
 
     int id=omp_get_thread_num();
   for(int p=0;p<num;p++) {
 
     int id=omp_get_thread_num();

+    uchar fX[h2];

     uchar RM2[h2];
     uchar RM2[h2];

+    

     for(int a=0;a<h2;a++) {
       RM2[a]=RM1[id*h2+a];           //Warning according to the size of h2, we can be outsize of Sbox1[a]
     }
 
     for(int a=0;a<h2;a++) {
       RM2[a]=RM1[id*h2+a];           //Warning according to the size of h2, we can be outsize of Sbox1[a]
     }
 

-  
+
+    uchar X[h2];
+
+
+    for(int a=0;a<h2;a+=4) {
+     X[a]=Sbox1[a&0xFF];           //Warning according to the size of h2, we can be outsize of Sbox1[a]
+     X[a+1]=Sbox1[(a+1)&0xFF];
+     X[a+2]=Sbox1[(a+2)&0xFF];
+     X[a+3]=Sbox1[(a+3)&0xFF];           
+   }
+
+    

     int offset=p*loc_len;
 
     
     for(int it=offset;it<offset+loc_len;it++) {
     int offset=p*loc_len;
 
     
     for(int it=offset;it<offset+loc_len;it++) {

-    int ind1,ind2;
+      int ind1,ind2;

 
 
 
 

-    //cout<<id<<" "<<it<<endl;
+      //cout<<id<<" "<<it<<endl;

     
     

-    if(enc) {
-      ind1=it*h2;
-      ind2=Pbox[it]*h2;
-    }
-    else {
-      ind2=it*h2;
-      ind1=Pbox[it]*h2;
-    }
+      if(enc) {
+       ind1=it*h2;
+       ind2=Pbox[it]*h2;
+      }
+      else {
+       ind2=it*h2;
+       ind1=Pbox[it]*h2;
+      }

   
 
   
 

-    /*for(int a=0;a<h2;a+=4){
-      fX[a]=RM1[X[a]];
-      fX[a+1]=RM1[X[a+1]];
-      fX[a+2]=RM1[X[a+2]];
-      fX[a+3]=RM1[X[a+3]];
-      }*/
-
-    for(int a=0;a<h2;a+=4){
-      fX[id*h2+a]=X[a];
-      fX[id*h2+a+1]=X[a+1];
-      fX[id*h2+a+2]=X[a+2];
-      fX[id*h2+a+3]=X[a+3];
-    }
-
-
+      for(int a=0;a<h2;a+=4) {
+       X[a]=X[Sbox1[a]];
+       X[a+1]=X[Sbox1[a+1]];
+       X[a+2]=X[Sbox1[a+2]];
+       X[a+3]=X[Sbox1[a+3]];
+      }

 
 

+     
+      
+      for(int a=0;a<h2;a+=4) {
+       fX[a]=X[a]^RM2[a];
+       fX[a+1]=X[a+1]^RM2[a+1];
+       fX[a+2]=X[a+2]^RM2[a+2];
+       fX[a+3]=X[a+3]^RM2[a+3];
+      }

     
     

-    /*   if(it<513) {
-      for(int a=0;a<h2;a++)
-       printf("%d ",fX[a]);
-      printf("\n");
-      }*/
-    
-    *(int*)&fX[id*h2+0]^=it;
-
-    /* if(it<513) {
-      for(int a=0;a<h2;a++)
-       printf("%d ",fX[a]);
-      printf("\n");
-      }*/
-
-
-  

     
     

-    /*for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=fX[id*h2+a]^RM1[id*h2+a];
-      fX[id*h2+a+1]=fX[id*h2+a+1]^RM1[id*h2+a+1];
-      fX[id*h2+a+2]=fX[id*h2+a+2]^RM1[id*h2+a+2];
-      fX[id*h2+a+3]=fX[id*h2+a+3]^RM1[id*h2+a+3];
-      }*/
-
-
-

     
     

-    for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=Sbox2[fX[id*h2+a]];
-      fX[id*h2+a+1]=Sbox2[fX[id*h2+a+1]];
-      fX[id*h2+a+2]=Sbox2[fX[id*h2+a+2]];
-      fX[id*h2+a+3]=Sbox2[fX[id*h2+a+3]];
-    }
+      /*for(int a=0;a<h2;a+=4) {
+       fX[a]=Sbox2[fX[a]];
+       fX[a+1]=Sbox2[fX[a+1]];
+       fX[a+2]=Sbox2[fX[a+2]];
+       fX[a+3]=Sbox2[fX[a+3]];
+       }*/

 
 //    rotate(RM1, &RM2[id*h2], h2, Pbox[it]%h2);
     
 
 //    rotate(RM1, &RM2[id*h2], h2, Pbox[it]%h2);
     

-    for(int a=0;a<h2;a+=4) {
-      RM2[a]=RM2[PboxRM[a]];
-      RM2[a+1]=RM2[PboxRM[a+1]];
-      RM2[a+2]=RM2[PboxRM[a+2]];
-      RM2[a+3]=RM2[PboxRM[a+3]];
-    }

     
     

-    for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=fX[id*h2+a]^RM2[id*h2+a];
-      fX[id*h2+a+1]=fX[id*h2+a+1]^RM2[id*h2+a+1];
-      fX[id*h2+a+2]=fX[id*h2+a+2]^RM2[id*h2+a+2];
-      fX[id*h2+a+3]=fX[id*h2+a+3]^RM2[id*h2+a+3];
-    }
+    
+     

 
 
     
 
 
     

-     for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=fX[id*h2+a]^seq_in[ind2+a];
-      fX[id*h2+a+1]=fX[id*h2+a+1]^seq_in[ind2+a+1];
-      fX[id*h2+a+2]=fX[id*h2+a+2]^seq_in[ind2+a+2];
-      fX[id*h2+a+3]=fX[id*h2+a+3]^seq_in[ind2+a+3];
-    }
+      for(int a=0;a<h2;a+=4) {
+       fX[a]=fX[a]^seq_in[ind2+a];
+       fX[a+1]=fX[a+1]^seq_in[ind2+a+1];
+       fX[a+2]=fX[a+2]^seq_in[ind2+a+2];
+       fX[a+3]=fX[a+3]^seq_in[ind2+a+3];
+      }

 
  
 
  

-    for(int a=0;a<h2;a+=4) {
-      seq_out[ind1+a]=fX[id*h2+a];
-      seq_out[ind1+a+1]=fX[id*h2+a+1];
-      seq_out[ind1+a+2]=fX[id*h2+a+2];
-      seq_out[ind1+a+3]=fX[id*h2+a+3];
-    }
-    
-    /*for(int a=0;a<h2;a+=4) {
-      RM1[id*h2+a]=RM1[id*h2+PboxRM[a]];
-      RM1[id*h2+a+1]=RM1[id*h2+PboxRM[a+1]];
-      RM1[id*h2+a+2]=RM1[id*h2+PboxRM[a+2]];
-      RM1[id*h2+a+3]=RM1[id*h2+PboxRM[a+3]];
-    }
-    */
+      for(int a=0;a<h2;a+=4) {
+       seq_out[ind1+a]=fX[a];
+       seq_out[ind1+a+1]=fX[a+1];
+       seq_out[ind1+a+2]=fX[a+2];
+       seq_out[ind1+a+3]=fX[a+3];
+      }
+
+
+
+      for(int a=0;a<h2;a+=4) {
+       RM2[a]=RM2[PboxRM[a]];
+       RM2[a+1]=RM2[PboxRM[a+1]];
+       RM2[a+2]=RM2[PboxRM[a+2]];
+       RM2[a+3]=RM2[PboxRM[a+3]];
+      }
+
+      
+      /*for(int a=0;a<h2;a+=4) {
+       RM1[id*h2+a]=RM1[id*h2+PboxRM[a]];
+       RM1[id*h2+a+1]=RM1[id*h2+PboxRM[a+1]];
+       RM1[id*h2+a+2]=RM1[id*h2+PboxRM[a+2]];
+       RM1[id*h2+a+3]=RM1[id*h2+PboxRM[a+3]];
+       }
+      */

 
     
 
     

-  }
+    }

   }
 
 }
   }
 
 }


@@ -297,146 +274,177 @@ void encrypt(uchar* seq_in, uchar *seq_out, int len,uchar* RM1,int *Pbox, int *P
 
 
 /*  uchar *X=new uchar[h2];
 
 
 /*  uchar *X=new uchar[h2];

-  uchar *fX=new uchar[h2];
-  unsigned int *lX=(unsigned int*)X;
-  unsigned int *lseq_in=(unsigned int*)seq_in;
+    uchar *fX=new uchar[h2];
+    unsigned int *lX=(unsigned int*)X;
+    unsigned int *lseq_in=(unsigned int*)seq_in;

 */
 */

-  uchar X[h2*num];
-  uchar fX[h2*num];
-  uchar RM2[h2*num];
+

 //  unsigned int *lX=(unsigned int*)X;
 //  unsigned int *lseq_in=(unsigned int*)seq_in;
 //  unsigned int *lX=(unsigned int*)X;
 //  unsigned int *lseq_in=(unsigned int*)seq_in;

+
+
+  int loc_len=len/num;

   
   

-#pragma omp parallel for 
-  for(int it=0;it<len;it++) {
-    int ind1=it*h2;
-    int ind2=Pbox[it]*h2;
+#pragma omp parallel for
+  for(int p=0;p<num;p++) {
+

     int id=omp_get_thread_num();
     int id=omp_get_thread_num();

- 
-    for(int a=0;a<h2;a+=4) {
-      X[id*h2+a]=seq_in[ind2+a];
-      X[id*h2+a+1]=seq_in[ind2+a+1];
-      X[id*h2+a+2]=seq_in[ind2+a+2];
-      X[id*h2+a+3]=seq_in[ind2+a+3];
-    }

 
 

-    for(int a=0;a<h2;a+=4){
-      fX[id*h2+a]=Sbox1[X[id*h2+a]];
-      fX[id*h2+a+1]=Sbox1[X[id*h2+a+1]];
-      fX[id*h2+a+2]=Sbox1[X[id*h2+a+2]];
-      fX[id*h2+a+3]=Sbox1[X[id*h2+a+3]];
+    uchar fX[h2];
+
+    
+    uchar RM2[h2];
+    for(int a=0;a<h2;a++) {
+      RM2[a]=RM1[id*h2+a];           

     }
 
     }
 

+  
+    int offset=p*loc_len;

 
 

-/*    for(int a=0;a<h2;a+=4) {
-      fX[a]=fX[a]^RM1[a];
-      fX[a+1]=fX[a+1]^RM1[a+1];
-      fX[a+2]=fX[a+2]^RM1[a+2];
-      fX[a+3]=fX[a+3]^RM1[a+3];
-      }*/
-    rotate(RM1, &RM2[id*h2], h2, Pbox[it]%h2);
-    for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=fX[id*h2+a]^RM2[id*h2+a];
-      fX[id*h2+a+1]=fX[id*h2+a+1]^RM2[id*h2+a+1];
-      fX[id*h2+a+2]=fX[id*h2+a+2]^RM2[id*h2+a+2];
-      fX[id*h2+a+3]=fX[id*h2+a+3]^RM2[id*h2+a+3];
-    }
+    
+    for(int it=offset;it<offset+loc_len;it++) {

 
 

-    for(int a=0;a<h2;a+=4) {
-      seq_out[ind1+a]=Sbox2[fX[id*h2+a]];
-      seq_out[ind1+a+1]=Sbox2[fX[id*h2+a+1]];
-      seq_out[ind1+a+2]=Sbox2[fX[id*h2+a+2]];
-      seq_out[ind1+a+3]=Sbox2[fX[id*h2+a+3]];
-    }
+  

 
 

-    /* for(int a=0;a<h2;a+=4) {
-      RM1[a]=RM1[PboxRM[a]];
-      RM1[a+1]=RM1[PboxRM[a+1]];
-      RM1[a+2]=RM1[PboxRM[a+2]];
-      RM1[a+3]=RM1[PboxRM[a+3]];
+      int ind1=it*h2;
+      int ind2=Pbox[it]*h2;

 
 

-      }*/
+ 
+      for(int a=0;a<h2;a+=4) {
+       fX[a]=seq_in[ind2+a];
+       fX[a+1]=seq_in[ind2+a+1];
+       fX[a+2]=seq_in[ind2+a+2];
+       fX[a+3]=seq_in[ind2+a+3];
+      }

 
 
 
 

-  }
+      for(int a=0;a<h2;a+=4){
+       fX[a]=Sbox1[fX[a]];
+       fX[a+1]=Sbox1[fX[a+1]];
+       fX[a+2]=Sbox1[fX[a+2]];
+       fX[a+3]=Sbox1[fX[a+3]];
+      }

 
 
 
 

+      for(int a=0;a<h2;a+=4) {
+       fX[a]=fX[a]^RM2[a];
+       fX[a+1]=fX[a+1]^RM2[a+1];
+       fX[a+2]=fX[a+2]^RM2[a+2];
+       fX[a+3]=fX[a+3]^RM2[a+3];
+      }
+
+
+
+
+      for(int a=0;a<h2;a+=4) {
+       seq_out[ind1+a]=Sbox2[fX[a]];
+       seq_out[ind1+a+1]=Sbox2[fX[a+1]];
+       seq_out[ind1+a+2]=Sbox2[fX[a+2]];
+       seq_out[ind1+a+3]=Sbox2[fX[a+3]];
+      }

 
 
 
 

+    for(int a=0;a<h2;a+=4) {
+       RM2[a]=RM2[PboxRM[a]];
+       RM2[a+1]=RM2[PboxRM[a+1]];
+       RM2[a+2]=RM2[PboxRM[a+2]];
+       RM2[a+3]=RM2[PboxRM[a+3]];
+      }
+      
+
+
+
+    }
+    
+  }

 }
 
 
 }
 
 

+
+

 template<int h2>
 void decrypt(uchar* seq_in, uchar *seq_out, int len,uchar* RM1,int *Pbox, int *PboxRM, uchar *Inv_Sbox1, uchar *Inv_Sbox2, int debug, int num) {
 
 
   /*uchar *fX=new uchar[h2];
 template<int h2>
 void decrypt(uchar* seq_in, uchar *seq_out, int len,uchar* RM1,int *Pbox, int *PboxRM, uchar *Inv_Sbox1, uchar *Inv_Sbox2, int debug, int num) {
 
 
   /*uchar *fX=new uchar[h2];

-  uchar *Inv_Sbox1=new uchar[256];
-  uchar *Inv_Sbox2=new uchar[256];
+    uchar *Inv_Sbox1=new uchar[256];
+    uchar *Inv_Sbox2=new uchar[256];

   */
   */

-  uchar fX[h2*num];
-  uchar RM2[h2*num];

 
 

-#pragma omp parallel for 
-  for(int it=0;it<len;it++) {

 
 

-    int ind1=it*h2;
-    int ind2=Pbox[it]*h2;
+
+
+
+  int loc_len=len/num;
+  
+#pragma omp parallel for
+  for(int p=0;p<num;p++) {
+

     int id=omp_get_thread_num();
     int id=omp_get_thread_num();

+    uchar RM2[h2];
+    for(int a=0;a<h2;a++) {
+      RM2[a]=RM1[id*h2+a];           //Warning according to the size of h2, we can be outsize of Sbox1[a]
+    }
+    uchar fX[h2];
+  
+    int offset=p*loc_len;

 
 

+    
+    for(int it=offset;it<offset+loc_len;it++) {

 
 
 
 

-    for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=seq_in[ind1+a];
-      fX[id*h2+a+1]=seq_in[ind1+a+1];
-      fX[id*h2+a+2]=seq_in[ind1+a+2];
-      fX[id*h2+a+3]=seq_in[ind1+a+3];
+      int ind1=it*h2;
+      int ind2=Pbox[it]*h2;
+
+
+
+
+      for(int a=0;a<h2;a+=4) {
+       fX[a]=seq_in[ind1+a];
+       fX[a+1]=seq_in[ind1+a+1];
+       fX[a+2]=seq_in[ind1+a+2];
+       fX[a+3]=seq_in[ind1+a+3];

            
            

-    }
-    for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=Inv_Sbox2[fX[id*h2+a]];
-      fX[id*h2+a+1]=Inv_Sbox2[fX[id*h2+a+1]];
-      fX[id*h2+a+2]=Inv_Sbox2[fX[id*h2+a+2]];
-      fX[id*h2+a+3]=Inv_Sbox2[fX[id*h2+a+3]];
-    }
-    /*for(int a=0;a<h2;a+=4) {
-      fX[a]=fX[a]^RM1[a];
-      fX[a+1]=fX[a+1]^RM1[a+1];
-      fX[a+2]=fX[a+2]^RM1[a+2];
-      fX[a+3]=fX[a+3]^RM1[a+3];
+      }
+      for(int a=0;a<h2;a+=4) {
+       fX[a]=Inv_Sbox2[fX[a]];
+       fX[a+1]=Inv_Sbox2[fX[a+1]];
+       fX[a+2]=Inv_Sbox2[fX[a+2]];
+       fX[a+3]=Inv_Sbox2[fX[a+3]];

       }
 
       }
 

-    for(int a=0;a<h2;a+=4) {
-      RM1[a]=RM1[PboxRM[a]];
-      RM1[a+1]=RM1[PboxRM[a+1]];
-      RM1[a+2]=RM1[PboxRM[a+2]];
-      RM1[a+3]=RM1[PboxRM[a+3]];
-      }*/

 
 
 
 

-    rotate(RM1, &RM2[id*h2], h2, Pbox[it]%h2);
-    for(int a=0;a<h2;a+=4) {
-      fX[id*h2+a]=fX[id*h2+a]^RM2[id*h2+a];
-      fX[id*h2+a+1]=fX[id*h2+a+1]^RM2[id*h2+a+1];
-      fX[id*h2+a+2]=fX[id*h2+a+2]^RM2[id*h2+a+2];
-      fX[id*h2+a+3]=fX[id*h2+a+3]^RM2[id*h2+a+3];      
-    }
+      for(int a=0;a<h2;a+=4) {
+       fX[a]=fX[a]^RM2[a];
+       fX[a+1]=fX[a+1]^RM2[a+1];
+       fX[a+2]=fX[a+2]^RM2[a+2];
+       fX[a+3]=fX[a+3]^RM2[a+3];      
+      }

 
     
 
     

-    for(int a=0;a<h2;a+=4) {
-      seq_out[ind2+a]=Inv_Sbox1[fX[id*h2+a]];
-      seq_out[ind2+a+1]=Inv_Sbox1[fX[id*h2+a+1]];
-      seq_out[ind2+a+2]=Inv_Sbox1[fX[id*h2+a+2]];
-      seq_out[ind2+a+3]=Inv_Sbox1[fX[id*h2+a+3]];
-    }
+      for(int a=0;a<h2;a+=4) {
+       seq_out[ind2+a]=Inv_Sbox1[fX[a]];
+       seq_out[ind2+a+1]=Inv_Sbox1[fX[a+1]];
+       seq_out[ind2+a+2]=Inv_Sbox1[fX[a+2]];
+       seq_out[ind2+a+3]=Inv_Sbox1[fX[a+3]];
+      }
+      for(int a=0;a<h2;a+=4) {
+       RM2[a]=RM2[PboxRM[a]];
+       RM2[a+1]=RM2[PboxRM[a+1]];
+       RM2[a+2]=RM2[PboxRM[a+2]];
+       RM2[a+3]=RM2[PboxRM[a+3]];
+       }

 
      
 
      

-  }
+    }

 
 
 
 

+  }

 }
 
 
 }
 
 

+

 int main(int argc, char** argv) {
 
 
 int main(int argc, char** argv) {
 
 


@@ -455,10 +463,10 @@ int main(int argc, char** argv) {
   }
 
 /*  printf("nb times %d\n",nb_test);
   }
 
 /*  printf("nb times %d\n",nb_test);

-  printf("ctr %d\n",ctr);
-  printf("h %d\n",h);
-  printf("lena %d\n",lena);
-  printf("size_buf %d\n",size_buf);
+    printf("ctr %d\n",ctr);
+    printf("h %d\n",h);
+    printf("lena %d\n",lena);
+    printf("size_buf %d\n",size_buf);

 */
   int h2=h*h;
   
 */
   int h2=h*h;
   


@@ -574,12 +582,6 @@ int main(int argc, char** argv) {
   int num=omp_get_max_threads();
   cout<<"num "<<num<<endl;
   
   int num=omp_get_max_threads();
   cout<<"num "<<num<<endl;
   

-  uchar RM1[num*(h * h)];
-  uchar RM2[num*(h * h)];
-  prga(sc, num*(h * h), RM1);
-  for(int i=0;i<num*h2;i++) {
-    RM2[i]=RM1[i];
-  }

   
 
 
   
 
 


@@ -587,10 +589,30 @@ int main(int argc, char** argv) {
   
     
   
   
     
   

-  uchar keyp[16];
-  for (int i = 48; i < 64; i++)
-    keyp[i-48] = DK[i];

 
 

+
+  uchar RM1[num*(h * h)];
+  uchar RM2[num*(h * h)];
+  /*for(int i=0;i<num;i++) {
+    
+    rc4key(&DK[48+i*16], sc, 16);
+    prga(sc, h2, &RM1[h2*i]);
+    for(int a=0;a<h2;a++) {
+      cout<<(int)RM1[h2*i+a]<<" ";
+    }
+    cout<<endl<<endl;
+  }*/
+
+  rc4key(&DK[48], sc, 16);
+  prga(sc, h2*num, RM1);
+
+  rc4key(&DK[64], sc, 16);
+  prga(sc, h2, RM2);
+
+
+
+
+  

 //  cout<<len<<endl;
   int *Pbox=new int[len];
 
 //  cout<<len<<endl;
   int *Pbox=new int[len];
 


@@ -598,15 +620,14 @@ int main(int argc, char** argv) {
   
   int *PboxRM=new int[h2];
 
   
   int *PboxRM=new int[h2];
 

-  rc4keyperm(keyp, len, rp, Pbox, 16);
+  rc4keyperm(&DK[48+16*num], len, rp, Pbox, 16);

 
 

-//  printf("len %d\n",len);
-  for(int i=0;i<len;i++) {
-//    printf("%d \n",Pbox[i]);
-  }
-  

 
 

+  rc4keyperm(RM2, h2, rp, PboxRM, h2);

 
 

+  for(int i=0;i<num*h2;i++) {
+    RM2[i]=RM1[i];
+  }

  
   double time=0;
   double t=TimeStart();
  
   double time=0;
   double t=TimeStart();