-
+typedef __uint64_t mylong;
typedef unsigned char uchar;
return x;
}
-ulong xorshift64(ulong t)
+
+mylong xorseed;
+
+mylong xorshift64()
{
/* Algorithm "xor" from p. 4 of Marsaglia, "Xorshift RNGs" */
- ulong x = t;
+ mylong x = xorseed;
x ^= x >> 12; // a
x ^= x << 25; // b
x ^= x >> 27; // c
- return x;
+ return xorseed=x;
}
+/*
__uint128_t g_lehmer64_state;
inline uint64_t splitmix64_stateless(uint64_t index) {
return g_lehmer64_state >> 64;
}
-
+*/
template<int h>
-void encrypt_ecb_prng(uchar* seq_in, uchar *seq_out, int len,uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, ulong myrand, int debug) {
+void encrypt_ecb_prng(uchar* seq_in, uchar *seq_out, int len,uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, mylong myrand, int debug) {
uchar X[h];
uchar Y[h];
uchar fX[h];
uchar gY[h];
- uchar RM1[h];
- uchar RM2[h];
+ uchar *RM1=&RM[0];
+ uchar *RM2=&RM[h];
uchar tmp[h];
- ulong *rm1=(ulong*)RM1;
- ulong *rm2=(ulong*)RM2;
+ mylong *rm1=(mylong*)RM1;
+ mylong *rm2=(mylong*)RM2;
+
+
for(int it=0;it<len/2;it++) {
int ind1=Pbox[it]*h;
for(int a=0;a<(h>>3);a++) {
- myrand=lehmer64();
+ myrand=xorshift64();
rm1[a]=myrand;
- myrand=lehmer64();
+ myrand=xorshift64();
rm2[a]=myrand;
}
template<int h>
-void decrypt_ecb_prng(uchar* seq_in, uchar *seq_out, int len, uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *Inv_Sbox1, uchar *Inv_Sbox2, ulong myrand, int debug) {
+void decrypt_ecb_prng(uchar* seq_in, uchar *seq_out, int len, uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *Inv_Sbox1, uchar *Inv_Sbox2, mylong myrand, int debug) {
uchar invfX[h];
uchar invgY[h];
uchar RM1[h];
uchar RM2[h];
uchar tmp[h];
- ulong *rm1=(ulong*)RM1;
- ulong *rm2=(ulong*)RM2;
+ mylong *rm1=(mylong*)RM1;
+ mylong *rm2=(mylong*)RM2;
+
for(int it=0;it<len/2;it++) {
int ind1=Pbox[it]*h;
for(int a=0;a<(h>>3);a++) {
- myrand=lehmer64();
+ myrand=xorshift64();
rm1[a]=myrand;
- myrand=lehmer64();
+ myrand=xorshift64();
rm2[a]=myrand;
}
template<int h>
-void encrypt_ecb_rm(uchar* seq_in, uchar *seq_out, int len,uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, ulong myrand, int debug) {
+void encrypt_ecb_rm(uchar* seq_in, uchar *seq_out, int len,uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *IV, int debug) {
uchar X[h];
uchar Y[h];
uchar fX[h];
uchar gY[h];
- uchar RM1[h];
- uchar RM2[h];
+ uchar IV1[h];
+ uchar IV2[h];
+ uchar *RM1;
+ uchar *RM2;
uchar tmp[h];
- ulong *rm1=(ulong*)RM1;
- ulong *rm2=(ulong*)RM2;
+
+
+
+ for(int a=0;a<h;a+=4) {
+ IV1[a]=IV[a];
+ IV1[a+1]=IV[a+1];
+ IV1[a+2]=IV[a+2];
+ IV1[a+3]=IV[a+3];
+ }
+
+
+ for(int a=0;a<h;a+=4) {
+ IV2[a]=IV[h+a];
+ IV2[a+1]=IV[h+a+1];
+ IV2[a+2]=IV[h+a+2];
+ IV2[a+3]=IV[h+a+3];
+
+ }
+
+
+
+
+
for(int it=0;it<len/2;it++) {
int ind1=Pbox[it]*h;
int ind2=Pbox[it+len/2]*h;
-
- for(int a=0;a<(h>>3);a++) {
- myrand=lehmer64();
- rm1[a]=myrand;
- myrand=lehmer64();
- rm2[a]=myrand;
- }
+ RM1=&RM[PboxSRM[it]*h];
+ RM2=&RM[h*h+PboxSRM[len/2-it]*h];
+
for(int a=0;a<h;a+=4) {
template<int h>
-void decrypt_ecb_rm(uchar* seq_in, uchar *seq_out, int len, uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *Inv_Sbox1, uchar *Inv_Sbox2, ulong myrand, int debug) {
+void decrypt_ecb_rm(uchar* seq_in, uchar *seq_out, int len, uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *Inv_Sbox1, uchar *Inv_Sbox2, uchar *IV, int debug) {
uchar invfX[h];
uchar invgY[h];
uchar fX[h];
uchar gY[h];
- uchar RM1[h];
- uchar RM2[h];
+ uchar IV1[h];
+ uchar IV2[h];
+ uchar *RM1;
+ uchar *RM2;
uchar tmp[h];
- ulong *rm1=(ulong*)RM1;
- ulong *rm2=(ulong*)RM2;
+
+
+
+ for(int a=0;a<h;a+=4) {
+ IV1[a]=IV[a];
+ IV1[a+1]=IV[a+1];
+ IV1[a+2]=IV[a+2];
+ IV1[a+3]=IV[a+3];
+ }
+
+
+ for(int a=0;a<h;a+=4) {
+ IV2[a]=IV[h+a];
+ IV2[a+1]=IV[h+a+1];
+ IV2[a+2]=IV[h+a+2];
+ IV2[a+3]=IV[h+a+3];
+
+ }
+
+
for(int it=0;it<len/2;it++) {
int ind1=Pbox[it]*h;
int ind2=Pbox[it+len/2]*h;
- for(int a=0;a<(h>>3);a++) {
- myrand=lehmer64();
- rm1[a]=myrand;
- myrand=lehmer64();
- rm2[a]=myrand;
- }
+ RM1=&RM[PboxSRM[it]*h];
+ RM2=&RM[h*h+PboxSRM[len/2-it]*h];
uchar gY[h];
uchar IV1[h];
uchar IV2[h];
- uchar RM1[h];
- uchar RM2[h];
+ uchar *RM1=&RM[0];
+ uchar *RM2=&RM[h];
uchar tmp[h];
- ulong *rm1=(ulong*)RM1;
- ulong *rm2=(ulong*)RM2;
+ mylong *rm1=(mylong*)RM1;
+ mylong *rm2=(mylong*)RM2;
+
+
+
for(int it=0;it<len/2;it++) {
int ind1=Pbox[it]*h;
for(int a=0;a<(h>>3);a++) {
- myrand=lehmer64();
+ myrand=xorshift64();
rm1[a]=myrand;
- myrand=lehmer64();
+ myrand=xorshift64();
rm2[a]=myrand;
}
uchar gY[h];
uchar IV1[h];
uchar IV2[h];
- uchar RM1[h];
- uchar RM2[h];
+ uchar *RM1=&RM[0];
+ uchar *RM2=&RM[h];
+// uchar RM1[h];
+// uchar RM2[h];
uchar tmp[h];
- ulong *rm1=(ulong*)RM1;
- ulong *rm2=(ulong*)RM2;
+ mylong *rm1=(mylong*)RM1;
+ mylong *rm2=(mylong*)RM2;
for(int it=0;it<len/2;it++) {
int ind1=Pbox[it]*h;
for(int a=0;a<(h>>3);a++) {
- myrand=lehmer64();
+ myrand=xorshift64();
rm1[a]=myrand;
- myrand=lehmer64();
+ myrand=xorshift64();
rm2[a]=myrand;
}
uchar RM[h*h*2+256];
uchar IV[2*h];
- ulong myrand=0;
-
+ mylong myrand=0;
+
double time_encrypt=0;
double time_decrypt=0;
inverse_tables(Sbox2,256,Inv_Sbox2);
-
- lehmer64_seed(myrand);
+ xorseed=myrand;
+// lehmer64_seed(myrand);
time_encrypt=0;
t=TimeStart();
if(cbcrm)
encrypt_cbc_rm<4>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<4>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<4>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<4>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
}
if(cbcrm)
encrypt_cbc_rm<8>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<8>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<8>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<8>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
}
if(cbcrm)
encrypt_cbc_rm<16>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<16>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<16>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<16>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
}
if(cbcrm)
encrypt_cbc_rm<32>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<32>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<32>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<32>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
}
if(cbcrm)
encrypt_cbc_rm<64>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<64>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<64>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<64>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
if(cbcrm)
encrypt_cbc_rm<128>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<128>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<128>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<128>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
if(cbcrm)
encrypt_cbc_rm<256>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbrm)
- encrypt_ecb_rm<256>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
+ encrypt_ecb_rm<256>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,IV,0);
if(ecbprng)
encrypt_ecb_prng<256>(seq, seq2,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,myrand,0);
store_RGB_pixmap("lena2.ppm", data_R, data_G, data_B, width, height);
}
-
- lehmer64_seed(myrand);
+
+ xorseed=myrand;
+ // lehmer64_seed(myrand);
time_decrypt=0;
t=TimeStart();
switch(h) {
if(cbcrm)
decrypt_cbc_rm<4>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<4>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<4>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<4>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
if(cbcrm)
decrypt_cbc_rm<8>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<8>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<8>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<8>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
if(cbcrm)
decrypt_cbc_rm<16>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<16>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<16>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<16>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
if(cbcrm)
decrypt_cbc_rm<32>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<32>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<32>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<32>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
if(cbcrm)
decrypt_cbc_rm<64>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<64>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<64>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<64>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
if(cbcrm)
decrypt_cbc_rm<128>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<128>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<128>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<128>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
if(cbcrm)
decrypt_cbc_rm<256>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbrm)
- decrypt_ecb_rm<256>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
+ decrypt_ecb_rm<256>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,IV,0);
if(ecbprng)
decrypt_ecb_prng<256>(seq2,seq,len,RM,Pbox,PboxSRM,Sbox1,Sbox2,Inv_Sbox1,Inv_Sbox2,myrand,0);
}
}
-
return 0;
}