BZOJ 1006 [HNOI2008] 神奇的国度(简单弦图的染色)
题目大意
K 国是一个热衷三角形的国度,连人的交往也只喜欢三角原则。他们认为三角关系:即 AB 相互认识,BC 相互认识,CA 相互认识,是简洁高效的。为了巩固三角关系,K 国禁止四边关系,五边关系等等的存在。所谓 N 边关系,是指 N 个人 A1 A2 ... An 之间仅存在 N 对认识关系:(A1 A2) (A2 A3) ... (An A1),而没有其它认识关系,比如四边关系指 A B C D 四个人 AB,BC,CD,DA 相互认识,而 AC,BD 不认识。全民比赛时,为了防止做弊,规定任意一对相互认识的人不得在一队,国王相知道,最少可以分多少支队。
Input
第一行两个整数 N,M。1<=N<=10000,1<=M<=1000000。表示有 N 个人,M 对认识关系.。接下来 M 行每行输入一对朋友
Output
输出一个整数,最少可以分多少队
做法分析
根据提题意,不难看出,所有的人构成的关系图是一个弦图(长度超过 3 的环中必有一条弦),求出它的完美性消除序列,根据完美消除序列逆序贪心的染色,最终所用的色数就是本题的答案
完美消除序列的求法:
使用陈丹琦讲述的 MCS 法,可以在 o(n+m) 的时间复杂度中求出一个图的完美消除序列,并在 o(n+m) 的时间复杂度下判断这个完美消除序列是否合法,不过,我不知道怎么在 o(n+m) 的时间复杂度下求出这个完美消除序列,或许是利用桶优化吧,我写了个堆优化的,时间复杂度也能接受
求完美消除序列的 MCS 法是倒着解的,也就是先求序列的第 n 个再求序列的第 n-1 个
每次选取图中具有最大标号的点作为完美消除序列中对应位置的点,并用这个点更新所有和他邻接的不再序列中的点的标号值
对于一个弦图的染色,用完美消除序列可以很好的解决,按照完美消除序列中的点倒着给图中的点贪心的然尽可能小的颜色
最终,一定能够用最小的颜色数量给图中的所有点染色
本题可以先求出来这个弦图的完美消除序列,由于一定是一个弦图,所以序列一定合法,直接根据消除序列染色就行
更多的和弦图区间图相关的知识,请看陈丹琦的PPT:弦图与区间图
参考代码
1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 #include <vector> 5 #include <queue> 6 7 using namespace std; 8 9 const int N=10006; 10 11 vector <int> arc[N]; 12 int n, m, R[N], SA[N], label[N]; 13 14 priority_queue <pair<int, int> > heap; 15 void Construct() { 16 fill(R+1, R+1+n, -1); 17 fill(label+1, label+1+n, 0); 18 for(int i=1; i<=n; i++) heap.push(make_pair(0, i)); 19 for(int cnt=n; cnt>=1; ) { 20 int id=heap.top().second; 21 heap.pop(); 22 if(R[id]!=-1) continue; 23 SA[cnt]=id, R[id]=cnt--; 24 for(int i=0, len=(int)arc[id].size(); i<len; i++) { 25 int u=arc[id][i]; 26 if(R[u]!=-1) continue; 27 label[u]++; 28 heap.push(make_pair(label[u], u)); 29 } 30 } 31 } 32 33 void Color(int u) { 34 for(int i=0, len=(int)arc[u].size(); i<len; i++) { 35 int v=arc[u][i]; 36 if(label[v]==-1) continue; 37 R[label[v]]=u; 38 } 39 for(int i=1; label[u]==-1; i++) if(R[i]!=u) label[u]=i; 40 } 41 42 int Color_Graph() { 43 fill(label+1, label+1+n, -1); 44 fill(R+1, R+1+n, -1); 45 for(int i=n; i>0; i--) Color(SA[i]); 46 int ans=0; 47 for(int i=1; i<=n; i++) ans=max(ans, label[i]); 48 return ans; 49 } 50 51 int main() { 52 scanf("%d%d", &n, &m); 53 for(int i=1; i<=n; i++) arc[i].clear(); 54 for(int i=0, a, b; i<m; i++) { 55 scanf("%d%d", &a, &b); 56 arc[a].push_back(b); 57 arc[b].push_back(a); 58 } 59 Construct(); 60 printf("%d\n", Color_Graph()); 61 return 0; 62 }
题目链接 & AC 通道