1、密钥产生：

  1. #include <bits/stdc++.h>
  2. using namespace std;
  4. int permChoiceOne[] = {           //PC-1 置换选择1
  5.    57, 49, 41, 33, 25, 17, 9 ,
  6.    1 , 58, 50, 42, 34, 26, 18,
  7.    10, 2 , 59, 51, 43, 35, 27,
  8.    19, 11, 3 , 60, 52, 44, 36,
  9.    63, 55, 47, 39, 31, 23, 15,
  10.    7 , 62, 54, 46, 38, 30, 22,
  11.    14, 6 , 61, 53, 45, 37, 29,
  12.    21, 13, 5 , 28, 20, 12, 4 };
  14. int permChoiceTwo[] = {             //PC-2 置换选择2
  15.    14, 17, 11, 24, 1 , 5 , 3 , 28,
  16.    15, 6 , 21, 10, 23, 19, 12, 4 ,
  17.    26, 8 , 16, 7 , 27, 20, 13, 2 ,
  18.    41, 52, 31, 37, 47, 55, 30, 40,
  19.    51, 45, 33, 48, 44, 49, 39, 56,
  20.    34, 53, 46, 42, 50, 36, 29, 32 };
  23. //左移次数确定
  24. int leftShiftTable[] = { 1, 1, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 1 };
  26. string rotateSubKey(string s, int rot)            // rot是左移旋转次数
  27. {
  28.     return s.substr(rot, s.length() - rot) + s.substr(0, rot); //截取字符
  29. }
  31. string firstPermute(string input)
  32. {
  33.     string res = "";
  34.     for (int i = 0; i < 56; i++)
  35.     {
  36.         res += input[permChoiceOne[i] - 1];
  37.     }
  38.     return res;
  39. }
  41. string secondPermute(string input)
  42. {
  43.     string res = "";
  44.     for (int i = 0; i < 48; i++)
  45.     {
  46.         res += input[permChoiceTwo[i] - 1];
  47.     }
  48.     return res;
  49. }
  51. void genKeys(string left, string right)
  52. {
  53.     ofstream fout;
  54.     fout.open("keygen.txt");        
  56.     for (int i = 0; i < 16; i++)
  57.     {
  58.         left = rotateSubKey(left, leftShiftTable[i]);
  59.         right = rotateSubKey(right, leftShiftTable[i]);
  61.         string key = secondPermute(left + right);
  63.         cout << "key " << i + 1 << " \t: " << key << endl; 
  64.         fout << key << endl;          // 保存到文件
  65.     }
  66. }
  68. int main()
  69. {
  70.     unsigned long long hexkey;
  71.     cout << "\n 请输入64bit的数字(16进制,共16个数字) : ";
  72.     cin >> hex >> hexkey; 
  74.     string key = bitset<64>(hexkey).to_string(); // 16 -> 2
  75.     cout << "二进制密钥 \t: " << key << endl;
  77.     key = firstPermute(key);
  78.     cout << "PC-1 key (k+) \t: " << key << endl;
  80.     cout << "\nSubKeys: " << endl;
  81.     genKeys(key.substr(0, 28), key.substr(28, 28));
  83.     cout << endl << endl;
  84. }

2、S盒变换

  1. #include <bits/stdc++.h>
  2. using namespace std;
  4. int expPermute[] = {        //扩展置换E
  5.     32, 1 , 2 , 3 , 4 , 5 ,
  6.     4 , 5 , 6 , 7 , 8 , 9 ,
  7.     8 , 9 , 10, 11, 12, 13,
  8.     12, 13, 14, 15, 16, 17,
  9.     16, 17, 18, 19, 20, 21,
  10.     20, 21, 22, 23, 24, 25,
  11.     24, 25, 26, 27, 28, 29,
  12.     28, 29, 30, 31, 32, 1 };
  14. string expansionPermute(string input)
  15. {
  16.     string res = "";
  17.     for(int i=0; i<48; i++)
  18.     {
  19.         res += input[expPermute[i]-1];
  20.     }
  21.     return res;
  22. }
  24. string XOR(string input1, string input2)
  25. {
  26.     string res = "";
  27.     for(int i=0; i<input1.length(); i++)
  28.     {
  29.         res += (input1[i] == input2[i]) ? "0" : "1";
  30.     }
  31.     return res;
  32. }
  34. int main()
  35. {
  36.     int i; 
  37.     unsigned long long hexInput;
  38.     string Ki; 
  39.     ifstream fin;
  41.     cout << "\n输入第几轮(i) : ";
  42.     cin >> i;
  44.     cout << "Enter 64-bit (i-1)th round output in hex: " ;
  45.     cin >> hex >> hexInput;
  46.     string input = bitset<64>(hexInput).to_string();
  48.     fin.open("keygen.txt");
  49.     for(int j=1; j<=i; j++)
  50.     {
  51.         fin >> Ki;
  52.     }
  55.     if(Ki.length() == 0)
  56.     {
  57.         cout << "\n 错误 \n" << endl;
  58.         exit(1);
  59.     }
  61.     cout << "\n64-bit 二进制输入 = " << input << endl ;
  62.     cout << " (Ki) = " << Ki << endl ;
  64.     string Ri_1 = input.substr(32,32); 
  65.     cout << "\n Ri_1 = " << Ri_1 << endl;
  67.     string R48 = expansionPermute(Ri_1);
  68.     cout << "Ri_1 变换后 = " << R48 << endl;
  70.     string sBoxInput = XOR(R48, Ki);
  71.     cout << "\n S变换后 : " << sBoxInput << endl << endl;
  72. }

3、Ri

  1. #include <bits/stdc++.h>
  2. using namespace std;
  4. unsigned int sBoxes[8][4][16] = {
  5.     //S1
  6.    {14,4,13,1,2,15,11,8,3,10,6,12,5,9,0,7,
  7.     0,15,7,4,14,2,13,1,10,6,12,11,9,5,3,8,
  8.     4,1,14,8,13,6,2,11,15,12,9,7,3,10,5,0,
  9.     15,12,8,2,4,9,1,7,5,11,3,14,10,0,6,13},
  10.     //S2
  11.     {15,1,8,14,6,11,3,4,9,7,2,13,12,0,5,10,
  12.     3,13,4,7,15,2,8,14,12,0,1,10,6,9,11,5,
  13.     0,14,7,11,10,4,13,1,5,8,12,6,9,3,2,15,
  14.     13,8,10,1,3,15,4,2,11,6,7,12,0,5,14,9},
  15.     //S3
  16.     {10,0,9,14,6,3,15,5,1,13,12,7,11,4,2,8,
  17.     13,7,0,9,3,4,6,10,2,8,5,14,12,11,15,1,
  18.     13,6,4,9,8,15,3,0,11,1,2,12,5,10,14,7,
  19.     1,10,13,0,6,9,8,7,4,15,14,3,11,5,2,12},
  20.     //S4
  21.     {7,13,14,3,0,6,9,10,1,2,8,5,11,12,4,15,
  22.     13,8,11,5,6,15,0,3,4,7,2,12,1,10,14,9,
  23.     10,6,9,0,12,11,7,13,15,1,3,14,5,2,8,4,
  24.     3,15,0,6,10,1,13,8,9,4,5,11,12,7,2,14},
  25.     //S5
  26.     {2,12,4,1,7,10,11,6,8,5,3,15,13,0,14,9,
  27.     14,11,2,12,4,7,13,1,5,0,15,10,3,9,8,6,
  28.     4,2,1,11,10,13,7,8,15,9,12,5,6,3,0,14,
  29.     11,8,12,7,1,14,2,13,6,15,0,9,10,4,5,3},
  30.     //S6
  31.     {12,1,10,15,9,2,6,8,0,13,3,4,14,7,5,11,
  32.     10,15,4,2,7,12,9,5,6,1,13,14,0,11,3,8,
  33.     9,14,15,5,2,8,12,3,7,0,4,10,1,13,11,6,
  34.     4,3,2,12,9,5,15,10,11,14,1,7,6,0,8,13,},
  35.     //S7
  36.     {4,11,2,14,15,0,8,13,3,12,9,7,5,10,6,1,
  37.     13,0,11,7,4,9,1,10,14,3,5,12,2,15,8,6,
  38.     1,4,11,13,12,3,7,14,10,15,6,8,0,5,9,2,
  39.     6,11,13,8,1,4,10,7,9,5,0,15,14,2,3,12},
  40.     //S8
  41.     {13,2,8,4,6,15,11,1,10,9,3,14,5,0,12,7,
  42.     1,15,13,8,10,3,7,4,12,5,6,11,0,14,9,2,
  43.     7,11,4,1,9,12,14,2,0,6,10,13,15,3,5,8,
  44.     2,1,14,7,4,10,8,13,15,12,9,0,3,5,6,11}
  45. };
  47. int permTable[] = {               //置换函数P
  48.     16, 7 , 20, 21, 29, 12, 28, 17,
  49.     1 , 15, 23, 26, 5 , 18, 31, 10,
  50.     2 , 8 , 24, 14, 32, 27, 3 , 9 ,
  51.     19, 13, 30, 6 , 22, 11, 4 , 25 };
  53. string substitution(string input)
  54. {
  55.     string res = ""; 
  56.     for(int i=0; i<8; i++)
  57.     {
  58.         string sInput = input.substr(6*i, 6) ;
  59.         int row = bitset<2>( sInput.substr(0,1) + sInput.substr(5,1) ).to_ulong() ;
  60.         int col = bitset<4>( sInput.substr(1,4) ).to_ulong() ;
  61.         res += bitset<4>(sBoxes[i][row][col]).to_string() ;
  63.     }
  64.     return res;
  65. }
  67. string permute(string input)
  68. {
  69.     string res = "";
  70.     for(int i=0; i<32; i++)
  71.     {
  72.         res += input[permTable[i]-1];
  73.     }
  74.     return res;
  75. }
  77. string XOR(string input1, string input2)
  78. {
  79.     string res = "";
  80.     for(int i=0; i<input1.length(); i++)
  81.     {
  82.         res += (input1[i] == input2[i]) ? "0" : "1";
  83.     }
  84.     return res;
  85. }
  87. int main()
  88. {
  89.     unsigned long long hexSBoxInput, hexInput;
  91.     cout << "\n  输入十六进制的S-box的48位输入（12bit）     : " ;
  92.     cin >> hex >> hexSBoxInput;
  94.     cout << "输入 十六进制的 64位（i-1）的输出(16-digits) : " ;
  95.     cin >> hex >> hexInput;
  97.     string sBoxinput = bitset<48>(hexSBoxInput).to_string();
  98.     cout << "\n S-Box Input       : " << sBoxinput << endl;
  100.     string input = bitset<64>(hexInput).to_string();
  101.     cout << " Round(i-1) output : " << input << endl;
  103.     string Li_1 = input.substr(0,32);
  104.     cout << "\n  Li_1              : " << Li_1 << endl;
  106.     string sBoxOutput = substitution(sBoxinput);
  107.     cout << "\n  S-Box output    = " << sBoxOutput << endl;
  109.     string P = permute(sBoxOutput);
  110.     cout << " P置换输出 = " << P << endl;
  112.     string Ri = XOR(P, Li_1); // P is permuted string and Li_1 is left half output of (i-1)th round
  113.     cout << "\n (Ri) = " << Ri << endl << endl;
  114. }