一种量子安全多方计算系统

1.本发明涉及量子密码学技术领域，尤其涉及一种量子安全多方计算系统。


背景技术：

2.随着信息技术的快速发展，对信息处理任务的需求日益增加，数据操作变得不可或缺。然而，数据通常包含大量机密信息。因此，如何在不泄露私人信息的情况下处理这些数据成为一个重要问题。在此背景下，出现了安全多方计算。安全多方计算是一种密码学原语，它允许分布式参与者协同计算任意函数，同时保护参与者输入数据的隐私。目前，安全多方计算主要由公钥密码系统实现，包括rsa方案和elgamal，它们都基于数学计算的复杂性。然而，随着量子计算的发展，shor算法、grover算法和其他量子算法使得在短时间内破解上述经典密码系统成为可能，经典信道下的smc协议的安全受到挑战。
3.现阶段引入量子计算执行安全多方计算，比较成熟的有两方求和或求积协议，而缺少多方求和或求积的协议。更进一步的，量子安全多方计算几何是量子安全多方计算的一个分支，主要包含与几何相关的问题，包括两点之间的距离、两个多边形的交点，以及确定几何对象(如安全点、线和平面)的相对位置。目前，量子安全多方计算几何的研究主要集中在量子安全双方距离计算上，并且几乎所有关于量子安全双方距离计算的研究都基于量子私有查询，这需要双方共同建立一个数据库，然后通过重复查询数据库来获得距离。然而，该方法在准确性和效率方面存在缺点。首先，数据库查询的结果是近似值而不是精确值，但在某些情况下，这种不精确是需要被严格避免的。其次，当实现高维量子安全双方距离计算协议时，数据库的大小呈指数级增长，导致更多的资源消耗和更低的效率。
4.因此，亟需一种新的量子安全多方计算系统，满足多种安全多方计算要求。


技术实现要素：

5.鉴于此，本发明实施例提供了一种量子安全多方计算系统，以消除或改善现有技术中存在的一个或更多个缺陷，解决现有技术缺少量子多方求和求积、距离计算、多面体体积计算和圆圆相交面积计算能力的问题。
6.本发明提供了一种量子安全多方计算系统，所述系统包括多方参与者，所述参与者执行以下操作：
7.执行量子多方求和操作：对于第一设定数量个第一参与者，两两随机配对，按照基于shamir秘密共享协议的第一预设量子双方求和方法，由多个第三方协助对各第一参与者所持有的第一秘密执行量子安全双方求和，各第三方保留本地的所有测量结果并计算总和后广播，对所有广播的结果计算总和得到各第一参与者所持有的所述第一秘密的和；其中，所述第一设定数量为偶数；
8.执行量子多方求积操作：对于第二设定数量个第二参与者，两两随机配对，按照基于shamir秘密共享协议的第二预设量子双方求积方法，由多个第三方协助对各第二参与者所持有的第二秘密执行量子安全双方求积，各第三方保留本地的所有测量结果并计算总和
后广播，对所有广播的结果计算总和得到各第二参与者所持有的所述第二秘密的积；其中，所述第二设定数量为偶数；
9.执行量子安全双方距离计算：对随机的第三参与者所持有的第一横坐标和第一纵坐标，以及随机的第四参与者所持有的第二横坐标和第二纵坐标；所述第三参与者计算所述第一横坐标和所述第一纵坐标的平方和得到第一数值秘密，将所述第一横坐标乘以负二得到第二数值秘密，将所述第一纵坐标乘以负二得到第三数值秘密；所述第四参与者计算所述第二横坐标和所述第二纵坐标的平方和得到第四数值秘密，将所述第二横坐标作为第五数值秘密，将所述第二纵坐标作为第六数值秘密；按照基于shamir秘密共享协议的所述第一预设量子双方求和方法，由多个第三方协助对所述第一数值秘密和所述第四数值秘密执行量子双方求和，并由各第三方保留获得的第一测量值；按照基于shamir秘密共享协议的所述第二预设量子双方求积方法，由多个第三方协助对所述第二数值秘密和所述第五数值秘密执行量子双方求积，各第三方保留获得的第二测量值；按照基于shamir秘密共享协议的所述第二预设量子双方求积方法，由多个第三方协助对所述第三数值秘密和所述第六数值秘密执行量子双方求积，各第三方保留获得的第三测量值；各第三方将本地的所述第一测量值、所述第二测量值和所述第三测量值求和并广播，对所有广播的结果求和得到所述第三参与者与所述第四参与者之间的平方距离，对所述平方距离开根得到所述第三参与者与所述第四参与者之间的距离。
10.在一些实施例中，所述系统执行量子多方求和操作过程中，当所述第一设定数量为奇数时，添加一个第一虚拟参与者参与各第一参与者的两两随机配对，所述第一虚拟参与者持有的所述第一秘密设为0。
11.在一些实施例中，所述系统执行量子多方求积操作过程中，当所述第二设定数量为奇数时，添加一个第二虚拟参与者参与各第二参与者的两两随机配对，所述第二虚拟参与者持有的所述第二秘密设为1。
12.在一些实施例中，基于shamir秘密共享协议的所述第一预设量子双方求和方法，包括：
13.由持有秘密a的一方随机构建多项式f(y)＝a+c1y+c2y2+
…
+c
t-1yt-1
mod d，并基于shamir秘密共享协议计算经典份额f(yi)，i＝1，2
…
，n；由持有秘密b的一方随机构建多项式并基于shamir秘密共享协议计算经典份额g(yi)，i＝1，2
…
，n；将f(yi)与g(yi)在一组第三方tp1，tp2，
…
tpn之间进行分配；
14.对于参与恢复的第三方tpk，k＝1，2，
…
，t，计算h(yk)＝f(yk)+g(yk)mod d，保守秘密h(yk)，并计算h(yk)的份额的投影
15.由tp1制备t个单粒子|0》1，|0》2，
…
，|0》
t
，对第一个粒子|0》1施加量子离散傅里叶变换操作，得到的结果为由tp1通过在粒子1和粒子k，k＝2，3，
…
，t之间分别执行t-1次纠缠门操作，生成纠缠态并将粒子k分发给tpk，k＝1，2，
…
，t；
16.每个tpk对分发到的粒子k做量子离散傅里叶变换操作，再执行广义泡利算子
量子态|φ2》变为
17.每个tpk测量自己的粒子|mk+ak》并宣布mk+ak，各tpk将测量结果相加，计算总和得到秘密a和秘密b的和为：
18.在一些实施例中，基于shamir秘密共享协议的第二预设量子双方求积方法，包括：
19.由持有秘密a的一方随机构建多项式f(y)＝a+c1y+c2y2+
…
+c
t-1yt-1
modd，并基于shamir秘密共享协议计算经典份额f(yi)，i＝1，2
…
，n；；c和d为参数，y为变量；由持有秘密b的一方随机构建多项式并基于shamir秘密共享协议计算经典份额g(yi)，i＝1，2
…
，n；c和d为参数，y为变量；将f(yi)与g(yi)在一组第三方tp1，tp2，
…
tpn之间进行分配；
20.对于每个tpi，i＝1，2，
…
，n，计算h(yi)＝f(yi)
×
g(yi)mod d并使用一个随机多项式zi(x)＝h
′
(yi)+β1x+β2x2+
…
+β
t-1
x
t-1
mod d在n个tp之间计算共享；使用范德蒙德矩阵表示t，对于参与恢复的第三方tpk，每个tpk仅知晓其对应的tk，并计算份额投影为：
21.由tp1制备t个单粒子|0》1，|0》2，
…
，|0》
t
，对第一个粒子|0》1施加量子离散傅里叶变换操作，得到的结果为由tp1通过在粒子1和粒子k，k＝2,3，
…
，t之间分别执行t-1次纠缠门操作，生成纠缠态并将粒子k分发给tpk，k＝1，2，
…
，t；
22.每个tpk对分发到的粒子k做量子离散傅里叶变换操作，再执行广义泡利算子量子态|φ2》变为
23.每个tpk测量自己的粒子|mk+bk》并宣布mk+bk，各tpk将测量结果相加，计算总和得到秘密a和秘密b的积为：
24.在一些实施例中，所述系统还用于执行量子安全四面体体积计算，包括：
25.对于单个四面体由4个参与者坐标作为顶点，记为：a0(x0，y0，z0)，a1(x1，y1，z1)，a2(x2，y2，z2)和a3(x3，y3，z3)；
26.采用三阶行列式计算单个四面体的体积：
[0027][0028]
记行列式：
[0029][0030][0031][0032][0033]
对于行列式d0,d1,d2和d3中的每个项，由每个项所涉及坐标对应的三个参与者执行所述量子多方求积操作进行量子安全三方求积，得到各项的值，对各项求和得到d0,d1,d2和d3的和，并求解v(a0a1a2a3)。
[0034]
在一些实施例中，所述系统还用于量子安全多方多面体体积计算，由多个参与者在不暴露私有坐标的情况下计算由其坐标组成的多面体体积，包括：
[0035]
将所述多面体拆分为多个四面体，每个四面体由四个参与者的坐标作为顶点；
[0036]
每个四面体对应的四个参与者执行量子安全四面体体积计算，并累加得到所述多面体体积。
[0037]
在一些实施例中，所述系统还用于执行量子安全圆-圆相交面积计算，由两个参与者在不暴露私有圆坐标的情况下，计算两个私有圆的相交面积，包括：
[0038]
由第五参与者在其私有圆内随机选择第三设定数量个标记点；
[0039]
由所述第五参与者与第六参与者执行所述量子安全双方距离计算，以计算各标记点到所述第六参与者私有圆圆心的距离，并计算处于所述第六参与者私有圆范围内的标记
点数量；
[0040]
则所述第五参与者与所述第六参与者之间私有圆的相交面积s
∩
计算式为：
[0041][0042]
其中，k表示处于所述第六参与者私有圆范围内的标记点数量，n表示所述第三设定数量，sc表示所述第五参与者的私有圆的面积。
[0043]
在一些实施例中，由两个参与者在不暴露私有圆坐标的情况下，计算两个私有圆的相交面积，包括：
[0044]
由所述第五参与者在其私有圆内随机选择n个标记点{p1，p2，
…
，pn}；
[0045]
由所述第六参与者构造一组随机数rn(n＝1，2
…
，n)；
[0046]
由所述第五参与者与第六参与者执行所述量子安全双方距离计算，以计算各标记点pi到所述第六参与者私有圆的圆心o1的平方距离参与协助的每个第三方tp保留关于份额不公布，i＝1，2，
…
，n；
[0047]
基于shamir秘密共享协议的第二预设量子双方求积方法，所述第五参与者与所述第六参与者在rn(n＝1，2，
…
，n)与常数1之间进行n次量子双方求积，每个第三方tp获得r1，r2，
…
，rn的份额；
[0048]
tpi将和ri的份额对应相加，以获得的份额；由tpi意外的其他tp公布份额测量结果，则tp1获得
[0049]
由所述第六参与者与tp1在第三方tp0的协助下适用量子百万富翁协议对r
12
+r1和r
12
+r2和
…
，r
12
+rn和分别进行比较，并由第三方tp0计算所述第六参与者拥有较大元素的个数，记为k，k为处于所述第六参与者私有圆范围内的标记点数量；
[0050]
由第三方tp0将k发送至所述第五参与者，并由所述第五参与者计算所述第五参与者与所述第六参与者之间私有圆的相交面积s
∩
计算式为：
[0051][0052]
其中，k表示处于所述第六参与者私有圆范围内的标记点数量，n表示所述第三设定数量，sc表示所述第五参与者的私有圆的面积。
[0053]
在一些实施例中，由所述第五参与者计算所述第五参与者与所述第六参与者之间私有圆的相交面积之后，还包括：
[0054]
由所述第五参与者将所述相交面积发送至所述第六参与者。
[0055]
本发明的有益效果至少是：
[0056]
本发明所述量子安全多方计算系统，基于对shamir秘密共享协议两方求和求积方案的延伸，借助几个半诚实的第三方，为多方参与者提供了量子多方求和求积的执行方案，以及量子安全双方距离计算方案，消除了对数据库的依赖，解决中间值泄露的问题，在准确性和效率方面都有所提高，极大减少了多方量子安全计算所需的资源。
[0057]
进一步的，所述系统提供对以多方参与者坐标作为顶点的多面体求体积的计算方案，通过将多面体拆解为多个四面体，引入代数对四面体体积的计算进行表达，将对四面体求体积的问题转化为多方求积和求和问题，实现了高拓展能力的量子安全多面体体积计
算。
[0058]
进一步的，所述系统提供对两方参与者私有圆计算圆圆相交面积的方案，通过执行量子安全双方距离计算，获取一方参与者私有圆内多个随机标记点到另一方参与者圆心的距离，并判断是否在另一方参与者私有圆范围内。在标记点数量足够多的情况下，通过处于另一方参与者私有圆范围内的标记点的数量与标记点总数的比，估计圆圆相交面积与另一方参与者私有圆面积的比，实现对双方私有圆相交面积的计算。
[0059]
本发明的附加优点、目的，以及特征将在下面的描述中将部分地加以阐述，且将对于本领域普通技术人员在研究下文后部分地变得明显，或者可以根据本发明的实践而获知。本发明的目的和其它优点可以通过在说明书以及附图中具体指出的结构实现到并获得。
[0060]
本领域技术人员将会理解的是，能够用本发明实现的目的和优点不限于以上具体所述，并且根据以下详细说明将更清楚地理解本发明能够实现的上述和其他目的。
附图说明
[0061]
此处所说明的附图用来提供对本发明的进一步理解，构成本技术的一部分，并不构成对本发明的限定。在附图中：
[0062]
图1为本发明一实施例所述量子安全多方计算系统的结构示意图。
[0063]
图2为本发明另一实施例所述量子安全多方计算系统执行量子安全圆-圆相交面积的示意图。
具体实施方式
[0064]
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
[0065]
在此，还需要说明的是，为了避免因不必要的细节而模糊了本发明，在附图中仅仅示出了与根据本发明的方案密切相关的结构和/或处理步骤，而省略了与本发明关系不大的其他细节。
[0066]
应该强调，术语“包括/包含”在本文使用时指特征、要素、步骤或组件的存在，但并不排除一个或更多个其它特征、要素、步骤或组件的存在或附加。
[0067]
在此，还需要说明的是，如果没有特殊说明，术语“连接”在本文不仅可以指直接连接，也可以表示存在中间物的间接连接。
[0068]
在下文中，将参考附图描述本发明的实施例。在附图中，相同的附图标记代表相同或类似的部件，或者相同或类似的步骤。
[0069]
量子安全多方计算是将量子力学的一些基本原理应用于安全多方计算的协议设计，使得在完成安全多方计算功能的同时保证协议能够抵抗量子计算的攻击且具有更优越的安全性能。这种技术用于解决一些隐私保护的问题，例如在不泄露参与方的隐私输入和输出的前提下让分布式参与方合作计算任意函数。现有技术中存在基于shamir的秘密共享提出了两方求和与求积协议，但是只限于两方参与者主体，无法满足多方参与者同时进行量子安全多方求和或求积的需求。同时，现有技术中，缺少关于距离计算、多面体面积计算
和圆圆相交面积计算的协议。本技术提供一种量子安全多方计算系统，在量子安全两方求和求积计算协议的基础上进行拓展，提供了量子安全多方求和求积计算协议以及量子安全多方几何计算的协议。
[0070]
具体的，如图1所示，本技术提供一种量子安全多方计算系统，所述系统包括多方参与者，本技术中所述的参与者是秘密的持有者，各参与者采用能够存储和执行计算机程序的电子设备作为客户端进行数据通信。本实施例中，所述参与者执行以下量子多方求和、量子多方求积以及量子安全双方距离计算操作，具体说明如下：
[0071]
1)由多方参与者共同执行量子多方求和操作：对于第一设定数量个第一参与者，两两随机配对，按照基于shamir秘密共享协议的第一预设量子双方求和方法，由多个第三方协助对各第一参与者所持有的第一秘密执行量子安全双方求和，各第三方保留本地的所有测量结果并计算总和后广播，对所有广播的结果计算总和得到各第一参与者所持有的第一秘密的和；其中，第一设定数量为偶数。
[0072]
在一些实施例中，当第一设定数量为奇数时，添加一个第一虚拟参与者参与各第一参与者的两两随机配对，第一虚拟参与者持有的第一秘密设为0。这样可以保证求和过程不会影响最终的结果。
[0073]
本实施例提供的量子多方求和方案，是基于量子安全双方求和进行拓展得到的。
[0074]
其中，通过对多方参与者进行两两配对，分别执行基于shamir秘密共享协议的第一预设量子双方求和方法，包括步骤s101～s105：
[0075]
步骤s101：由持有秘密a的一方随机构建多项式f(y)＝a+c1y+c2y2+
…
+c
t-1yt-1
mod d，并基于shamir秘密共享协议计算经典份额f(yi)，i＝1，2
…
，n；c和d为参数，y为变量；由持有秘密b的一方随机构建多项式并基于shamir秘密共享协议计算经典份额g(yi)，i＝1，2
…
，n；c和d为参数，y为变量；将f(yi)与g(yi)在一组第三方tp1，tp2，
…
tpn之间进行分配；
[0076]
步骤s102：对于参与恢复的第三方tpk，k＝1，2，
…
，t，计算h(yk)＝f(yk)+g(yk)mod d，保守秘密h(yk)，并计算h(yk)的份额的投影
[0077]
步骤s103：由tp1制备t个单粒子|0》1，|0》2，
…
，|0》
t
，对第一个粒子|0》1施加量子离散傅里叶变换操作，得到的结果为由tp1通过在粒子1和粒子k，k＝2，3，
…
，t之间分别执行t-1次纠缠门操作，生成纠缠态并将粒子k分发给tpk，k＝1，2，
…
，t。
[0078]
步骤s104：每个tpk对分发到的粒子k做量子离散傅里叶变换操作，再执行广义泡利算子量子态|φ2》变为
[0079]
步骤s105：每个tpk测量自己的粒子|mk+ak》并宣布mk+ak，各tpk将测量结果相加，计算总和得到秘密a和秘密b的和为：
[0080]
结合步骤s101～s105，多方参与者执行量子安全求和的步骤包括步骤s201～s202：
[0081]
步骤s201：要求参与量子安全多方求和的参与者的第一设定数量为偶数，参与者两两随机配对，每对参与者执行如步骤s101～s105的量子安全双方求和协议，每个协助的第三方tp在协助执行后，保留自己的测量结果暂不公布。
[0082]
步骤s202：在所有参与者两两执行了双方量子求和协议后，由参与恢复的第三方tpk，k＝1，2，
…
，t，计算其测量结果的总和并广播；由所有参与恢复的tp将其结果相加，则可以得到多方参与者所持有秘密的和：
[0083]
在步骤s201～s202中，如果参与者的第一设定数量为奇数，则添加一个第一虚拟参与者，第一虚拟参与者持有的第一秘密设为0，以保证能够完成参与者的两两配对，并且不影响最终的求和结果。
[0084]
2)执行量子多方求积操作：对于第二设定数量个第二参与者，两两随机配对，按照基于shamir秘密共享协议的第二预设量子双方求积方法，由多个第三方协助对各第二参与者所持有的第二秘密执行量子安全双方求积，各第三方保留本地的所有测量结果并计算总和后广播，对所有广播的结果计算总和得到各第二参与者所持有的所述第二秘密的积；其中，第二设定数量为偶数。
[0085]
在一些实施例中，当第二设定数量为奇数时，添加一个第二虚拟参与者参与各第二参与者的两两随机配对，第二虚拟参与者持有的第二秘密设为1。
[0086]
其中，通过对多方参与者进行两两配对，基于shamir秘密共享协议的第二预设量子双方求积方法，包括步骤s301～s305：
[0087]
步骤s301：由持有秘密a的一方随机构建多项式f(y)＝a+c1y+c2y2+
…
+c
t-1yt-1
mod d，并基于shamir秘密共享协议计算经典份额f(yi)，i＝1，2
…
，n；由持有秘密b的一方随机构建多项式并基于shamir秘密共享协议计算经典份额g(yi)，i＝1，2
…
，n；将f(yi)与g(yi)在一组第三方tp1，tp2，
…
tpn之间进行分配。
[0088]
步骤s302：对于每个tpi，i＝1，2，
…
，n，计算h
′
(yi)＝f(yi)
×
g(yi)mod d并使用一个随机多项式zi(x)＝h
′
(yi)+β1x+β2x2+
…
+β
t-1
x
t-1
mod d在n个tp之间计算共享；使用范德蒙德矩阵表示t，对于参与恢复的第三方tpk，每个tpk仅知晓其对应的tk，并计算份额投影为：
[0089]
步骤s303：由tp1制备t个单粒子|0》1，|0》2，
…
，|0》
t
，对第一个粒子|0》1施加量子离散傅里叶变换操作，得到的结果为由tp1通过在粒子1和粒子k，k＝2，3，
…
，t之间分别执行t-1次纠缠门操作，生成纠缠态并将粒子k分发给tpk，k＝1，2，
…
，t。
[0090]
步骤s304：每个tpk对分发到的粒子k做量子离散傅里叶变换操作，再执行广义泡利算子量子态|φ2》变为
[0091]
步骤s305：每个tpk测量自己的粒子|mk+bk》并宣布mk+bk，各tpk将测量结果相加，计算总和得到秘密a和秘密b的积为：
[0092]
结合步骤s301～s305，多方参与者执行量子安全求积的步骤包括步骤s401～s402：
[0093]
步骤s401：要求参与量子安全多方求积的参与者的第二设定数量为偶数，参与者两两随机配对，每对参与者执行如步骤s301～s305的量子安全双方求积协议，每个协助的第三方tp在协助执行后，保留自己的测量结果暂不公布。
[0094]
步骤s402：在所有参与者两两执行了双方量子求积协议后，由参与恢复的第三方tpk，k＝1，2，
…
，t，计算测量结果的总和并广播。由所有参与恢复的tp将其结果相加，则可以得到多方参与者所持有秘密的积为：
[0095]
步骤s401～s402中，如果参与者的第二设定数量为奇数，则添加一个第二虚拟参与者参与各第二参与者的两两随机配对，第二虚拟参与者持有的第二秘密设为1。使得保证能够完成参与者的两两配对的基础上，不影响最终的求积结果。
[0096]
3)执行量子安全双方距离计算：
[0097]
对随机的第三参与者所持有的第一横坐标和第一纵坐标，以及随机的第四参与者所持有的第二横坐标和第二纵坐标。
[0098]
第三参与者计算第一横坐标和第一纵坐标的平方和得到第一数值秘密，将第一横坐标乘以负二得到第二数值秘密，将第一纵坐标乘以负二得到第三数值秘密。
[0099]
第四参与者计算第二横坐标和第二纵坐标的平方和得到第四数值秘密，将第二横坐标作为第五数值秘密，将第二纵坐标作为第六数值秘密；按照基于shamir秘密共享协议的第一预设量子双方求和方法，由多个第三方协助对第一数值秘密和第四数值秘密执行量子双方求和，并由各第三方保留获得的第一测量值；按照基于shamir秘密共享协议的第二预设量子双方求积方法，由多个第三方协助对第二数值秘密和第五数值秘密执行量子双方求积，各第三方保留获得的第二测量值；按照基于shamir秘密共享协议的第二预设量子双方求积方法，由多个第三方协助对第三数值秘密和第六数值秘密执行量子双方求积，各第三方保留获得的第三测量值；各第三方将本地的第一测量值、第二测量值和第三测量值求和并广播，对所有广播的结果求和得到第三参与者与第四参与者之间的平方距离，对平方距离开根得到第三参与者与第四参与者之间的距离。
[0100]
具体的，对于输入，两方参与者a和b所持有的坐标秘密分别为a(x1，y1)和b(x2，y2)，输出的距离可以表达为：为了便于计算，可以优先计算平方距离：|ab|2＝x
12
+x
22-2x1x2+y
12
+y
22-2y1y2。参与者a和b可以基于步骤s101～s105的量子两方求和协议以及步骤s301～s305的量子双方求积协议计算获得|ab|2。具体的，步骤如下s501～s505：
[0101]
步骤501：参与者a首先独自计算(x
12
+y
12
)，类似地，b计算(x
22
+y
22
)。
[0102]
步骤502：参与者a和b分别用(x
12
+y
12
)和(x
22
+y
22
)执行量子双方求和协议。每tpk，k＝1，2，
…
，t都获得测量值m
k1
+a
k1
。
[0103]
步骤503：参与者a和b使用-2x1和x2执行量子双方求积协议，tpk获得m
k2
+a
k2
。
[0104]
步骤504：参与者a和b使用-2y1和y2执行量子双方求积协议，tpk获得m
k3
+a
k3
。
[0105]
步骤505：所有tpk将三个测量结果相加并公布，然后计算其结果的总和：
得到|ab|2。
[0106]
在一些实施例中，所述系统还用于执行量子安全四面体体积计算，包括步骤s601～s604：
[0107]
步骤s601：对于单个四面体由4个参与者坐标作为顶点，记为：a0(x0,y0,z0)，a1(x1,y1,z1)，a2(x2,y2,z2)和a3(x3,y3,z3)。
[0108]
步骤s602：采用三阶行列式计算单个四面体的体积：
[0109][0110]
步骤s603：记行列式：
[0111][0112][0113][0114][0115]
步骤s604：对于行列式d0,d1,d2和d3中的每个项，由每个项所涉及坐标对应的三个参与者执行步骤s401～s402中的量子多方求积操作进行量子安全三方求积，得到各项的值，对各项求和得到d0，d1，d2和d3的和，并求解v(a0a1a2a3)。
[0116]
其中，对于d0，可以由三名参与者使用6次三方求积分别计算x1y2z3，x3y1z2，x2y3z1，-x3y2z1，-x2y1z3，-x1y3z2，同理计算d1，d2和d3，所有tp发布他们的测量结果，然后相加所有测量结果得到d0+d1+d2+d3，再将结果除6可得v(a
0a
1a2a3)。
[0117]
在一些实施例中，所述系统还用于量子安全多方多面体体积计算，由多个参与者在不暴露私有坐标的情况下计算由其坐标组成的多面体体积，包括步骤s701～s702：
[0118]
步骤s701：将多面体拆分为多个四面体，每个四面体由四个参与者的坐标作为顶点。
[0119]
步骤s702：每个四面体对应的四个参与者执行量子安全四面体体积计算，并累加得到所多面体体积。
[0120]
在一些实施例中，所述系统还用于执行量子安全圆-圆相交面积计算，由两个参与者在不暴露私有圆坐标的情况下，计算两个私有圆的相交面积，包括步骤s801～s803：
[0121]
步骤s801：参照图2，由第五参与者在其私有圆内随机选择第三设定数量个标记点。
[0122]
步骤s802：由第五参与者与第六参与者执行量子安全双方距离计算，以计算各标记点到第六参与者私有圆圆心的距离，并计算处于第六参与者私有圆范围内的标记点数量。
[0123]
步骤s803：则第五参与者与第六参与者之间私有圆的相交面积s
∩
计算式为：
[0124][0125]
其中，k表示处于第六参与者私有圆范围内的标记点数量，n表示第三设定数量，sc表示第五参与者的私有圆的面积。随着n的增加，获得的结果将更接近实际值。
[0126]
在一些实施例中，由两个参与者在不暴露私有圆坐标的情况下，计算两个私有圆的相交面积，包括步骤s901～s907：
[0127]
其中，第五参与者的私有圆记为c2(o2，r2)，第六参与者的私有圆记为c1(o1，r1)。
[0128]
步骤s901：由第五参与者在其私有圆内随机选择n个标记点{p1，p2，
…
，pn}。
[0129]
步骤s902：由第六参与者构造一组随机数rn(n＝1，2，
…
，n)。
[0130]
步骤s903：由第五参与者与第六参与者执行量子安全双方距离计算，以计算各标记点pi到第六参与者私有圆的圆心o1的平方距离参与协助的每个第三方tp保留关于份额不公布，i＝1，2，
…
，n。
[0131]
步骤s904：基于shamir秘密共享协议的第二预设量子双方求积方法，第五参与者与第六参与者在rn(n＝1，2，
…
，n)与常数1之间进行n次量子双方求积，每个第三方tp获得r1，r2，
…
，rn的份额。
[0132]
步骤s905：tpi将和ri的份额对应相加，以获得的份额；由tpi意外的其他tp公布份额测量结果，则tp1获得
[0133]
步骤s906：由第六参与者与tp1在第三方tp0的协助下适用量子百万富翁协议对r
12
+r1和r
12
+r2和
…
，r
12
+rn和分别进行比较，并由第三方tp0计算第六参与者拥有较大元素的个数，记为k，k为处于第六参与者私有圆范围内的标记点数量。
[0134]
步骤s907：由第三方tp0将k发送至第五参与者，并由第五参与者计算第五参与者与第六参与者之间私有圆的相交面积s
∩
计算式为：
[0135]
[0136]
其中，k表示处于第六参与者私有圆范围内的标记点数量，n表示第三设定数量，sc表示第五参与者的私有圆的面积。
[0137]
在一些实施例中，由第五参与者计算第五参与者与第六参与者之间私有圆的相交面积之后，还包括：由第五参与者将相交面积发送至第六参与者。
[0138]
本技术中所述量子安全多方计算系统，在量子安全两方求和求积计算协议的基础上进行拓展，提供了量子安全多方求和求积计算协议以及量子安全多方几何计算的协议，通过引入多个第三方进行协助，能够省去对数据库的依赖，在准确性和效率方面都有所提高，当实现高维量子安全双方距离计算协议时，所需的资源也少了很多。就量子安全双方距离计算协议而言，现有技术主要基于量子私有查询的方式进行，估算值与数据库大小正相关，而本技术是基于求和求积协议执行的，通信复杂度更低，且无条件精准。
[0139]
综上所述，本发明所述量子安全多方计算系统，基于对shamir秘密共享协议两方求和求积方案的延伸，借助几个半诚实的第三方，为多方参与者提供了量子多方求和求积的执行方案，以及量子安全双方距离计算方案，消除了对数据库的依赖，解决中间值泄露的问题，在准确性和效率方面都有所提高，极大减少了多方量子安全计算所需的资源。
[0140]
进一步的，所述系统提供对以多方参与者坐标作为顶点的多面体求体积的计算方案，通过将多面体拆解为多个四面体，引入代数对四面体体积的计算进行表达，将对四面体求体积的问题转化为多方求积和求和问题，实现了高拓展能力的量子安全多面体体积计算。
[0141]
进一步的，所述系统提供对两方参与者私有圆计算圆圆相交面积的方案，通过执行量子安全双方距离计算，获取一方参与者私有圆内多个随机标记点到另一方参与者圆心的距离，并判断是否在另一方参与者私有圆范围内。在标记点数量足够多的情况下，通过处于另一方参与者私有圆范围内的标记点的数量与标记点总数的比，估计圆圆相交面积与另一方参与者私有圆面积的比，实现对双方私有圆相交面积的计算。
[0142]
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时以实现前述边缘计算服务器部署方法的步骤。该计算机可读存储介质可以是有形存储介质，诸如随机存储器(ram)、内存、只读存储器(rom)、电可编程rom、电可擦除可编程rom、寄存器、软盘、硬盘、可移动存储盘、cd-rom、或技术领域内所公知的任意其它形式的存储介质。
[0143]
本领域普通技术人员应该可以明白，结合本文中所公开的实施方式描述的各示例性的组成部分、系统和方法，能够以硬件、软件或者二者的结合来实现。具体究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。当以硬件方式实现时，其可以例如是电子电路、专用集成电路(asic)、适当的固件、插件、功能卡等等。当以软件方式实现时，本发明的元素是被用于执行所需任务的程序或者代码段。程序或者代码段可以存储在机器可读介质中，或者通过载波中携带的数据信号在传输介质或者通信链路上传送。
[0144]
需要明确的是，本发明并不局限于上文所描述并在图中示出的特定配置和处理。为了简明起见，这里省略了对已知方法的详细描述。在上述实施例中，描述和示出了若干具体的步骤作为示例。但是，本发明的方法过程并不限于所描述和示出的具体步骤，本领域的
技术人员可以在领会本发明的精神后，作出各种改变、修改和添加，或者改变步骤之间的顺序。
[0145]
本发明中，针对一个实施方式描述和/或例示的特征，可以在一个或更多个其它实施方式中以相同方式或以类似方式使用，和/或与其他实施方式的特征相结合或代替其他实施方式的特征。
[0146]
以上所述仅为本发明的优选实施例，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种量子安全多方计算系统，其特征在于，所述系统包括多方参与者，所述参与者执行以下操作：执行量子多方求和操作：对于第一设定数量个第一参与者，两两随机配对，按照基于shamir秘密共享协议的第一预设量子双方求和方法，由多个第三方协助对各第一参与者所持有的第一秘密执行量子安全双方求和，各第三方保留本地的所有测量结果并计算总和后广播，对所有广播的结果计算总和得到各第一参与者所持有的所述第一秘密的和；其中，所述第一设定数量为偶数；执行量子多方求积操作：对于第二设定数量个第二参与者，两两随机配对，按照基于shamir秘密共享协议的第二预设量子双方求积方法，由多个第三方协助对各第二参与者所持有的第二秘密执行量子安全双方求积，各第三方保留本地的所有测量结果并计算总和后广播，对所有广播的结果计算总和得到各第二参与者所持有的所述第二秘密的积；其中，所述第二设定数量为偶数；执行量子安全双方距离计算：对随机的第三参与者所持有的第一横坐标和第一纵坐标，以及随机的第四参与者所持有的第二横坐标和第二纵坐标；所述第三参与者计算所述第一横坐标和所述第一纵坐标的平方和得到第一数值秘密，将所述第一横坐标乘以负二得到第二数值秘密，将所述第一纵坐标乘以负二得到第三数值秘密；所述第四参与者计算所述第二横坐标和所述第二纵坐标的平方和得到第四数值秘密，将所述第二横坐标作为第五数值秘密，将所述第二纵坐标作为第六数值秘密；按照基于shamir秘密共享协议的所述第一预设量子双方求和方法，由多个第三方协助对所述第一数值秘密和所述第四数值秘密执行量子双方求和，并由各第三方保留获得的第一测量值；按照基于shamir秘密共享协议的所述第二预设量子双方求积方法，由多个第三方协助对所述第二数值秘密和所述第五数值秘密执行量子双方求积，各第三方保留获得的第二测量值；按照基于shamir秘密共享协议的所述第二预设量子双方求积方法，由多个第三方协助对所述第三数值秘密和所述第六数值秘密执行量子双方求积，各第三方保留获得的第三测量值；各第三方将本地的所述第一测量值、所述第二测量值和所述第三测量值求和并广播，对所有广播的结果求和得到所述第三参与者与所述第四参与者之间的平方距离，对所述平方距离开根得到所述第三参与者与所述第四参与者之间的距离。2.根据权利要求1所述的量子安全多方计算系统，其特征在于，所述系统执行量子多方求和操作过程中，当所述第一设定数量为奇数时，添加一个第一虚拟参与者参与各第一参与者的两两随机配对，所述第一虚拟参与者持有的所述第一秘密设为0。3.根据权利要求1所述的量子安全多方计算系统，其特征在于，所述系统执行量子多方求积操作过程中，当所述第二设定数量为奇数时，添加一个第二虚拟参与者参与各第二参与者的两两随机配对，所述第二虚拟参与者持有的所述第二秘密设为1。4.根据权利要求1所述的量子安全多方计算系统，其特征在于，基于shamir秘密共享协议的所述第一预设量子双方求和方法，包括：由持有秘密a的一方随机构建多项式f(y)＝a+c1y+c2y2+
…
+c
t-1
y
t-1
mod d，并基于shamir秘密共享协议计算经典份额f(y
i
)，i＝1，2
…
，n；c和d为参数，y为变量；由持有秘密b的一方随机构建多项式并基于shamir秘密共享协议计算经典份额g(y
i
)，i＝1，2
…
，n；c和d为参数，y为变量；将f(y
i
)与g(y
i
)在一组第三方tp1,
tp2,
…
tp
n
之间进行分配；对于参与恢复的第三方tp
k
，k＝1，2，...，t，计算h(y
k
)＝f(y
k
)+g(y
k
)mod d，保守秘密h(y
k
)，并计算h(y
k
)的份额的投影由tp1制备t个单粒子对第一个粒子|0>1施加量子离散傅里叶变换操作，得到的结果为由tp1通过在粒子1和粒子k，k＝2，3，...，t之间分别执行t-1次纠缠门操作，生成纠缠态并将粒子k分发给tp
k
，k＝1，2，...，t；每个tp
k
对分发到的粒子k做量子离散傅里叶变换操作，再执行广义泡利算子量子态|φ2>变为每个tp
k
测量自己的粒子|m
k
+a
k
>并宣布m
k
+a
k
，各tp
k
将测量结果相加，计算总和得到秘密a和秘密b的和为：5.根据权利要求4所述的量子安全多方计算系统，其特征在于，基于shamir秘密共享协议的第二预设量子双方求积方法，包括：由持有秘密a的一方随机构建多项式f(y)＝a+c1y+c2y2+
…
+c
t-1
y
t-1
mod d，并基于shamir秘密共享协议计算经典份额f(y
i
)，i＝1，2
…
，n；由持有秘密b的一方随机构建多项式并基于shamir秘密共享协议计算经典份额g(y
i
)，i＝1，2
…
，n；将f(y
i
)与g(y
i
)在一组第三方tp1,tp2,
…
tp
n
之间进行分配；对于每个tp
i
，i＝1，2，...，n，计算h(y
i
)＝f(y
i
)
×
g(y
i
)mod d并使用一个随机多项式z
i
(x)＝h
′
(y
i
)+β1x+β2x2+
…
+β
t-1
x
t-1
mod d在n个tp之间计算共享；使用范德蒙德矩阵表示t，对于参与恢复的第三方tp
k
，每个tp
k
仅知晓其对应的t
k
，并计算份额投影为：由tp1制备t个单粒子|0>1，|0>2，...，|0>
t
，对第一个粒子|0>1施加量子离散傅里叶变换操作，得到的结果为由tp1通过在粒子1和粒子k，k＝2,3，...，t之间分别执行t-1次纠缠门操作，生成纠缠态并将粒子k分发给tp
k
，k＝1，2，...，t；每个tp
k
对分发到的粒子k做量子离散傅里叶变换操作，再执行广义泡利算子量子态|φ2>变为每个tp
k
测量自己的粒子|m
k
+b
k
>并宣布m
k
+b
k
，各tp
k
将测量结果相加，计算总和得到秘密a和秘密b的积为：6.根据权利要求5所述的量子安全多方计算系统，其特征在于，所述系统还用于执行量子安全四面体体积计算，包括：对于单个四面体由4个参与者坐标作为顶点，记为：a0(x0,y0,z0)，a1(x1，y1，z1)，a2(x2，
y2，z2)和a3(x3，y3，z3)；采用三阶行列式计算单个四面体的体积：记行列式：记行列式：记行列式：记行列式：对于行列式d0，d1，d2和d3中的每个项，由每个项所涉及坐标对应的三个参与者执行所述量子多方求积操作进行量子安全三方求积，得到各项的值，对各项求和得到d0，d1，d2和d3的和，并求解v(a0a1a2a3)。7.根据权利要求6所述的量子安全多方计算系统，其特征在于，所述系统还用于量子安全多方多面体体积计算，由多个参与者在不暴露私有坐标的情况下计算由其坐标组成的多面体体积，包括：将所述多面体拆分为多个四面体，每个四面体由四个参与者的坐标作为顶点；每个四面体对应的四个参与者执行量子安全四面体体积计算，并累加得到所述多面体体积。8.根据权利要求5所述的量子安全多方计算系统，其特征在于，所述系统还用于执行量子安全圆-圆相交面积计算，由两个参与者在不暴露私有圆坐标的情况下，计算两个私有圆的相交面积，包括：
由第五参与者在其私有圆内随机选择第三设定数量个标记点；由所述第五参与者与第六参与者执行所述量子安全双方距离计算，以计算各标记点到所述第六参与者私有圆圆心的距离，并计算处于所述第六参与者私有圆范围内的标记点数量；则所述第五参与者与所述第六参与者之间私有圆的相交面积s
∩
计算式为：其中，k表示处于所述第六参与者私有圆范围内的标记点数量，n表示所述第三设定数量，s
c
表示所述第五参与者的私有圆的面积。9.根据权利要求8所述的量子安全多方计算系统，其特征在于，由两个参与者在不暴露私有圆坐标的情况下，计算两个私有圆的相交面积，包括：由所述第五参与者在其私有圆内随机选择n个标记点{p1，p2，...，p
n
}；由所述第六参与者构造一组随机数r
n
(n＝1，2，...，n)；由所述第五参与者与第六参与者执行所述量子安全双方距离计算，以计算各标记点pi到所述第六参与者私有圆的圆心o1的平方距离参与协助的每个第三方tp保留关于份额不公布，i＝1，2，...，n；基于shamir秘密共享协议的第二预设量子双方求积方法，所述第五参与者与所述第六参与者在r
n
(n＝1,2，...，n)与常数1之间进行n次量子双方求积，每个第三方tp获得r1,r2，...，r
n
的份额；tp
i
将和r
i
的份额对应相加，以获得的份额；由tp
i
意外的其他tp公布份额测量结果，则tp1获得由所述第六参与者与tp1在第三方tp0的协助下适用量子百万富翁协议对r
12
+r1和r
12
+r2和和分别进行比较，并由第三方tp0计算所述第六参与者拥有较大元素的个数，记为k，k为处于所述第六参与者私有圆范围内的标记点数量；由第三方tp0将k发送至所述第五参与者，并由所述第五参与者计算所述第五参与者与所述第六参与者之间私有圆的相交面积s
∩
计算式为：其中，k表示处于所述第六参与者私有圆范围内的标记点数量，n表示所述第三设定数量，s
c
表示所述第五参与者的私有圆的面积。10.根据权利要求9所述的量子安全多方计算系统，其特征在于，由所述第五参与者计算所述第五参与者与所述第六参与者之间私有圆的相交面积之后，还包括：由所述第五参与者将所述相交面积发送至所述第六参与者。

技术总结
本发明提供一种量子安全多方计算系统，基于对Shamir秘密共享协议两方求和求积方案的延伸，借助几个半诚实的第三方，为多方参与者提供了量子多方求和求积的执行方案，以及量子安全双方距离计算方案，消除了对数据库的依赖，解决中间值泄露的问题，在准确性和效率方面都有所提高，极大减少了多方量子安全计算所需的资源。提供对以多方参与者坐标作为顶点的多面体求体积的计算方案，通过将多面体拆解为多个四面体，引入代数对四面体体积进行表达，将对四面体求体积的问题转化为多方求积和求和问题；以及，提供对两方参与者私有圆计算圆圆相交面积的方案，通过执行量子安全双方距离计算，根据随机标记点落入对方私有圆的数量计算双方私有圆相交面积。算双方私有圆相交面积。算双方私有圆相交面积。


技术研发人员：窦钊 王艺霏 刘钊乾 彭海朋 杨榆 毕经国 赖裕平
受保护的技术使用者：北京邮电大学
技术研发日：2023.07.04
技术公布日：2023/10/11