如何在Zimpl中定义有向图的α弧集函数及相关约束？

嘿，我来帮你搞定这个Zimpl里的集合定义问题！其实Zimpl本身就很擅长处理这类集合筛选的逻辑，不用写复杂的函数，用集合推导或者宏就能轻松实现你要的α(W, W')。

第一步：先声明基础集合与变量 #

首先得把你的顶点集、弧集和变量先定义好，这是基础：

// 定义总顶点集V，以及它的两个子集V₁、V₂
set V;
set V1 subsetof V;
set V2 subsetof V;

// 定义有向弧集A，每条弧是V中顶点的有序对(u,v)
set A subsetof V cross V;

// 定义每条弧对应的变量x_{(u,v)}，这里假设是非负整数，你可以根据需求调整类型
var x{A} integer >= 0;

第二步：直接用集合推导表示α(W, W') #

α(W, W')本质是弧集A中所有起点在W、终点在W'的弧，在Zimpl里可以直接用集合推导式在表达式中筛选，不需要单独定义函数。比如你需要的α(V₁, V₂)可以直接写成：

{(u,v) in A | u in V1, v in V2}

对应的约束就可以直接写：

// 约束：从V₁到V₂的弧变量之和≥2
constraint MinFlowV1ToV2:
    sum{(u,v) in A where u in V1 and v in V2} x[(u,v)] >= 2;

第三步：用宏封装重复逻辑（可选） #

如果你的问题里需要多次用到α(W, W')这个逻辑（比如针对不同的W和W'子集），可以定义一个宏来封装，让代码更简洁易读：

// 定义宏alpha，接收两个集合参数W和W'，返回对应的弧集合
macro alpha(W, W') = { (u,v) in A | u in W, v in W' };

之后调用宏就能直接生成目标弧集，约束可以写成：

constraint MinFlowV1ToV2:
    sum{(u,v) in alpha(V1, V2)} x[(u,v)] >= 2;

额外提示 #

• 如果V₁和V₂是不相交的（题目里V=V₁∪V₂，通常默认不相交），可以在声明子集时加上disjoint关键字：set V1 subsetof V disjoint from V2;，这样能避免逻辑错误。
• Zimpl的集合推导支持多种筛选条件，你可以根据需求扩展α的定义，比如加上弧的权重限制之类的。

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