#include <fstream>
#include <iostream>
using namespace std;
#include <stdio.h>
#include <stdlib.h>
typedef unsigned char uint8_t;
typedef unsigned short uint16_t;
typedef unsigned int uint32_t;
typedef unsigned long long uint64_t;
typedef struct {
uint16_t hi16;
uint32_t lo32;
}uint48_t;
#define BLOCK 8
#define NUMROUND 16
uint8_t IPTable[] = {
58, 50, 42, 34, 26, 18, 10, 2,
60, 52, 44, 36, 28, 20, 12, 4,
62, 54, 46, 38, 30, 22, 14, 6,
64, 56, 48, 40, 32, 24, 16, 8,
57, 49, 41, 33, 25, 17, 9, 1,
59, 51, 43, 35, 27, 19, 11, 3,
61, 53, 45, 37, 29, 21, 13, 5,
63, 55, 47, 39, 31, 23, 15, 7
};
uint8_t InvertIPTable[] = {
40, 8, 48, 16, 56, 24, 64, 32,
39, 7, 47, 15, 55, 23, 63, 31,
38, 6, 46, 14, 54, 22, 62, 30,
37, 5, 45, 13, 53, 21, 61, 29,
36, 4, 44, 12, 52, 20, 60, 28,
35, 3, 43, 11, 51, 19, 59, 27,
34, 2, 42, 10, 50, 18, 58, 26,
33, 1, 41, 9, 49, 17, 57, 25
};
uint8_t ExpanTable[] = {
32, 1, 2, 3, 4, 5,
4, 5, 6, 7, 8, 9,
8, 9, 10, 11, 12, 13,
12, 13, 14, 15, 16, 17,
16, 17, 18, 19, 20, 21,
20, 21, 22, 23, 24, 25,
24, 25, 26, 27, 28, 29,
28, 29, 30, 31, 32, 1
};
uint8_t SBoxTable[8][4][16] = {
{
{14, 4, 13, 1, 2, 15, 11, 8, 3, 10, 6, 12, 5, 9, 0, 7} ,
{ 0, 15, 7, 4, 14, 2, 13, 1, 10, 6, 12, 11, 9, 5, 3, 8} ,
{ 4, 1, 14, 8, 13, 6, 2, 11, 15, 12, 9, 7, 3, 10, 5, 0} ,
{15, 12, 8, 2, 4, 9, 1, 7, 5, 11, 3, 14, 10, 0, 6, 13}
},
{
{15, 1, 8, 14, 6, 11, 3, 4, 9, 7, 2, 13, 12, 0, 5, 10},
{ 3, 13, 4, 7, 15, 2, 8, 14, 12, 0, 1, 10, 6, 9, 11, 5},
{ 0, 14, 7, 11, 10, 4, 13, 1, 5, 8, 12, 6, 9, 3, 2, 15},
{13, 8, 10, 1, 3, 15, 4, 2, 11, 6, 7, 12, 0, 5, 14, 9}
},
{
{10, 0, 9, 14, 6, 3, 15, 5, 1, 13, 12, 7, 11, 4, 2, 8},
{13, 7, 0, 9, 3, 4, 6, 10, 2, 8, 5, 14, 12, 11, 15, 1},
{13, 6, 4, 9, 8, 15, 3, 0, 11, 1, 2, 12, 5, 10, 14, 7},
{ 1, 10, 13, 0, 6, 9, 8, 7, 4, 15, 14, 3, 11, 5, 2, 12}
},
{
{ 7, 13, 14, 3, 0, 6, 9, 10, 1, 2, 8, 5, 11, 12, 4, 15},
{13, 8, 11, 5, 6, 15, 0, 3, 4, 7, 2, 12, 1, 10, 14, 9},
{10, 6, 9, 0, 12, 11, 7, 13, 15, 1, 3, 14, 5, 2, 8, 4},
{ 3, 15, 0, 6, 10, 1, 13, 8, 9, 4, 5, 11, 12, 7, 2, 14}
},
{
{ 2, 12, 4, 1, 7, 10, 11, 6, 8, 5, 3, 15, 13, 0, 14, 9},
{14, 11, 2, 12, 4, 7, 13, 1, 5, 0, 15, 10, 3, 9, 8, 6},
{ 4, 2, 1, 11, 10, 13, 7, 8, 15, 9, 12, 5, 6, 3, 0, 14},
{11, 8, 12, 7, 1, 14, 2, 13, 6, 15, 0, 9, 10, 4, 5, 3}
},
{
{12, 1, 10, 15, 9, 2, 6, 8, 0, 13, 3, 4, 14, 7, 5, 11},
{10, 15, 4, 2, 7, 12, 9, 5, 6, 1, 13, 14, 0, 11, 3, 8},
{ 9, 14, 15, 5, 2, 8, 12, 3, 7, 0, 4, 10, 1, 13, 11, 6},
{ 4, 3, 2, 12, 9, 5, 15, 10, 11, 14, 1, 7, 6, 0, 8, 13}
},
{
{ 4, 11, 2, 14, 15, 0, 8, 13, 3, 12, 9, 7, 5, 10, 6, 1},
{13, 0, 11, 7, 4, 9, 1, 10, 14, 3, 5, 12, 2, 15, 8, 6},
{ 1, 4, 11, 13, 12, 3, 7, 14, 10, 15, 6, 8, 0, 5, 9, 2},
{ 6, 11, 13, 8, 1, 4, 10, 7, 9, 5, 0, 15, 14, 2, 3, 12}
},
{
{13, 2, 8, 4, 6, 15, 11, 1, 10, 9, 3, 14, 5, 0, 12, 7},
{ 1, 15, 13, 8, 10, 3, 7, 4, 12, 5, 6, 11, 0, 14, 9, 2},
{ 7, 11, 4, 1, 9, 12, 14, 2, 0, 6, 10, 13, 15, 3, 5, 8},
{ 2, 1, 14, 7, 4, 10, 8, 13, 15, 12, 9, 0, 3, 5, 6, 11}
}
};
uint8_t PBoxTable[] = {
16, 7, 20, 21,
29, 12, 28, 17,
1, 15, 23, 26,
5, 18, 31, 10,
2, 8, 24, 14,
32, 27, 3, 9,
19, 13, 30, 6,
22, 11, 4, 25
};
uint8_t PC1[2][28] = {
{57, 49, 41, 33, 25, 17, 9,
1, 58, 50, 42, 34, 26, 18,
10, 2, 59, 51, 43, 35, 27,
19, 11, 3, 60, 52, 44, 36
},
{63, 55, 47, 39, 31, 23, 15,
7, 62, 54, 46, 38, 30, 22,
14, 6, 61, 53, 45, 37, 29,
21, 13, 5, 28, 20, 12, 4
}
};
uint8_t PC2[] = {
14, 17, 11, 24, 1, 5,
3, 28, 15, 6, 21, 10,
23, 19, 12, 4, 26, 8,
16, 7, 27, 20, 13, 2,
41, 52, 31, 37, 47, 55,
30, 40, 51, 45, 33, 48,
44, 49, 39, 56, 34, 53,
46, 42, 50, 36, 29, 32
};
uint8_t ShiftTable[] = {
1, 1, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 1
};
uint48_t SubKeyTable[16];
uint64_t DESEncipher(uint64_t in, uint64_t key);
uint64_t DESDecipher(uint64_t in, uint64_t key);
uint64_t DESCore(uint64_t in, uint48_t *SubKeyTable);
uint8_t BitOnPos(uint64_t in, uint8_t pos, uint8_t numOfBit);
uint64_t IP(uint64_t in);
uint64_t FP(uint64_t in);
void KeyInvert(uint48_t *subkeyTable);
uint32_t LeftRoundShift(uint32_t in, uint8_t len, uint32_t shift);
uint32_t PC1Permut(uint64_t in, uint32_t table);
uint64_t PC2Permut(uint32_t cn, uint32_t dn);
void KeySchedule(uint64_t key, uint48_t *subkeyTable);
uint64_t OneRound(uint64_t in, uint48_t subKey);
uint32_t F(uint32_t in, uint48_t subKey);
uint48_t Expansion(uint32_t in);
uint48_t XORSubKey(uint48_t in, uint48_t subKey);
uint32_t SBox(uint48_t in);
uint8_t SubSBox(uint8_t in, uint8_t tableIndex);
uint32_t PBox(uint32_t in);
uint64_t Swap64(uint64_t in);
int main(){
ifstream plaintext("plaintext", ios::binary);
ofstream ciphertext("ciphertext", ios::binary);
uint64_t block = 0;
uint64_t encipherBlock= 0;
uint64_t decipherBlock= 0;
uint64_t key = 0x5B5A57676A56676ELLU;
char buf[BLOCK];
plaintext.read(buf, BLOCK);
for(int i=0; i<BLOCK; i++){
block <<= 8;
block += buf[i];
}
block = 0x0123456789abcdefLLU;
key = 0x123456789abcdefLLU;
cout<<"plain block is: "<<hex<<block<<endl;
encipherBlock = DESEncipher(block, key);
cout<<"cipher block is: "<<hex<<encipherBlock<<endl;
decipherBlock = DESDecipher(encipherBlock, key);
cout<<"Decipher block result: "<<hex<<decipherBlock<<endl;
KeySchedule(key, SubKeyTable);
uint64_t a = IP(block);
cout<<"a: "<<hex<<a<<endl;
uint32_t l0 = a >> 32;
uint32_t r0 = a;
cout<<"l0: "<<l0<<endl;
cout<<"r0: "<<r0<<endl;
uint48_t e = Expansion(r0);
cout<<"e: "<<hex<<e.hi16<<e.lo32<<endl;
uint64_t temp = 0;
temp = e.hi16;
temp = (temp<<32) + e.lo32;
uint64_t suba = temp ^ 0xB02679B49A5LLU;
cout<<"A: "<<suba<<endl;
uint48_t ins;
ins.hi16 = suba >> 32;
ins.lo32 = suba;
uint32_t B;
B = SBox(ins);
cout<<"B: "<<B<<endl;
uint32_t p;
p = PBox(B);
cout<<"p: "<<p<<endl;
uint64_t one = 0;
cout<<"key0: "<<SubKeyTable[0].hi16<<SubKeyTable[0].lo32<<endl;
one = OneRound(a, SubKeyTable[0]);
cout<<one<<endl;
plaintext.close();
ciphertext.close();
return 0;
}
uint64_t DESEncipher(uint64_t in, uint64_t key){
KeySchedule(key, SubKeyTable);
return DESCore(in, SubKeyTable);
}
void KeyInvert(uint48_t *subkeyTable){
uint32_t lo, hi;
for(int i=0; i<8; i++){
lo = subkeyTable[i].lo32;
hi = subkeyTable[i].hi16;
subkeyTable[i].lo32 = subkeyTable[15-i].lo32;
subkeyTable[i].hi16 = subkeyTable[15-i].hi16;
subkeyTable[15-i].hi16 = hi;
subkeyTable[15-i].lo32 = lo;
}
}
uint64_t DESDecipher(uint64_t in, uint64_t key){
KeySchedule(key, SubKeyTable);
KeyInvert(SubKeyTable);
return DESCore(in, SubKeyTable);
}
uint64_t DESCore(uint64_t in, uint48_t *SubKeyTable){
uint64_t res;
uint48_t subKey;
res = IP(in);
for(int i=0; i<NUMROUND; i++){
subKey = SubKeyTable[i];
res = OneRound(res, subKey);
}
res = Swap64(res);
res = FP(res);
return res;
}
uint8_t BitOnPos(uint64_t in, uint8_t pos, uint8_t numOfBit){
uint64_t temp = 0;
uint64_t hi_indicator = 1;
uint64