弦图

弦图是一种特殊的图,很多在一般图上的 NP-Hard 问题在弦图上都有优秀的线性时间复杂度算法.

一些定义与性质

子图:点集和边集均为原图点集和边集子集的图.

导出子图(诱导子图):点集为原图点集子集,边集为所有满足 两个端点均在选定点集中 的图.

团:完全子图.

极大团:不是其他团子图的图.

最大团:点数最大的团.

团数:最大团的点数,记为 ω(G)\omega(G).

最小染色:用最少的颜色给点染色使得所有边连接的两点颜色不同.

色数:最小染色的颜色数,记为 χ(G)\chi(G).

最大独立集:最大的点集使得点集中任意两点都没有边直接相连.该集合的大小记为 α(G)\alpha(G).

最小团覆盖:用最少的团覆盖所有的点.使用团的数量记为 κ(G)\kappa(G).

弦:连接环中不相邻两点的边.

弦图:任意长度大于 33 的环都有一个弦的图称为弦图.

Lemma 1:团数 ω(G)≤χ(G)\omega(G)\le \chi(G) 色数

证明:考虑单独对最大团的导出子图进行染色,至少需要 ω(G)\omega(G) 种颜色.

Lemma 2:最大独立集数 α(G)≤κ(G)\alpha(G)\le \kappa(G) 最小团覆盖数

证明:每个团中至多选择一个点.

Lemma 3:弦图的任意导出子图一定是弦图.

证明:如果弦图有导出子图不是弦图,说明在这个导出子图上存在大于 33 的无弦环,则无论原图如何(怎么加边)都不会使得原图是弦图,矛盾.

Lemma 4:弦图的任意导出子图一定不可能是一个点数大于 33 的环.

证明:一个点数大于 33 的环不是弦图,用以上定理即可.

弦图的判定

问题描述

给定一个无向图,判断其是否为弦图.

点割集

对于图 GG 上的两点 u,vu,v,定义这两点间的 点割集 为满足删除这一集合后,u,vu,v 两点之间不连通.如果关于 u,vu,v 两点间的一个点割集的任意子集都不是点割集,则称这个点割集为 极小点割集.

Lemma 5:图关于 u,vu,v 的极小点割集将原图分成了若干个连通块,设包含 uu 的连通块为 V1V_1,包含 vv 的连通块为 V2V_2,则对于极小点割集上的任意一点 aa,N(a)N(a) 一定包含 V1V_1 和 V2V_2 中的点.

证明:若 N(a)N(a) 只包含 V1V_1 或 V2V_2 中的至多一个连通块中的点,从点割集中删去 aa 点,仍不连通,则原点割集不是最小点割集.

Lemma 6:弦图上任意两点间的极小点割集的导出子图一定为一个团.

证明:极小点割集大小 ≤1\le 1 时,导出子图一定为一个团.

否则,设极小点割集上有两点为 x,yx,y,由 Lemma 5 得,N(x)N(x) 中有 V1,V2V_1,V_2 中的点,设为 x1,x2x_1,x_2,同样的,设 y1,y2y_1,y_2,注意,可能有 x1=y1,x2=y2x_1=y_1,x_2=y_2.

由于 V1,V2V_1,V_2 均为连通块,则在 x1,y1x_1,y_1 和 x2,y2x_2,y_2 两个点对之间存在最短路径.设 x,yx,y 在 V1,V2V_1,V_2 内部的最短路为 x−x1∼y1−y,x−x2∼y2−yx-x_1\sim y_1-y,x-x_2\sim y_2-y,则图上存在一个环 x−x1∼y1−y−y2∼x2−xx-x_1\sim y_1-y-y_2\sim x_2-x,该环的大小一定 ≥4\ge 4,根据弦图的定义,此时该环上一定存在一条弦.

若这条弦连接了 V1,V2V_1,V_2 两个连通块,则点集不是点割集.若这条弦连接了单个连通块内部的两个点或一个连通块内部的一个点和一个点割集上的点,都不满足最短路的性质.所以这条弦只能连接 x,yx,y 两点.

由此,可证弦图中每个极小点割集中的两点都有边直接相连,故性质得证.

单纯点

设 N(x)N(x) 表示与点 xx 相邻的点集.若点集 {x}+N(x)\{x\}+N(x) 的导出子图为一个团,则称点 xx 为单纯点.

Lemma 7:任何一个弦图都至少有一个单纯点,不是完全图的弦图至少有两个不相邻的单纯点.

证明:数学归纳法.单独考虑每一连通块.

归纳基底:当图与完全图同构时,图上任意一点都是单纯点.当图的点数 ≤3\le 3 时,引理成立.

若图上的点数 ≥4\ge 4 且图不为完全图,可知必然存在 u,vu,v 使得 (u,v)∉E(u,v)\notin E.设 II 是图关于 u,vu,v 的极小点割集.设 A,BA,B 分别是删去 II 后的导出子图上 u,vu,v 所在的连通块.由于问题的对称性,我们只考虑 AA 一侧的情况,设 L=A+IL=A+I.若 LL 为完全图,则 uu 为单纯点;若不是,因为 LL 是原图的导出子图,一定也是弦图,所以有两个不相邻的单纯点,因为 II 是一个团,其上两点都相邻,所以 AA 中一定有一个单纯点.该单纯点扩展到全图也为单纯点.

由于每次将整个图分成若干个连通块证明,大小一定减小,且都满足性质,故归纳成立.

完美消除序列

令 n=∣V∣n=|V|,完美消除序列 v1,v2,…,vnv_1,v_2,\ldots ,v_n 为 1,2,…,n1,2,\ldots ,n 的一个排列,满足 viv_i 在 {vi,vi+1,…,vn}\{v_i,v_{i+1},\ldots ,v_n\} 的导出子图中为单纯点.

Lemma 8:一个无向图是弦图当且仅当其有一个完全消除序列.

充分性:点数为 11 的弦图有完全消除序列.由 Lemma 3 和 Lemma 7,点数为 nn 的弦图的完美消除序列可以由点数为 n−1n-1 的弦图的完美消除序列加上一个单纯点得到.

必要性:假设有无向图存在结点数 >3>3 的环且拥有完美消除序列,设在完美消除序列中出现的第一个环上的点为 vv,设 vv 在环上与 v1,v2v_1,v_2 相连,则有完美消除序列的性质即单纯点的定义可得 v1,v2v_1,v_2 直接有边相连,矛盾.

朴素算法

每次找到一个 单纯点 vv,加入到完美消除序列中.

将点 vv 与其相邻的边从图上删除.

重复以上过程,若所有点都被删除,则原图是弦图且求得了一个完美消除序列;若图上不存在单纯点,则原图不是弦图.

时间复杂度 O(n4)O(n^4).

MCS 算法

最大势算法(Maximum Cardinality Search)是一种可以在 O(n+m)O(n+m) 的时间复杂度内求出无向图的完美消除序列的方法.

逆序给结点编号,即按从 nn 到 11 的顺序给点标号.

设 labelxlabel_x 表示第 xx 个点与多少个已经标号的点相邻,每次选择 labellabel 值最大的未标号结点进行标号.

用链表维护对于每个 ii,满足 labelx=ilabel_x=i 的 xx.

由于每条边对 ∑i=1nlabeli\sum_{i=1}^n label_i 的贡献最多是 22,时间复杂度 O(n+m)O(n+m).

正确性证明:

设 α(x)\alpha(x) 为 xx 在这个序列中的位置. 我们需要证明对于任何一个弦图,算法求出的序列一定是一个完美消除序列,即在序列中位于某个点后面且与这个点相连的所有点两两相连.

Lemma 9:考虑三个点 u,v,wu,v,w 满足 α(u)<α(v)<α(w)\alpha(u)<\alpha(v)<\alpha(w),如果 uwuw 相连,vwvw 不相连,则 ww 只给 uu 的 labellabel 贡献,不给 vv 贡献.为了让 vv 比 uu 先加入序列,需要一个 xx 满足 α(v)<α(x)\alpha(v)<\alpha(x) 且 vxvx 相连,uxux 不相连,即 xx 只给 vv 贡献而不给 uu 贡献.

Lemma 10:任意一个弦图一定不存在一个序列 v0,v1,…,vk(k≥2)v_0,v_1,\dots,v_k(k\ge 2) 满足下列性质:

  1. vivjv_iv_j 相连当且仅当 ∣i−j∣=1|i-j|=1.
  2. α(v0)>α(vi)(i∈[1,k])\alpha(v_0)>\alpha(v_i)(i\in[1,k]).
  3. 存在 i∈[1,k−1]i\in[1,k-1],满足 α(vi)<α(vi+1)<⋯<α(vk)\alpha(v_i)<\alpha(v_{i+1})<\dots<\alpha(v_k) 且 α(vi)<α(vi−1)<⋯<α(v1)<α(vk)<α(v0)\alpha(v_i)<\alpha(v_{i-1})<\dots<\alpha(v_1)<\alpha(v_k)<\alpha(v_0).

证明:

由于 α(v1)<α(vk)<α(v0)\alpha(v_1)<\alpha(v_k)<\alpha(v_0),且 v1v0v_1v_0 相连,vkv0v_kv_0 不相连,所以由性质一,存在 xx 满足 α(vk)<α(x)\alpha(v_k)<\alpha(x) 且 vkxv_kx 相连,v1xv_1x 不相连.

考虑最小的 j∈(1,k]j\in(1,k] 满足 vjxv_jx 相连,我们可以推出 v0xv_0x 不相连,否则 v0v1⋯vjxv_0v_1\cdots v_jx 构成了一个长度 ≥4\ge 4 且无弦的环.

如果 x<v0x<v_0,则 v0,v1,…,vj,xv_0,v_1,\dots,v_j,x 也是一个满足性质的序列;如果 v0<xv_0<x 则 x,vj,…,v1,v0x,v_j,\dots,v_1,v_0 也是一个满足性质的序列.

在上面的推导中,我们扩大了 min⁡(v0,vk)\min(v_0,v_k),于是一直推下去,一定会产生矛盾.

Theorem 1:对于任何一个弦图,最大势算法求出的序列一定是一个完美消除序列.

证明:考虑任意三个点 u,v,wu,v,w 满足 α(u)<α(v)<α(w)\alpha(u)<\alpha(v)<\alpha(w),我们需要证明若 uvuv 相连,uwuw 相连,则 vwvw 一定相连.

考虑反证法,假设不相连,那么 w,u,vw,u,v 就是一个满足 Lemma 10 中性质的序列,我们证明了这样序列不存在,所以矛盾,vwvw 相连.

参考代码:

while (cur) {
  p[cur] = h[nww];
  rnk[p[cur]] = cur;
  h[nww] = nxt[h[nww]];
  lst[h[nww]] = 0;
  lst[p[cur]] = nxt[p[cur]] = 0;
  tf[p[cur]] = true;
  for (vector<int>::iterator it = G[p[cur]].begin(); it != G[p[cur]].end();
       it++)
    if (!tf[*it]) {
      if (h[deg[*it]] == *it) h[deg[*it]] = nxt[*it];
      nxt[lst[*it]] = nxt[*it];
      lst[nxt[*it]] = lst[*it];
      lst[*it] = nxt[*it] = 0;
      deg[*it]++;
      nxt[*it] = h[deg[*it]];
      lst[h[deg[*it]]] = *it;
      h[deg[*it]] = *it;
    }
  cur--;
  if (h[nww + 1]) nww++;
  while (nww && !h[nww]) nww--;
}

如果此时原图是弦图,此时求出的就是完美消除序列;但是由于原图可能不是弦图,此时求出的一定不是完美消除序列,所以问题转化为 判断求出的序列是否是原图的完美消除序列.

判断一个序列是否是完美消除序列

朴素算法

根据定义,依次判断完美消除序列 vv 上 {vi,vi+1,…,vn}\{v_i,v_{i+1},\ldots ,v_n\} 中与 viv_i 相邻的点是否构成了一个团.时间复杂度 O(nm)O(nm).

优化后的算法

根据完美消除序列的定义,设 viv_i 在 vi,vi+1,…,vn{v_i,v_{i+1},\ldots , v_n} 中相邻的点从小到大为 {vc1,vc2,…,vck}\{v_{c_1},v_{c_2},\ldots ,v_{c_k} \},则只需判断 vc1v_{c_1} 与其他点是否直接连通即可.时间复杂度 O(n+m)O(n+m).

jud = true;
for (int i = 1; i <= n; i++) {
  cur = 0;
  for (vector<int>::iterator it = G[p[i]].begin(); it != G[p[i]].end(); it++)
    if (rnk[p[i]] < rnk[*it]) {
      s[++cur] = *it;
      if (rnk[s[cur]] < rnk[s[1]]) swap(s[1], s[cur]);
    }
  for (int j = 2; j <= cur; j++)
    if (!st[s[1]].count(s[j])) {
      jud = false;
      break;
    }
}
if (!jud)
  printf("Imperfect\n");
else
  printf("Perfect\n");

至此,弦图判定问题 可以在 O(n+m)O(n+m) 的时间复杂度内解决.

弦图的极大团

令 N(x)N(x) 为满足与 xx 直接有边相连且在完美消除序列上的 xx 之后的序列.则弦图的极大团一定为 {x}+N(x)\{x\}+N(x).

证明:考虑弦图的一个极大团 VV,其中的点在完美消除序列中出现的第一个点 xx,一定有 V⊆{x}+N(x)V\subseteq \{x\}+N(x),又因为 VV 是极大团,所以 V={x}+N(x)V=\{x\}+N(x).

弦图最多有 nn 个极大团.求出弦图的每个极大团,可以判断每个 {x}+N(x)\{x\}+N(x) 是否为极大团.

设 A={x}+N(x),B={y}+N(y)A=\{x\}+N(x),B=\{y\}+N(y),若 A⫋BA\subsetneqq B,则 AA 不是极大团.此时在完美消除序列上显然有 yy 在 xx 前.

设 nxtxnxt_x 表示 N(x)N(x) 中在完美消除序列上最靠前的点,y∗y* 表示所有满足 A⊆BA\subseteq B 的 yy 中的最靠后的点.此时必然有 nxty∗=xnxt_{y*}=x,否则 y∗y* 不是最靠后的,令 y∗=nxty∗y*=nxt_{y*} 仍然满足条件.

A⫋BA\subsetneqq B 当且仅当 ∣A∣+1≤∣B∣|A|+1\le |B|.

问题转化为判断是否存在 yy,满足 nxty=xnxt_y=x 且 ∣N(x)∣+1≤∣N(y)∣|N(x)|+1\le |N(y)|.时间复杂度 O(n+m)O(n+m).

for (int i = 1; i <= n; i++) {
  cur = 0;
  for (vector<int>::iterator it = G[p[i]].begin(); it != G[p[i]].end(); it++)
    if (rnk[p[i]] < rnk[*it]) {
      s[++cur] = *it;
      if (rnk[s[cur]] < rnk[s[1]]) swap(s[1], s[cur]);
    }
  fst[p[i]] = s[1];
  N[p[i]] = cur;
}
for (int i = 1; i <= n; i++) {
  if (!vis[p[i]]) ans++;
  if (N[p[i]] >= N[fst[p[i]]] + 1) vis[fst[p[i]]] = true;
}

弦图的色数/弦图的团数

一种构造方法:按完美消除序列从后往前依次给每个点染色,给每个点染上可以染的最小颜色.时间复杂度 O(m+n)O(m+n).

正确性证明:设以上方法使用了 tt 种颜色,则 t≥χ(G)t\ge \chi(G).由于团上每个点都是不同的颜色,所以 t=ω(G)t=\omega(G),由 Lemma 1,t=ω(G)≤χ(G)t=\omega(G)\le \chi(G).综上,可得 t=χ(G)=ω(G)t=\chi(G)=\omega(G).

无需染色方案,只需求出弦图的色数/团数时,可以取 ∣{x}+N(x)∣|\{x\}+N(x)| 的最大值得到.

for (int i = 1; i <= n; i++) ans = max(ans, deg[i] + 1);

弦图的最大独立集/最小团覆盖

最大独立集:完美消除序列从前往后,选择所有没有与已经选择的点有直接连边的点.

最小团覆盖:设最大独立集为 {v1,v2,…,vt}\{v_1,v_2,\ldots ,v_t\},则团的集合 {{v1+N(v1)},{v2+N(v2)},…,{vt+N(vt)}}\{\{v_1+N(v_1)\},\{v_2+N(v_2)\},\ldots ,\{v_t+N(v_t)\} \} 为图的最小团覆盖.时间复杂度均为 O(n+m)O(n+m).

正确性证明:设以上方案独立集数和团覆盖数为 tt,由定义得 t≤α(G),t≥κ(G)t\le \alpha(G),t\ge \kappa(G),由 Lemma 2 得,α(G)≤κ(G)\alpha(G)\le \kappa(G),所以 t=α(G)=κ(G)t=\alpha(G)=\kappa(G).

for (int i = 1; i <= n; i++)
  if (!vis[p[i]]) {
    ans++;
    for (vector<int>::iterator it = G[p[i]].begin(); it != G[p[i]].end(); it++)
      vis[*it] = true;
  }

习题

SPOJ FISHNET - Fishing Net

P3196[HNOI2008] 神奇的国度

P3852[TJOI2007] 小朋友

参考资料

弦图相关

2009 WC 讲稿

弦图总结 - 租酥雨

R. E. Tarjan and M. Yannakakis, Simple linear-time algorithms to test chordality of graphs,test acyclicity of hypergraphs,and selectively reduce acyclic hypergraphs, SIAM J. Comput., 13 (1984), pp. 566–579.

本文内容来自 OI Wiki《弦图》;作者与修订历史见原文。AIOJ 对排版、站内链接和资源路径做了格式适配。

内容依 CC BY-SA 4.0 与 SATA 许可共享;附加条款以原文声明为准。感谢 OI Wiki 作者与贡献者;请按 SATA 约定前往 OI Wiki 仓库 Star 支持。