DES算法详解与源码

网上关于DES算法的讲述有很多,大致思路一致。但是很多细节的处理上没有交代清楚,源码质量也参差不齐,为此也花了很多时间研究了一下,现在把完整思路和源码整理如下。

 

1. DES算法简介:

  

DES算法为密码体制中的对称密码体制,又被称为美国数据加密标准,是1972年美国IBM公司研制的对称密码体制加密算法。 明文按64位进行分组,密钥长64位,密钥事实上是56位参与DES运算(第8、16、24、32、40、48、56、64位是校验位, 使得每个密钥都有奇数个1)分组后的明文组和56位的密钥按位替代或交换的方法形成密文组的加密方法。

 

2. DES步骤:

  1.首先输入一个8字节的密钥M8。将8字节的密钥转换成成64位的二进制位M64。其中第8,16,24,32,40,48,56,64位是校验位, 使得每个密钥都有奇数个1。所以密钥事实上是56位,储存时可以先储存64位,在下一步骤中的PC1置换中可以扣除。

    这里注意高低位问题:比如'L'的ASCII为75,转换为二进制为'01001011',在这里二进制从左到右为从高到低位,在密钥转换时,则从低位向高位储存

    如:'L'的ASCII为75,二进制为'01001011',转换至数组内为:char a = {0,1,0,0,1,0,1,1};

    网上各源代码对此处的处理有所不同,这种处理办法能够获得正确的解析结果。

  2. 将64位M64经过PC1置换转换为56位(扣除奇偶校验位)。

    对这56位密钥进行如下表的换位。

    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,

       表的意思是第57位移到第1位,第49位移到第2位,...... 以此类推。

  3. 变换后得到56位数据M56,将它分成两部分,C[0][28], D[0][28]。

  4. 计算16个子密钥。计算方法C[i][28] D[i][28]为对前一个C[i-1][28], D[i-1][28]做循环左移操作。16次的左移位数如下表:

      1,  1,  2,  2,  2,  2,  2,  2,  1,   2,  2,   2,   2,   2,   2,  1    (左移位数)

    在这里的16次循环中,每一次循环中,左移后,将C[i][28],D[i][28]重新凑成56位后,对这56位进行一次PC2置换,得到48位的子密钥sub[i][48]。16次循环后就有16个子密钥。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,

    表的意思是第14位移到第1位,第17位移到第2位,以此类推。在此过程中,发现第9,18,22,25, 35,38,43,54位丢弃。

*********************************************************************

  5.得到子密钥以后,就是对明文进行加密和解密操作了,加密和解密操作具有对称性,因此介绍加密操作:

  6.输入明文,要求达到8字节,若没有达到8字节则进行扩展,扩展至8字节(补0或补'\0')。

  7.将8字节的明文转换成成64位的二进制位。

  8. 对64bit的明文输入进行换位变换。换位表如下:

    58, 50, 12, 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

   表的意思就是第一次变换时,第58位移到第1位,第50位移到第2位,...... 依此类推。得到64位数据data[64]。

  9.接下来将进入循环处理阶段,总共进行16次循环,每次循环遵循如下操作:

  for(int i=0;i<16;i++)

  {

    10.将64位数据data[64]前后分成两块L[i][32], R[i][32]。

    11.对R[i][32]进行扩展变换成48位数,方法如下, 记为E(R[i][32])

      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,

    12.将E(R[i][32])与sub[i][48](由密钥处理获得的子密钥,见步骤4)作异或运算,得到48位数B[48]。

    13.将48位数顺序分成8份,6位一份,B[48]->B[8][6]。

    14.使用8个S盒对分成8份的B[8][6]进行替换,步骤如下:

        a) 取出B[j][6]的第1位和第6位连成一个2位二进制数m,范围(0-3),转换成十进制数后作为第j个S盒的行标

          如:假设B[0]为'010011',则取出第一位0和第六位1,组成01,01转换为十进制为1,因此m代表了第0个S盒中的第2行(数组下标从0开始)

        b) 取出B[j][6]的第2、3、4、5位连成一个4位二进制数n,范围(0-15),转换成十进制数后作为第j个S盒的列标

          如:假设B[0]为'010011',则取出'1001',转换成十进制为9,因此n代表了第0个S盒中的第10列(数组下标从0开始)

        c) 有了行标和列表则进入第j个S盒中查找对应的数字

          接上例,第0个S盒中第二行第十列获得数字:6

        d) 将获得的数字转换成二进制,替换48位的B中的前4位。

          接上例,6转换为二进制为0110,这样B的前四位就为0110,循环后从第5位继续。

        e) 经过8次循环,每次替换4位,最终B的前32位则为我们新获得的加密数据B[32]

    15.对B[32]进行如下变换得到B'[32]:

         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,

    16.B'[32]与L[i][32]作异或运算,把结果覆盖data[64]的左半部分,若i!=15,则交换data[64]的左右部分(即异或运算结果作为下一次循环的右边,原来的右边作为下      一次运算结果的左边)。若i=15,则不进行交换。

  }

  17. 结束循环后,将16次处理过后的data[64]进行一次逆初始置换,这是对第8步的逆变换

     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

  这样就得到了最终的密文了QUQ~

 

  解密过程同加密过程完全相同,不同的是对子密钥的异或顺序不同,加密过程 i 从0-15与子密钥配对,解密过程中则从15-0与子密钥进行异或。

  

源代码见链接:

http://pan.baidu.com/s/1c0D1RAo

 

5.1.cpp为字符char与8位二进制char之间的转换

 

5.2&5.3.cpp为DES的加密程序,其中输出中间过程如下:

a. 密钥对应密钥的64位比特(8X8形式输出,使用5.1中的函数进行转换);

b. 经过置换选择1后,分别输出的C0(28位)和D0(28位)对应的比特;

c. 输出第5次迭代时,产生的子密钥K5所对应的48位比特;

d. 原文对应的64位2进制比特序列;

e. 64比特序列经过初始置换IP后,得到新的64位比特;

f. 乘积变换需要进行16次迭代,输出经过5次迭代后,得到的64位比特;

g. 输出完成乘积变换(16次迭代)后得到的64位比特;

h. 输出加密后的密文的64位比特;

 

5.4.cpp是针对一段指定的密文进行解密的程序。

 

5.5.cpp,将给定的文本文件“DES.txt”进行加密,密钥依旧是“KEYMIYAO”,把密文保存在另一个文本文件(“DES_Encrypt.txt”)中,然后对这个文件进行解密,得到的原文保存到文本文件(“DES_Final.txt”)中。其中包括了字符不足时的扩展。将此文件加入vs中后要自己创建以上提到的3个txt资源文件。

 

添加txt文本的方法如下:

  VC++6.0中,点击菜单栏 文件 - 新建

    在弹出的对话框中选择文本文件,然后键入指定的文件名。然后将其拖入Resource File文件夹。

    

    

    也可以右键Resource File文件夹,直接添加。

  VS2012中可以直接右键Resource File文件夹。新建项-添加,在实用工具内选择txt文件即可。

 

 

DESc版本.cpp,是DES的c语言实现,效率高,写的好。别看我我是从信息安全专业那里扒来的捂脸QUQ。

 

在低版本的VS中可能会出现for循环中不允许声明int的问题,自行修改int变量提至函数开头。

有空再来补充3DES算法了··

 

转载于:https://www.cnblogs.com/xiaoxiaff/p/3960302.html

C代码 //初始置换表IP int IP_Table[64] = { 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, 56,48,40,32,24,16,8,0, 58,50,42,34,26,18,10,2, 60,52,44,36,28,20,12,4, 62,54,46,38,30,22,14,6}; //逆初始置换表IP^-1 int IP_1_Table[64] = {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, 32,0,40,8,48,16,56,24}; //扩充置换表E int E_Table[48] = {31, 0, 1, 2, 3, 4, 3, 4, 5, 6, 7, 8, 7, 8,9,10,11,12, 11,12,13,14,15,16, 15,16,17,18,19,20, 19,20,21,22,23,24, 23,24,25,26,27,28, 27,28,29,30,31, 0}; //置换函数P int P_Table[32] = {15,6,19,20,28,11,27,16, 0,14,22,25,4,17,30,9, 1,7,23,13,31,26,2,8, 18,12,29,5,21,10,3,24}; //S盒 int S[8][4][16] =//S1 {{{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}}, //S2 {{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}}, //S3 {{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}}, //S4 {{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}}, //S5 {{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}}, //S6 {{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}}, //S7 {{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}}, //S8 {{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}}}; //置换选择1 int PC_1[56] = {56,48,40,32,24,16,8, 0,57,49,41,33,25,17, 9,1,58,50,42,34,26, 18,10,2,59,51,43,35, 62,54,46,38,30,22,14, 6,61,53,45,37,29,21, 13,5,60,52,44,36,28, 20,12,4,27,19,11,3}; //置换选择2 int PC_2[48] = {13,16,10,23,0,4,2,27, 14,5,20,9,22,18,11,3, 25,7,15,6,26,19,12,1, 40,51,30,36,46,54,29,39, 50,44,32,46,43,48,38,55, 33,52,45,41,49,35,28,31}; //对左移次数的规定 int MOVE_TIMES[16] = {1,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1}; //初始置换表IP int IP_Table[64] = { 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, 56,48,40,32,24,16,8,0, 58,50,42,34,26,18,10,2, 60,52,44,36,28,20,12,4, 62,54,46,38,30,22,14,6}; //逆初始置换表IP^-1 int IP_1_Table[64] = {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, 32,0,40,8,48,16,56,24}; //扩充置换表E int E_Table[48] = {31, 0, 1, 2, 3, 4, 3, 4, 5, 6, 7, 8, 7, 8,9,10,11,12, 11,12,13,14,15,16, 15,16,17,18,19,20, 19,20,21,22,23,24, 23,24,25,26,27,28, 27,28,29,30,31, 0}; //置换函数P int P_Table[32] = {15,6,19,20,28,11,27,16, 0,14,22,25,4,17,30,9, 1,7,23,13,31,26,2,8, 18,12,29,5,21,10,3,24}; //S盒 int S[8][4][16] =//S1 {{{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}}, //S2 {{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}}, //S3 {{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}}, //S4 {{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}}, //S5 {{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}}, //S6 {{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}}, //S7 {{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}}, //S8 {{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}}}; //置换选择1 int PC_1[56] = {56,48,40,32,24,16,8, 0,57,49,41,33,25,17, 9,1,58,50,42,34,26, 18,10,2,59,51,43,35, 62,54,46,38,30,22,14, 6,61,53,45,37,29,21, 13,5,60,52,44,36,28, 20,12,4,27,19,11,3}; //置换选择2 int PC_2[48] = {13,16,10,23,0,4,2,27, 14,5,20,9,22,18,11,3, 25,7,15,6,26,19,12,1, 40,51,30,36,46,54,29,39, 50,44,32,46,43,48,38,55, 33,52,45,41,49,35,28,31}; //对左移次数的规定 int MOVE_TIMES[16] = {1,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1}; 二、模块化。 对面向过程的程序,模块化是否清晰是至关重要的。 下面是函数的声明: C代码 int ByteToBit(ElemType ch,ElemType bit[8]); int BitToByte(ElemType bit[8],ElemType *ch); int Char8ToBit64(ElemType ch[8],ElemType bit[64]); int Bit64ToChar8(ElemType bit[64],ElemType ch[8]); int DES_MakeSubKeys(ElemType key[64],ElemType subKeys[16][48]); int DES_PC1_Transform(ElemType key[64], ElemType tempbts[56]); int DES_PC2_Transform(ElemType key[56], ElemType tempbts[48]); int DES_ROL(ElemType data[56], int time); int DES_IP_Transform(ElemType data[64]); int DES_IP_1_Transform(ElemType data[64]); int DES_E_Transform(ElemType data[48]); int DES_P_Transform(ElemType data[32]); int DES_SBOX(ElemType data[48]); int DES_XOR(ElemType R[48], ElemType L[48],int count); int DES_Swap(ElemType left[32],ElemType right[32]); int DES_EncryptBlock(ElemType plainBlock[8], ElemType subKeys[16][48], ElemType cipherBlock[8]); int DES_DecryptBlock(ElemType cipherBlock[8], ElemType subKeys[16][48], ElemType plainBlock[8]); int DES_Encrypt(char *plainFile, char *keyStr,char *cipherFile); int DES_Decrypt(char *cipherFile, char *keyStr,char *plainFile); int ByteToBit(ElemType ch,ElemType bit[8]); int BitToByte(ElemType bit[8],ElemType *ch); int Char8ToBit64(ElemType ch[8],ElemType bit[64]); int Bit64ToChar8(ElemType bit[64],ElemType ch[8]); int DES_MakeSubKeys(ElemType key[64],ElemType subKeys[16][48]); int DES_PC1_Transform(ElemType key[64], ElemType tempbts[56]); int DES_PC2_Transform(ElemType key[56], ElemType tempbts[48]); int DES_ROL(ElemType data[56], int time); int DES_IP_Transform(ElemType data[64]); int DES_IP_1_Transform(ElemType data[64]); int DES_E_Transform(ElemType data[48]); int DES_P_Transform(ElemType data[32]); int DES_SBOX(ElemType data[48]); int DES_XOR(ElemType R[48], ElemType L[48],int count); int DES_Swap(ElemType left[32],ElemType right[32]); int DES_EncryptBlock(ElemType plainBlock[8], ElemType subKeys[16][48], ElemType cipherBlock[8]); int DES_DecryptBlock(ElemType cipherBlock[8], ElemType subKeys[16][48], ElemType plainBlock[8]); int DES_Encrypt(char *plainFile, char *keyStr,char *cipherFile); int DES_Decrypt(char *cipherFile, char *keyStr,char *plainFile); 其实,模块化与速度也是一对矛盾,因为了解函数运行机制的人就知道,我们的计算机在运行某个函数时,是要用栈来保存入口状态的,在运行结束后又要恢复现场,这些操作势必会影像系统性能,但我们不能将所有代码写在Main函数里,虽然那样做我们的加密算法效率又会大增,但是那种代码未免太过于丑陋不堪。因此,为了帅,还是牺牲一下性能吧。 三、实现。 代码里能用移位操作都尽量用了移位操作,能用逻辑运算符的都用了逻辑运算符。 详细的行注相信你可以看懂吧。有问题可以M我。 C代码 //字节转换成二进制 int ByteToBit(ElemType ch, ElemType bit[8]){ int cnt; for(cnt = 0;cnt >cnt)&1; } return 0; } //二进制转换成字节 int BitToByte(ElemType bit[8],ElemType *ch){ int cnt; for(cnt = 0;cnt < 8; cnt++){ *ch |= *(bit + cnt)<<cnt; } return 0; } //将长度为8的字符串转为二进制位串 int Char8ToBit64(ElemType ch[8],ElemType bit[64]){ int cnt; for(cnt = 0; cnt < 8; cnt++){ ByteToBit(*(ch+cnt),bit+(cnt<<3)); } return 0; } //将二进制位串转为长度为8的字符串 int Bit64ToChar8(ElemType bit[64],ElemType ch[8]){ int cnt; memset(ch,0,8); for(cnt = 0; cnt < 8; cnt++){ BitToByte(bit+(cnt<<3),ch+cnt); } return 0; } //生成子密钥 int DES_MakeSubKeys(ElemType key[64],ElemType subKeys[16][48]){ ElemType temp[56]; int cnt; DES_PC1_Transform(key,temp);//PC1置换 for(cnt = 0; cnt < 16; cnt++){//16轮跌代,产生16个子密钥 DES_ROL(temp,MOVE_TIMES[cnt]);//循环左移 DES_PC2_Transform(temp,subKeys[cnt]);//PC2置换,产生子密钥 } return 0; } //密钥置换1 int DES_PC1_Transform(ElemType key[64], ElemType tempbts[56]){ int cnt; for(cnt = 0; cnt < 56; cnt++){ tempbts[cnt] = key[PC_1[cnt]]; } return 0; } //密钥置换2 int DES_PC2_Transform(ElemType key[56], ElemType tempbts[48]){ int cnt; for(cnt = 0; cnt < 48; cnt++){ tempbts[cnt] = key[PC_2[cnt]]; } return 0; } //循环左移 int DES_ROL(ElemType data[56], int time){ ElemType temp[56]; //保存将要循环移动到右边的位 memcpy(temp,data,time); memcpy(temp+time,data+28,time); //前28位移动 memcpy(data,data+time,28-time); memcpy(data+28-time,temp,time); //后28位移动 memcpy(data+28,data+28+time,28-time); memcpy(data+56-time,temp+time,time); return 0; } //IP置换 int DES_IP_Transform(ElemType data[64]){ int cnt; ElemType temp[64]; for(cnt = 0; cnt < 64; cnt++){ temp[cnt] = data[IP_Table[cnt]]; } memcpy(data,temp,64); return 0; } //IP逆置换 int DES_IP_1_Transform(ElemType data[64]){ int cnt; ElemType temp[64]; for(cnt = 0; cnt < 64; cnt++){ temp[cnt] = data[IP_1_Table[cnt]]; } memcpy(data,temp,64); return 0; } //扩展置换 int DES_E_Transform(ElemType data[48]){ int cnt; ElemType temp[48]; for(cnt = 0; cnt < 48; cnt++){ temp[cnt] = data[E_Table[cnt]]; } memcpy(data,temp,48); return 0; } //P置换 int DES_P_Transform(ElemType data[32]){ int cnt; ElemType temp[32]; for(cnt = 0; cnt < 32; cnt++){ temp[cnt] = data[P_Table[cnt]]; } memcpy(data,temp,32); return 0; } //异或 int DES_XOR(ElemType R[48], ElemType L[48] ,int count){ int cnt; for(cnt = 0; cnt < count; cnt++){ R[cnt] ^= L[cnt]; } return 0; } //S盒置换 int DES_SBOX(ElemType data[48]){ int cnt; int line,row,output; int cur1,cur2; for(cnt = 0; cnt < 8; cnt++){ cur1 = cnt*6; cur2 = cnt<<2; //计算在S盒中的行与列 line = (data[cur1]<<1) + data[cur1+5]; row = (data[cur1+1]<<3) + (data[cur1+2]<<2) + (data[cur1+3]<>3; data[cur2+1] = (output&0X04)>>2; data[cur2+2] = (output&0X02)>>1; data[cur2+3] = output&0x01; } return 0; } //交换 int DES_Swap(ElemType left[32], ElemType right[32]){ ElemType temp[32]; memcpy(temp,left,32); memcpy(left,right,32); memcpy(right,temp,32); return 0; } //加密单个分组 int DES_EncryptBlock(ElemType plainBlock[8], ElemType subKeys[16][48], ElemType cipherBlock[8]){ ElemType plainBits[64]; ElemType copyRight[48]; int cnt; Char8ToBit64(plainBlock,plainBits); //初始置换(IP置换) DES_IP_Transform(plainBits); //16轮迭代 for(cnt = 0; cnt = 0; cnt--){ memcpy(copyRight,cipherBits+32,32); //将右半部分进行扩展置换,从32位扩展到48位 DES_E_Transform(copyRight); //将右半部分与子密钥进行异或操作 DES_XOR(copyRight,subKeys[cnt],48); //异或结果进入S盒,输出32位结果 DES_SBOX(copyRight); //P置换 DES_P_Transform(copyRight); //将明文左半部分与右半部分进行异或 DES_XOR(cipherBits,copyRight,32); if(cnt != 0){ //最终完成左右部的交换 DES_Swap(cipherBits,cipherBits+32); } } //逆初始置换(IP^1置换) DES_IP_1_Transform(cipherBits); Bit64ToChar8(cipherBits,plainBlock); return 0; } //加密文件 int DES_Encrypt(char *plainFile, char *keyStr,char *cipherFile){ FILE *plain,*cipher; int count; ElemType plainBlock[8],cipherBlock[8],keyBlock[8]; ElemType bKey[64]; ElemType subKeys[16][48]; if((plain = fopen(plainFile,"rb")) == NULL){ return PLAIN_FILE_OPEN_ERROR; } if((cipher = fopen(cipherFile,"wb")) == NULL){ return CIPHER_FILE_OPEN_ERROR; } //设置密钥 memcpy(keyBlock,keyStr,8); //将密钥转换为二进制流 Char8ToBit64(keyBlock,bKey); //生成子密钥 DES_MakeSubKeys(bKey,subKeys); while(!feof(plain)){ //每次读8个字节,并返回成功读取的字节数 if((count = fread(plainBlock,sizeof(char),8,plain)) == 8){ DES_EncryptBlock(plainBlock,subKeys,cipherBlock); fwrite(cipherBlock,sizeof(char),8,cipher); } } if(count){ //填充 memset(plainBlock + count,'\0',7 - count); //最后一个字符保存包括最后一个字符在内的所填充的字符数量 plainBlock[7] = 8 - count; DES_EncryptBlock(plainBlock,subKeys,cipherBlock); fwrite(cipherBlock,sizeof(char),8,cipher); } fclose(plain); fclose(cipher); return OK; } //解密文件 int DES_Decrypt(char *cipherFile, char *keyStr,char *plainFile){ FILE *plain, *cipher; int count,times = 0; long fileLen; ElemType plainBlock[8],cipherBlock[8],keyBlock[8]; ElemType bKey[64]; ElemType subKeys[16][48]; if((cipher = fopen(cipherFile,"rb")) == NULL){ return CIPHER_FILE_OPEN_ERROR; } if((plain = fopen(plainFile,"wb")) == NULL){ return PLAIN_FILE_OPEN_ERROR; } //设置密钥 memcpy(keyBlock,keyStr,8); //将密钥转换为二进制流 Char8ToBit64(keyBlock,bKey); //生成子密钥 DES_MakeSubKeys(bKey,subKeys); //取文件长度 fseek(cipher,0,SEEK_END); //将文件指针置尾 fileLen = ftell(cipher); //取文件指针当前位置 rewind(cipher); //将文件指针重指向文件头 while(1){ //密文的字节数一定是8的整数倍 fread(cipherBlock,sizeof(char),8,cipher); DES_DecryptBlock(cipherBlock,subKeys,plainBlock); times += 8; if(times < fileLen){ fwrite(plainBlock,sizeof(char),8,plain); } else{ break; } } //判断末尾是否被填充 if(plainBlock[7] < 8){ for(count = 8 - plainBlock[7]; count < 7; count++){ if(plainBlock[count] != '\0'){ break; } } } if(count == 7){//有填充 fwrite(plainBlock,sizeof(char),8 - plainBlock[7],plain); } else{//无填充 fwrite(plainBlock,sizeof(char),8,plain); } fclose(plain); fclose(cipher); return OK; } //字节转换成二进制 int ByteToBit(ElemType ch, ElemType bit[8]){ int cnt; for(cnt = 0;cnt >cnt)&1; } return 0; } //二进制转换成字节 int BitToByte(ElemType bit[8],ElemType *ch){ int cnt; for(cnt = 0;cnt < 8; cnt++){ *ch |= *(bit + cnt)<<cnt; } return 0; } //将长度为8的字符串转为二进制位串 int Char8ToBit64(ElemType ch[8],ElemType bit[64]){ int cnt; for(cnt = 0; cnt < 8; cnt++){ ByteToBit(*(ch+cnt),bit+(cnt<<3)); } return 0; } //将二进制位串转为长度为8的字符串 int Bit64ToChar8(ElemType bit[64],ElemType ch[8]){ int cnt; memset(ch,0,8); for(cnt = 0; cnt < 8; cnt++){ BitToByte(bit+(cnt<<3),ch+cnt); } return 0; } //生成子密钥 int DES_MakeSubKeys(ElemType key[64],ElemType subKeys[16][48]){ ElemType temp[56]; int cnt; DES_PC1_Transform(key,temp);//PC1置换 for(cnt = 0; cnt < 16; cnt++){//16轮跌代,产生16个子密钥 DES_ROL(temp,MOVE_TIMES[cnt]);//循环左移 DES_PC2_Transform(temp,subKeys[cnt]);//PC2置换,产生子密钥 } return 0; } //密钥置换1 int DES_PC1_Transform(ElemType key[64], ElemType tempbts[56]){ int cnt; for(cnt = 0; cnt < 56; cnt++){ tempbts[cnt] = key[PC_1[cnt]]; } return 0; } //密钥置换2 int DES_PC2_Transform(ElemType key[56], ElemType tempbts[48]){ int cnt; for(cnt = 0; cnt < 48; cnt++){ tempbts[cnt] = key[PC_2[cnt]]; } return 0; } //循环左移 int DES_ROL(ElemType data[56], int time){ ElemType temp[56]; //保存将要循环移动到右边的位 memcpy(temp,data,time); memcpy(temp+time,data+28,time); //前28位移动 memcpy(data,data+time,28-time); memcpy(data+28-time,temp,time); //后28位移动 memcpy(data+28,data+28+time,28-time); memcpy(data+56-time,temp+time,time); return 0; } //IP置换 int DES_IP_Transform(ElemType data[64]){ int cnt; ElemType temp[64]; for(cnt = 0; cnt < 64; cnt++){ temp[cnt] = data[IP_Table[cnt]]; } memcpy(data,temp,64); return 0; } //IP逆置换 int DES_IP_1_Transform(ElemType data[64]){ int cnt; ElemType temp[64]; for(cnt = 0; cnt < 64; cnt++){ temp[cnt] = data[IP_1_Table[cnt]]; } memcpy(data,temp,64); return 0; } //扩展置换 int DES_E_Transform(ElemType data[48]){ int cnt; ElemType temp[48]; for(cnt = 0; cnt < 48; cnt++){ temp[cnt] = data[E_Table[cnt]]; } memcpy(data,temp,48); return 0; } //P置换 int DES_P_Transform(ElemType data[32]){ int cnt; ElemType temp[32]; for(cnt = 0; cnt < 32; cnt++){ temp[cnt] = data[P_Table[cnt]]; } memcpy(data,temp,32); return 0; } //异或 int DES_XOR(ElemType R[48], ElemType L[48] ,int count){ int cnt; for(cnt = 0; cnt < count; cnt++){ R[cnt] ^= L[cnt]; } return 0; } //S盒置换 int DES_SBOX(ElemType data[48]){ int cnt; int line,row,output; int cur1,cur2; for(cnt = 0; cnt < 8; cnt++){ cur1 = cnt*6; cur2 = cnt<<2; //计算在S盒中的行与列 line = (data[cur1]<<1) + data[cur1+5]; row = (data[cur1+1]<<3) + (data[cur1+2]<<2) + (data[cur1+3]<>3; data[cur2+1] = (output&0X04)>>2; data[cur2+2] = (output&0X02)>>1; data[cur2+3] = output&0x01; } return 0; } //交换 int DES_Swap(ElemType left[32], ElemType right[32]){ ElemType temp[32]; memcpy(temp,left,32); memcpy(left,right,32); memcpy(right,temp,32); return 0; } //加密单个分组 int DES_EncryptBlock(ElemType plainBlock[8], ElemType subKeys[16][48], ElemType cipherBlock[8]){ ElemType plainBits[64]; ElemType copyRight[48]; int cnt; Char8ToBit64(plainBlock,plainBits); //初始置换(IP置换) DES_IP_Transform(plainBits); //16轮迭代 for(cnt = 0; cnt = 0; cnt--){ memcpy(copyRight,cipherBits+32,32); //将右半部分进行扩展置换,从32位扩展到48位 DES_E_Transform(copyRight); //将右半部分与子密钥进行异或操作 DES_XOR(copyRight,subKeys[cnt],48); //异或结果进入S盒,输出32位结果 DES_SBOX(copyRight); //P置换 DES_P_Transform(copyRight); //将明文左半部分与右半部分进行异或 DES_XOR(cipherBits,copyRight,32); if(cnt != 0){ //最终完成左右部的交换 DES_Swap(cipherBits,cipherBits+32); } } //逆初始置换(IP^1置换) DES_IP_1_Transform(cipherBits); Bit64ToChar8(cipherBits,plainBlock); return 0; } //加密文件 int DES_Encrypt(char *plainFile, char *keyStr,char *cipherFile){ FILE *plain,*cipher; int count; ElemType plainBlock[8],cipherBlock[8],keyBlock[8]; ElemType bKey[64]; ElemType subKeys[16][48]; if((plain = fopen(plainFile,"rb")) == NULL){ return PLAIN_FILE_OPEN_ERROR; } if((cipher = fopen(cipherFile,"wb")) == NULL){ return CIPHER_FILE_OPEN_ERROR; } //设置密钥 memcpy(keyBlock,keyStr,8); //将密钥转换为二进制流 Char8ToBit64(keyBlock,bKey); //生成子密钥 DES_MakeSubKeys(bKey,subKeys); while(!feof(plain)){ //每次读8个字节,并返回成功读取的字节数 if((count = fread(plainBlock,sizeof(char),8,plain)) == 8){ DES_EncryptBlock(plainBlock,subKeys,cipherBlock); fwrite(cipherBlock,sizeof(char),8,cipher); } } if(count){ //填充 memset(plainBlock + count,'\0',7 - count); //最后一个字符保存包括最后一个字符在内的所填充的字符数量 plainBlock[7] = 8 - count; DES_EncryptBlock(plainBlock,subKeys,cipherBlock); fwrite(cipherBlock,sizeof(char),8,cipher); } fclose(plain); fclose(cipher); return OK; } //解密文件 int DES_Decrypt(char *cipherFile, char *keyStr,char *plainFile){ FILE *plain, *cipher; int count,times = 0; long fileLen; ElemType plainBlock[8],cipherBlock[8],keyBlock[8]; ElemType bKey[64]; ElemType subKeys[16][48]; if((cipher = fopen(cipherFile,"rb")) == NULL){ return CIPHER_FILE_OPEN_ERROR; } if((plain = fopen(plainFile,"wb")) == NULL){ return PLAIN_FILE_OPEN_ERROR; } //设置密钥 memcpy(keyBlock,keyStr,8); //将密钥转换为二进制流 Char8ToBit64(keyBlock,bKey); //生成子密钥 DES_MakeSubKeys(bKey,subKeys); //取文件长度 fseek(cipher,0,SEEK_END); //将文件指针置尾 fileLen = ftell(cipher); //取文件指针当前位置 rewind(cipher); //将文件指针重指向文件头 while(1){ //密文的字节数一定是8的整数倍 fread(cipherBlock,sizeof(char),8,cipher); DES_DecryptBlock(cipherBlock,subKeys,plainBlock); times += 8; if(times < fileLen){ fwrite(plainBlock,sizeof(char),8,plain); } else{ break; } } //判断末尾是否被填充 if(plainBlock[7] < 8){ for(count = 8 - plainBlock[7]; count < 7; count++){ if(plainBlock[count] != '\0'){ break; } } } if(count == 7){//有填充 fwrite(plainBlock,sizeof(char),8 - plainBlock[7],plain); } else{//无填充 fwrite(plainBlock,sizeof(char),8,plain); } fclose(plain); fclose(cipher); return OK; } 最后,写一个简单的main函数来检验它: C代码 int main() { clock_t a,b; a = clock(); DES_Encrypt("1.txt","key.txt","2.txt"); b = clock(); printf("加密消耗%d毫秒\n",b-a); system("pause"); a = clock(); DES_Decrypt("2.txt","key.txt","3.txt"); b = clock(); printf("解密消耗%d毫秒\n",b-a); getchar(); return 0; }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值