]> AND Private Git Repository - Cipher_code.git/blob - OneRoundIoT/OneRound/one_round_hash.cpp
Logo AND Algorithmique Numérique Distribuée

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