本发明涉及一种新型密钥生成的shield分组密码实现方法、装置及可读存储介质。
背景技术：
近年来，世界超级计算机的发展是突飞猛进的，它的计算能力在近年来一直都是呈现指数增长，密码系统及密码算法本身的安全性收到了严重的冲击。利用强大的计算能力进行最基本的暴力攻击方法破解密码的密钥。密钥长度决定了可能的密钥数量，因此提高密钥长度，是增加了破解难度，对于des被破解之后，采用3des方案，目的是通过增加密钥长度来提高密码算法的安全性，而aes密码算法，在从事高机密信息加密传输，则使用密钥长度为256位的aes算法。
分组密码可用于数据加解密、消息认证码和哈希函数的构造，因此在信息安全保障机制中扮演着非常重要的角色。目前，绝大多数密码系统在理论上遵循柯克霍夫原则，在柯克霍夫原则提出即使密码系统的任何细节已为人悉知，只要密钥未被泄露，它也应该是安全的，这表明一个密码系统的安全性仅依赖于密钥的安全性，并且原则当中提到，密钥必须易于沟通和记忆，而不需写下，且双方可以容易的改变密钥。
在国内分组密码算法设计竞赛当中，设计目标是安全性高、可扩展性好、适应性强的分组密码，以满足多个行业领域对分组密码算法的应用需求，要求分组密码算法的分组长度与密钥长度至少分别为128位。目前分组密码算法，具有分组长度、密钥长度分别为128位的密码算法，例如aes密码算法、sm4密码算法以及一些比较知名的轻量级分组密码skinny与simon等，这些密码加密函数采用spn结构与feistel结构，密钥扩展函数也同样采用spn结构或feistel结构。而这种分组密码算法的密钥扩展函数存在问题：一是密钥长度过长，不利用保存，保存的需要花费大量的资源；二是产生的轮密钥之间有相互关联关系，暴露一些信息，容易遭受相关密钥攻击，参见文献：徐林宏,郭建胜,崔竞一,李明明.piccolo算法的相关密钥-不可能差分攻击[j].软件学报,2019,30(8):2349-2361。另外，分组长度为128位的密码算法的s盒替换一般采用8×8的s盒，造成大量资源的开销，参见文献：王沁,梁静,齐悦.一种有效缩减aes算法s盒面积的组合逻辑优化设计[j].电子学报,2010,38(4):205-208，并且s盒是固定不动，容易遭受侧信道等攻击。
技术实现要素：
本发明提供了本发明涉及一种新型密钥生成的shield分组密码实现方法、装置及可读存储介质，其目的在于，克服现有技术中的分组密码算法密钥扩展函数存在不易保存、易受攻击以及8×8的固定s盒占用资源占用多等问题。
本发明的技术方案如下：
一方面，一种新型密钥生成的shield分组密码实现方法，包括以下步骤：
步骤1：将128位明文/密文加载至寄存器，进行加/解密运算；
步骤2：将所述待加/解密数据按以下步骤进行32轮迭代轮运算操作；
若输入的待加密数据进行加密运算操作，则在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待加密数据进行轮密钥加变换、s盒替换变换、行移位变换及列混合变换，以列混合变换后得到的数据作为下一轮的待加密数据，完成32轮迭代运算操作后，再进行一次轮密钥加变换，得到密文数据；
若输入的待解密数据进行解密运算操作，先将输入的待解密数据进行一次轮密钥加逆变换，然后在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待解密数据进行列混合逆变换、行移位逆变换、s盒替换逆变换及轮密钥加逆变换，以轮密钥加逆变换后得到的数据作为下一轮的待解密数据，完成32轮迭代运算操作后，得到明文数据；
所述列混合逆变换和列混合变换、行移位逆变换和行移位变换、s盒替换逆变换和s盒替换变换以及轮密钥加逆变换和轮密钥加变换均互为逆运算。
进一步地，所述s盒替换变换操作具体方法如下，且s盒为具有对合性质的s盒sb0与s盒sb1，在加密与解密过程中，s盒sb0与s盒sb1元素不变：
128位需要进行s盒替换变换的中间数据被分成32个4位数据单元a1、a2、…、a31、a32，随着加/解密轮数变化，数据单元分别进行s盒sb0与s盒sb1递增/递减替换变换；
在加密过程的第r轮s盒替换变换中，前32-r 1个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1；：
在解密过程的第r轮s盒替换变换中，前r个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1，1≤r≤32。
加密第1轮(32个4位数据单元进行sb0替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb0(a32)
加密第2轮(第1至第31个4位数据单元进行sb0替换变换，第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb1(a32)
加密第3轮(第1至第30个4位数据单元进行sb0替换变换，第31至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb1(a31)||sb1(a32)
…
加密第31轮(第1至第2个4位数据单元进行sb0替换变换，第3个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
加密第32轮(第1个4位数据单元进行sb0替换变换，第2个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb1(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)解密第1轮(第1个4位数据单元进行sb0替换变换，第2个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb1(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
解密第2轮(第1至第2个4位数据单元进行sb0替换变换，第3个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
解密第3轮(第1至第3个4位数据单元进行sb0替换变换，第4个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
…
解密第31轮(第1至第31个4位数据单元进行sb0替换变换，第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb1(a32)
解密第32轮(32个4位数据单元进行sb0替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb0(a32)。
进一步地，在所述轮密钥加变换与轮密钥加逆变换操作中使用的轮密钥构造过程如下：
输入原始密钥key，将原始密钥key从高位到低位分为三个key1、key2与key3子密钥，将key1子密钥进行一次md5哈希函数变换得到第1轮至第16轮中，轮密钥加变换中使用的128位轮密钥rkey1；
将key2子密钥进行一次md5哈希函数变换得到第16轮至32轮中，轮密钥加变换中使用的128位轮密钥rkey2；
将key3子密钥进行一次md5哈希函数变换得到最后一次轮密钥加变换中使用的128位轮密钥rkey3。
轮密钥的构造，是将分得的子密钥通过md5哈希函数进行运算得到轮密钥；
进一步地，所述轮密钥加变换和轮密钥加逆变换的操作过程如下：
轮密钥加变换操作是将明文或每一轮中间值state与轮密钥rkey以及轮常数rc进行异或运算；
轮密钥加变换如下：
轮密钥加逆变换使用轮密钥及轮常数的顺序与轮密钥加变换相反，轮密钥加逆变换如下：
其中，state表示进行轮密钥加变换操作的数据，rc[i]表示轮常数rc的第i位，rc[33-i]表示轮常数rc的第33-i位。
一方面，一种新型密钥生成的shield分组密码实现装置，包括：
数据加载模块：将128位明文/密文加载至寄存器，进行加/解密运算；
轮运算操作模块：将加载到寄存器的所述待加/解密数据进行32轮迭代轮运算操作；
若输入的待加密数据进行加密运算操作，则在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待加密数据进行轮密钥加变换单元、s盒替换变换单元、行移位变换单元及列混合变换单元，以列混合变换单元得到的数据作为下一轮的待加密数据，完成32轮迭代运算操作后，再进行一次轮密钥加变换，得到密文数据；
若输入的待解密数据进行解密运算操作，先将输入的待解密数据进行一次轮密钥加逆变换，然后在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待解密数据进行列混合逆变换单元、行移位逆变换单元、s盒替换逆变换单元及轮密钥加逆变换单元，以轮密钥加逆变换单元得到的数据作为下一轮的待解密数据，完成32轮迭代运算操作后，得到明文数据；
所述列混合逆变换单元和列混合变换单元、行移位逆变换单元和行移位变换单元、s盒替换逆变换单元和s盒替换变换单元，以及轮密钥加逆变换单元和轮密钥加变换单元均互为逆运算单元。
进一步地，s盒替换变换单元的具体操作过程如下，且s盒为具有对合性质的s盒sb0与s盒sb1，在加密与解密过程中，s盒sb0与s盒sb1元素不变：
128位需要进行s盒替换变换的中间数据被分成32个4位数据单元a1、a2、…、a31、a32，随着加/解密轮数变化，数据单元分别进行s盒sb0与s盒sb1递增/递减替换变换；
在加密过程的第r轮s盒替换变换中，前32-r 1个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1；
在解密过程的第r轮s盒替换变换中，前r个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1，1≤r≤32。
进一步地，在所述轮密钥加变换单元与轮密钥加逆变换单元操作中使用的轮密钥构造过程如下：
输入原始密钥key，将原始密钥key从高位到低位分为三个key1、key2与key3子密钥，将key1子密钥进行一次md5哈希函数变换得到第1轮至第16轮中，轮密钥加变换中使用的128位轮密钥rkey1；
将key2子密钥进行一次md5哈希函数变换得到第16轮至32轮中，轮密钥加变换中使用的128位轮密钥rkey2；
将key3子密钥进行一次md5哈希函数变换得到最后一次轮密钥加变换中使用的128位轮密钥rkey3。
轮密钥的构造，是将分得的子密钥通过md5哈希函数进行运算得到轮密钥；
进一步地，所述轮密钥加变换单元和轮密钥加逆变换单元的操作过程如下：
轮密钥加变换单元是用于将明文或每一轮中间值state与轮密钥rkey以及轮常数rc进行异或运算；
轮密钥加变换如下：
轮密钥加逆变换使用轮密钥及轮常数的顺序与轮密钥加变换相反，轮密钥加逆变换如下：
其中，state表示进行轮密钥加变换操作的数据，rc[i]表示轮常数rc的第i位，rc[33-i]表示轮常数rc的第33-i位。
另一方面，一种可读存储介质，包括计算机程序指令，其特征在于：所述计算机程序指令被处理终端执行时使所述处理终端执行一种新型密钥生成的shield分组密码实现方法。
有益效果
本发明提供了一种新型密钥生成的shield分组密码实现方法、装置及可读存储介质，本发明技术方案中的方法从加密算法出发，其算法设计原则是安全，其次是速度、算法复杂度、串并行。由于shield算法当中，采用了md5作为密钥产生算法，md5在运算过程4个函数分别是需要16轮操作，在实现相当于需要64次运算操作，导致需要花费大量运算操作时间，为了满足加密函数与密钥产生操作实现时间同步运算操作，将轮函数当中轮密钥加变换的轮密钥每16轮进行一次md5运算，更新其密钥，利用哈希函数数字指纹特性，使得轮函数当中密钥空间变得很大，大大增加了密钥空间，从而提高算法的安全性，极大增加密码分析难度，从而使得目前的攻击方法很难做到攻击shield算法，特别是密码算法当中相关密钥攻击。
在shield算法当中，我们原始(主)密钥不再需要大于等于分组长度，不过原始(主)密钥长度长一些，组合方式复杂一些，对于算法本身安全而言是非常有利的，但原始(主)密钥甚至可以做到只有6位密钥，都可以达到算法高安全的应用需求。从而减少密钥管理中的密钥安全存储压力，节省密钥存储资源。
shield算法的s盒替换变换采用两个4×4的s盒进行128位数据替换变换，将数据分成32个4位单元，则随着加/解密轮数变化，数据单元分别进行对应s盒递增/递减替换操作，相比固定一个8×8的s盒，可以节省了大量资源，并且这种活动性s盒变换，可以抵抗侧信道等攻击方法。
综上所述，本发明的技术方案所述的一种新型密钥生成的shield分组密码实现方法在保证性能的前提下，大大提高算法的安全性，能有效地抵抗现有攻击方法。
附图说明
图1为本发明实例所述方法的加密过程示意图；
图2为本发明实例所述加密方法对应的解密过程示意图；
图3为md5算法4轮变换运算图；
图4为md5算法核心变换运算图；
图5为行移位变换运算关系图；
图6为行移位逆变换运算关系图。
具体实施方式
下面结合附图和实施例对本发明做进一步的说明。
一种新型密钥生成的shield分组密码实现方法，在本实例中，shield分组密码的实现方法采用spn结构，分组长度为128位，密钥长度不固定，迭代32轮。shield加密运算流程如图1所示，包含密钥生成函数与加密轮函数，加密轮函数运算中包含轮密钥加变换(addroundkey)，s盒替换变换(subcells)，行移位变换(shiftrows)，列混合变换(mixcolumns)四个模块。解密运算流程如图2，算法解密轮运算包含，列混合逆变换(invmixcolumns)，行移位逆变换(invshiftrows)，s盒替换逆变换(invsubcells)，轮密钥加逆变换(addroundkey)四个模块。
shield加密描述如下算法1，参考图1。
算法1：shield加密
输入:plaintext,key；
输出:ciphertxet；
其中：key为原始密钥，rkey1、rkey2及rkey3为密钥生成函数生成的轮密钥。
密钥生成函数(keyschedule)：输入原始密钥key，将原始密钥key分为三分key1、key2与key3子密钥，key1子密钥进行一次md5哈希函数变换得到第1轮至第16轮，16次轮密钥加变换128位加密轮密钥rkey1；key2子密钥进行一次md5哈希函数变换得到第16轮至32轮，16次轮密钥加变换128位加密轮密钥rkey2，key3子密钥进行一次md5哈希函数变换得到最后一次(第33次)轮密钥加变换128位轮密钥rkey3。
md5哈希函数是输入的处理信息分组为512位，且每个keyi(1≤i≤3)信息又被划分为16个32位子分组，经过算法处理后，md5哈希函数的输出信息为4个32位子分组，将这4个分组联级得到128位哈希值消息作为轮密钥rkeyi(1≤i≤3)输出。
对于一个任意长度的keyi(1≤i≤3)作为md5哈希函数输入，计算输出得到128位哈希值作为轮密钥rkeyi(1≤i≤3)，那么对输入子密钥的处理包括如下几个步骤：
(1)第一步：附加填充位。keyi(1≤i≤3)被补充使其位数b满足b＝448mod512，也就是说补充上去的位数，使之得到448，补充的内容第一位1以外，其他位补0。(例如，keyi为56位，则填充392位，使其变成448。)
(2)第二步：附加消息长度。附上64位以表示keyi(1≤i≤3)的长度位数，并将其加进填充后的消息末尾。
(3)第三步：初始化算法参数。128位的缓冲区用来存放中间结果及最后结果，128位分为32位(a，b，c与d)4个寄存器，它们的初始值为：a＝0x01234567，b＝0x89abcdef，c＝0xfedcba98，d＝0x76543210。
(4)第四步：对输入的数据消息做运算。运算算法包含4轮“循环”的压缩函数，该4轮压缩函数结构相似，但每一轮各有一个逻辑函数，分别表示为：ff、gg、hh及ii，如图3所示，每个逻辑函数具体表示为：
ff(a,b,c,d,m[k],s,t[i])<<<s:a←b ((a f(b,c,d) m[k] t[i])<<<s)
gg(a,b,c,d,m[k],s,t[i])<<<s:a←b ((a g(b,c,d) m[k] t[i])<<<s)
hh(a,b,c,d,m[k],s,t[i])<<<s:a←b ((a h(b,c,d) m[k] t[i])<<<s)
hh(a,b,c,d,m[k],s,t[i])<<<s:a←b ((a i(b,c,d) m[k] t[i])<<<s)
其中<<<是表示循环左移操作，t[i]为4294967296×abs(sin(i))的32位整数部分，i的单位弧度。
每一循环由对缓存a，b，c与d的16步操作组成，每一步骤详细操作如图4中部分所示。其中r在1至4轮分别代以f，g，h，i。第四次循环的输出加到第一次循环的四个缓存输入上作为处理下一分组的输入。
(5)第五步：输出轮密钥结果。在产生输出128位的哈希值作为轮密钥，输出以高位字节a开始，低位字节d结束。
轮密钥加变换(addroundkey)：将128位明文或每一轮中间值state与通过md5算法生成的128位轮密钥rkey0、rkey1及rkey2以及轮常数rc进行异或运算，其中轮常数与数据最后8位进行异或运算：
当1≤i≤16时，轮密钥加变换：
当17≤i≤32时，轮密钥加变换：
当i为33时，轮密钥加变换：
为了节省算法资源的开销，采用led密码算法的线性移位寄存器的方式产生33个轮常数，轮常数排列数组为rc[33]＝{0x01，0x03，0x07，0x0f，0x1f，0x3e，0x3d，0x3b，0x37，0x2f，0x1e，0x3c，0x39，0x33，0x27，0x0e，0x1d，0x3a，0x35，0x2b，0x16，0x2c，0x18，0x30，0x21，0x02，0x05，0x0b，0x17，0x2e，0x1c，0x38，0x31}。
s盒替换变换(subcells)：s盒变换是算法的非线性组件，加/解密s盒采用midori密码算法的sb0与sb1，该sb0与sb1是具有对合性质的两个s盒，从而加密s盒替换变换sb0与sb1元素表示为：sb0[16]＝{0xc，0xa，0xd，0x3，0xe，0xb，0xf，0x7，0x8，0x9，0x1，0x5，0x0，0x2，0x4，0x6}，sb1[16]＝{0x1，0x0，0x5，0x3，0xe，0x2，0xf，0x7，0xd，0xa，0x9，0xb，0xc，0x8，0x4，0x6}。128位加密中间值state分成32个4位数据单元a1、a2、…、a31、a32，则随着加密轮数变化，数据单元分别进行sb0递减与sb1递增替换变换，运算关系如下：
加密第1轮(32个4位数据单元进行sb0替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb0(a32)
加密第2轮(第1至第31个4位数据单元进行sb0替换变换，第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb1(a32)
加密第3轮(第1至第30个4位数据单元进行sb0替换变换，第31至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb1(a31)||sb1(a32)
…
加密第31轮(第1至第2个4位数据单元进行sb0替换变换，第3个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
加密第32轮(第1个4位数据单元进行sb0替换变换，第2个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb1(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
行移位变换：对于16个单元组成4×4矩阵，矩阵的每一行向左循环不同的单元移量，第零行单元移量保持不变，第一行单元移量为循环左移一个单元，第二行单元移量为循环左移两个单元，第三行单元移量为循环左移3个单元，行移位变换运算关系如图5所示。
列混合变换：采用aes密码的列混淆变换矩阵m，m矩阵组成元素表示如下，其中数据以16进制表示：
列混合变换运算是列混合变换矩阵m与state中16个单元组成4×4矩阵在有限域gf(28)上的乘法变换，变换公式，其中数据以16进制表示：
shield解密算法描述如下算法2，参考图2。
输入:ciphertxet,key；
输出:plaintext；
shield解密模块使用了加密运算变换中四种逆变换(轮密钥加逆变换、行移位逆变换、列混淆逆变换以及s盒替换逆变换)，其中s盒替换满足对合性质，从而加密与解密当中，s盒不变，还是s盒sb0与s盒sb1，以加密运算相反的顺序对密文进行解密。
轮密钥加逆变换使用轮密钥以及轮常数的顺序与轮密钥加变换相反，轮密钥加逆变换如下：
当i为33时，轮密钥加变换：
当32≥i≥17时，轮密钥加变换：
当16≥i≥1时，轮密钥加变换：
s盒替换变换(invsubcells)：sb0与sb1是具有对合性质的两个s盒，从而解密与加密同为一个s盒元素。128位解密中间值state分成32个4位数据单元a1、a2、…、a31、a32，则随着解密轮数变化，数据单元分别进行sb0递增与sb1递减替换变换，运算关系如下：
解密第1轮(第1个4位数据单元进行sb0替换变换，第2个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb1(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
解密第2轮(第1至第2个4位数据单元进行sb0替换变换，第3个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb1(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
解密第3轮(第1至第3个4位数据单元进行sb0替换变换，第4个至第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb1(a30)||sb1(a31)||sb1(a32)
…
解密第31轮(第1至第31个4位数据单元进行sb0替换变换，第32个4位数据单元进行sb1替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb1(a32)
解密第32轮(32个4位数据单元进行sb0替换变换)：
state←sb0(a1)||sb0(a2)||sb0(a3)||…||sb0(a30)||sb0(a31)||sb0(a32)
行移位逆变换：对于16个单元组成4×4矩阵，矩阵的每一行向右循环不同的单元移量，第零行单元移量保持不变，第一行单元循环右移一个单元，第二行单元循环右移两个单元，第三行单元循环右移3个单元，行移位逆变换运算关系如图6所示。
列混合逆变换：采用aes密码的列混淆逆变换矩阵m-1，m-1矩阵组成元素表示如下，其中数据以16进制表示：
列混合变换运算是列混合变换矩阵m-1与state中16个单元组成4×4矩阵在有限域gf(28)上的乘法变换，变换公式，其中数据以16进制表示：
shield算法测试向量(数据以16进制表示)：
(1)明文：0000_0000_0000_0000_0000_0000_0000_0000
(1)密钥：012345(key1为01，key2为23，key3为45)
第一个轮密钥：96a3_be3c_f272_e017_046d_1d26_74a5_2bd3
第二个轮密钥：3769_3cfc_7480_49e4_5d87_b8c7_d8b9_aacd
第三个轮密钥：6c83_49cc_7260_ae62_e3b1_3968_31a8_398f
(1)密文：8273__ac00_4f64_9a86_dd6a_955c_c784_57d0
(2)明文：0000_0000_0000_0000_0000_0000_0000_0000
(2)密钥：012345(key1为0，key2为12，key3为345)
第一个轮密钥：cfcd_2084_95d5_65ef_66e7_dff9_f987_64da
第二个轮密钥：c20a_d4d7_6fe9_7759_aa27_a0c9_9bff_6710
第三个轮密钥：289d_ff07_669d_7a23_de0e_f88d_2f71_29e7
(2)密文：b22e__6064_2105_3da9_76c0_3067_69d8_47ea
(3)明文：ffff_ffff_ffff_ffff_ffff_ffff_ffff_ffff
(3)密钥：012345(key1为01，key2为23，key3为45)
第一个轮密钥：96a3_be3c_f272_e017_046d_1d26_74a5_2bd3
第二个轮密钥：3769_3cfc_7480_49e4_5d87_b8c7_d8b9_aacd
第三个轮密钥：6c83_49cc_7260_ae62_e3b1_3968_31a8_398f
(3)密文：d167_b8bc_fddd_cf4c_2df2_8a19_8a36_29b0
按照本发明所提出的实现方法，利用给出的测试向量，实现者对明文进行加密处理得到测试向量中所述的密文，对密文进行解密处理得到测试向量中所述的明文。
本发明所述的shield算法分组密码实现方法在软硬实现当中，应该与相比传统的aes、sm4等分组密码算法比较，占用资源相当，甚至还少一些。
一种新型密钥生成的shield分组密码实现装置，包括：
数据加载模块：将128位明文/密文加载至寄存器，进行加/解密运算；
轮运算操作模块：将加载到寄存器的所述待加/解密数据进行32轮迭代轮运算操作；
若输入的待加密数据进行加密运算操作，则在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待加密数据进行轮密钥加变换单元、s盒替换变换单元、行移位变换单元及列混合变换单元，以列混合变换单元得到的数据作为下一轮的待加密数据，完成32轮迭代运算操作后，再进行一次轮密钥加变换，得到密文数据；
若输入的待解密数据进行解密运算操作，先将输入的待解密数据进行一次轮密钥加逆变换，然后在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待解密数据进行列混合逆变换单元、行移位逆变换单元、s盒替换逆变换单元及轮密钥加逆变换单元，以轮密钥加逆变换单元得到的数据作为下一轮的待解密数据，完成32轮迭代运算操作后，得到明文数据；
所述列混合逆变换单元和列混合变换单元、行移位逆变换单元和行移位变换单元、s盒替换逆变换单元和s盒替换变换单元，以及轮密钥加逆变换单元和轮密钥加变换单元均互为逆运算单元。
s盒替换变换单元的具体操作过程如下，且s盒为具有对合性质的s盒sb0与s盒sb1，在加密与解密过程中，s盒sb0与s盒sb1元素不变：
128位需要进行s盒替换变换的中间数据被分成32个4位数据单元a1、a2、…、a31、a32，随着加/解密轮数变化，数据单元分别进行s盒sb0与s盒sb1递增/递减替换变换；
在加密过程的第r轮s盒替换变换中，前32-r 1个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1；
在解密过程的第r轮s盒替换变换中，前r个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1，1≤r≤32。
在所述轮密钥加变换单元与轮密钥加逆变换单元操作中使用的轮密钥构造过程如下：
输入原始密钥key，将原始密钥key从高位到低位分为三个key1、key2与key3子密钥，将key1子密钥进行一次md5哈希函数变换得到第1轮至第16轮中，轮密钥加变换中使用的128位轮密钥rkey1；
将key2子密钥进行一次md5哈希函数变换得到第16轮至32轮中，轮密钥加变换中使用的128位轮密钥rkey2；
将key3子密钥进行一次md5哈希函数变换得到最后一次轮密钥加变换中使用的128位轮密钥rkey3。
轮密钥的构造，是将分得的子密钥通过md5哈希函数进行运算得到轮密钥；
所述轮密钥加变换单元和轮密钥加逆变换单元的操作过程如下：
轮密钥加变换单元是用于将明文或每一轮中间值state与轮密钥rkey以及轮常数rc进行异或运算；
轮密钥加变换如下：
轮密钥加逆变换使用轮密钥及轮常数的顺序与轮密钥加变换相反，轮密钥加逆变换如下：
其中，state表示进行轮密钥加变换操作的数据，rc[i]表示轮常数rc的第i位，rc[33-i]表示轮常数rc的第33-i位。
应当理解，本发明各个实施例中的功能单元模块可以集中在一个处理单元中，也可以是各个单元模块单独物理存在，也可以是两个或两个以上的单元模块集成在一个单元模块中，可以采用硬件或软件的形式来实现。
一种可读存储介质，包括计算机程序指令，其特征在于：所述计算机程序指令被处理终端执行时使所述处理终端执行一种新型密钥生成的shield分组密码实现方法，其有益效果参见方法部分的有益效果，在此不再赘述。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
1.一种新型密钥生成的shield分组密码实现方法，其特征在于，包括以下步骤：
步骤1：将128位明文/密文加载至寄存器，进行加/解密运算；
步骤2：将所述待加/解密数据按以下步骤进行32轮迭代轮运算操作；
若输入的待加密数据进行加密运算操作，则在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待加密数据进行轮密钥加变换、s盒替换变换、行移位变换及列混合变换，以列混合变换后得到的数据作为下一轮的待加密数据，完成32轮迭代运算操作后，再进行一次轮密钥加变换，得到密文数据；
若输入的待解密数据进行解密运算操作，先将输入的待解密数据进行一次轮密钥加逆变换，然后在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待解密数据进行列混合逆变换、行移位逆变换、s盒替换逆变换及轮密钥加逆变换，以轮密钥加逆变换后得到的数据作为下一轮的待解密数据，完成32轮迭代运算操作后，得到明文数据；
所述列混合逆变换和列混合变换、行移位逆变换和行移位变换、s盒替换逆变换和s盒替换变换以及轮密钥加逆变换和轮密钥加变换均互为逆运算。
2.根据权利要求1所述的一种新型密钥生成的shield分组密码实现方法，其特征在于，所述s盒替换变换操作具体方法如下，且s盒为具有对合性质的s盒sb0与s盒sb1，在加密与解密过程中，s盒sb0与s盒sb1元素不变：
128位需要进行s盒替换变换的中间数据被分成32个4位数据单元a1、a2、…、a31、a32，随着加/解密轮数变化，数据单元分别进行s盒sb0与s盒sb1递增/递减替换变换；
在加密过程的第r轮s盒替换变换中，前32-r 1个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1；：
在解密过程的第r轮s盒替换变换中，前r个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1，1≤r≤32。
3.根据权利要求1所述的一种新型密钥生成的shield分组密码实现方法，其特征在于，在所述轮密钥加变换与轮密钥加逆变换操作中使用的轮密钥构造过程如下：
输入原始密钥key，将原始密钥key从高位到低位分为三个key1、key2与key3子密钥，将key1子密钥进行一次md5哈希函数变换得到第1轮至第16轮中，轮密钥加变换中使用的128位轮密钥rkey1；
将key2子密钥进行一次md5哈希函数变换得到第16轮至32轮中，轮密钥加变换中使用的128位轮密钥rkey2；
将key3子密钥进行一次md5哈希函数变换得到最后一次轮密钥加变换中使用的128位轮密钥rkey3。
4.根据权利要求3所述的方法，其特征在于，所述轮密钥加变换和轮密钥加逆变换的操作过程如下：
轮密钥加变换操作是将明文或每一轮中间值state与轮密钥rkey以及轮常数rc进行异或运算；
轮密钥加变换如下：
轮密钥加逆变换使用轮密钥及轮常数的顺序与轮密钥加变换相反，轮密钥加逆变换如下：
其中，state表示进行轮密钥加变换操作的数据，rc[i]表示轮常数rc的第i位，rc[33-i]表示轮常数rc的第33-i位。
5.一种新型密钥生成的shield分组密码实现装置，其特征在于，包括：
数据加载模块：将128位明文/密文加载至寄存器，进行加/解密运算；
轮运算操作模块：将加载到寄存器的所述待加/解密数据进行32轮迭代轮运算操作；
若输入的待加密数据进行加密运算操作，则在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待加密数据进行轮密钥加变换单元、s盒替换变换单元、行移位变换单元及列混合变换单元，以列混合变换单元得到的数据作为下一轮的待加密数据，完成32轮迭代运算操作后，再进行一次轮密钥加变换，得到密文数据；
若输入的待解密数据进行解密运算操作，先将输入的待解密数据进行一次轮密钥加逆变换，然后在第1轮至第32轮迭代轮运算中的每一轮轮运算中依次对每一轮输入的待解密数据进行列混合逆变换单元、行移位逆变换单元、s盒替换逆变换单元及轮密钥加逆变换单元，以轮密钥加逆变换单元得到的数据作为下一轮的待解密数据，完成32轮迭代运算操作后，得到明文数据；
所述列混合逆变换单元和列混合变换单元、行移位逆变换单元和行移位变换单元、s盒替换逆变换单元和s盒替换变换单元，以及轮密钥加逆变换单元和轮密钥加变换单元均互为逆运算单元。
6.根据权利要求5所述的装置，其特征在于，s盒替换变换单元的具体操作过程如下，且s盒为具有对合性质的s盒sb0与s盒sb1，在加密与解密过程中，s盒sb0与s盒sb1元素不变：
128位需要进行s盒替换变换的中间数据被分成32个4位数据单元a1、a2、…、a31、a32，随着加/解密轮数变化，数据单元分别进行s盒sb0与s盒sb1递增/递减替换变换；
在加密过程的第r轮s盒替换变换中，前32-r 1个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1；
在解密过程的第r轮s盒替换变换中，前r个数据单元进行s盒sb0替换操作，剩余的数据单元进行s盒sb1，1≤r≤32。
7.根据权利要求5所述的装置，其特征在于，在所述轮密钥加变换单元与轮密钥加逆变换单元操作中使用的轮密钥构造过程如下：
输入原始密钥key，将原始密钥key从高位到低位分为三个key1、key2与key3子密钥，将key1子密钥进行一次md5哈希函数变换得到第1轮至第16轮中，轮密钥加变换中使用的128位轮密钥rkey1；
将key2子密钥进行一次md5哈希函数变换得到第16轮至32轮中，轮密钥加变换中使用的128位轮密钥rkey2；
将key3子密钥进行一次md5哈希函数变换得到最后一次轮密钥加变换中使用的128位轮密钥rkey3。
8.根据权利要求7所述的装置，其特征在于，所述轮密钥加变换单元和轮密钥加逆变换单元的操作过程如下：
轮密钥加变换单元是用于将明文或每一轮中间值state与轮密钥rkey以及轮常数rc进行异或运算；
轮密钥加变换如下：
轮密钥加逆变换使用轮密钥及轮常数的顺序与轮密钥加变换相反，轮密钥加逆变换如下：
其中，state表示进行轮密钥加变换操作的数据，rc[i]表示轮常数rc的第i位，rc[33-i]表示轮常数rc的第33-i位。
9.一种可读存储介质，包括计算机程序指令，其特征在于：所述计算机程序指令被处理终端执行时使所述处理终端执行权利要求1至4任一项所述的方法。
技术总结