跳转至

xVGAE: A Hierarchical Variational Graph Autoencoder for Exchangeable Graphs

会议: NeurIPS 2026
arXiv: 2610.30815
OpenReview: https://openreview.net/forum?id=yfPdA2AWpR
论文: https://arxiv.org/abs/2610.30815
领域: 图学习
关键词: 图生成模型, 交换图, 变分图自编码器, Aldous-Hoover定理, 图元(Graphon)

一句话总结

基于 Aldous-Hoover 交换图表示定理构建双层隐变量变分图自编码器 xVGAE,通过图级隐变量驱动的高斯混合节点先验与双线性图元解码器,实现了单步前向极速且高质量的无向图生成。

研究背景与动机

图生成模型在分子药物发现、社交网络建模及生物蛋白拓扑等领域具有基础性作用。然而,现有图生成方法在面对节点规模达数百个的较稠密图时面临严峻的计算瓶颈:以 GraphRNN、GRAN 为代表的自回归模型需要按序列逐步预测邻接矩阵,破坏了图的节点置换不变性且推理开销随节点数二次膨胀;以 DiGress、SparseDiff、EDGE 等为代表的扩散模型虽保持了置换等变性,但往往需要在成百上千个去噪步长中反复迭代更新所有候选边的显式隐变量,且其隐状态本质上是带噪邻接矩阵,缺乏可解释的语义压缩表征。

这种困境的根源在于:现有模型缺乏一个能够同时保证置换不变性、避免边级显式状态存储并具有完备统计表达能力的紧凑隐变量结构。传统的随机块模型(SBM)与潜在空间模型(LSM)虽然具有天然的置换不变性,但其参数化假设过于僵化;而早期变分图自编码器(VGAE)的非灵活内积解码器也仅适用于局部表征学习,难以生成具有全局拓扑异质性的大图。

Aldous-Hoover 定理从理论上证明了任意联合可交换无向图的生成过程均可充分解耦为:一个全局图级隐变量、一组条件独立同分布的节点级隐变量,以及一个映射隐变量到边连接概率的对称图元(Graphon)函数。本文的核心 idea 是将经典 Aldous-Hoover 交换图表示定理作为深度图生成的归纳偏置,构建分层变分图自编码器 xVGAE,利用轻量全局扩散先验与条件高斯混合节点先验解耦拓扑结构,仅通过双线性图元解码器实现无边级隐变量的极速单步生成。

方法详解

整体框架

xVGAE 包含推断模型(编码器)与生成模型(解码器)。推断阶段,编码器输入图邻接矩阵与拉普拉斯特征,分别推断全局后验分布 \(q_\phi(z_0 \mid A)\) 以及依赖于全局隐变量的节点级后验分布 \(q_\phi(z_i \mid z_0, A)\);生成阶段,首先从低维全局隐扩散先验中采样图级变量 \(z_0\),随后以此为条件从高斯混合先验中独立同分布采样 \(n\) 个节点变量 \(z_i\),最后通过条件双线性解码器一次前向并行计算所有节点对的边概率。

%%{init: {'flowchart': {'rankSpacing': 24, 'nodeSpacing': 28, 'padding': 6, 'wrappingWidth': 400}}}%%
flowchart TD
    A["输入图邻接矩阵 A 与特征"] --> B["分层编码推断<br/>全局后验与条件节点后验"]
    B --> C["分层隐空间先验采样<br/>全局扩散先验与高斯混合节点先验"]
    C --> D["分块双线性图元解码器<br/>FiLM 调制与对称交互矩阵预测"]
    D --> E["生成图邻接概率矩阵 A_hat"]

关键设计

1. 分层编码推断:置换不变与等变的分层后验

针对传统图自编码器无法对整图语义与局部角色进行解耦的痛点,xVGAE 的推断模型遵循生成层次严格解耦为全局因子与节点局部因子: $\(q_\phi(z_0, z_{1:n} \mid A) = q_\phi(z_0 \mid A) \prod_{i=1}^n q_\phi(z_i \mid A, z_0)\)$ 编码器首先通过图拉普拉斯矩阵的前 \(d\) 个特征向量构造初始节点特征,再使用多层集合注意力块(Set Attention Blocks, SAB)捕获图结构依赖。全局编码器利用特征值通过 FiLM 调制 SAB 层,并经由多头注意力池化(PMA)输出置换不变的全局高斯分布参数;局部编码器则在采样出 \(z_0\) 后进行条件调制,输出满足置换等变性的各节点高斯分布参数。这种解耦避免了为每条边引入二次规模的显式隐变量。

2. 分层隐空间先验采样:全局扩散先验与高斯混合节点先验

为了刻画复杂网络中图级结构异质性与节点间的多模态社区聚集特性,xVGAE 在隐空间建立了以全局变量为中介的双层先验结构。对于节点局部先验 \(p_\theta(z_i \mid z_0)\),采用以 \(z_0\) 为条件的 \(K\) 模态高斯混合模型,由 MLP 预测各高斯分量的权重、均值和对角方差,使得节点隐表征天然契合混合隶属社区结构;而对于图级全局先验 \(p_\theta(z_0)\),鉴于全局变量承载了整图全局拓扑的调节职责,模型引入了一个针对低维向量 \(z_0\) 的去噪扩散概率模型(DDPM)。由于扩散过程只在低维全局隐空间执行(维度仅为 32),相比在完整图邻接矩阵或边空间做去噪扩散,其采样速度实现了成百上千倍的飞跃。

3. 分块双线性图元解码器:参数化连续点式 Graphon

根据 Aldous-Hoover 定理,边连接概率函数需关于节点变量严格对称。解码器将图元函数 \(\omega_\theta(z_0, z_i, z_j)\) 参数化为条件分块双线性交互模型: $\(s_{ij} = \sum_{k=1}^{K_h} \alpha_k(z_0) h_i^{(k)\top} W_k(z_0) h_j^{(k)} + b_0(z_0) + c(h_i) + c(h_j), \quad \omega_\theta(z_0, h_i, h_j) = \sigma(s_{ij})\)$ 其中节点特征 \(h_i\) 由共享 MLP 在 \(z_0\) 的 FiLM 特征调制下映射得到(或可选地经过 SAB 细化);交互矩阵 \(W_k(z_0) = W_k^{\text{base}} + W_k^{\text{mod}}(z_0)\) 严格对称,由超网络以低秩分解形式预测。各注意力头针对不相交的特征切块独立运算后加权汇总,既严格保证了节点置换的对称性,又自然推广并涵盖了经典 SBM、MMSBM 与内积图模型。

损失函数 / 训练策略

训练采用两阶段解耦策略:第一阶段固定 \(p(z_0) = \mathcal{N}(0, I)\),端到端优化双层变分下界(ELBO): $\(\mathcal{L}(\theta,\phi) = \mathbb{E}_{q_\phi}\left[\log p_\theta(A \mid z_0, z_{1:n})\right] - \beta_0 \text{KL}\left(q_\phi(z_0 \mid A) \parallel p(z_0)\right) - \beta \mathbb{E}_{q_\phi(z_0 \mid A)}\left[\text{KL}\left(q_\phi(z_{1:n} \mid A, z_0) \parallel p_\theta(z_{1:n} \mid z_0)\right)\right]\)$ 其中重构损失计算邻接矩阵上三角元素交叉熵,局部 KL 散度采用蒙特卡洛无偏估计,避免显式采样高斯混合离散分配;\(\beta\) 采用余弦退火策略从 \(10^{-4}\) 逐步增至 1。第二阶段冻结自编码器,使用编码的训练集 \(z_0\) 作为真实目标,以标准均方误差损失训练低维全局潜扩散网络 \(\mathcal{L}_{\text{DM}} = \mathbb{E}[\|\epsilon - \epsilon_\theta(z_0^{(t)}, t)\|_2^2]\)。

实验关键数据

主实验

评估在统一计算预算(单张 48GB GPU,训练限时 12 小时,采样限时 10 小时)下进行,采用偏置高斯 TV-MMD\(^2\) 比率(除以训练-测试经验基线,接近 1.0 为优)、WL 子树核 MMD\(^2\)、分类器双样本检验(C2ST\(_\Delta \in [0, 0.5]\),越接近 0 越优)以及批处理单图生成时间(s/G)。

数据集 模型 度分布 MMD\(^2\) ↓ 聚类系数 MMD\(^2\) ↓ 谱分布 MMD\(^2\) ↓ WL MMD\(^2\) ↓ C2ST\(_\Delta\) ↓ 生成速度 (s/G) ↓
MMSBM (稠密社区) GRAN 2.07 1.07 1.24 1.35 0.037 4.15
EDGE 1.86 2.32 64.3 1.27 0.037 0.16
DiGress 74.3 2.78 846 26.4 0.356 18.2
ILE OOM OOM OOM OOM OOM OOM
xVGAE-MLP (本文) 2.27 1.05 0.44 1.28 0.032 <0.02
xVGAE-SAB (本文) 2.53 1.01 1.86 1.34 0.062 <0.02
CiteSeer-Ego (引文网络) GraphRNN 3.74 26.8 16.2 5.11 0.464 0.14
DiGress 10.9 4.13 3.57 4.41 0.217 7.05
EDGE 18.0 6.78 14.5 7.08 0.148 0.011
ILE 21.7 5.85 5.90 8.98 0.299 55.7
xVGAE-SAB (本文) 8.30 5.20 4.97 3.85 0.102 <0.01
LFR (幂律社区) DiGress 10.8 3.61 5.22 8.64 0.103 33.0
EDGE 3.88 3.16 3.78 3.65 0.060 0.05
xVGAE-SAB (本文) 19.5 2.50 1.17 8.38 0.165 <0.01

消融实验

针对全局扩散先验与解码器注意力机制进行消融分析:

数据集 / 变体 聚类 MMD\(^2\) ↓ 谱分布 MMD\(^2\) ↓ WL MMD\(^2\) ↓ C2ST\(_\Delta\) ↓ 说明
MMSBM xVGAE-MLP 1.05 0.44 1.28 0.032 完整模型(带 \(z_0\) 扩散先验)
MMSBM w/o 扩散 (\(z_0 \sim \mathcal{N}(0, I)\)) 1.03 0.96 1.47 0.040 去除扩散先验后,各分布统计指标出现轻微退化
MMSBM xVGAE-SAB 1.01 1.86 1.34 0.062 引入自注意力细化,在纯块模型上反而轻微过拟合
LFR xVGAE-SAB 2.50 1.17 8.38 0.165 包含自注意力细化,有效拟合幂律度分布与复杂混合率
LFR xVGAE-MLP 3.55 3.28 9.41 0.153 去除自注意力后,聚类与谱指标明显变差
LFR w/o 扩散先验 5.39 1.46 7.70 0.190 失去全局扩散先验调控后,全局可分性显著降低

关键发现

  • 生成效率呈压倒性优势:由于 xVGAE 将扩散过程限制在 32 维的全局向量 \(z_0\),节点与边概率矩阵可并行单步生成,生成耗时均小于 0.02 秒/图;相比 DiGress 加速超 3,000 倍,相比基于高阶消息传递的 ILE 加速超 5,000 倍,且彻底杜绝了 ILE 在稠密图上的显存溢出(OOM)崩溃。
  • 解码器结构的归纳偏置适配性:在纯块结构清晰的 MMSBM 稠密图上,点式独立 MLP 解码器已足以达到最优性能(谱 MMD\(^2\) 达 0.44);而在度分布呈重尾幂律的 LFR 网络中,SAB 集合注意力解码器通过捕获高阶依赖显著提升了聚类与谱特征的拟合质量。
  • 无监督表征的可解释性:PCA 投影分析表明,全局变量 \(z_0\) 自发编码了社区总数、混合率、平均度等宏观属性;而局部变量 \(z_i\) 则精准对应了节点的社区归属及在 Ego 网络中距核心节点的跳数距离。

亮点与洞察

  • 扎实的理论奠基:以联合可交换阵列的 Aldous-Hoover 定理为指导原则设计深度生成架构,突破了早期 VGAE 无法表达复杂宏观混合拓扑的理论壁垒,为图隐空间解耦提供了严格而优雅的数学范式。
  • 将扩散限定在紧凑隐流形:完全规避了在节点对 \(n \times n\) 边空间运行昂贵多步扩散去噪的传统模式,使得生成兼具扩散模型的高生成质量与单步自编码器的高速并发优势。

局限与展望

  • 极端链状不规则几何结构建模受限:在空间几何链状主导的真实生物蛋白图(Protein D&D)上,xVGAE 的聚类和谱指标落后于自回归 GRAN 与局部迭代扩展 ILE,表明将高度各向异性的局部不规则三维折叠结构强行压缩为点式隐向量仍存在保真度瓶颈。
  • 未来拓展方向:当前设计仅专注于无向无权简单图,理论框架可自然扩展至部分可交换阵列(处理有向图、二分图),并可通过拓展图元解码维度支持边权重与多维节点属性。

相关工作与启发

  • vs. VGAE / GraphVAE: 经典 VGAE 依赖不可学习的内积解码器且缺少全局隐变量,无法建模图间拓扑漂移;GraphVAE 严重依赖图匹配,难以扩展到超百节点的图。xVGAE 通过全局-局部双层分层与分块双线性解码器彻底解决了可扩展性与表达力矛盾。
  • vs. 图扩散模型 (DiGress / EDGE / LGD): 现有图扩散模型显式维护边级别的高维带噪变量,采样极慢且训练耗费巨大显存。xVGAE 仅在低维全局空间运行扩散,解码单步完成,推理吞吐量提升了数个数量级。

评分

  • 新颖性: ⭐⭐⭐⭐⭐ 首次以 Aldous-Hoover 定理系统性指导百节点级深度图生成并成功结合隐扩散。
  • 实验充分度: ⭐⭐⭐⭐⭐ 覆盖稠密、幂律、引文与生物四大结构体系,全面对比 12 类基线并在统一算力预算下完成严谨复现。
  • 写作质量: ⭐⭐⭐⭐⭐ 统计理论阐述透彻,实验消融与可解释性隐空间可视化论证详实。
  • 价值: ⭐⭐⭐⭐⭐ 为解决图生成模型面临的效率瓶颈提供了兼具统计理论根基与极高工程实用性的通用范式。