云计算环境下使用全同态加密进行安全统计分析的方法

文档序号:25044034发布日期:2021-05-14 11:37阅读:117来源:国知局
云计算环境下使用全同态加密进行安全统计分析的方法

1.本发明涉及云计算信息安全的技术领域,更具体地,涉及一种云计算环境下使用全同态加密进行安全统计分析的方法。


背景技术:

2.近年来,随着海量数据的涌入和时代的发展,各类业务都需要更多的数据存储空间和更强的数据处理能力,在数据处理方面,单一的计算引擎已经满足不了用户的需求,云计算和大数据技术应运而生,随着云计算和大数据技术的快速发展,云服务器的存储空间和计算能力使得云中进行数据挖掘和统计处理成为可能,但是在云服务模式下如何保证用户的隐私安全成为其推广和应用中面临的首要问题。
3.为了防止数据泄露事件的发生,对数据加密后再存储在云平台上是公认最有效的方法,但这带来了一个新的问题:加密后的密文毫无特征,云平台无法对其进行有效处理,丧失了它在大数据处理方面的优势,而一旦将数据解密又极易造成隐私数据泄露,甚至带来严重的灾难性后果,在云端进行数据挖掘和统计分析的过程中,既要保证数据的隐私性,又要保证其准确性,全同态加密技术作为解决这一问题的关键手段,是一种可以对密文进行操作但仍可以恢复明文的加密算法,如2019年2月12日中国专利(公布号:cn109327304a)中公开了一种云计算中实现隐私保护的轻量级同态加密方法,在不解密的情况下实现运算,有效解决了云环境下数据上传服务器端,云端不可信的问题,但在保证隐私的同时,无法对数据进行统计分析,因此不能确定计算结果的准确性,无法保证数据计算精度。


技术实现要素:

4.为解决在计算外包时,云端不可信的情况下如何进行统计计算,保证数据计算精度的问题,本发明提出一种云计算环境下使用全同态加密进行安全统计分析的方法,在保证数据隐私的同时,对数据进行统计分析,从而确定计算结果的准确性,保证数据的计算精度。
5.为了达到上述技术效果,本发明的技术方案如下:
6.一种云计算环境下使用全同态加密进行安全统计分析的方法,至少包括:
7.s1.明确参与安全统计分析的用户客户端;
8.s2.在本地生成公钥pk和私钥sk并保存,将公钥pk告知云服务器,私钥sk告知私有服务器;
9.s3.利用公钥pk对用户客户端的统计数据进行同态加密,加密后上传至云服务器;
10.s4.云服务器通过全同态加密的乘法协议、加法协议及整数比较协议,在加密域下进行用户客户端的统计量分析、统计结果检验及一元线性回归拟合;
11.s5.云服务器将加密的统计量分析及统计结果检验返回至本地,本地利用私钥sk对统计量分析、统计结果检验及一元线性回归拟合结果进行解密。
12.在本技术方案中,首先基于本地、用户客户端与云服务器构成的三方系统,明确参
与安全统计分析的用户客户端,然后从隐私保护的角度出发,基于全同态加密的乘法协议、加法协议及整数比较协议,对用户客户端的统计数据进行同态加密,加密后上传至不受信任的云服务器,结合统计学方法,在云端环境下,直接对密文状态下的数据进行计算,在不泄露隐私的情况下完成计算,得到相应的统计量分析、统计结果检验及一元线性回归拟合,实现了真正的计算外包,而且基于全同态加密技术可以保证计算精度。
13.优选地,步骤s2所述在本地生成公钥pk和私钥sk的过程为:
14.s21.设置正整数q及t,其中,q表示密文模用于约减密文多项式的系数;t表示明文模用于约减明文多项式的系数;q>>t;
15.s22.令[q/t],定义r=z/(x
n
+1),r
q
=z
q
[x]/(x
n
+1)以及r
t
=z
t
[x]/(x
n
+1);其中,z表示具有整数系数的多项式,z
q
[x]表示一组具有整数系数的多项式集,整数系数范围为[0,q

1),χ表示一个离散的高斯误差分布;
[0016]
s23.令e

χ
n
,其中,s为私钥sk,a、e均视为r
q
中的元素,n总为2的幂次方,且n个系数独立于给定的分布采样,u
k
表示在z∩[

k/2,k/2)上的随机分布,公钥pk表示为:
[0017][0018]
其中,[.]
q
表示约减密文多项式的系数。
[0019]
在此,t通常为质数。
[0020]
优选地,步骤s3中所述利用公钥pk对用户客户端的统计数据进行同态加密的过程包括:
[0021]
设用户客户端的统计数据矩阵为x,统计数据矩阵x中的任意一个明文元素设为m,使用公钥pk对统计数据矩阵x中的每一个明文元素m进行加密,将s、a、e
i
视为r
q
中的元素,e1,e2←
χ
n
;创建的加密密文表示为:
[0022]
ct=a
·
pk+(m+e1,e2)。
[0023]
优选地,步骤s4所述的全同态加密的乘法协议为:
[0024]
s411.令q=q2,设s
’←
s2,在r
q
上均匀选取元素a’,在高斯分布χ上选取误差e’,将执行密钥设置为其中,b’=[

a’s+e’+qs’]
q

[0025]
s412.将明文m1及明文m2分别同态加密为密文ct1=(c0,c1)及密文ct2=(d0,d1),其中,c0、c1为密文ct1的两个成员,d0、d1为密文ct2的两个成员;并发送给云服务器;
[0026]
s413.密文ct1=(c0,c1)及密文ct2=(d0,d1)在云服务器中进行乘法运算ct1×
ct2,得到(x0,x1,x2)=([c0d0]
q
,[c0d1+c1d0]
q
,[c1d1]
q
),其中x0,x1,x2为乘积密文的三个成员;
[0027]
s414.云服务器将乘积结果返回本地,本地解密后得到明文域的乘积结果。
[0028]
优选地,步骤s4所述的全同态加密的加法协议为:
[0029]
s401.将明文m1及明文m2分别同态加密为密文ct1=(c0,c1)及密文ct2=(d0,d1),其中,c0、c1为密文ct1的两个成员,d0、d1为密文ct2的两个成员;并发送给云服务器;
[0030]
s402.密文ct1=(c0,c1)及密文ct2=(d0,d1)在云服务器中进行加法运算,即:
[0031]
s403.云服务器得到明文m1及明文m2的加密域求和结果ct
sum
=[ct1+ct2]
q
并返回本地,本地解密后得到明文域的加法结果。
[0032]
优选地,步骤s4所述的全同态加密的整数比较协议为:
[0033]
s421.由一位加密x与一位加密y,得到([[x]],[[y]]);
[0034]
s422.判断x是否小于y,若是,返回[[1]];否则,不变;
[0035]
s423.对于b=(x
l
‑1<y
l
‑1)∨((x
l
‑1=y
l
‑1)∧(x
l
‑2<y
l
‑2)).....∨((x
l
‑1=y
l
‑1)(x
l
‑2=y
l
‑2)....(x
l
=y
l
)(x0<y0)),l表示下标参数,满足(x
i
<y
i
)=(1

x
i
)(y
i
),且
[0036]
优选地,步骤s4所述用户客户端的统计量分析包括:均值分析、方差分析及中位数分析,统计结果检验包括:
[0037]
t

检验,检验用户客户端统计数据的总体方差未知、正态分布或近似正态分布单样本的均值是否与已知的总体均值相等;
[0038]
2*2列联表χ2检验,检验用户客户端两个及两个以上的统计数据样本率以及两个分类变量的关联性;
[0039]
单因素方差分析及两样本wilcoxon

秩和检验,所述两样本wilcoxon

秩和检验与t检验相对应。
[0040]
在此,2*2列联表χ2检验根本思想在于比较理论频数和实际频数的吻合程度或拟合优度问题,单因素方差分析是将所有测量值间的总变异按照其变异的来源分解为多个部份,然后进行比较,评价由某种因素所引起的变异是否具有统计学意义,两样本wilcoxon

秩和检验用于不满足正态分布的小样本中,检验两独立样本之间的差异是否显著。
[0041]
优选地,步骤s4所述的一元线性回归拟合为用户客户端的统计数据中两个变量之间相关关系的拟合,设自变量为x,因变量为y,自变量x与因变量y之间的一元线性回归方程为:
[0042]
y=β0+β1x+ε
[0043]
其中,β0和β1均表示未知参数,β0为回归常数,β1为回归系数;ε表示其他随机因素的影响;
[0044]
利用最小二乘准则进行参数估计:
[0045][0046]
其中,y
i
表示第i次观测,因变量的观测值;表示第i次观测,因变量的预测值。
[0047]
优选地,自变量x与因变量y之间的一元线性回归方程的斜率为:
[0048][0049]
其中,表示自变量x的均值,表示自变量y的均值;表示解决估计量,表达式为:
[0050]
[0051]
其中,w表示自变量x或因变量y的总个数,x
i
表示第i个自变量;y
i
表示第i个因变量。
[0052]
优选地,步骤s5所述本地利用私钥sk对统计量分析、统计结果检验及一元线性回归拟合结果进行解密的过程为:
[0053]
设公钥pk加密后的密文为ct=(c0,c1),利用私钥s=sk对密文ct=(c0,c1)进行解密,得到解密明文m:
[0054][0055]
其中,表示舍入到最近的整数,为取整操作。
[0056]
与现有技术相比,本发明技术方案的有益效果是:
[0057]
本发明提出一种云计算环境下使用全同态加密进行安全统计分析的方法,首先基于本地、用户客户端与云服务器构成的三方系统,明确参与安全统计分析的用户客户端,然后从隐私保护的角度出发,基于全同态加密的乘法协议、加法协议及整数比较协议,对用户客户端的统计数据进行同态加密,加密后上传至不受信任的云服务器,在云端环境下直接对密文状态下的数据进行计算,在不泄露隐私的情况下完成计算,得到相应的统计量分析、统计结果检验及一元线性回归拟合,实现了真正的计算外包,而且基于全同态加密技术保证了计算精度。
附图说明
[0058]
图1表示本发明实施例中提出的云计算环境下使用全同态加密进行安全统计分析的方法。
具体实施方式
[0059]
附图仅用于示例性说明,不能理解为对本专利的限制;
[0060]
为了更好地说明本实施例,附图某些部位会有省略、放大或缩小,并不代表实际尺寸;
[0061]
对于本领域技术人员来说,附图中某些公知内容说明可能省略是可以理解的。
[0062]
下面结合附图和实施例对本发明的技术方案做进一步的说明。
[0063]
实施例1
[0064]
如图1所示的云计算环境下使用全同态加密进行安全统计分析的方法的流程示意图,参见图1,包括:
[0065]
s1.明确参与安全统计分析的用户客户端;
[0066]
s2.在本地生成公钥pk和私钥sk并保存,将公钥pk告知云服务器,私钥sk告知私有服务器;
[0067]
在本地生成公钥pk和私钥sk的过程为:
[0068]
s21.设置正整数q及t,其中,q表示密文模用于约减密文多项式的系数;t表示明文模用于约减明文多项式的系数;q>>t;
[0069]
s22.令[q/t],定义r=z/(x
n
+1),r
q
=z
q
[x]/(x
n
+1)以及r
t
=z
t
[x]/(x
n
+1);其中,z表示具有整数系数的多项式,z
q
[x]表示一组具有整数系数的多项式集,整数系数范围为
[0,q

1),χ表示一个离散的高斯误差分布;
[0070]
s23.令e

χ
n
,其中,s为私钥sk,a、e均视为r
q
中的元素,n总为2的幂次方,且n个系数独立于给定的分布采样,u
k
表示在z∩[

k/2,k/2)上的随机分布,公钥pk表示为:
[0071][0072]
其中,[.]
q
表示约减密文多项式的系数。
[0073]
在此,t通常为质数。
[0074]
s3.利用公钥pk对用户客户端的统计数据进行同态加密,加密后上传至云服务器;
[0075]
利用公钥pk对用户客户端的统计数据进行同态加密的过程包括:
[0076]
设用户客户端的统计数据矩阵为x,统计数据矩阵x中的任意一个明文元素设为m,使用公钥pk对统计数据矩阵x中的每一个明文元素m进行加密,将s、a、e
i
视为r
q
中的元素,e1,e2←
χ
n
;创建的加密密文表示为:
[0077]
ct=a
·
pk+(m+e1,e2)。
[0078]
具体的,以用户客户端alice为例,设其统计数据矩阵x具体表示a,表达式为:
[0079][0080]
对a中的每一个元素m进行加密,随机选定一个正整数q,质数t,满足q>>t,则元素m加密后的密文ct,对于每一个元素加密后,上传至云服务器。
[0081]
s4.云服务器通过全同态加密的乘法协议、加法协议及整数比较协议,在加密域下进行用户客户端的统计量分析、统计结果检验及一元线性回归拟合;
[0082]
用户客户端的统计量分析包括:均值分析、方差分析及中位数分析,统计结果检验包括:
[0083]
t

检验,检验用户客户端统计数据的总体方差未知、正态分布或近似正态分布单样本的均值是否与已知的总体均值相等;包括:
[0084]
1)假设检验:
[0085]
h0:样本均值与总体均值相等,h1:样本均值与总体均值不相等
[0086]
2)构造t检验统计量:
[0087]
对于要检验样本均值是否等于总体均值的双侧检验,若则拒绝原假设,认为样本值与均值不等,否则不拒绝原假设;对于要检验样本均值是否比总体均值大的单侧检验,若则拒绝原假设,认为样本值不大于总体均值,否则不拒绝原假设;对于要检验样本均值是否比总体均值小的单侧检验,若则拒绝原假设,认为样本值不小于总体均值,否则不拒绝原假设。
[0088]
2*2列联表χ2检验,检验用户客户端两个及两个以上的统计数据样本率以及两个分类变量的关联性;
[0089]
单因素方差分析及两样本wilcoxon

秩和检验,所述两样本wilcoxon

秩和检验与t检验相对应。
[0090]
单因素方差分析包括:
[0091]
总方差(total variation):全部测量值x
ij
与总均数间μ的差异,与总均数间μ的差异,自由度为dt=n

1;
[0092]
组间方差(between group variation):各组的均数μ
i
与总均数μ间的差异,自由度为db=m

1;
[0093]
组内方差(within group variation):每组的每个测量值x
ij
与该组均数μ
i
的差异,自由为其中,m表示观察组的个数,n
i
表示每一组的观察数,n表示总的观察数;构造f统计量进行检验,当给定显著性水平为α时,f的拒绝域为
[0094]
两样本wilcoxon

秩和检验:也叫mann

whitney u检验,是一种非参数检验方法,与参数检验中的t检验相对应,主要是用于不满足正态分布的小样本,检验两独立样本x和y之间的差异是否显著,包括:
[0095]
建立假设检验:
[0096]
h0:x、y的中位数相等h1:x、y的中位数不相等
[0097]
将两个样本数据x和y混合并由小到大进行等级排列,最小的数据秩次编为1,最大的数据秩次编为n1+n2;
[0098]
令r
i
和r
j
分别为样本x
i
和样本点y
j
在混合样本中的秩,w
x
和w
y
分别为r
i
和r
j
的总和(秩和):
[0099][0100]
令w为w
x
和w
y
中较小者,即w=min(w
x
,w
y
)称为wilcoxon秩和统计量;把w值与秩和检验表中某α显著性水平下的临界值相比较,如果w1<w<w2,则两样本差异不显著;如果w不等于w1或w大于等于w2,则表明两样本差异显著。
[0101]
具体的,步骤s4包括:
[0102]
s41.设用户客户端的列向量均值其中:
[0103][0104][0105]
[0106][0107]
s42.设用户客户端的列向量方差σ2=(σ
12
,σ
22
,...,σ
n2
),其中:
[0108][0109][0110][0111][0112]
s43.对用户客户端alice的数据矩阵b求解中位数:
[0113]
b=(x1,x2,

,x
n
)
[0114]
s431.对b中的每一个元素m进行加密,同时将每一个元素m转化为二进制再按位加密,随机选定一个正整数q,质数t,满足q>>t,则元素m加密后的密文ct满足:
[0115]
ct=a
·
pk+(m+e1,e2),
[0116]
其中,e1,e2←
χ
n
,pk为公钥,对于每一个元素加密及二进制按位加密后,上传至云服务器;
[0117]
s432.使用全同态加密整数比较协议,比较相邻的元素。如果x1>x2,则二者交换;
[0118]
s433.对每一对相邻元素做同样的工作,从开始第一对到结尾的最后一对。这步做完后,最后的元素会是最大的数;
[0119]
s434.针对所有的元素重复步骤s433,除了最后一个;
[0120]
s435.持续每次对越来越少的元素重复上面的步骤,直到没有任何一对数字需要比较;
[0121]
s436.上述排序完成后,根据中位数的定义,若b中的元素个数n为偶数,则中位数为(x
n/2
+x
n/2+1
)的和,解密后的平均值;若b中的元素个数为奇数,则中位数为x
(n+1)/2
解密后的值。
[0122]
s44.对用户客户端alice的数据矩阵c进行t

检验,c=(x1,x2,

,x
n
);
[0123]
s441.对c中的每一个元素m进行加密,随机选定一个正整数q,质数t,满足q>>t,则元素m加密后的密文ct满足:
[0124]
ct=a
·
pk+(m+e1,e2),
[0125]
其中,e1,e2←
χ
n
,pk为公钥,对于每一个元素加密后,上传至云服务器;
[0126]
s442.基本假设(总体服从正态分布,μ已知)
[0127]
h0:样本均值与总体均值μ相等,h1:样本均值与总体均值μ相等;
[0128]
s443.根据具体实施步骤s41、s42,求得样本均值标准差σ
x
,解密后,构造t检验统计量:
[0129]
s444.若则拒绝原假设,认为样本均值与总体均值μ不等;反之样本均值与总体均值μ相等。
[0130]
s45.对表1中2*2列联表进行χ2检验:
[0131]
表1
[0132]
组别\属性wz总和xaba+bycdc+d总和a+cb+da+b+c+d
[0133]
s451.基本假设为:
[0134]
h0:观察频数与期望频数没有差别;h1:观察频数与期望频数有差别。
[0135]
s452.对表1每一个元素m进行全同态加密,得到矩阵,
[0136][0137]
各元素m加密后的密文ct满足:
[0138]
ct=a
·
pk+(m+e1,e2),
[0139]
其中,e1,e2←
χ
n
,pk为公钥,对于每一个元素加密后,上传至云服务器;
[0140]
s453.据全同态加密的加法、乘法协议,求解=(a
enc
*d
enc

b
enc
*c
enc
)2*(a+b+c+d)
enc
,q=((a+b)
enc
*(c+d)
enc
*(a+c)
enc
*(b+d)
enc
),解密后得到卡方统计量为:
[0141][0142]
s454.给定显著性水平α,若则拒绝h0,认为“属性”和“组别”有关;反之,则无关。
[0143]
s46.假设某个因素有s个分组,对应的观测值分别为:
[0144][0145][0146]
进行单因素方差分析;
[0147]
s461.基本假设
[0148]
h0:μ1=μ2=

=μ
n h1:μ1、μ2…
μ
n
不全相等
[0149]
s462.对样本数据中的每一个元素m进行加密,随机选定一个正整数q,质数t,满足q>>t,则元素m加密后的密文ct满足
[0150]
ct=a
·
pk+(m+e1,e2)。
[0151]
其中,e1,e2←
χ
n
,pk为公钥,对于每一个元素加密后,上传至云服务器;
[0152]
s463.根据具体实施步骤s41、s42,求得组内平方和ssw与组间平方和ssb,如下,
[0153]
组内平方和:自由度为n

s
[0154]
组间平方和:自由度为s
‑1[0155]
s464.解密后,构造f统计量,
[0156][0157]
s465.若f>f1‑
α
(s

1,n

s),则拒绝原假设,认为μ1、μ2…
μ
n
不全相等。
[0158]
s47.wilcoxon

秩和检验考察两样本总体x=(x1,x2,...,x
n
)和y=(y1,y2,...,y
n
)所代表的总体的中位数是否一样。
[0159]
s471.基本假设:h0:m
x
=m
y
h1:m
x
、m
y
不相等
[0160]
s472.对样本x的每一个数都标记为类标签0,对样本y的每一个数都标记为类标签1,并将其混合起来。对样本数据x、y中的每一个元素m进行加密,随机选定一个正整数q,质数t,满足q>>t,则元素m加密后的密文ct满足
[0161]
ct=a
·
pk+(m+e1,e2)。
[0162]
其中,e1,e2←
χ
n
,pk为公钥,对于每一个元素加密及二进制按位加密后,上传至云服务器;
[0163]
s473.根据全同态加密整数比较协议,参照具体实施步骤s43中的排序原理,将x、y混合样本中的2n个数进行排序编秩;
[0164]
s474.根据标签指示计算样本x
i
和样本点y
j
在混合样本中的秩r
i
和r
j
,令w
x
和w
y
分别为r
i
和r
j
的总和(秩和):
[0165][0166]
s475.令w为w
x
和w
y
中较小者,即w=min(w
x
,w
y
)称为wilcoxon秩和统计量;
[0167]
s476.把w值与秩和检验表中某α显著性水平下的临界值相比较,如果w1<w<w2,则两样本差异不显著;如果w不等于w1或w大于等于w2,则表明两样本差异显著。
[0168]
s48.假设从总体中获取了n组观察值(x1,y1),(x2,y2),

,(x
n
,y
n
),使用一元线性回归方程拟合,假设方程为:
[0169]
y=β0+β1x+ε
[0170]
s481.对样本数据中的每一个元素m进行加密,随机选定一个正整数q,质数t,满足q>>t,则元素m加密后的密文ct满足
[0171]
ct=a
·
pk+(m+e1,e2)。
[0172]
其中,e1,e2←
χ
n
,pk为公钥,对于每一个元素加密后,上传至云服务器;
[0173]
s482.求得回归方程的斜率与截距估计量:
[0174]
s483.解密后得到一元线性回归方程
[0175]
全同态加密加法协议包括的过程为:
[0176]
s401.将明文m1及明文m2分别同态加密为密文ct1=(c0,c1)及密文ct2=(d0,d1),其中,c0、c1为密文ct1的两个成员,d0、d1为密文ct2的两个成员,并发送给云服务器;
[0177]
s402.云服务器中进行加法运算
[0178]
s403.云服务器得到明文m1及明文m2的加密域求和结果ct
sum
=[ct1+ct2]
q
并返回本地,本地解密后得到明文域的加法结果。
[0179]
全同态加密乘法协议包括的过程为:
[0180]
s411.令q=q2,设s
’←
s2,在r
q
上均匀选取元素a’,在高斯分布χ上选取误差e’,将执行密钥设置为其中,b’=[

a’s+e’+qs’]
q

[0181]
s412.将明文m1及明文m2分别同态加密为密文ct1=(c0,c1)及密文ct2=(d0,d1),其中,c0、c1为密文ct1的两个成员,d0、d1为密文ct2的两个成员;并发送给云服务器;
[0182]
s413.密文ct1=(c0,c1)及密文ct2=(d0,d1)在云服务器中进行乘法运算ct1×
ct2,得到(x0,x1,x2)=([c0d0]
q
,[c0d1+c1d0]
q
,[c1d1]
q
),其中x0,x1,x2为乘积密文的三个成员;
[0183]
s414.云服务器将乘积结果返回本地,本地解密后得到明文域的乘积结果。
[0184]
全同态加密整数比较协议包括的过程为:
[0185]
s421.由一位加密x与一位加密y,得到([[x]],[[y]]);
[0186]
s422.判断x是否小于y,若是,返回[[1]];否则,不变;
[0187]
s423.对于b=(x
l
‑1<y
l
‑1)∨((x
l
‑1=y
l
‑1)∧(x
l
‑2<y
l
‑2)).....∨((x
l
‑1=y
l
‑1)(x
l
‑2=y
l
‑2)....(x
l
=y
l
)(x0<y0)),l表示下标参数,满足(x
i
<y
i
)=(1

x
i
)(y
i
),且
[0188]
s5.云服务器将加密的统计量分析及统计结果检验返回至本地,本地利用私钥sk对统计量分析、统计结果检验及一元线性回归拟合结果进行解密。具体为:
[0189]
设公钥加密后的密文为ct=(c0,c1),利用私钥s=sk对密文ct=(c0,c1)进行解密,得到解密明文m:
[0190][0191]
其中,表示舍入到最近的整数,为取整操作。
[0192]
附图中描述位置关系的用于仅用于示例性说明,不能理解为对本专利的限制;显然,本发明的上述实施例仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
当前第1页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1