本发明提出一种基于明文相似矩阵的对称全同态加密方法，用于解决现有对称全同态加密效率低的技术问题。实现步骤为：用户根据要求生成两个长度相等的大素数，并根据生成的大素数构建剩余类环，根据剩余类环构建一般线性群，计算同态计算公钥和对称密钥，使用对称密钥对明文矩阵的相似矩阵进行加密，使用对称密钥对密文矩阵进行解密；云服务器使用同态计算公钥对密文矩阵进行同态计算；用户使用对称密钥对同态密文矩阵进行解密。本发明密钥选取和加密过程简单，将明文矩阵随机隐藏，提高加密算法的安全性，在密文计算过程中没有引入噪声，可以根据需求对密文矩阵进行任意计算。本发明可用于云计算、大数据环境等对重要数据的全程密态保护。
1.一种基于明文相似矩阵的对称全同态加密方法，其特征在于包括如下步骤：(1)参数生成：用户根据安全要求随机生成两个长度相等的大素数p和q；
(2)用户构建剩余类环和一般线性群，实现步骤为:(2.1)用户构建关于大素数p的剩余类环 和关于大素数q的剩余类环(2.2)用户利用剩余类环 构建模p意义下的所有n阶可逆矩阵构成的一般线性群同时利用剩余类环 构建模q意义下的所有n阶可逆矩阵构成的一般线性群(3)用户计算同态计算公钥和对称密钥，实现步骤为：(3.1)用户计算RSA模数N，N＝pq，并将RSA模数N作为同态计算公钥；
(3.2)用户在一般线性群 中随机选取一个n阶可逆矩阵A，并求取可逆矩阵A模p意义下的逆矩阵A-1，A、A-1和p构成对称密钥sk，sk＝(A,A-1,p)；
(4)用户对明文矩阵进行加密，实现步骤为：(4.1)用户根据需要在一般线性群 中选取k个n阶可逆的明文矩阵并分别计算k个明文矩阵Mk的相似矩阵 其中，k≥2， 是明文矩阵Mk中第i行第j列的元素， 是相似矩阵Bk中第i行第j列的元素；
(4.2)用户在一般线性群 中随机选取k个n阶矩阵 并根据中国剩余定理计算k个密文矩阵
其中， 是矩阵Dk中第i行第j列的元素， 是密文矩阵Ck中第i行第j列的元素且对任意的i,j＝1,2,…,n满足同余式方程组：(5)用户对密文矩阵进行解密，实现步骤为：(5.1)用户通过同余式 构建k个n阶矩阵(5.2)用户对k个n阶矩阵Bk左乘矩阵A-1，右乘矩阵A，得到k个明文矩阵Mk≡A-1BkA(mod p)；
(6)云服务器对密文矩阵进行同态计算，实现步骤为：(6.1)云服务器根据用户需要，从k个密文矩阵Ck中选取明文矩阵M1、M2、M1′和M2′对应的密文矩阵C1、C2、C1′和C2′，其中明文矩阵M1、M2、M1′和M2′相同或不同；
(6.2)云服务器使用同态计算公钥N，对密文矩阵C1和C2进行同态加法计算，同时对密文矩阵C1′和C2′进行同态乘法计算，得到同态加法密文矩阵C+和同态乘法密文矩阵C×；
(7)用户对同态加法密文矩阵C+和同态乘法密文矩阵C×进行解密，得到对应明文矩阵M+和M×，实现步骤为：(7.1)用户通过同态加法密文矩阵 和同余式(c+)ij≡(b+)ij(mod p)，构建n阶矩阵
(7.2)用户通过同态乘法密文矩阵 和同余式(c×)ij≡(b×)ij(mod p)，构建n阶矩阵
(7.3)用户对n阶矩阵B+左乘矩阵A-1，右乘矩阵A，得到明文矩阵M+≡A-1B+A(mod p)，其中M+＝M1+M2；
(7.4)用户对n阶矩阵B×左乘矩阵A ，右乘矩阵A，得到明文矩阵M×≡A B×A(mod p)，其中M×＝M1′×M2′。
2.根据权利要求1所述的基于明文相似矩阵的对称全同态加密方法，其特征在于，步骤(6.2)中所述的云服务器使用同态计算公钥N，对密文矩阵C1和C2进行同态加法计算，同时对密文矩阵C1′和C2′进行同态乘法计算，计算公式分别为：对密文矩阵C1和C2进行同态加法计算的公式为：C+≡C1+C2(mod N)，
对密文矩阵C1′和C2′进行同态乘法计算的公式为：C×≡C1′×C2′(mod N)。
基于明文相似矩阵的对称全同态加密方法
技术领域
[0001] 本发明属于信息安全领域，涉及一种对称的全同态加密方法，具体涉及一种基于明文相似矩阵的对称全同态加密方法，可应用于云计算、大数据环境等对重要数据的全程密态保护，在对密文数据不进行解密的情况下完成对明文数据的计算。
背景技术
[0002] 随着互联网的发展，尤其是云计算概念的诞生，人们在加密数据搜索与处理等方面的需求日益增加。但是对于大数据的处理，用户就必须委托给第三方(云)来进行操作；用户在云端存储的数据可能包含一些敏感信息，所以在将数据存储到云端之前必须对数据进行加密保护；然而，明文数据一旦被加密，明文数据结构会发生变化，一些在明文数据上的操作将不再适用于密文数据。
[0003] 同态加密是一种新型的加密算法，对密文数据进行同态运算得到一个输出，将这一输出进行解密，其结果与用同一方法处理未加密的原始数据得到的输出结果相等，这种加密方法的优点是可以直接对密文数据进行计算，而不用等到解密后在进行处理，但缺陷是只能对密文数据进行有限次的同态运算或只能对密文数据进行单一代数运算(加或乘)。然而，单运算同态密码和浅同态密码同态计算功能受限，例如，Paillier密码体制、Damgard-Jurik密码体制等只能满足加法同态；RSA密码体制、ElGamal等只能满足乘法同态；BGN密码体制满足任意次加法同态和一次的乘法同态。
[0004] 全同态加密是在同态加密的基础上提出来的，全同态加密可以根据用户要求，允许第三方对密文数据进行任意计算，且对密文数据的计算结果进行解密即得到对应明文数据的计算结果，而在整个处理过程中无需对数据进行解密；其意义在于，从根本上解决重要数据的隐私保护和密文数据的计算相互冲突的瓶颈问题。全同态加密方法的基本步骤包括：参数生成，密钥选取，对明文消息加密，对密文消息解密，对密文消息进行同态计算。
[0005] 2009年Gentry设计了第一个基于格的全同态密码体制，实现了历史性的跨越。然而，一些全同态密码方案由于加密方案中密文存在噪声，所以当密文计算到一定程度，其噪声将超过上限，对这样的密文解密将可能失败。还有一些对称全同态密码方案，由于过程繁琐导致效率低下。例如，2017年Khalil Hariss、Hassan Noura和Abed Ellatif Samhat在论文“Fully Enhanced Homomorphic Encryption algorithm of MORE approach for real world applications”(Journal of Information Security andApplications.2017.Pages 233-242)中提出了一种改进的用于随机化和加密的矩阵运算的对称全同态加密算法。该方法的具体步骤是：两个用户协商一对秘密参数：一个私密钥、一个初始向量IV，使用哈希算法生成64-bits的动态密钥DK；使用动态置换密钥DKp生成一个置换盒π＝[pi]1≤i≤N，将其用于输入的明文消息；基于动态扩散密钥DKd和流密码算法(例如：RC4),两个用户共享一个秘密序列s，使用重构函数将先前的秘密序列s变换为w个尺寸为 的小矩阵，使用矩阵生成方程构建包含w个矩阵和其逆矩阵的矩阵密钥库；创建另一个置换盒Δ＝[δi]1≤i≤H，Δ的尺寸为H，使用动态选择密钥DKs选择一个名为k的密钥块，k经过置换盒Δ的作用，生成序列 从序列 中选出下标为δk的密钥，这样就
生成了一对密钥矩阵 这一对密钥矩阵就作为相应矩阵块的动态加密密钥；使用动态选择算法从矩阵密钥库中选择动态加密密钥，用动态的选择密钥对明文消息进行分块加密；解密过程是加密过程的一个逆过程，用户拥有动态密钥DK和初始向量IV就可以生成一切秘密参数，用MORE算法进行解密。该方案可以实现对称的全同态加密，但是不足之处在于密钥选取过程繁琐，在加密过程中需要对明文消息进行分块加密，导致效率低。
发明内容
[0006] 本发明的目的是针对上述现有技术的缺陷，提出了一种基于明文相似矩阵的对称全同态加密方法，用于解决现有对称全同态加密效率低的技术问题。
[0007] 为实现上述目的，本发明采取的技术方案，包括如下步骤：
[0008] (1)参数生成：用户根据安全要求随机生成两个长度相等的大素数p和q；
[0009] (2)用户构建剩余类环和一般线性群，实现步骤为:
[0010] (2.1)用户构建关于大素数p的剩余类环 和关于大素数q的剩余类环
[0011] (2.2)用户利用剩余类环 构建模p意义下的所有n阶可逆矩阵构成的一般线性群 同时利用剩余类环 构建模q意义下的所有n阶可逆矩阵构成的一般线性群
[0012] (3)用户计算同态计算公钥和对称密钥，实现步骤为：
[0013] (3.1)用户计算RSA模数N，N＝pq，并将RSA模数N作为同态计算公钥；
[0014] (3.2)用户在一般线性群 中随机选取一个n阶可逆矩阵A，并求取可逆矩阵A模p意义下的逆矩阵A-1，A、A-1和p构成对称密钥sk，sk＝(A,A-1,p)；
[0015] (4)用户对明文矩阵进行加密，实现步骤为：
[0016] (4.1)用户根据需要在一般线性群 中选取k个n阶可逆的明文矩阵并分别计算k个明文矩阵Mk的相似矩阵 其中，k≥2，
是明文矩阵Mk中第i行第j列的元素， 是相似矩阵Bk中第i行第j列的元素；
[0017] (4.2)用户在一般线性群 中随机选取k个n阶矩阵 并根据中国剩余定理计算k个密文矩阵
[0018]
[0019] 其中， 是矩阵Dk中第i行第j列的元素， 是密文矩阵Ck中第i行第j列的元素且对任意的i,j＝1,2,…,n满足同余式方程组：
[0020]
[0021] (5)用户对密文矩阵进行解密，实现步骤为：
[0022] (5.1)用户通过同余式 构建k个n阶矩阵
[0023] (5.2)用户对k个n阶矩阵Bk左乘矩阵A-1，右乘矩阵A，得到k个明文矩阵Mk≡A-1BkA(modp)；
[0024] (6)云服务器对密文矩阵进行同态计算，实现步骤为：
[0025] (6.1)云服务器根据用户需要，从k个密文矩阵Ck中选取明文矩阵M1、M2、M1′和M2′对应的密文矩阵C1、C2、C1′和C2′，其中明文矩阵M1、M2、M1′和M2′相同或不同；
[0026] (6.2)云服务器使用同态计算公钥N，对密文矩阵C1和C2进行同态加法计算，同时对密文矩阵C1′和C2′进行同态乘法计算，得到同态加法密文矩阵C+和同态乘法密文矩阵C×；
[0027] (7)用户对同态加法密文矩阵C+和同态乘法密文矩阵C×进行解密，得到对应明文矩阵M+和M×。
[0028] 本发明与现有技术相比，具有以下优点：
[0029] 1.本发明利用随机生成的大素数构建剩余类环，根据剩余类环构建一般线性群，在构建的一般线性群中进行密钥选取，并且使用相似矩阵对明文矩阵进行加密，过程简单，与现有技术相比，提高了效率。
[0030] 2.本发明在加密过程中对明文矩阵的相似矩阵进行加密，在解密时，通过乘矩阵的逆矩阵，结果满足矩阵乘自身的逆矩阵等于单位矩阵，从而计算过程不会引入噪声，支持对密文矩阵进行任意次的同态计算，满足用户的需求。
[0031] 3.本发明在加密过程中对明文矩阵的相似矩阵进行加密，使得明文矩阵被其相似矩阵随机隐藏，提高了加密方法的安全性。
附图说明
[0032] 图1为本发明的实现流程图。
具体实施方式
[0033] 以下结合附图和具体实施例，对本发明作进一步详细说明。
[0034] 步骤1)参数生成：用户根据安全要求随机生成两个长度相等的大素数p和q，取p和q的长度均为1024比特；
[0035] 步骤2)用户构建剩余类环和一般线性群，实现步骤为:
[0036] 步骤2.1)用户构建关于大素数p的剩余类环 和关于大素数q的剩余类环
[0037] 步骤2.2)用户利用剩余类环 构建模p意义下的所有n阶可逆矩阵构成的一般线性群 利用剩余类环 构建模q意义下的所有n阶可逆矩阵构成的一般线性群
[0038] 步骤3)用户计算同态计算公钥和对称密钥，实现步骤为：
[0039] 步骤3.1)用户计算RSA模数N，N＝pq，并将RSA模数N作为同态计算公钥；
[0040] 步骤3.2)用户在一般线性群 中随机选取一个n阶可逆矩阵A，并求取可逆矩阵A模p意义下的逆矩阵A-1，A、A-1和p构成对称密钥sk，sk＝(A,A-1,p)；
[0041] 步骤4)用户对明文矩阵进行加密，实现步骤为：
[0042] 步骤4.1)用户根据需要在一般线性群 中选取k个n阶可逆的明文矩阵并分别计算k个明文矩阵Mk的相似矩阵 其中，k≥
2， 是明文矩阵Mk中第i行第j列的元素， 是相似矩阵Bk中第i行第j列的元素；
[0043] 步骤4.2)用户在一般线性群 中随机选取k个n阶矩阵 并根据中国剩余定理计算k个密文矩阵
[0044]
[0045] 其中， 是矩阵Dk中第i行第j列的元素， 是密文矩阵Ck中第i行第j列的元素且对任意的i,j＝1,2,…,n满足同余式方程组：
[0046]
[0047] 步骤5)用户对密文矩阵进行解密，实现步骤为：
[0048] 步骤5.1)用户通过同余式 构建k个n阶矩阵
[0049] 步骤5.2)用户对k个n阶矩阵Bk左乘矩阵A-1，右乘矩阵A，得到k个明文矩阵Mk≡A-1
BkA(modp)；
[0050] 步骤6)云服务器对密文矩阵进行同态计算，实现步骤为：
[0051] 步骤6.1)云服务器根据用户需要，从k个密文矩阵Ck中选取明文矩阵M1、M2、M1′和M2′对应的密文矩阵C1、C2、C1′和C2′，其中明文矩阵M1、M2、M1′和M2′相同或不同；
[0052] 步骤6.2)云服务器使用同态计算公钥N，对密文矩阵C1和C2进行同态加法计算，同时对密文矩阵C1′和C2′进行同态乘法计算，得到同态加法密文矩阵C+和同态乘法密文矩阵C×，对密文矩阵C1和C2进行同态加法计算的公式为：
[0053] C+≡C1+C2(modN)，
[0054] 对密文矩阵C1′和C2′进行同态乘法计算的公式为：
[0055] C×≡C1′×C2′(modN)；
[0056] 步骤7)用户对同态加法密文矩阵C+和同态乘法密文矩阵C×进行解密，得到对应明文矩阵M+和M×，实现步骤为：
[0057] 步骤7.1)用户通过同态加法密文矩阵 和同余式(c+)ij≡(b+)ij(mod p)，构建n阶矩阵
[0058] 步骤7.2)用户通过同态乘法密文矩阵 和同余式(c×)ij≡(b×)ij(mod p)，构建n阶矩阵
[0059] 步骤7.3)用户对n阶矩阵B+左乘矩阵A-1，右乘矩阵A，得到明文矩阵M+≡A-1B+A(modp)，其中M+＝M1+M2；
[0060] 步骤7.4)用户对n阶矩阵B×左乘矩阵A-1，右乘矩阵A，得到明文矩阵M×≡A-1B×A(modp)，其中M×＝M1′×M2′。
