bzoj1006 [HNOI2008]神奇的国度

最新推荐文章于 2024-11-05 17:17:10 发布

anmob303090

最新推荐文章于 2024-11-05 17:17:10 发布

阅读量89

点赞数

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/StupidBoy/p/4586812.html

版权

1006: [HNOI2008]神奇的国度

Time Limit: 20 Sec Memory Limit: 162 MB
Submit: 2304 Solved: 1043

Description

K国是一个热衷三角形的国度,连人的交往也只喜欢三角原则.他们认为三角关系:即AB相互认识,BC相互认识,CA相互认识,是简洁高效的.为了巩固三角关系,K国禁止四边关系,五边关系等等的存在.所谓N边关系,是指N个人 A1A2...An之间仅存在N对认识关系:(A1A2)(A2A3)...(AnA1),而没有其它认识关系.比如四边关系指ABCD四个人 AB,BC,CD,DA相互认识,而AC,BD不认识.全民比赛时,为了防止做弊，规定任意一对相互认识的人不得在一队，国王相知道，最少可以分多少支队。

Input

第一行两个整数N，M。1<=N<=10000,1<=M<=1000000.表示有N个人，M对认识关系. 接下来M行每行输入一对朋友

Output

输出一个整数，最少可以分多少队

Sample Input

4 5
1 2
1 4
2 4
2 3
3 4

Sample Output

HINT

一种方案(1,3)(2)(4)

Source

题目大意：给出一个图，此图中不存在点数>3的环（这里的环指内部没边相连的环），问给此图的点染色，使得相邻的点颜色不同，所需要的最少颜色种类

分析：

首先此题有论文：

http://wenku.baidu.com/link?url=Vz-Cjrq24tvblbk0l1I4IPYYfZ82jkGpzGBoNBDhi6qSCJko9hSTcrKGWrunO91TlU_eei4KBIdw1Y9U9bTnZOUJCPau8qTmUzcIzik_nQ3

简述一下：（省略一大堆证明之类的东东）

首先，像题目中那样的图（不存在点数>3的环）就做弦图，

一些点v1, v2, v3, ....., vn以及这些点之间相连的边叫这些点的诱导子图，

一些点的诱导子图如果是完全图，那就叫做团，

一个图中与一个点 i 相连的点为N( i )，

一个图中，如果有一个点 i 与N( i )的诱导子图是完全图的话，那么这个点 i 叫做单纯点

弦图中有样东西叫完美消除数列，这个数列是这样的：数列中第 i 个点在数列第 i - n 个点的诱导子图中是单纯点

怎么求完美消除序列：

我们定义Label[i]表示第 i 个点与多少已标记的点相连，一开始所有的点都没有标记，

然后选择一个Label最大的点 i ，标记这个点，并更新其他点的label，

按照标记顺序得到的数列便叫完美消除序列

这个可以用堆（优先队列）做，复杂度 nlogn，

也可以优化至O(n+m)，建n-1个链表List，分别表示Label为x的有哪些（存入List[x]），

High记录最大的Label，

选择时：每次选择Label最大的点，从list[high]中从前往后选（注意忽略已经在完美消除序列中的点），选择到后，在list中删除此点，并加入完美消除序列

更新时：每次更新Label，比如从x改为x+1，那么不需要删除List[x]中的那个点，在List[x+1]的后面加入即可

复杂度O(n+m)

当然，本程序使用了优先队列

根据这个完美消除序列染色，每次染色都尽量染编号小(就是说尽量染以前的颜色，并保证不与相邻的已经染色的点的颜色相同)的颜色，得到的颜色数就是答案，即最小染色数

综上所述，本题得解：

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <cstdlib>
 4 #include <cmath>
 5 #include <deque>
 6 #include <vector>
 7 #include <queue>
 8 #include <iostream>
 9 #include <algorithm>
10 #include <map>
11 #include <set>
12 #include <ctime>
13 using namespace std;
14 typedef long long LL;
15 #define For(i, s, t) for(int i = (s); i <= (t); i++)
16 #define Ford(i, s, t) for(int i = (s); i >= (t); i--)
17 #define MIT (2147483647)
18 #define INF (1000000001)
19 #define MLL (1000000000000000001LL)
20 #define sz(x) ((bnt) (x).size())
21 #define clr(x, y) memset(x, y, sizeof(x))
22 #define puf push_front
23 #define pub push_back
24 #define pof pop_front
25 #define pob pop_back
26 #define ft first
27 #define sd second
28 #define mk make_pair
29 inline void SetIO(string Name) {
30     string Input = Name+".in",
31     Output = Name+".out";
32     freopen(Input.c_str(), "r", stdin),
33     freopen(Output.c_str(), "w", stdout);
34 }
35 
36 const int N = 10010, M = 1000010;
37 typedef pair<int, int> II;
38 int First[N], To[M*2], Next[M*2], Tot;
39 int n, m;
40 int Label[N], High, Que[N];
41 priority_queue<II> S;
42 bool Visit[N];
43 int Color[N], Draw[N];
44 
45 inline void Insert(int u, int v) {
46     Tot++;
47     To[Tot] = v, Next[Tot] = First[u];
48     First[u] = Tot;
49 }
50 
51 inline void Input() {
52     scanf("%d%d", &n, &m);
53     For(i, 1, m) {
54         int u, v;
55         scanf("%d%d", &u, &v);
56         Insert(u, v), Insert(v, u);
57     }
58 }
59 
60 inline void Solve() {
61     High = 0, clr(Label, 0), clr(Visit, 0);
62     For(i, 1, n) S.push(II(0, i));
63     For(t, 1, n) {
64         II It = S.top();
65         while(Visit[It.sd]) {
66             S.pop();
67             It = S.top();
68         }
69         Visit[It.sd] = 1, Que[t] = It.sd;
70         for(int v, Tab = First[It.sd]; Tab; Tab = Next[Tab])
71             if(!Visit[v = To[Tab]]) {
72                 Label[v]++;
73                 S.push(II(Label[v], v));
74             }
75     }
76     
77     int Ans = 0;
78     For(i, 1, n) Color[i] = n+1;
79     clr(Draw, 0);
80     For(t, 1, n) {
81         int x = Que[t];
82         for(int Tab = First[x]; Tab; Tab = Next[Tab])
83             Draw[Color[To[Tab]]] = t;
84         Color[x] = 1;
85         while(Draw[Color[x]] == t) Color[x]++;
86         Ans = max(Ans, Color[x]);
87     }
88     
89     printf("%d\n", Ans);
90 }
91 
92 int main() {
93     SetIO("1006");
94     Input();
95     Solve();
96     return 0;
97 }