OneRoundIoT/OneRound/one_round_hash.cpp

   1 //gcc pixmap_io.c  -c
   2 //g++ -O3 one_round_test.cpp pixmap_io.o  -o one_round_test -std=c++11
   3
   4 //
   5
   6
   7 #include <iostream>
   8 #include <list>
   9 #include<math.h>
  10 #include<stdlib.h>
  11 #include<stdio.h>
  12 #include<string.h>
  13 #include <fstream>
  14 #include <sys/time.h>
  15
  16 /*#include <cryptopp/hex.h>
  17 #include <cryptopp/sha.h>
  18 #include <cryptopp/osrng.h>
  19 #include <cryptopp/secblock.h>
  20 */
  21
  22
  23 extern "C" {
  24   int load_RGB_pixmap(char *filename, int *width, int *height, unsigned char**R_data, unsigned char**G_data, unsigned char**B_data);
  25   void store_RGB_pixmap(char *filename, unsigned char *R_data, unsigned char *G_data, unsigned char *B_data, int width, int height);
  26 }
  27
  28
  29 //using namespace CryptoPP;
  30 using namespace std;
  31
  32
  33 int key_size=256;
  34 int nb_test=1;
  35 int ctr=0;
  36
  37
  38
  39
  40
  41
  42
  43 typedef unsigned char   uchar;
  44
  45
  46 double TimeStart()
  47 {
  48   struct timeval tstart;
  49   gettimeofday(&tstart,0);
  50   return( (double) (tstart.tv_sec + tstart.tv_usec*1e-6) );
  51 }
  52
  53 double TimeStop(double t)
  54 {
  55   struct timeval tend;
  56
  57   gettimeofday(&tend,0);
  58   t = (double) (tend.tv_sec + tend.tv_usec*1e-6) - t;
  59   return (t);
  60 }
  61
  62
  63
  64
  65
  66
  67 void inverse_tables(uchar *tab, int size_tab,uchar *inv_perm_tabs) {
  68
  69   for(int i=0;i<size_tab;i++) {
  70     inv_perm_tabs[tab[i]] = i;
  71   }
  72
  73 }
  74
  75 void inverse_tables_int(int *tab, int size_tab,int *inv_perm_tabs) {
  76
  77   for(int i=0;i<size_tab;i++) {
  78     inv_perm_tabs[tab[i]] = i;
  79   }
  80
  81 }
  82
  83
  84
  85 void rc4key(uchar *key, uchar *sc, int size_DK) {
  86
  87   for(int i=0;i<256;i++) {
  88     sc[i]=i;
  89   }
  90
  91
  92   uchar j0 = 0;
  93   for(int i0=0; i0<256; i0++) {
  94     j0 = (j0 + sc[i0] + key[i0%size_DK] )&0xFF;
  95     uchar tmp = sc[i0];
  96     sc[i0] = sc[j0 ];
  97     sc[j0] = tmp;
  98   }
  99 }
 100
 101
 102
 103 void rc4keyperm(uchar *key,int len, int rp,int *sc, int size_DK) {
 104
 105   //sc=1:len;
 106
 107
 108
 109   for (int i=0;i<len;i++) {
 110     sc[i]=i;
 111   }
 112   for (int it = 0; it < rp; it++) {
 113     int j0 = 1;
 114     for(int i0 = 0; i0<len; i0++) {
 115       j0 = (j0 + sc[i0] + sc[j0] + key[i0%size_DK] )% len;
 116       int tmp = sc[i0];
 117       sc[i0] = sc[j0];
 118       sc[j0] = tmp;
 119     }
 120
 121   }
 122 }
 123
 124 void prga(uchar *sc, int ldata, uchar *r) {
 125   uchar i0=0;
 126   uchar j0=0;
 127
 128   for (int it=0; it<ldata; it++) {
 129     i0 = ((i0+1)&0xFE); //%255);
 130     j0 = (j0 + sc[i0])&0xFF;
 131     uchar tmp = sc[i0];
 132     sc[i0] = sc[j0];
 133     sc[j0] = tmp;
 134     r[it]=sc[(sc[i0]+sc[j0])&0xFF];
 135   }
 136 }
 137
 138 inline uchar  circ(uchar x,int n) {return (x << n) | (x >> (8 - n));}
 139
 140 void diff(uchar *Y, uchar *X, int h) {
 141
 142   if(h==4) {
 143     Y[0] = X[1]^X[2]^X[3];
 144     Y[1] = X[0]^X[2]^X[3];
 145     Y[2] = X[0]^X[1]^X[3];
 146     Y[3] = X[0]^X[1]^X[2];
 147   }
 148   else if(h==8) {
 149     /*Y[0] = X[0]^X[2]^X[3]^X[5]^X[6]^X[7];
 150     Y[1] = X[0]^X[1]^X[3]^X[4]^X[6]^X[7];
 151     Y[2] = X[0]^X[1]^X[2]^X[4]^X[5]^X[7];
 152     Y[3] = X[1]^X[2]^X[3]^X[4]^X[5]^X[6];
 153     Y[4] = X[0]^X[1]^X[5]^X[6]^X[7];
 154     Y[5] = X[1]^X[2]^X[4]^X[6]^X[7];
 155     Y[6] = X[2]^X[3]^X[4]^X[5]^X[7];
 156     Y[7] = X[0]^X[3]^X[4]^X[5]^X[6];
 157     */
 158
 159
 160     uchar X0=X[0];
 161     uchar X1=X[1];
 162     uchar X2=X[2];
 163     uchar X3=X[3];
 164     uchar X4=X[4];
 165     uchar X5=X[5];
 166     uchar X6=X[6];
 167     uchar X7=X[7];
 168
 169     Y[0] = X0^X2^X3^X5^X6^X7;
 170     Y[1] = X0^X1^X3^X4^X6^X7;
 171     Y[2] = X0^X1^X2^X4^X5^X7;
 172     Y[3] = X1^X2^X3^X4^X5^X6;
 173     Y[4] = X0^X1^X5^X6^X7;
 174     Y[5] = X1^X2^X4^X6^X7;
 175     Y[6] = X2^X3^X4^X5^X7;
 176     Y[7] = X0^X3^X4^X5^X6;
 177
 178
 179   }
 180   else if(h==16) {
 181     /*
 182     Y[0] = X[3] ^ X[4] ^ X[6] ^ X[8] ^ X[9] ^ X[13] ^ X[14];
 183   Y[1] = X[2] ^ X[5] ^ X[7] ^ X[8] ^ X[9] ^ X[12] ^ X[15];
 184   Y[2] = X[1] ^ X[4] ^ X[6] ^ X[10] ^ X[11] ^ X[12] ^ X[15];
 185   Y[3] = X[0] ^ X[5] ^ X[7] ^ X[10] ^ X[11] ^ X[13] ^ X[14];
 186   Y[4] = X[0] ^ X[2] ^ X[5] ^ X[8] ^ X[11] ^ X[14] ^ X[15];
 187   Y[5] = X[1] ^ X[3] ^ X[4] ^ X[9] ^ X[10] ^ X[14] ^ X[15];
 188   Y[6] = X[0] ^ X[2] ^ X[7] ^ X[9] ^ X[10] ^ X[12] ^ X[13];
 189   Y[7] = X[1] ^ X[3] ^ X[6] ^ X[8] ^ X[11] ^ X[12] ^ X[13];
 190   Y[8] = X[0] ^ X[1] ^ X[4] ^ X[7] ^ X[10] ^ X[13] ^ X[15];
 191   Y[9] = X[0] ^ X[1] ^ X[5] ^ X[6] ^ X[11] ^ X[12] ^ X[14];
 192   Y[10] = X[2] ^ X[3] ^ X[5] ^ X[6] ^ X[8] ^ X[13] ^ X[15];
 193   Y[11] = X[2] ^ X[3] ^ X[4] ^ X[7] ^ X[9] ^ X[12] ^ X[14];
 194   Y[12] = X[1] ^ X[2] ^ X[6] ^ X[7] ^ X[9] ^ X[11] ^ X[12];
 195   Y[13] = X[0] ^ X[3] ^ X[6] ^ X[7] ^ X[8] ^ X[10] ^ X[13];
 196   Y[14] = X[0] ^ X[3] ^ X[4] ^ X[5] ^ X[9] ^ X[11] ^ X[14];
 197   Y[15] = X[1] ^ X[2] ^ X[4] ^ X[5] ^ X[8] ^ X[10] ^ X[15];
 198     */
 199
 200     uchar X0=X[0];
 201     uchar X1=X[1];
 202     uchar X2=X[2];
 203     uchar X3=X[3];
 204     uchar X4=X[4];
 205     uchar X5=X[5];
 206     uchar X6=X[6];
 207     uchar X7=X[7];
 208     uchar X8=X[8];
 209     uchar X9=X[9];
 210     uchar X10=X[10];
 211     uchar X11=X[11];
 212     uchar X12=X[12];
 213     uchar X13=X[13];
 214     uchar X14=X[14];
 215     uchar X15=X[15];
 216
 217
 218   Y[0] = X3 ^ X4 ^ X6 ^ X8 ^ X9 ^ X13 ^ X14;
 219   Y[1] = X2 ^ X5 ^ X7 ^ X8 ^ X9 ^ X12 ^ X15;
 220   Y[2] = X1 ^ X4 ^ X6 ^ X10 ^ X11 ^ X12 ^ X15;
 221   Y[3] = X0 ^ X5 ^ X7 ^ X10 ^ X11 ^ X13 ^ X14;
 222   Y[4] = X0 ^ X2 ^ X5 ^ X8 ^ X11 ^ X14 ^ X15;
 223   Y[5] = X1 ^ X3 ^ X4 ^ X9 ^ X10 ^ X14 ^ X15;
 224   Y[6] = X0 ^ X2 ^ X7 ^ X9 ^ X10 ^ X12 ^ X13;
 225   Y[7] = X1 ^ X3 ^ X6 ^ X8 ^ X11 ^ X12 ^ X13;
 226   Y[8] = X0 ^ X1 ^ X4 ^ X7 ^ X10 ^ X13 ^ X15;
 227   Y[9] = X0 ^ X1 ^ X5 ^ X6 ^ X11 ^ X12 ^ X14;
 228   Y[10] = X2 ^ X3 ^ X5 ^ X6 ^ X8 ^ X13 ^ X15;
 229   Y[11] = X2 ^ X3 ^ X4 ^ X7 ^ X9 ^ X12 ^ X14;
 230   Y[12] = X1 ^ X2 ^ X6 ^ X7 ^ X9 ^ X11 ^ X12;
 231   Y[13] = X0 ^ X3 ^ X6 ^ X7 ^ X8 ^ X10 ^ X13;
 232   Y[14] = X0 ^ X3 ^ X4 ^ X5 ^ X9 ^ X11 ^ X14;
 233   Y[15] = X1 ^ X2 ^ X4 ^ X5 ^ X8 ^ X10 ^ X15;
 234
 235
 236
 237   }
 238   else if(h==32) {
 239
 240
 241     Y[0]=X[0]^X[1]^X[2]^X[3]^X[4]^X[7]^X[8]^X[10]^X[12]^X[15]^X[16]^X[17]^X[18]^X[20]^X[21]^X[24]^X[25]^X[28]^X[30];
 242     Y[1]=X[0]^ X[1]^X[2]^X[3]^X[5]^X[6]^X[9]^X[11]^X[13]^X[14]^X[16]^X[17]^X[19]^X[20]^X[21]^ X[24]^X[25]^X[29]^X[31];
 243     Y[2]=X[0]^X[1]^X[2]^X[3]^X[5]^X[6]^X[8]^X[10]^X[13]^X[14]^X[16]^X[18]^X[19]^X[22]^X[23]^X[26]^X[27]^X[28]^X[30];
 244     Y[3]=X[0]^X[1]^X[2]^X[3]^X[4]^X[7]^X[9]^X[11]^X[12]^X[15]^X[17]^X[18]^X[19]^X[22]^X[23]^X[26]^X[27]^X[29]^X[31];
 245     Y[4]=X[0]^X[3]^X[5]^X[6]^X[7]^X[10]^X[11]^ X[12]^X[13]^X[14]^ X[15]^X[16]^X[19]^X[21]^X[23]^  X[25]^X[27]^X[30]^X[31];
 246     Y[5]=X[1]^X[2]^X[4]^X[6]^X[7]^X[10]^X[11]^X[12]^X[13]^X[14]^X[16 ]^X[17]^X[18]^X[20]^X[22]^X[24]^X[26]^X[30]^X[31];
 247     Y[6]=X[1]^X[2]^X[4]^X[5]^X[7]^X[8]^X[9]^X[12]^X[13]^X[14]^X[15]^ X[17]^X[18]^X[21]^X[23]^X[25]^X[27]^X[28]^X[29];
 248     Y[7]=X[0]^X[3]^X[4]^X[5]^X[6]^X[9 ]^X[9]^X[12]^X[13]^X[14]^X[15]^X[16]^X[19]^X[20]^X[22]^X[24]^X[26]^X[28]^X[29];
 249     Y[8]=X[0]^X[2]^X[6]^X[7]^X[8]^X[10]^X[11]^X[14]^X[15]^X[16]^X[18]^X[21]^X[22]^X[25]^X[26];
 250     Y[9]=X[1]^ X[3]^X[6]^X[7]^X[9]^X[10]^X[11]^X[14]^X[15]^X[17]^X[19]^X[20]^X[23]^X[24]^X[27];
 251     Y[10]=X[0]^X[2]^X[4]^X[5]^X[8]^X[9]^X[10]^X[12]^X[13]^X[16]^X[18]^X[20]^X[23]^ X[24]^X[27];
 252     Y[11]=X[1]^X[3]^X[4]^X[5]^X[8]^X[9]^X[11]^X[12]^X[13]^X[17]^X[19]^X[21]^X[22]^X[25]^X[26];
 253     Y[12]=X[0]^X[3]^X[4]^X[5]^X[6]^X[7]^X[10]^X[11]^X[13]^X[14]^X[15]^X[16]^X[19]^X[21]^X[23]^X[25]^X[27]^X[30]^X[31];
 254     Y[13]=X[1]^X[2]^X[4]^X[5]^X[6]^X[7]^X[10]^X[11]^X[12]^X[14]^X[15]^X[17]^ X[18]^X[20]^X[22]^X[24]^X[26]^X[30]^X[31];
 255     Y[14]=X[1]^X[2]^X[4]^X[5]^X[6]^X[7]^X[8]^X[9]^X[12]^X[13]^X[15]^X[17]^X[18]^X[21]^X[23]^X[25]^X[27]^X[28]^X[29];
 256     Y[15]=X[0]^X[3]^X[4]^X[5]^X[6]^X[7]^X[8]^X[9]^X[12]^X[13]^X[14]^X[16]^X[19]^X[20]^X[22]^ X[24]^X[26]^X[28]^X[29];
 257     Y[16]=X[0]^X[1]^X[2]^X[4]^X[8 ]^X[8]^X[10]^X[13 ]^X[15]^X[16]^X[17]^X[18]^X[19]^X[20]^X[21]^X[24]^X[25]^X[28]^X[30];
 258     Y[17]=X[0]^X[1]^X[3]^X[5]^X[6]^X[9]^X[11]^X[13]^X[14]^X[16]^X[17]^X[18]^X[19]^X[20]^X[21]^X[24]^X[25]^X[29]^X[31];
 259     Y[18]=X[0]^X[2]^X[3]^X[5]^X[6]^X[8]^X[10]^X[13]^X[14]^X[16]^X[17]^X[18]^X[19]^X[22]^X[23]^X[26]^X[27]^X[28]^X[30];
 260     Y[19]=X[1]^X[2]^X[3]^X[4]^X[7]^X[9]^X[11]^X[12]^X[15]^X[16]^X[17]^X[18]^X[19]^X[22]^X[23]^X[26]^X[28 ]^X[29]^X[31];
 261     Y[20]=X[0]^X[1]^X[5]^X[7]^X[10 ]^X[10]^X[13]^X[15]^X[16]^X[17]^X[20]^X[21]^X[23]^X[29]^X[30];
 262     Y[21]=X[0]^X[1]^X[4]^X[6]^X[8]^X[11]^X[12]^X[14]^X[16]^X[17]^X[20]^X[21]^X[22]^X[28]^X[31];
 263     Y[22]=X[2]^X[3]^X[5]^X[7]^X[8]^X[11]^X[13]^X[15]^X[18]^X[19]^X[21]^X[22]^X[23]^X[28]^X[31];
 264     Y[23]=X[2]^X[3]^X[4]^X[6]^X[9]^X[10]^X[12]^X[14]^ X[18]^X[19]^X[20]^X[22]^X[23]^X[29]^X[30];
 265     Y[24]=X[0]^X[1]^X[5]^X[7]^X[9]^X[10]^X[13]^X[15]^X[16]^X[17]^X[24]^X[25]^X[27]^X[29]^X[30];
 266     Y[25]=X[0]^X[1]^X[4]^X[6]^X[8]^X[11]^X[12]^X[14]^X[16]^X[17]^X[24]^X[25]^X[26]^X[28]^X[31];
 267     Y[26]=X[2]^X[3]^X[5]^X[7]^X[8]^X[11]^X[13]^X[15]^X[18]^X[19]^X[25]^X[26]^X[27]^X[28]^ X[31];
 268     Y[27]=X[2]^X[3]^X[4]^X[6]^X[9]^X[10]^X[12]^X[14]^X[18]^X[19]^X[24]^X[26]^X[27]^X[29]^X[30];
 269     Y[28]=X[0]^X[2]^X[6]^X[7]^X[14]^X[15]^X[16]^X[18]^X[21]^X[22]^X[25]^X[26]^X[28]^X[30]^X[31];
 270     Y[29]=X[2]^X[3]^X[6]^X[7]^X[14]^X[15]^X[17]^X[19]^X[20]^X[23]^X[24]^X[27]^X[29]^X[30]^X[31];
 271     Y[30]=X[1]^X[2]^X[4]^X[5]^X[12]^X[13]^X[16]^X[18]^X[20]^X[23]^X[24]^X[27]^X[28]^X[29]^X[30];
 272     Y[31]=X[2]^X[3]^X[4]^X[5]^X[12]^X[13]^X[17]^X[19]^X[21]^X[22]^X[25]^X[26]^X[28]^X[29]^X[31];
 273   }
 274
 275 }
 276
 277
 278
 279
 280 //the proposed hash function, which is based on DSD structure. Sensitivity is ensured by employing the binary diffusion
 281
 282 void hash_DSD_BIN(uchar* seq_in, uchar* RM1,int len, int *PboxRM, uchar *Sbox1,  int h) {
 283
 284
 285   // Goal: Calculate the hash value
 286   // Output: RM (hash value)
 287
 288 //  uchar *X=new uchar[h2];
 289 //  uchar *fX=new uchar[h2];
 290   uchar X[h];
 291   uchar fX[h];
 292   uchar fX2[h];
 293   int ind2;
 294
 295
 296   for(int it=0;it<len;it++) {
 297
 298     ind2=it*h;
 299
 300     // Mix with dynamic RM
 301
 302     for(int a=0;a<h;a+=4) {
 303       fX[a]=RM1[a]^seq_in[ind2+a];
 304       fX[a+1]=RM1[a+1]^seq_in[ind2+a+1];
 305       fX[a+2]=RM1[a+2]^seq_in[ind2+a+2];
 306       fX[a+3]=RM1[a+3]^seq_in[ind2+a+3];
 307     }
 308
 309     // First Diffusion Operation
 310     diff(fX2,fX,h);
 311
 312
 313     // Substitution  Operation
 314     for(int a=0;a<h;a+=4) {
 315       fX[a]=Sbox1[fX2[a]];           //Warning according to the size of h2, we can be outsize of Sbox1[a]
 316       fX[a+1]=Sbox1[fX2[a+1]];
 317       fX[a+2]=Sbox1[fX2[a+2]];
 318       fX[a+3]=Sbox1[fX2[a+3]];
 319     }
 320
 321     // Second Diffusion Operation
 322
 323     /*for(int a=0;a<h;a+=4) {
 324       fX2[a]=circ(X[a],3);
 325       fX2[a+1]=circ(X[a+1],5);
 326       fX2[a+2]=circ(X[a+2],1);
 327       fX2[a+3]=circ(X[a+3],6);
 328       }*/
 329
 330     diff(fX2,fX,h);
 331
 332 // update RM and mix it with hashed block
 333     for(int a=0;a<h;a+=4) {
 334       RM1[a]=fX2[a]^RM1[PboxRM[a]];
 335       RM1[a+1]=fX2[a+1]^RM1[PboxRM[a+1]];
 336       RM1[a+2]=fX2[a+2]^RM1[PboxRM[a+2]];
 337       RM1[a+3]=fX2[a+3]^RM1[PboxRM[a+3]];
 338       }
 339 /*
 340      for(int a=0;a<h;a+=4) {
 341       RM1[a]=RM1[PboxRM[a]];
 342       RM1[a+1]=RM1[PboxRM[a+1]];
 343       RM1[a+2]=RM1[PboxRM[a+2]];
 344       RM1[a+3]=RM1[PboxRM[a+3]];
 345     }
 346
 347      for(int a=0;a<h;a+=4) {
 348       RM1[a]=fX2[a]^RM1[a];
 349       RM1[a+1]=fX2[a+1]^RM1[a+1];
 350       RM1[a+2]=fX2[a+2]^RM1[a+2];
 351       RM1[a+3]=fX2[a+3]^RM1[a+3];
 352     }
 353 */
 354
 355
 356
 357   }
 358
 359
 360 }
 361
 362
 363
 364
 365
 366
 367
 368
 369
 370
 371 int main(int argc, char** argv) {
 372
 373
 374   int h=16;
 375   int lena=0;
 376   int size_buf=1;
 377   int change=0;
 378
 379
 380   for(int i=1; i<argc; i++){
 381     if(strncmp(argv[i],"nb",2)==0)    nb_test = atoi(&(argv[i][2]));    //nb of test
 382     if(strncmp(argv[i],"ctr",3)==0) ctr = atoi(&(argv[i][3]));          //CTR ? 1  otherwise CBC like
 383     if(strncmp(argv[i],"h",1)==0) h = atoi(&(argv[i][1]));          //size of block
 384     if(strncmp(argv[i],"sizebuf",7)==0) size_buf = atoi(&(argv[i][7]));          //SIZE of the buffer
 385     if(strncmp(argv[i],"lena",4)==0) lena = atoi(&(argv[i][4]));          //Use Lena or buffer
 386     if(strncmp(argv[i],"c",1)==0) change = atoi(&(argv[i][1]));          //Use Lena or buffer
 387   }
 388
 389
 390   cout<<size_buf<<endl;
 391
 392   int seed=time(NULL);
 393 //  cout<<seed<<endl;
 394   srand48(12);
 395
 396   uchar Secretkey[key_size];
 397
 398   uchar counter[key_size];
 399
 400   for(int i=0;i<key_size;i++) {
 401     Secretkey[i]=lrand48()&0xFF;
 402     counter[i]=lrand48()&0xFF;
 403   }
 404
 405
 406   int size = 64;
 407   uchar DK[size];
 408
 409
 410
 411
 412   int width;
 413   int height;
 414
 415   uchar *data_R, *data_G, *data_B;
 416   int imsize;
 417   uchar *buffer;
 418
 419   if(lena==1) {
 420     load_RGB_pixmap("lena.ppm", &width, &height, &data_R, &data_G, &data_B);
 421     imsize=width*height*3;
 422 //  load_RGB_pixmap("No_ecb_mode_picture.ppm", &width, &height, &data_R, &data_G, &data_B);
 423   }
 424   else {
 425     imsize=size_buf;
 426     buffer=new uchar[imsize];
 427     for(int i=0;i<imsize;i++) {
 428       buffer[i]=lrand48();
 429     }
 430   }
 431
 432
 433
 434
 435   uchar* seq= new uchar[imsize];
 436   uchar* seq2= new uchar[imsize];
 437
 438   int oneD;
 439   if(lena) {
 440     oneD=width*height;
 441     for(int i=0;i<oneD;i++) {
 442       seq[i]=data_R[i];
 443       seq[oneD+i]=data_G[i];
 444       seq[2*oneD+i]=data_B[i];
 445     }
 446   }
 447   else {
 448     oneD=imsize;
 449     for(int i=0;i<oneD;i++) {
 450       seq[i]=buffer[i];
 451     }
 452   }
 453
 454
 455   if(change) {
 456
 457     seq[4]++;
 458   }
 459   if(change==2) {
 460
 461     seq[9]++;
 462   }
 463
 464
 465
 466
 467   int total_len=imsize;
 468   int rp=1;
 469   int len= total_len/h;
 470   cout<<len<<endl;
 471
 472
 473   uchar *mix=new uchar[256];
 474
 475
 476
 477
 478   for (int i = 0; i < 256 ; i++) {
 479     mix[i]=Secretkey[i]^counter[i];
 480   }
 481
 482
 483 //  cout<<"hash "<<endl;
 484   for (int i = 0; i < 64 ; i++) {
 485 //    DK[i]=digest[i];
 486     DK[i]=mix[i];
 487   }
 488
 489
 490
 491
 492   int *PboxRM=new int[h];
 493   uchar Sbox1[256];
 494   uchar sc[256];
 495   uchar RM1[h];
 496   uchar RM2[h];
 497
 498
 499
 500
 501   double time=0;
 502   double t=TimeStart();
 503   rc4key(DK, Sbox1, 8);
 504
 505
 506   rc4key(&DK[8], sc, 8);
 507
 508
 509   prga(sc, h, RM1);
 510
 511
 512
 513   rc4keyperm(&DK[16], h, rp, PboxRM, 8);
 514
 515   time+=TimeStop(t);
 516   cout<<"Time initializaton "<<time<<endl;
 517
 518
 519
 520
 521
 522
 523
 524
 525
 526
 527
 528   for(int i=0;i<h;i++){
 529     RM2[i]=RM1[i];
 530   }
 531
 532   cout<<"imsize "<<imsize<<endl;
 533
 534 /*  for(int i=0;i<imsize;i++){
 535     cout<<(int)seq[i]<<" ";
 536   }
 537   cout<<endl;
 538 */
 539
 540   time=0;
 541   t=TimeStart();
 542   for(int i=0;i<nb_test;i++)
 543   {
 544     hash_DSD_BIN(seq, RM1,len,PboxRM,Sbox1,h);
 545   }
 546
 547
 548
 549
 550   time+=TimeStop(t);
 551   cout<<"Hash Time  "<<time<<endl;
 552   cout<<(double)imsize*nb_test/time<<"\t";
 553
 554   for(int i=0;i<h;i++){
 555     cout<<(int)RM1[i]<<" ";
 556   }
 557   cout<<endl;
 558
 559
 560
 561   return 0;
 562 }