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 }
BZOJ 1006

 

题目链接 & AC 通道

 

BZOJ 1006 [HNOI2008] 神奇的国度

 

 

 

posted @ 2013-08-24 19:23  jianzhang.zj  阅读(3567)  评论(0编辑  收藏  举报