带连接范围约束的二部网络保度随机化实现方法咨询

你将该任务建模为二部网络保度随机化问题的思路完全正确——任务的两个核心约束（每个主题对应论文总数恒定、每个出版年份对应论文总数恒定），本质就是二部图左右两侧节点的度值需要严格保持不变，遇到的工具适配问题并非建模逻辑错误导致，而是原生工具的能力边界不覆盖「边连接存在范围约束」的场景。

• 该方法的边生成逻辑是基于度序列的无配对约束随机匹配，完全不校验两个节点是否属于历史合法连接范围，因此必然会生成类似「topic3连接1998年」这类越界边
• 其输出的随机化结果和全局年份置换无本质差异，本质是因为无约束配置模型生成的是「所有满足度约束的二部图空间内的均匀随机样本」，未纳入你定义的topic-年份合法连接边界，不符合你要求的约束场景

目前有成熟的、符合零模型统计要求的实现路径，核心是使用受限边交换法在合法的二部图空间内做随机游走，从原始真实数据构成的二部图出发迭代生成随机样本，不需要从零开发复杂逻辑，甚至可以直接基于networkx原生功能改造实现：

1. 基础准备
   • 先基于原始论文数据构建初始二部图：左侧节点为所有研究主题topic，右侧节点为所有出版年份pub_year，每篇论文对应一条topic与pub_year的连接边。此时每个左侧节点的度就是对应主题的总论文数，每个右侧节点的度就是对应年份的总论文数，完全匹配你的保度要求
   • 预先构建合法连接映射：记录每个topic允许连接的所有pub_year（即你梳理的各主题历史年份覆盖区间），作为边合法性的校验依据
2. 核心随机化逻辑
   改造networkx原生的double_edge_swap边交换流程，把原来的无约束交换替换为带合法性校验的版本：
   • 每次随机选中现有二部图中的两条边（t1-y1、t2-y2），尝试将边重连为（t1-y2、t2-y1）
   • 仅当重连后的两条边都属于合法连接集合（即t1允许连接y2、t2允许连接y1）时，才接受这次交换
   • 重复执行该交换操作足够多次（一般达到总边数的10-20倍迭代量即可保证样本的随机无偏性），最终得到的二部图就是满足所有约束的随机化结果
     这种边交换的实现方式是文献计量领域带约束零模型构建的通用标准方法，统计性质经过大量研究验证，采样得到的是所有满足约束的二部图空间内的无偏均匀样本。
3. 效果保障
   该方法的全程操作不会改变左右两侧节点的度值，同时所有边都严格落在topic对应的年份范围内，完全规避越界问题，同时满足两类核心保度要求。

• 不要尝试直接使用带权重的配置模型：权重逻辑不会限制边的生成范围，仍然会产出越界连接
• 不要采用「topic内部独立打乱年份」的实现方式：这类操作无法保证每个出版年份的总论文数恒定，会破坏你要求的第二类核心约束

内容的提问来源于stack exchange，提问作者Xiancheng LI