论文解读（SEP）《Structural Entropy Guided Graph Hierarchical Pooling》

论文信息

论文标题：Structural Entropy Guided Graph Hierarchical Pooling
论文作者：Junran Wu, Xueyuan Chen, Ke Xu, Shangzhe Li
论文来源：2022,ICLR
论文地址：download 
论文代码：download

1 Introduction 

　　之前池化的方法大多是 固定比例的池化 或者 基于逐步池化设计的。本文针对上述问题，设计了 SEP 模型，具体来说，在不分配特定于层的压缩配额的情况下，设计了一种全局优化算法来一次生成针对池化的簇分配矩阵。

　　然后，本文举例说明了以往的方法在重建环图和网格合成图中的局部结构损伤。除 SEP 外，还进一步设计了两个分类模型，SEP-G 和 SEP-N，分别用于图分类和节点分类。结果表明，SEP 在图分类基准上优于最先进的图池方法，并在节点分类上获得了更好的性能。

　　池化有两种代表性的工作：

• • 基于节点丢弃的池化策略：TopKPool、SAGPool、ASAP

• • 基于节点聚类的层次池化：DiffPool

　　对于第一种池化策略：由于不考虑池化过程中的消息聚合，不可避免的造成了信息丢失；对于第二种：虽然使用聚类思想可以有效减少信息丢失，但是此类方法目前还不是很完善。上述两种池化策略由于需要人工指定节点压缩比例，所以均存在对局部结构信息的破坏。

　　本文贡献：

• • 我们在以往的分层池工作中发现了阻碍 GNN 发展的两个关键问题，包括由于固定的压缩配额和逐步池设计而造成的局部结构损伤和次优问题。
  • 通过引入结构信息理论，我们提出了一种新的分层池化方法，称为SEP，以解决所揭示的问题。
  • 我们在图重建、图分类和节点分类任务上广泛地验证了SEP，与SOTA分层池化方法相比，它们的性能更好。

2 Proposed Method

　　整体框架图示：

　　

2.1 Preliminaries

Graph neural networks

　　　　Hi+1=Re⁡LU(D~−12A~D~−12HiWi)" role="presentation">Hi+1=ReLU(D~−12A~D~−12HiWi)$H_{i+1}=\mathrm{Re}LU\left({\tilde{\mathbf{D}}}^{-\frac{1}{2}}\tilde{\mathbf{A}}{\tilde{\mathbf{D}}}^{-\frac{1}{2}}{H}_i{W}_i\right)$

Hierarchical pooling

　　假设 Si∈Rni+1×ni" role="presentation">Si∈Rni+1×ni${\mathbf{S}}_i\in {\mathbb{R}}^{n_{i+1}\times n_i}$ 代表第 i" role="presentation">i$i$ 层的聚类分配矩阵，那么池化后的新的邻接矩阵和节点特征矩阵由下一个方程式计算：

　　　　Ai+1=SiAiSi⊤;Pi+1=SiHi" role="presentation">Ai+1=SiAiS⊤i;Pi+1=SiHi${\mathbf{A}}_{i+1}={\mathbf{S}}_i{\mathbf{A}}_i{\mathbf{S}}_i^{\mathrm{\top }};{\mathbf{P}}_{i+1}={\mathbf{S}}_i{H}_i$

　　其中：

• • Ai∈Rni×ni" role="presentation">Ai∈Rni×ni${\mathbf{A}}_i\in {\mathbb{R}}^{n_i\times n_i}$ 代表第 i" role="presentation">i$i$ 层的邻接矩阵；
  • Hi∈Rni×h" role="presentation">Hi∈Rni×h$H_i\in {\mathbb{R}}^{n_i\times h}$ 代表第 i" role="presentation">i$i$ 层由 GNN 导出的特征矩阵；
  • Pi+1" role="presentation">Pi+1${\mathbf{P}}_{i+1}$ 接收节点隐藏特征 Hi" role="presentation">Hi$H_i$ 并合并这些特征，以改进新图中 ni+1" role="presentation">ni+1$n_{i+1}$ 簇的初始表示；

2.2 Cluster Assignments via Structural Entropy Minimization

　　结构熵用来解码图的层次结构，通过结构熵最小化，可以将图的层次结构解码为相应的编码树，使噪声或随机变化产生的干扰最小化。本文认为一种有效的结构熵最小化算法可以揭示层次池化层之间的联系，并消除图中的噪声结构。

　　编码树 T" role="presentation">T$T$ 上 G" role="presentation">G$G$ 的结构熵的形式化方程可以写为：

　　　　HT(G)=−∑vt∈Tgvtvol⁡(V)log⁡vol⁡(vt)vol⁡(vt+)(3)" role="presentation">HT(G)=−∑vt∈Tgvtvol(V)logvol(vt)vol(v+t)(3)${\mathcal{H}}^T(G)=-\sum _{v_t\in T}\frac{g_{v_t}}{\mathrm{vol}(\mathcal{V})}\log \frac{\mathrm{vol}\left(v_t\right)}{\mathrm{vol}\left(v_t^{+}\right)}(3)$

　　其中，

• • vt" role="presentation">vt$v_t$ 是编码树 T" role="presentation">T$T$ 的一个非根节点，同时可以被认为是一个根据编码树 T" role="presentation">T$T$ 划分的节点子集 ⊂V" role="presentation">⊂V$\subset \mathcal{V}$；
  • vt+" role="presentation">v+t$v_t^{+}$ 是 vt" role="presentation">vt$v_t$ 的父节点；
  • gvt" role="presentation">gvt$g_{v_t}$ 是指在 vt" role="presentation">vt$v_t$ 的叶节点划分中具有端点的边数；
  • vol⁡(V)" role="presentation">vol(V)$\mathrm{vol}(\mathcal{V})$ 代表 V" role="presentation">V$\mathcal{V}$ 中的叶子节点度的和；
  • vol⁡(vt)" role="presentation">vol(vt)$\mathrm{vol}\left(v_t\right)$ 代表 vt" role="presentation">vt$v_t$ 节点度；

　　编码树是图的一种自然的层次划分，为了使结构熵最小化 H(G)=min∀T{HT(G)}" role="presentation">H(G)=min∀T{HT(G)}$\mathcal{H}(G)=\underset{\mathrm{\forall }T}{min}\{{\mathcal{H}}^T(G)\}$，本文建立了不同层之间的连接。此外，图中的局部结构将被保留，因为不需要分配图层中特定的节点压缩配额。

　　首选具有一定高度的自然编码树，因为大多数任务只需要几次固定时间的图粗化。在此背景下，提出了 G" role="presentation">G$G$ 的 k" role="presentation">k$k$ 维结构熵来解码具有固定高度 k" role="presentation">k$k$ 的最优编码树：

　　　　H(k)(G)=min∀T:Height⁡(T)=k{HT(G)}(4)" role="presentation">H(k)(G)=min∀T:Height(T)=k{HT(G)}(4)${\mathcal{H}}^{(k)}(G)=\underset{\mathrm{\forall }T:\mathrm{Height}(T)=k}{\text{min}}\{{\mathcal{H}}^T(G)\}(4)$

　　本文在 k" role="presentation">k$k$ 维结构熵的指导下，研究了对具有一定高度 k" role="presentation">k$k$ 的编码树进行解密的求解方法。首先，我们定义了三个函数。

　　Definition 3.1. Let  T  be any coding tree for graph  G=(V,E)" role="presentation">G=(V,E)$G=(\mathcal{V},\mathcal{E})$, vr" role="presentation">vr$v_r$  is the root node of  T" role="presentation">T$T$  and  V" role="presentation">V$\mathcal{V}$  are the leaf nodes of  T" role="presentation">T$T$ . Given any two nodes  (vi,vj)" role="presentation">(vi,vj)$(v_i,v_j)$  in  T" role="presentation">T$T$ , in which  vi∈vr" role="presentation">vi∈vr$v_i\in v_r$.children and  vj∈vr" role="presentation">vj∈vr$v_j\in v_r$.children.

　　Define a function MERGET⁡(vi,vj)" role="presentation">MERGET(vi,vj)${\mathrm{MERGE}}_T(v_i,v_j)$ for T" role="presentation">T$T$ to insert a new node vε" role="presentation">vε$v_{\epsilon }$ between vr" role="presentation">vr$v_r$ and (vi,vj)" role="presentation">(vi,vj)$(v_i,v_j)$ :

　　　　vε⋅ children ←vi(5)vε⋅ children ←vj(6)vr. children ←vε(7)" role="presentation">vε⋅ children ←vi(5)vε⋅ children ←vj(6)vr. children ←vε(7)$\begin{array}{l}{v}_{\epsilon }\cdot \text{ children }\leftarrow v_i(5)\\ v_{\epsilon }\cdot \text{ children }\leftarrow v_j(6)\\ v_r.\text{ children }\leftarrow v_{\epsilon }(7)\end{array}$

　　Definition 3.2. Following the setting in Definition 3.1, given any two nodes (vi,vj)" role="presentation">(vi,vj)$(v_i,v_j)$ , in which vi∈vj.children" role="presentation">vi∈vj.children$v_i\in v_j.children$.

　　Define a function REMOVET⁡(vi)" role="presentation">REMOVET(vi)${\mathrm{REMOVE}}_T\left(v_i\right)$ for T" role="presentation">T$T$ to remove vi" role="presentation">vi$v_i$ from T" role="presentation">T$T$ and merge vi.children" role="presentation">vi.children$v_i.children$ to vj.chileren" role="presentation">vj.chileren$v_j.chileren$:

　　　　vj.children ←vi.children (8)" role="presentation">vj.children ←vi.children (8)$v_j\text{.children }\leftarrow v_i\text{.children }(8)$

　　Definition 3.3. Following the setting in Definition 3.1, given any two nodes (vi,vj)" role="presentation">(vi,vj)$(v_i,v_j)$ , in which vi∈vj.children" role="presentation">vi∈vj.children$v_i\in v_j.children$ and ∣Heigth(vj)−Height(vi)∣>1" role="presentation">∣Heigth(vj)−Height(vi)∣>1$\mid Heigth\left(v_j\right)-Height\left(v_i\right)\mid >1$ .

　　Define a function FILL⁡(vi,vj)" role="presentation">FILL(vi,vj)$\mathrm{FILL}(v_i,v_j)$ for T" role="presentation">T$T$ to insert a new node vε" role="presentation">vε$v_{\epsilon }$ between vi" role="presentation">vi$v_i$ and vj" role="presentation">vj$v_j$ :

　　　　vε.children ←vi;(9)vj.children ←vε;(10)" role="presentation">vε.children ←vi;(9)vj.children ←vε;(10)$\begin{array}{l}{v}_{\epsilon }\text{.children }\leftarrow v_i;(9)\\ v_j\text{.children }\leftarrow v_{\epsilon };(10)\end{array}$

　　使用上述 3 个 Definition 构造一个 高度为 k" role="presentation">k$k$ 的编码树：

　　

Algorithm 1

　　Stage 1: Bottom to top construction —— 从上到下 构造一棵二进制编码树，将根的两个节点和并乘一个新的划分，最大限度的减少结构熵；

　　Stage 2: Compress T" role="presentation">T$T$ to the certain height k" role="presentation">k$k$ —— 为满足固定数量的图粗化，需要通过删除节点来压缩之前的全高二进制编码树。每次从 T" role="presentation">T$T$ 中选择一个内部节点，这使得 T" role="presentation">T$T$ 在去除该节点后的结构熵变化最小。

　　Stage 3: Fill T to avoid cross-layer links —— 在第二阶段结束时，已经在结构熵的指导下得到了一个一定高度为 k" role="presentation">k$k$ 的编码树。但是，由于跨层链接，节点在下一层中没有直接后继，这将导致在基于这种编码树实现分层池时节点丢失。因此，我们需要执行第三个阶段，以确保层间信息传输的完整性，也不需要干扰 G" role="presentation">G$G$ 在编码树 T" role="presentation">T$T$ 上的结构熵。

　　最终，G" role="presentation">G$G$ 的编码树 T" role="presentation">T$T$ 被获得，其中 T=(VT,ET)" role="presentation">T=(VT,ET)$T=({\mathcal{V}}^T,{\mathcal{E}}^T)$，VT=(V0T,…,VkT)" role="presentation">VT=(VT0,…,VTk)${\mathcal{V}}^T=({\mathcal{V}}_0^T,\dots ,{\mathcal{V}}_k^T)$ 和V0T=V" role="presentation">VT0=V${\mathcal{V}}_0^T=\mathcal{V}$。此外，从 ET" role="presentation">ET${\mathcal{E}}^T$ 中也可以得到簇分配矩阵，即 S=(S1,…,Sk)" role="presentation">S=(S1,…,Sk)$\mathbb{S}=({\mathbf{S}}_1,\dots ,{\mathbf{S}}_k)$。

　　Proposition 3.4. Let T be a coding tree after the second stage of Algorithm 1, and given two adjacent nodes (vi,vj)inT" role="presentation">(vi,vj)inT$(v_i,v_j)inT$ , in which vi∈vj.children" role="presentation">vi∈vj.children$v_i\in v_j.children$ and ∣Heigth⁡(vj)−Height⁡(vi)∣>1" role="presentation">∣Heigth(vj)−Height(vi)∣>1$\mid Heig\mathrm{th}\left(v_j\right)-\mathrm{Height}\left(v_i\right)\mid >1$ . Then, HT(G)=HTFILL(vi,vj)(G)" role="presentation">HT(G)=HTFILL(vi,vj)(G)${\mathcal{H}}^T(G)={\mathcal{H}}^{T_{FILL(v_i,v_j)}}(G)$.

　　Proof：

　　

2.3 Graph Neural Network for Graph Classification

　　图分类框架如下：

　　

　　图表示如下：

　　　　hG=Concat⁡(Readout⁡(SEPi(GCNi(Hi,Ai),Si))∣∀i∈{1,2,3})(14)" role="presentation">hG=Concat(Readout(SEPi(GCNi(Hi,Ai),Si))∣∀i∈{1,2,3})(14)$h_G=\mathrm{Concat}(\mathrm{Readout}\left(SE{P}_i(GC{N}_i(H_i,{\mathbf{A}}_i),{\mathbf{S}}_i)\right)\mid \mathrm{\forall }i\in \{1,2,3\})(14)$

　　Proposition 3.5. Given a permutation matrix  P∈{0,1}n×n" role="presentation">P∈{0,1}n×n$\mathcal{P}\in \{0,1{\}}^{n\times n}$ , then  SEP⁡(A,H)=SEP⁡(PAP⊤,PH)" role="presentation">SEP(A,H)=SEP(PAP⊤,PH)$\mathrm{SEP}(A,H)=\mathrm{SEP}(\mathcal{P}A{\mathcal{P}}^{\mathrm{\top }},\mathcal{P}H)$  (i.e., SEP is permutation invariant).

2.4 Graph Neural Network for Node Classification

　　节点分类框架如下：

　　

　　除了将图转换为高级表示的功能外，我们还可以采用相同的矩阵 Si" role="presentation">Si$S_i$ 来展开压缩图表示 Hi" role="presentation">Hi$H_i$，并将 Ai" role="presentation">Ai$A_i$ 结构到原始空间：

　　　　Ai+1=Si⊤AiSi;Pi+1=Si⊤Hi(15)" role="presentation">Ai+1=S⊤iAiSi;Pi+1=S⊤iHi(15)${\mathbf{A}}_{i+1}={\mathbf{S}}_i^{\mathrm{\top }}{\mathbf{A}}_i{\mathbf{S}}_i;{\mathbf{P}}_{i+1}={\mathbf{S}}_i^{\mathrm{\top }}{H}_i(15)$

3 Experiments

3.1 Graph Reconstruction

　　在介绍 GNN 中两个主要分类任务的实验之前，我们首先采用图重构实验，对池化层保留的结构信息进行量化，直接揭示以往的分层池化方法对图的局部结构造成的损害。

Confifiguration

　　训练一个自动编码器，用池化层和非池化层重建输入图。学习目标是最小化输入特征 X" role="presentation">X$\mathbf{X}$ 和输出特征 Xr" role="presentation">Xr${\mathbf{X}}^r$ 之间的均方误差(MSE)，即 min‖X−Xr‖2" role="presentation">min∥X−Xr∥2$min{\Vert \mathbf{X}-{\mathbf{X}}^r\Vert }^2$。对于配置，我们使用合成图，包括一个环和一个网格，其中输入特征是二维欧几里得空间中节点的坐标。详细的实验配置和模型描述见附录A.1。

Reconstruction results

　　Figure 4 展示了使用原始图和使用不同池化方法重构的图：

　　

　　我们选择 TopKPool 来表示节点丢弃设计的方法，选择 MinCutPool 来表示节点聚类设计的方法。我们首先注意到 TopKPool 的重建结果，其中原始图的基本形状甚至无法识别。这证实了节点下降方法会导致严重的信息缺失，并意味着节点下降方法在图分类中的性能较差。另一方面，MinCutPool 确实保留了原始图形的基本形状。但是，我们仍然可以看到环的边缘和网格中心的显著扭曲，这代表了环和网格的关键结构，这验证了图中的局部结构会被人为指定的节点压缩配额破坏的假设。相反地，SEP 几乎重建了环，并保留了网格中心的基本结构，这表明我们的池化方法获得了原始图的关键结构信息。

3.2 Graph Classifification

图分类结果：

　　

Visualization case study

　　

Visualization case study

聚类可视化

　　

3.3 Node Classification

　　

Node classifification accuracies with different network depths

　　

 

4 Conclusions

　　在本文中，我们开发了一个优化算法来解决现有的分层池方法的几个局限性。特别是，在结构熵最小化的指导下，我们的池化方法SEP不仅可以捕获池化层之间的连通性，而且还可以解决由于节点压缩的超参数而导致的局部结构破坏的问题。在提出的SEP基础上，我们引入了两种学习模型，SEP-G和SEP-n，分别用于图分类和节点分类。实验结果表明，SEP-G在图分类方面取得了显著的改进，并且SEP-N在节点分类任务上取得了优于其他gnn的性能。结构熵和节点特征的结合是未来工作的一个有趣的方向。

• 论文信息
• 1 Introduction 
• 2 Proposed Method
•     2.1 Preliminaries
•     2.2 Cluster Assignments via Structural Entropy Minimization
•     2.3 Graph Neural Network for Graph Classification
•     2.4 Graph Neural Network for Node Classification
• 3 Experiments
•     3.1 Graph Reconstruction
•     3.2 Graph Classifification
•     3.3 Node Classification
• 4 Conclusions

__EOF__