1006: [HNOI2008]神奇的国度
Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 3494 Solved: 1592
[ Submit][ Status][ Discuss]
Description
K国是一个热衷三角形的国度,连人的交往也只喜欢三角原则.他们认为三角关系:即AB相互认识,BC相互认识,CA
相互认识,是简洁高效的.为了巩固三角关系,K国禁止四边关系,五边关系等等的存在.所谓N边关系,是指N个人 A1A2
...An之间仅存在N对认识关系:(A1A2)(A2A3)...(AnA1),而没有其它认识关系.比如四边关系指ABCD四个人 AB,BC,C
D,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
1 2
1 4
2 4
2 3
3 4
Sample Output
3
HINT
一种方案(1,3)(2)(4)
Source
基于图中的一个环,我们定义以下概念:
弦:连接一个简单环中不相邻的两个点的边。
弦图:任意长度大于3的简单环都至少包含一条弦的图。
单纯点:若点x和它的邻域组成的点集的生成子图为一个团,则称x是单纯点。
显然,基于题设的限制,给出的图一定是一张弦图
原题等价于,给每个点染色,要求相邻的点不同色
对这张弦图,求出一个完美消除序列,即一个点集序列,满足:
对于序列中任意一个后缀点集的生成子图该后缀的第一个点是单纯点
有个结论,在完美消除序列里,从后往前染色,每个点取当前能染的编号最小的颜色,这样就是最优了
然后完美消除序列可用最大势算法求解,从n~1给每个点标号
记label[i]:与i相邻的已经标号的点的数量
每次取label值最大的点给出当前标号
最后,完美消除序列的第k个点就是标号为k的点
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<bitset>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;
const int maxn = 1E4 + 10;
struct data{
int Num,w; data(){}
data(int Num,int w): Num(Num),w(w){}
bool operator < (const data &B) const {return w < B.w;}
};
typedef __gnu_pbds::priority_queue<data,less<data>,__gnu_pbds::pairing_heap_tag> Heap;
int n,m,cnt,label[maxn],Num[maxn],c[maxn],bo[maxn];
bool vis[maxn];
Heap Q;
vector <int> v[maxn];
Heap::point_iterator id[maxn];
int getint()
{
char ch = getchar(); int ret = 0;
while (ch < '0' || '9' < ch) ch = getchar();
while ('0' <= ch && ch <= '9')
ret = ret * 10 + ch - '0',ch = getchar();
return ret;
}
int main()
{
#ifdef DMC
freopen("DMC.txt","r",stdin);
#endif
n = getint(); m = getint();
while (m--)
{
int x = getint(),y = getint();
v[x].push_back(y);
v[y].push_back(x);
}
for (int i = 1; i <= n; i++) id[i] = Q.push(data(i,0));
for (int i = n; i; i--)
{
data K = Q.top(); Q.pop();
int k = K.Num; Num[i] = k; vis[k] = 1;
for (int j = 0; j < v[k].size(); j++)
{
int to = v[k][j];
if (vis[to]) continue; ++label[to];
Q.modify(id[to],data(to,label[to]));
}
}
for (int i = n; i; i--)
{
int k = Num[i]; ++cnt;
for (int j = 0; j < v[k].size(); j++)
bo[c[v[k][j]]] = cnt;
for (int j = 1; ; j++)
if (bo[j] != cnt) {c[k] = j; break;}
}
int Ans = 0;
for (int i = 1; i <= n; i++) Ans = max(Ans,c[i]);
cout << Ans << endl;
return 0;
}