【自然语言处理】【知识图谱】SEU:无监督、非神经网络实体对齐超越有监督图神经网络?

SEU:从对齐到分配----简单有效的无监督实体对齐
《From Alignment to Assignment:Frustratingly Simple Unsupervised Entity Alignment》

论文地址:https://arxiv.org/pdf/2109.02363.pdf

相关博客:
【自然语言处理】【知识图谱】利用属性、值、结构来实现实体对齐
【自然语言处理】【知识图谱】基于图匹配神经网络的跨语言知识图谱对齐
【自然语言处理】【知识图谱】使用属性嵌入实现知识图谱间的实体对齐
【自然语言处理】【知识图谱】用于实体对齐的多视角知识图谱嵌入
【自然语言处理】【知识图谱】MTransE:用于交叉知识对齐的多语言知识图谱嵌入
【自然语言处理】【知识图谱】SEU:无监督、非神经网络实体对齐超越有监督图神经网络?

一、简介

  • 知识图谱

    知识图谱已经被广泛用于推荐系统、问答系统等下游的任务中。近些年,由不同的组织使用不同的语言在不同领域构建了大量的知识图谱。这些交叉语言知识图谱通常具有独特的信息,但也存在着一些覆盖。如果能够集成这些交叉语言知识图谱,那么就能够为那些缺乏语言资源的用户提供更广阔的视角。因此,交叉语言知识图谱吸引了越来越多的注意力。

  • 传统实体对齐

    交叉语言实体对齐的目标是,发现知识图谱中的等价实体。传统的方法主要依靠词汇匹配和概率推理,但需要机器翻译系统来解决交叉语言问题。但是,现有的机器翻译系统在有效上下文信息的情况下,并不能实现高的准确率。特别是对于那些不相似的语言,例如 Chinese-English \text{Chinese-English} Chinese-English Japanese-English \text{Japanese-English} Japanese-English

  • 图神经网络实体对齐

    最近, GCN \text{GCN} GCN及其变体已经在各种图相关的应用中取得了state-of-the-art的结果。直觉上, GNN \text{GNN} GNN能够更好的捕获知识图谱的结构信息。一些基于 GNN \text{GNN} GNN的实体对齐方法确实明显地改善了公开数据集上的效果。所有基于 GNN \text{GNN} GNN的实体对齐方法都基于一个核心假设,即需要对齐的实体对具有相似的局部结构。然而,基于 GNN \text{GNN} GNN的方法在带来优秀结果的同时,也不可避免地继承了神经网络的缺点:

    • 糟糕的解释性

      集合了非线性运算和大量参数的 GNN \text{GNN} GNN方法难以被解释。因此,许多研究人员将 GNN \text{GNN} GNN当做是黑盒,并只专注于提高性能指标。这种风气导致很难判断新的设计是真的有效,还是仅仅在特定数据集上过拟合。近期的一些研究表明,一些先进的实体对齐方法在若干个公开数据集上被传统方法打败。

    • 效率低

      为了改进效果,新的实体对齐方法堆叠了大量的新技术,例如:图注意力网络、图匹配网络和联合训练等。因此,整个架构变的越来越复杂,导致时间和空间复杂度显著增加。

  • 本文贡献

    本文尝试在不损失准确率的情况下,消除现有实体对齐方法的冗余设计来增强解释性和效率。基于假设:源和目标 KG \text{KG} KG的结构和文本特征是同构的,本文成功将实体对齐问题转换为分配问题。分配问题是一个基础的且被很好研究的组合优化问题,可以通过 Hungarian \text{Hungarian} Hungarian算法或者 Sinkhorn \text{Sinkhorn} Sinkhorn算法求解。

    基于上面的发现,本文提出了使用神经网络的实体对齐方法 SEU \text{SEU} SEU。相比于基于 GNN \text{GNN} GNN的实体对齐方法, SEU \text{SEU} SEU仅保留了用于特征传播的基础图卷积操作,抛弃了复杂的神经网络。在公开数据集上的实验结果显示:

    • GPU \text{GPU} GPU上, SEU \text{SEU} SEU能在若干秒内完成;即使在 CPU \text{CPU} CPU上,也能在数十秒内完成;
    • SEU \text{SEU} SEU在几乎所有公开数据集上都超越了有监督的state-of-the-art方法;

    此外,本文还讨论了现有的实体对齐方法表现不佳的原因,以及交叉语言实体对齐中的神经网络的必要性。

二、任务定义

​ 知识图谱以三元组 ( h , r , t ) (h,r,t) (h,r,t)的形式存储了真实世界的知识。一个知识图谱可以被定义为 G = ( E , R , T ) G=(E,R,T) G=(E,R,T),其中分别表示实体集合、关系集合和三元组集合。给定一个源图谱 G s = ( E s , R s , T s ) G_s=(E_s,R_s,T_s) Gs=(Es,Rs,Ts)和目标图谱 G t = ( E t , R t , T t ) G_t=(E_t,R_t,T_t) Gt=(Et,Rt,Tt),实体对齐的目标是发现两个图谱监督的对应关系 P \textbf{P} P

三、相关工作

1. 分配问题

​ 分配问题是一个基础的且被很好研究的组合优化问题。一个简单的例子:为N个工人分配N个工作。假设每个工人都能够在规定时间段内完成每项工作,令 x i j x_{ij} xij表示第 i i i个工人被分配到第 j j j项工作收益。分配问题的目标是寻找最优的分配计划,使得所有工作的收益最大化。该问题的求解形式化的定义为
a r g    m a x P ∈ P N ⟨ P , X ⟩ F (1) \mathop{arg\;max}_{\textbf{P}\in\mathbb{P}_N}\langle \textbf{P},\textbf{X} \rangle_F \tag{1} argmaxPPNP,XF(1)
X ∈ R N × N \textbf{X}\in\mathbb{R}^{N\times N} XRN×N是收益矩阵; P \textbf{P} P是排列矩阵,表示分配计划。 P \textbf{P} P中的每行或者每列仅有一个1,其余均为0。 P N \mathbb{P}_N PN表示所有N维排列矩阵的集合。 ⟨ ⋅ ⟩ F \langle\cdot\rangle_F F表示 Frobenius \text{Frobenius} Frobenius内积。

本文采用 Hungarian \text{Hungarian} Hungarian算法和 Sinkhorn \text{Sinkhorn} Sinkhorn算法解决分配问题。

四、 SEU \text{SEU} SEU

1. 对齐问题的分配建模

​ 方法 SEU \text{SEU} SEU的输入为四个矩阵: A s ∈ R ∣ E s ∣ × ∣ E s ∣ \textbf{A}_s\in\mathbb{R}^{|E_s|\times |E_s|} AsREs×Es A t ∈ R ∣ E t ∣ × ∣ E t ∣ \textbf{A}_t\in\mathbb{R}^{|E_t|\times |E_t|} AtREt×Et表示源图谱 G s G_s Gs和目标图谱 G t G_t Gt的邻接矩阵。 H s ∈ R ∣ E s ∣ × d \textbf{H}_s\in\mathbb{R}^{|E_s|\times d} HsREs×d H t ∈ R ∣ E t ∣ × d \textbf{H}_t\in\mathbb{R}^{|E_t|\times d} HtREt×d表示实体的文本特征,该特征通过机器翻译系统或者交叉语言词嵌入向量映射至了统一语义空间。

​ 类似于分配问题,实体对齐也需要满足 1-to-1 \text{1-to-1} 1-to-1约束。令排列矩阵 P ∈ P ∣ E ∣ \textbf{P}\in\mathbb{P}_{|E|} PPE表示 G s G_s Gs G t G_t Gt间的实体对应, P i j = 1 \textbf{P}_{ij}=1 Pij=1表示 e i ∈ G s e_i\in G_s eiGs e j ∈ G t e_j\in G_t ejGt是相等的实体对。 SEU \text{SEU} SEU的目标是根据 { A s , A t , H s , H t } \{\textbf{A}_s,\textbf{A}_t,\textbf{H}_s,\textbf{H}_t\} {As,At,Hs,Ht}求解 P \textbf{P} P。考虑下面的理想情况:

  • A s \textbf{A}_s As A t \textbf{A}_t At是同构的,即根据 P \textbf{P} P重新排序 A s \textbf{A}_s As的节点能够转换为 A t \textbf{A}_t At
    PA s P − 1 = A t (2) \textbf{PA}_s\textbf{P}^{-1}=\textbf{A}_t \tag{2} PAsP1=At(2)

  • 等价实体的文本特征能够被翻译系统完美映射。因此, H s \textbf{H}_s Hs H t \textbf{H}_t Ht也可以根据实体对应的 P \textbf{P} P来被对齐,
    PH s = H t (3) \textbf{PH}_s=\textbf{H}_t \tag{3} PHs=Ht(3)

​ 通过合并等式 ( 2 ) (2) (2) ( 3 ) (3) (3),五元组 { A s , A t , H s , H t , P } \{\textbf{A}_s,\textbf{A}_t,\textbf{H}_s,\textbf{H}_t,\textbf{P}\} {As,At,Hs,Ht,P}的关系如下:
( PA s P − 1 ) l PH s = A t l H t ∀ l ∈ N ⇒ PA s l H s = A t l H t (4) (\textbf{PA}_s\textbf{P}^{-1})^l\textbf{PH}_s=\textbf{A}_t^l\textbf{H}_t\quad \forall l\in \mathbb{N} \\ \Rightarrow\quad\textbf{PA}_s^l\textbf{H}_s=\textbf{A}_t^l\textbf{H}_t \tag{4} (PAsP1)lPHs=AtlHtlNPAslHs=AtlHt(4)
​ 在一对一约束 P ∈ P ∣ E ∣ \textbf{P}\in\mathbb{P}_{|E|} PPE的条件下,等式 ( 4 ) (4) (4)中的 P \textbf{P} P可以通过最小化 Frobenius \text{Frobenius} Frobenius ∥ PA s l H s − A t l H t ∥ \parallel \textbf{PA}_s^l\textbf{H}_s-\textbf{A}_t^l\textbf{H}_t \parallel PAslHsAtlHt求解。理论上,对于任意深度 l ∈ N l\in\mathbb{N} lN P \textbf{P} P的解应该都相同的。然而,上面的推断都是基于理想同构条件的。但是,实际中 G s G_s Gs G t G_t Gt并不是严格同构的,且翻译系统并不能完美的将文本特征映射至统一语义空间。为了减少实际中噪音的影响, P \textbf{P} P应该拟合各种深度的 l l l。因此,本文提出了下面的等式来解决交叉语言实体对齐问题
arg min P ∈ P ∣ E ∣ ∑ l = 0 L ∥ PA s l H s − A t l H t ∥ F 2 (5) \mathop{\text{arg min}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L\parallel \textbf{PA}_s^l\textbf{H}_s-\textbf{A}_t^l\textbf{H}_t \parallel_F^2 \tag{5} arg minPPEl=0LPAslHsAtlHtF2(5)

定理1

等式 ( 5 ) (5) (5)等价于解决下面的分配问题
arg max P ∈ P ∣ E ∣ ⟨ P , ∑ l = 0 L A t l H t ( A s l H s ) T ⟩ F (6) \mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_{|E|}} \Big\langle \textbf{P},\sum_{l=0}^L\textbf{A}_t^l\textbf{H}_t(\textbf{A}_s^l\textbf{H}_s)^T \Big\rangle_F\tag{6} arg maxPPEP,l=0LAtlHt(AslHs)TF(6)
证明:

基于 Frobenius \text{Frobenius} Frobenius方式的性质 ∥ A-B ∥ F 2 = ∥ A ∥ F 2 + ∥ B ∥ F 2 − 2 ⟨ A,B ⟩ F \parallel \textbf{A-B} \parallel_F^2=\parallel\textbf{A}\parallel_F^2+\parallel\textbf{B}\parallel_F^2-2\langle\textbf{A,B}\rangle_F A-BF2=AF2+BF22A,BF,对等式 ( 5 ) (5) (5)进行如下推导
arg min P ∈ P ∣ E ∣ ∑ l = 0 L ∥ PA s l H s − A t l H t ∥ F 2 = arg min P ∈ P ∣ E ∣ ∑ l = 0 L ∥ PA s l H s ∥ F 2 + ∥ A t l H t ∥ F 2 − 2 ⟨ PA s l H s , A t l H t ⟩ (7) \begin{aligned} &\mathop{\text{arg min}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L\parallel \textbf{PA}_s^l\textbf{H}_s-\textbf{A}_t^l\textbf{H}_t \parallel_F^2 \\ =&\mathop{\text{arg min}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L\parallel\textbf{PA}_s^l\textbf{H}_s\parallel_F^2+\parallel\textbf{A}_t^l\textbf{H}_t\parallel_F^2-2\langle\textbf{PA}_s^l\textbf{H}_s,\textbf{A}_t^l\textbf{H}_t\rangle \end{aligned} \tag{7} =arg minPPEl=0LPAslHsAtlHtF2arg minPPEl=0LPAslHsF2+AtlHtF22PAslHs,AtlHt(7)
由于排列矩阵 P \textbf{P} P是正交的,所以 ∥ PA s l H s ∥ F 2 \parallel\textbf{PA}_s^l\textbf{H}_s\parallel_F^2 PAslHsF2 ∥ A t l H t ∥ F 2 \parallel\textbf{A}_t^l\textbf{H}_t\parallel_F^2 AtlHtF2是常数。因此,等式 ( 7 ) (7) (7)等价于
arg max P ∈ P ∣ E ∣ ∑ l = 0 L ⟨ PA s l H s , A t l H t ⟩ F (8) \mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L\langle \textbf{PA}_s^l\textbf{H}_s,\textbf{A}_t^l\textbf{H}_t \rangle_F \tag{8} arg maxPPEl=0LPAslHs,AtlHtF(8)
对于任意的实数矩阵 A \textbf{A} A B \textbf{B} B,有两个恒成立的等式: ⟨ A,B ⟩ F = Tr ( AB T ) \langle\textbf{A,B}\rangle_F=\text{Tr}(\textbf{AB}^T) A,BF=Tr(ABT) ⟨ A,B+C ⟩ F = ⟨ A,B ⟩ F + ⟨ A,C ⟩ F \langle \textbf{A,B+C} \rangle_F=\langle\textbf{A,B}\rangle_F+\langle\textbf{A,C}\rangle_F A,B+CF=A,BF+A,CF,其中 Tr ( X ) \text{Tr}(\textbf{X}) Tr(X)表示矩阵 X \textbf{X} X的迹。因此,定理1可以被证明:
arg max P ∈ P ∣ E ∣ ∑ l = 0 L ⟨ PA s l H s , A t l H t ⟩ F = arg max P ∈ P ∣ E ∣ ∑ l = 0 L Tr ( PA s l H s ( A t l H t ) T ) = arg max P ∈ P ∣ E ∣ ∑ l = 0 L ⟨ P , A t l H t ( A s l H s ) T ⟩ F = arg max P ∈ P ∣ E ∣ ⟨ P , ∑ l = 0 L A t l H t ( A s l H s ) ⟩ (9) \begin{aligned} &\mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L\langle\textbf{PA}_s^l\textbf{H}_s,\textbf{A}_t^l\textbf{H}_t\rangle_F \\ =&\mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L \text{Tr}(\textbf{PA}_s^l\textbf{H}_s(\textbf{A}_t^l\textbf{H}_t)^T)\\ =&\mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_{|E|}}\sum_{l=0}^L \langle\textbf{P},\textbf{A}_t^l\textbf{H}_t(\textbf{A}_s^l\textbf{H}_s)^T\rangle_F\\ =&\mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_{|E|}}\Big\langle \textbf{P},\sum_{l=0}^L\textbf{A}_t^l\textbf{H}_t(\textbf{A}_s^l\textbf{H}_s) \Big\rangle \end{aligned} \tag{9} ===arg maxPPEl=0LPAslHs,AtlHtFarg maxPPEl=0LTr(PAslHs(AtlHt)T)arg maxPPEl=0LP,AtlHt(AslHs)TFarg maxPPEP,l=0LAtlHt(AslHs)(9)
通过定理1,成功将实体对齐问题转换为分配问题。需要注意的是,实体规模 ∣ E s ∣ |E_s| Es ∣ E t ∣ |E_t| Et在实践中通常是不一致的,导致收益矩阵不是方阵。这种不平衡分配的问题能够轻易转换为平衡分配问题。假设 ∣ E s ∣ > ∣ E t ∣ |E_s|>|E_t| Es>Et,一个朴素的简化方法是使用0来填充收益矩阵,使其形状成为 R ∣ E s ∣ × ∣ E s ∣ \mathbb{R}^{|E_s|\times|E_s|} REs×Es。这种朴素的方法适合于 ∣ E s ∣ |E_s| Es ∣ E t ∣ |E_t| Et差距较小的数据集。对于差距较大的数据集,可以使用更高效的简化算法https://www.hpl.hp.com/techreports/2012/HPL-2012-40R1.pdf。

2. 求解分配问题

​ 解决分配问题的第一个多项式时间复杂度算法是 Hungarian \text{Hungarian} Hungarian算法,其原始的复杂度为 O ( n 4 ) O(n^4) O(n4)。后来的一些改进将时间复杂度降低为 O ( n 3 ) O(n^3) O(n3)

​ 除了 Hungarian \text{Hungarian} Hungarian算法,分配问题还可以看作是最优运输问题的特例。在最优运输问题中,分配计划 P \textbf{P} P可以是随机矩阵。基于 Sinkhorn \text{Sinkhorn} Sinkhorn操作, Cuturi \text{Cuturi} Cuturi等人提出了一个快速且完全并行的算法来解决最优运输问题。
S 0 ( X ) = e x p ( X ) S k ( X ) = N c ( N r ( S k − 1 ( X ) ) ) Sinkhorn ( X ) = l i m k → ∞ S k ( X ) \begin{aligned} S^0(\textbf{X})&=exp(\textbf{X}) \\ S^k(\textbf{X})&=\mathcal{N}_c(\mathcal{N}_r(S^{k-1}(\textbf{X}))) \\ \text{Sinkhorn}(\textbf{X})&=\mathop{lim}_{k\rightarrow\infty}S^k(\textbf{X}) \end{aligned} S0(X)Sk(X)Sinkhorn(X)=exp(X)=Nc(Nr(Sk1(X)))=limkSk(X)
其中, N r ( X ) = X ⊘ ( X 1 N 1 N T ) \mathcal{N}_r(\textbf{X})=\textbf{X}\oslash(\textbf{X}\textbf{1}_N\textbf{1}_N^T) Nr(X)=X(X1N1NT) N c = X ⊘ ( 1 N 1 N T X ) \mathcal{N}_c=\textbf{X}\oslash(\textbf{1}_N\textbf{1}_N^T\textbf{X}) Nc=X(1N1NTX)是矩阵行和列元素规范化操作, ⊘ \oslash 是element-wise除法。然后, Mena \text{Mena} Mena等人进一步证明了分配问题能够作为最优运输问题的特例,然后使用 Sinkhorn \text{Sinkhorn} Sinkhorn操作进行求解。
arg max P ∈ P N ⟨ P,X ⟩ F = l i m τ → 0 + Sinkhorn ( X / τ ) \begin{aligned} \mathop{\text{arg max}}_{\textbf{P}\in\mathbb{P}_N}\langle \textbf{P,X} \rangle_F \\ =\mathop{lim}_{\tau\rightarrow0^+}\text{Sinkhorn}(\textbf{X}/\tau) \end{aligned} arg maxPPNP,XF=limτ0+Sinkhorn(X/τ)
一般来说, Sinkhorn \text{Sinkhorn} Sinkhorn操作的时间复杂度是 O ( k n 2 ) O(kn^2) O(kn2)。由于迭代数目 k k k是有限的,实践中 Sinkhorn \text{Sinkhorn} Sinkhorn操作通常能够获得近似解。基于实践的经验,非常小的 k k k就能在实体对齐上实现好的表现。因此, Sinkhorn \text{Sinkhorn} Sinkhorn算法的实际时间复杂度为 O ( n 2 ) O(n^2) O(n2)

五、实现细节

1. 文本特征 H \textbf{H} H

  • Word-level

    先前的实体对齐工作中,最常使用的文本特征是词级别的实体名称向量。具体来说,这些方法会使用机器翻译系统或者交叉语言词嵌入将实体名称映射至统一语义空间,然后平均预训练实体名称向量来构造初始特征。为了公平比较,本文采用相同的实体名称和词向量。

  • Char-level

    由于专有名词之间存在广泛的矛盾以及词表的有限尺寸。单词级别的实体对齐方法遭受严重的 OOV \text{OOV} OOV问题。因此,许多实体对齐方法通过 char-CNN \text{char-CNN} char-CNN或者 name-BERT \text{name-BERT} name-BERT来使用字符级别的特征。为了保持 SEU \text{SEU} SEU的简单性和一致性,采用翻译后实体名称的字符bigram作为输入特征,而不使用复杂的神经网络。

2. 邻接矩阵 A \textbf{A} A

​ 上面所有的推论都是建立在邻接矩阵 A s \textbf{A}_s As A t \textbf{A}_t At是同构的。显然,令 D \textbf{D} D为邻接矩阵 A s / t \textbf{A}_{s/t} As/t的度矩阵,等概率随机游走矩阵 A r = D − 1 A s / t \textbf{A}_r=\textbf{D}^{-1}\textbf{A}_{s/t} Ar=D1As/t和对称归一化 Laplacian \text{Laplacian} Laplacian矩阵 A L = I − D − 1 / 2 A s / t D − 1 / 2 \textbf{A}_L=\textbf{I}-\textbf{D}^{-1/2}\textbf{A}_{s/t}\textbf{D}^{-1/2} AL=ID1/2As/tD1/2也是同构的。因此, A s / t \textbf{A}_{s/t} As/t被替换为 A r \textbf{A}_r Ar或者 A L \textbf{A}_L AL,本文顶点方法也是有效的。然而,上面的矩阵忽略了知识图谱中的关系,即所有类型的关系都同等重要。直觉上,低频关系应该有更高的权重,因为他们能够表示更加独特的信息。依据这个直觉,应用简单的策略来生成关系邻接矩阵 A r e l \textbf{A}_{rel} Arel,对应 a i j ∈ A r e l a_{ij}\in\textbf{A}_{rel} aijArel
a i j = ∑ r j ∈ R i , j ln ( ∣ T ∣ / ∣ T r j ∣ ) ∑ k ∈ N i ∑ r k ∈ R i , k ln ( ∣ T ∣ / ∣ T r k ∣ ) \textbf{a}_{ij}=\frac{\sum_{r_j\in R_{i,j}}\text{ln}(|T|/|T_{r_j}|)}{\sum_{k\in\mathcal{N}_i}\sum_{r_k\in R_{i,k}}\text{ln}(|T|/|T_{r_k}|)} aij=kNirkRi,kln(T/Trk)rjRi,jln(T/Trj)
其中, N i \mathcal{N}_i Ni表示实体 e i e_i ei的邻居集合, R i , j R_{i,j} Ri,j是实体 e i e_i ei e j e_j ej间的关系, ∣ T ∣ |T| T ∣ T r ∣ |T_r| Tr分别表示所有三元组的总数和包含关系 r r r的三元组数量。

六、实验结果

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-YMieznnI-1648347935623)(./图片/SEU实验.png)]

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

BQW_

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值