密码学优秀课程设计优质报告.docx

密码学
课程设计汇报

班级：信安09-2班
姓名：李明月
学号：08093755

目录
1 古典密码算法—凯撒密码4
1.1凯撒密码概述4
1.2算法原理及设计思想4
1.3主要算法分析4
1.4程序运行结果4
1.5密码安全性分析5
2 序列密码—RC45
2.1 RC4算法概述5
2.2 算法原理及设计思想5
2.3 程序主要算法分析6
2.4 程序运行结果7
2.5 算法分析7
3 分组密码算法8
3.1 DES加解密算法的实现8
3.1.1 DES算法概述8
3.1.2 算法原理及设计思想8
3.1.3 程序主要算法分析11
3.1.4 程序运行结果13
3.1.5 安全性分析14
3.2 AES加解密算法的实现14
3.2.1 AES算法概述15
3.2.2 算法原理及设计思想15
3.2.3 程序主要算法分析17
3.2.4 程序运行结果22
3.2.5 安全性分析22
4 HASH函数—MD5算法23
4.1 算法概述23
4.2 算法原理及设计思想23
4.3 程序主要算法分析26
4.4 程序运行结果28
4.5 安全性分析28

5 公钥密码算法---RSA29
5.1 算法概述29
5.2算法原理及设计思想29
5.2.1 算法描述—密钥生成29
5.2.2 算法描述—加密、解密29
5.2.3 原理29
5.3程序主要算法分析30
5.4程序运行结果31
5.5安全性分析31
6 设计体会32

一、 古典密码算法---凯撒密码
1.1凯撒密码概述
凯撒密码作为一个最为古老对称加密体制，在古罗马时候全部已经很流行，她基础思想是：经过把字母移动一定位数来实现加密和解密。比如，假如密钥是把明文字母位数向后移动三位，那么明文字母B就变成了密文E，依次类推，X将变成A，Y变成B，Z变成C，由此可见，位数就是凯撒密码加密和解密密钥。　
它是一个代换密码。听说凯撒是率先使用加密函古代将领之一，所以这种加密方法被称为凯撒密码。
　 在密码学中凯撒密码（或称凯撒加密、凯撒变换、变换加密）是一个最简单且最广为人知加密技术。它是一个替换加密技术，明文中全部字母全部在字母表上向后（或向前）根据一个固定数目进行偏移后被替换成密文。比如，当偏移量是3时候，全部字母A将被替换成D，B变成E，以这类推。这个加密方法是以凯撒名字命名，当年凯撒曾用此方法和其将军们进行联络。凯撒密码通常被作为其它更复杂加密方法中一个步骤，比如维吉尼亚密码。凯撒密码还在现代ROT13系统中被应用。不过和全部利用字母表进行替换加密技术一样，凯撒密码很轻易被破解，而且在实际应用中也无法确保通信安全。
1.2算法原理及设计思想
它是一个替换密码，经过将字母按次序推后起3位起到加密作用，如将字母A换作字母D，将字母B换作字母E。因听说恺撒是率先使用加密函古代将领之一，所以这种加密方法被称为恺撒密码。这是一个简单加密方法，这种密码密度是很低，只需简单地统计字频就能够破译。现今又叫“移位密码”，只不过移动为数不一定是3位而已。
密码术能够大致别分为两种，即易位和替换，当然也有二者结合更复杂方法。在易位中字母不变，位置改变；替换中字母改变，位置不变。
凯撒密码表就是用D代a，用E代b，……，用z代w，（注意！）用A带x，用B代y，C代z。这些替换规则也可用一张表格来表示（所以叫“密表”）。
1.3关键算法分析
//密码表定义
char a[26];
for(int i=0;i<26;i++) a[i]=char(65+i);
//明文转化为凯撒密码
for(int h=0;h<strlen(s);h++)
g[h]=int(s[h]);
l=((g[h]-65)+key)%26;
//凯撒密码转化为明文
for(int v=0;v<strlen(q);v++)
e[v]=int(q[v]);
b=((e[v]-65)-key1+26)%26;
1.4程序运行结果

1.5密码安全性分析
凯撒密码是没有密钥，即使没有密钥也能将它破解出来，因为凯撒移位密码只有25种密钥，最多就是将这25种可能性挨个检测一下能够了，这就是我们所说暴力破解法。也可在用软件破解，不过我提倡用人工。
推理方法：
1、对于有空格凯撒移位，单字母A和I是突破口，这无异相当于告诉了移动位数，这么很轻易就被破解了。所以，假如我们要用凯撒密码话一定要去掉空格加大破解难。
2、差数法
　　有空格时，而又没有单字母A和I时，这种方法很，假如我们令A=1,B=2,C=3......就是每个字母是字母第多个，经过移位后单词，每两相邻字母之间差值不变。如the差值为12,3（在这里我是用后面一个字母减前面一个字母，当然你也能够用后面一个字母减前面一个字母），移动后两个相邻字母差值也将会是1，2，3。对于没有空格恺撒破解起来就比有空格难部分，对于没有空格我们还要对密文进行分析，找出反复出现字母串，然后对字母串进行猜测，例，假如有3个字母串，出现次数比较高，我们就能够假设它为the因为3个字母串出现次最多就是the，当然这不是一成不变，这时应该就被破解了。
二、 序列密码—RC4
2.1 RC4算法概述
RC4加密算法是大名鼎鼎RSA三人组中头号人物Ron Rivest在1987年设计密钥长度可变流加密算法簇。之所以称其为簇，是因为其关键部分S-box长度可为任意，但通常为256字节。该算法速度能够达成DES加密10倍左右。RC4算法是一个在电子信息领域加密技术手段，用于无线通信网络，是一个电子密码，只有经过授权（缴纳对应费用）用户才能享受该服务。RC4算法原理很简单，包含初始化算法和伪随机子密码生成算法两大部分。
2.2算法原理及设计思想
RC4流密码是一个可变密钥长度、面向字节操作流密码。以随机置换为基础。广泛用于SSL/TLS标准当中。
RC4算法能够分为两个部分，第一是依据种子密钥，利用密钥调度算法对数据表S进行重新排列，第二部分是利用伪随机数生成算法，从已重新排列数据表S中取出一个字节。每取出一个字节，数据表S将发生改变。
RC4描述起来也很简单：用从1到256个字节(8-2048比特)可变长度密钥初始化一个256字节状态向量S，S元素标识为S[0],S[1],…，S[255]，从始至终置换后S包含从0-255全部8比特数。对于加密和解密，字节K由S中255个元素根据一定方法选出一个元素生成。每生成一个K值，S中元素个体就被重新置换一次。
2.3 程序关键算法分析
初始化S
开始时，S中值被置为按升序0-255，即S[0]=0，S[1]=1,…,S[255]=255。同时建立一个临时变量T。假如密钥K长度为256字节，则将K赋给T。不然若密钥长度为keylen字节，则将K值赋给T前keylen个元素，并循环反复用K值赋给T剩下元素，直到T全部元素被赋值。
然后用T产生S初始置换，从S[0]到S[255]，对每个S[i]，依据T[i]确定方案，将S[i]置换为S中另一字节。因为对S操作仅仅为交换，所以唯一改变就是置换。S仍然包含全部值0-255元素。
void rc4_setup(struct rc4_state *s, unsigned char *key, int length)
{
int i, j, k, *m, a;
s->x = 0;
s->y = 0;
m = s->m;
for(i = 0; i < 256; i++)
{
m[i] = i;
}
j = k = 0;
for(i = 0; i < 256; i++)
{
a = m[i];
j = (unsigned char)(j + a + key[k]);
m[i] = m[j]; m[j] = a;
if(++k >= length) k = 0;
}
}
密钥流生成
向量S一旦初始化完成，输入密钥就不再被使用。密钥流生成是从S[0]到S[255]，对每个S[i]，依据目前S值，将S[i]和S中另一字节置换。当S[255]完成置换后，操作继续反复从S[0]开始。
加密中，将k值和下一明文字节异或；解密中，将k值和下一密文字节异或。
void rc4_crypt(struct rc4_state *s, unsigned char *data, int length)
{
int i, x, y, *m, a, b;
x = s->x;
y = s->y;
m = s->m;
for(i = 0; i < length; i++)
{
x = (unsigned char)(x + 1); a = m[x];
y = (unsigned char)(y + a);
m[x] = b = m[y];
m[y] = a;
data[i] ^= m[(unsigned char)( a + b )];
}
s->x = x;
s->y = y;
}
2.4程序运行结果

2.5算法分析
RC4算法优点是：算法简单、高效，尤其适合软件实现，RC4是现在应用最广商密级序列密码，现在被用于SSL/TLS标准中。
因为RC4算法加密是采取xor，所以，一旦子密钥序列出现了反复，密文就有可能被破解。那么，RC4算法生成子密钥序列是否会出现反复呢？经过我测试，存在部分弱密钥，使得子密钥序列在不到100万字节内就发生了完全反复，假如是部分反复，则可能在不到10万字节内就能发生反复，所以，推荐在使用RC4算法时，必需对加密密钥进行测试，判定其是否为弱密钥。
但在就有以色列科学家指出RC4加密算法存在着漏洞，这可能对无线通信网络安全组成威胁。
　以色列魏茨曼研究所和美国思科企业研究者发觉，在使用“有线等效保密规则”（WEP）无线网络中，在特定情况下，大家能够逆转RC4算法加密过程，获取密钥，从而将己加密信息解密。实现这一过程并不复杂，只需要使用一台个人电脑对加密数据进行分析，经过多个小时时间就能够破译出信息全部内容。
　教授说，这并不表示全部使用RC4算法软件全部轻易泄密，但它意味着RC4算法并不像大家原先认为那样安全。这一发觉可能促进大家重新设计无线通信网络，而且使用新加密算法。
三、 分组密码
3.1 DES加解密算法实现
3.1.1 DES加解密算法概述
1977年1月，美国政府颁布：采纳IBM企业设计方案作为非机密数据正式数据加密标准（DES Data Encryption Standard）。DES算法在POS、ATM、磁卡及智能卡（IC卡）、加油站、高速公路收费站等领域被广泛应用，以此来实现关键数据保密，如信用卡持卡人PIN加密传输，IC卡和POS间双向认证、金融交易数据包MAC校验等，均用到DES算法。
DES算法入口参数有三个：Key、Data、Mode。其中Key为8个字节共64位，是DES算法工作密钥；Data也为8个字节64位，是要被加密或被解密数据；Mode为DES工作方法，有两种：加密或解密。
DES是一个分组密码算法，它使用56位密钥，以64位为单位对数据分组进行加密解密（密文和明文分组长度相同，均为64位），DES加密和解密使用同一密钥，DES保密性依靠于密钥。
DES加密过程可简单描述为三个阶段：

3.1.2算法原理及设计思想
DES算法把64位明文输入块变为64位密文输出块，它所使用密钥也是64位,其功效是把输入64位数据块按位重新组合，并把输出分为L0 、R0两部分，每部分各长32位，其置换规则见下表：
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位换到第一位，第50位换到第2位，……，依这类推，最终一位是原来第7位。 L0、R0则是换位输出后两部分，L0是输出左32位，R0 是右32位，例：设置换前输入值为D1D2D3……D64，则经过初始置换后结果为：L0=D550……D8；R0=D57D49...D7。
经过26次迭代运算后，得到L16、R16，将此作为输入，进行逆置换，即得到密文输出。逆置换恰好是初始置逆运算，比如，第1位经过初始置换后，处于第40位，而经过逆置换，又将第40位换回到第1位，其逆置换规则以下表所表示：
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,
放大换位表
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,
单纯换位表
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,
在f(Ri,Ki)算法描述图中，S1,S2...S8为选择函数，其功效是把6bit数据变为4bit数据。下面给出选择函数Si(i=1,2......8)功效表：
选择函数Si
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、子密钥Ki(48bit)生成算法
初始Key值为64位，但DES算法要求，其中第8、16、......64位是奇偶校验位，不参与DES运算。故Key 实际可用位数便只有56位。即：经过缩小选择换位表1变换后，Key 位数由64 位变成了56位，此56位分为C0、D0两部分，各28位，然后分别进行第1次循环左移，得到C1、D1，将C1（28位）、D1（28位）合并得到56位，再经过缩小选择换位2，从而便得到了密钥K0（48位）。依这类推，便可得到K1、K2、......、K15，不过需要注意是，16次循环左移对应左移位数要依据下述规则进行：
循环左移位数1,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1
以上介绍了DES算法加密过程。 DES算法解密过程是一样，区分仅仅在于第一次迭代时用子密钥K15，第二次K14、……，最终一次用K0，算法本身并没有任何改变。

密钥生成过程步骤图以下所表示：

3.1.3程序关键算法分析
（1）S盒功效经过下列函数来实现，将48位输入转换成32位输出以下所表示：
void S_func(bool Out[32],const bool In[48])//将48位转换成32位
{
int j,m,n;
//膨胀后比特串分为8组，每组6比特。
for(j=0;j<8;j++,In+=6,Out+=4)
{
m = (In[0]*2)+In[5];
n = (In[1]*8)+(In[2]*4)+(In[3]*2)+In[4];
ByteToBit(Out,&SBox[j][m][n],4);
}
}
（2）函数F包含扩展置换，和子密钥异或，S盒变换及P盒变换，输入为32位，产生48位中间结果，并最终产生32比特输出
void F_func(bool In[32],const bool Ki[48])
{
static bool MR[48];
Transform(MR,In,EC,48);
Xor(MR, Ki, 48);
//膨胀后比特串分为8组，每组6比特。各组经过各自S盒后，又变为4比特，合并后又成为32比特。
S_func(In, MR);
//该32比特经过P变换后，输出比特串才是32比特
f(Ri-1,Ki)
Transform(In, In, PP, 32);
}
（3）下面为子密钥生成函数，输入种子密钥首先经过PC-1置换，将奇偶校验位删除，且剩下56位密钥打乱重排然后再生成子密钥，具体过程以下所表示：
void SetKey(char key[8])//生成子密钥
{
int i;
static bool K[64], *KL = &K[0], *KR = &K[28];
ByteToBit(K,key,64); //转换为二进制
Transform(K,K,EP1,56); //64比特密钥K，经过EP1后，生成56比特串。
//生成16个子密钥
for(i=0;i<16;i++)
{
//循环左移，合并
RotateL(KL,28,LOOP[i]);
RotateL(KR,28,LOOP[i]);
Transform(SubKey[i],K,EP2,48);
}
}
（4）下面为加密函数：
void CDES::Encryption(char out[8],char In[8])//加密函数
{
ByteToBit(M,In,64); //转换为二进制
Transform(M,M,IP,64);
for(int i=0;i<16;i++)
{
memcpy(tmp,Ri,32);
F_func(Ri,SubKey[i]);
Xor(Ri,Li,32);//将所得结果和明文左32位进行异或
memcpy(Li,tmp,32);//将明文左右32位交换
}
Transform(M, M, LP, 64);

BitToByte(out, M, 64);
// return(out);
}
（5）下面为解密函数，实际上是加密一个逆运算：
void CDES::Decryption(char out[8],char In[8])//解密函数(加密逆过程)
{
ByteToBit(M,In,64); //转换为二进制
Transform(M,M,IP,64);
for(int i=15;i>=0;i--)
{
memcpy(tmp,Li,32);
F_func(Li,SubKey[i]);
Xor(Li,Ri,32);
memcpy(Ri,tmp,32);
}
Transform(M, M, LP, 64);
BitToByte(out, M, 64);
// return(out);
}
3.1.4程序运行结果为：

程序总步骤图以下所表示：

3.1.5安全性分析
对DES安全性关键争论：
（1）、对DESS盒、迭代次数、密钥长度等设计准则争议
（2）、DES存在着部分弱密钥和半弱密钥
（3）、DES56位密钥无法抵御穷举工具
对于DES算法能够利用互补性、弱密钥和半弱密钥、密钥搜索、差分分析和线性分析等方法进行攻击。对于DES密码也可使用穷举密钥攻击,n=256≈7×106,即使使用每秒种能够计算一百万个密钥大型计算机,也需要算106天才能求得所使用密钥,所以看来是很安全。不过密码教授Diffie和Hellman指出,假如设计一个一微秒能够核实一个密钥超大规模集成片,那么它在一天内能够核实8.64×1010个密钥。假如由一个百万个这么集成片组成专用机,她们当初估量:这种专用机造价约为两千万美元。在五年内分期偿还,平均天天约需付一万美元。因为用穷举法破译平均只需要计算半个密钥空间,所以取得解平均时间为半天。为确保DES安全性，又出现了2DES，三重DES等。
3.2 AES加解密算法实现

3.2.1 AES算法概述
AES加密算法即密码学中高级加密标准（Advanced Encryption Standard，AES），又称Rijndael加密法，是美国联邦政府采取一个区块加密标准。这个标准用来替换原先DES，已经被多方分析且广为全世界所使用。经过五年甄选步骤，高级加密标准由美国国家标准和技术研究院 （NIST）于11月26日公布于FIPS PUB 197，并在5月26日成为有效标准。
AES基础要求是，采取对称分组密码体制，密钥长度最少支持为128、192、256，分组长度128位，AES加密数据块大小最大是256bit，不过密钥大小在理论上没有上限。AES加密有很多轮反复和变换。大致步骤以下：1、密钥扩展（KeyExpansion），2、初始轮（Initial Round），3、反复轮（Rounds），每一轮又包含：SubBytes、ShiftRows、MixColumns、AddRoundKey，4、最终轮（Final Round），最终轮没有MixColumns。
3.2.2算法原理及设计思想
AES 算法基于排列和置换运算。排列是对数据重新进行安排，置换是将一个数据单元替换为另一个。AES 使用多个不一样方法来实施排列和置换运算。
AES 是一个迭代、对称密钥分组密码，它能够使用128、192 和 256 位密钥，而且用 128 位（16字节）分组加密和解密数据。和公共密钥密码使用密钥对不一样，对称密钥密码使用相同密钥加密和解密数据。经过分组密码返回加密数据位数和输入数据相同。迭代加密使用一个循环结构，在该循环中反复置换和替换输入数据。
（1）首先将明文以字节为单位进行处理，以128位分组、128位密钥为例。先将明文按字节分成列组，将明文前四字节组成一列，接下来4个字节组成第二列，后面字节依次组成第三列和第四列，则组成了一个4乘4矩阵。
（2）AES也是由基础变换单位“轮”数次迭代而成。AES轮变换由四个不一样变换组成：
1） 字节替换变换
非线性字节替换，单独处理每个字节：
求该字节在有限域GF(28)上乘法逆，"0"被映射为本身，即对于α∈GF(28)，求β∈GF(28)，
使得α·β=β·α=1mod(x8+x4+x2+x+1)。
对上一步求得乘法逆作仿射变换
yi=xi + x(i+4)mod8 + x(i+6)mod8 + x(i+7)mod8 + ci
(其中ci是6310即第i位）
2) 行移位变换
行移位变换完成基于行循环位移操作，变换方法： 

即行移位变换作用于行上，第0行不变，第1行循环左移1个字节，第2行循环左移2个字节，第3行循环左移3个字节。
3） 列混合变换（最终一轮中没有）
逐列混合，方法：
b(x) = (03·x3 + 01·x2 + 01·x + 02) · a(x) mod(x4 + 1)

矩阵表示形式：

4） 和子密钥异或
只是简单将密钥按位异或到一个状态上。每轮加密密钥按次序取自扩展密钥，扩展密钥是由初始密钥扩展而成。
密钥扩展
AES密钥扩展算法输入值是4字(16字节)，输出值是一个44字(176字节)一维线性数组，为初始轮密钥加阶段和其它10轮中每一轮提供4字轮秘密钥，输入密钥直接被复制到扩展密钥数组前四个字，然后每次用四个字填充扩展密钥数组余下部分
3.2.3程序关键算法分析
（1） 程序编写过程中严格根据AES算法实施过程，将用到参数及函数封装在AES类中，再进行调用，以下所表示：
class AES
{
public:
AES(unsigned char* key);
virtual ~AES();
unsigned char* Cipher(unsigned char* input);
unsigned char* InvCipher(unsigned char* input);
void* Cipher(void* input, int length=0);
void* InvCipher(void* input, int length);
private:
unsigned char Sbox[256];
unsigned char InvSbox[256];
unsigned char w[11][4][4];
void KeyExpansion(unsigned char* key, unsigned char w[][4][4]);
unsigned char FFmul(unsigned char a, unsigned char b);
void SubBytes(unsigned char state[][4]);
void ShiftRows(unsigned char state[][4]);
void MixColumns(unsigned char state[][4]);
void AddRoundKey(unsigned char state[][4], unsigned char k[][4]);
void InvSubBytes(unsigned char state[][4]);
void InvShiftRows(unsigned char state[][4]);
void InvMixColumns(unsigned char state[][4]);
};
（2）先将输入明文按列序组合成4*4矩阵，直接和第0组密钥（即输入密钥）相加（异或），作为轮加密输入
然后循环10次进行SubBytes、ShiftRows、MixColumns、AddRoundKey运算，最终恢复原序列
（3）需要注意是最终一轮并不进行MixColumns（列混淆变换）加密过程函数Cipher,它只有一个参数，为输入明文，函数返回值为加密以后密文，解密过程和加密过程类似。
unsigned char* AES::Cipher(unsigned char* input)
{
unsigned char state[4][4];
int i,r,c;
//将明文按字节分成列组
for(r=0; r<4; r++)
{
for(c=0; c<4 ;c++)
{
state[r][c] = input[c*4+r];
}
}
AddRoundKey(state,w[0]);
for(i=1; i<=10; i++)
{
SubBytes(state);//字节替换
ShiftRows(state);//行移位
if(i!=10)MixColumns(state);//列混合（最终一轮除外）
AddRoundKey(state,w[i]);//和子密钥异或
}
for(r=0; r<4; r++)
{
for(c=0; c<4 ;c++)
{
input[c*4+r] = state[r][c];
}
}
return input;
}
（4）下面是每一轮变换中四个小变换实现函数以下：
//字节替换，经过SBox表来实现；
void AES::SubBytes(unsigned char state[][4])
{
int r,c;
for(r=0; r<4; r++)
{
for(c=0; c<4; c++)
{
state[r][c] = Sbox[state[r][c]];
}
}
}
//行移位作用于行上，第0行不变，第1行循环左移1个字节，第2行循环左移2个字节，第3行循环左移3个字节。
void AES::ShiftRows(unsigned char state[][4])
{
unsigned char t[4];
int r,c;
for(r=1; r<4; r++)
{
for(c=0; c<4; c++)
{
t[c] = state[r][(c+r)%4];
}
for(c=0; c<4; c++)
{
state[r][c] = t[c];
}
}
}
//列混合 FFmul为有限域GF(28)上乘法，标准算法应该是循环8次（b和a每一位相乘，结果相加），但这里只用到最低2位，解密时用到逆列混淆也只用了低4位，所以在这里高4位运算是多出，只计算低4位。
void AES::MixColumns(unsigned char state[][4])
{
unsigned char t[4];
int r,c;
for(c=0; c< 4; c++)
{
for(r=0; r<4; r++)
{
t[r] = state[r][c];
}
for(r=0; r<4; r++)
{
state[r][c] = FFmul(0x02, t[r])
^ FFmul(0x03, t[(r+1)%4])
^ FFmul(0x01, t[(r+2)%4])
^ FFmul(0x01, t[(r+3)%4]);
}
}
}
//和子密钥异或
void AES::AddRoundKey(unsigned char state[][4], unsigned char k[][4])
{
int r,c;
for(c=0; c<4; c++)
{
for(r=0; r<4; r++)
{
state[r][c] ^= k[r][c];//异或运算
}
}
}
//密钥扩展

//将前一列即第n-1组第三列四个字节循环左移1个字节，并对每个字节进行字节替换变换SubBytes,将第一行(即第一个字节)和轮常量rc[n]相加，最终再和前一组该列相加
void AES::KeyExpansion(unsigned char* key, unsigned char w[][4][4])
{
Int i,j,r,c;
unsigned char rc[] = {0x01, 0x02, 0x04, 0x08, 0x10, 0x20, 0x40, 0x80, 0x1b, 0x36};
for(r=0; r<4; r++)
{
for(c=0; c<4; c++)
{
w[0][r][c] = key[r+c*4];
}
}
for(i=1; i<=10; i++)
{
for(j=0; j<4; j++)
{
unsigned char t[4];
for(r=0; r<4; r++)
{
t[r] = j ? w[i][r][j-1] : w[i-1][r][3];
}
if(j == 0)
{
unsigned char temp = t[0];
for(r=0; r<3; r++)
{
t[r] = Sbox[t[(r+1)%4]];
}
t[3] = Sbox[temp];
t[0