2

所以现在我正在使用 SAT 来解决最小顶点覆盖问题,这是我对图 G = {V,E} 的编码有 k 个顶点覆盖,下面是子句:

Let n = sizeof(V);

首先,至少一个顶点在顶点覆盖中:

For i in {1..k}
    Add clause (x<1,i> ∨ x<2,i> ∨ ··· ∨ x<n,i>);

那么,任何一个顶点都不能在顶点覆盖中出现两次:

For j in {1..n}
    For l and m in {1..k} with l < m
        Add clause (¬x<j,l> ∨ ¬x<j,m>) 

之后,顶点覆盖中的特定位置只能出现一个顶点:

For j in {1..k}
    For l and m in {1..n} with l < m
        Add clause (¬x<l,j> ∨ ¬x<m,j>) 

最后,顶点覆盖中的至少一个顶点应该来自边:

For i and j in each edge e from E
    Add clause (x<i,1> ∨ x<i,2> ∨ ... ∨ x<i,k> ∨ x<j,1> ∨ ... ∨ x<j,k>)

现在我可以通过使用这种编码获得最小的顶点覆盖,但是效率很差。我只能得到小于 20 个顶点的任何图的结果,否则只需要几分钟和几小时就能得到结果。我现在正在考虑将其从 SAT 进一步降低,也许是 3SAT。但看起来我不能简单地将所有子句从 nCNF 更改为 3CNF 以获得相同的结果。谁能帮我弄清楚下一步该怎么做?我需要全新的编码吗?

太感谢了。

顺便说一句,我正在使用 MiniSAT 作为求解器。

4

1 回答 1

2

首先,让我假设我在理解您的编码时遇到了一些麻烦,所以我将从头开始。这就是我解决问题的方式。

注意:我的示例基于SMT-LIB语法,可以使用MaxSMT求解器求解,例如z3optimathsat。但是,由于它不使用任何 SMT 功能,您实际上可以使用MaxSAT比赛中使用的标准 WCNF 格式编写相同的内容。在选择求解器来处理问题时,这将为您提供更多选择。我可能错了,但我推测MaxSAT求解器在这个特定问题上的表现可能优于MaxSMT求解器。


令 G = {V, E} 为图。

首先,为图中的每个顶点声明一个布尔变量:

(declare-fun vertex_1 () Bool)
...
(declare-fun vertex_K () Bool)

(任何没有边的顶点都应该被省略,因为它只会浪费时间。)

i其次,为连接顶点和顶点的图中的每条边声明一个布尔变量j(假设无向)

(declare-fun edge_i_j () Bool)
...

第三,断言edge_i_j必须覆盖每条边:

(assert edge_i_j)
...

第四,如果一条边edge_i_j被覆盖,那么顶点i或顶点j必须是true

(assert (=> edge_i_j (or vertex_i vertex_j)))
...

第五,对于每个vertex_i断言都带有一个应该是的软子句。如果不是这种情况,则对 的值进行惩罚:vertex_ifalse1cover

(assert-soft (not vertex_i) :weight 1 :id cover)

最后解决问题:

(set-option :config opt.maxsmt_engine=maxres) ; only for optimathsat
(minimize cover)                              ; only for optimathsat
(check-sat)
(get-objectives)
(get-model)

此时,可以使用任何高效的 MaxSAT/MaxSMT引擎(例如MaxRes)来获得一个覆盖所有边的模型,同时使用最少数量的顶点(如果有的话)。

于 2019-12-02T18:18:19.483 回答