【bzoj 1006】 神奇的国度 【HNOI2008】

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

Sample Output

3

HINT

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

这道题是最大势算法，可以用链表优化到O(n+m)或者用优先队列优化到O(mlogm)，下面是程序：

链表：

#include<stdio.h>
#include<queue>
#include<string.h>
#include<iostream>
using namespace std;
const int N=10005,M=1000005,INF=0x7f7f7f7f;
struct edge{
	int v,next;
}a[2*M];
queue<int>q[N];
int head[N],k,ans[N];
bool vis[N];
void read(int &s){
	s=0;
	char c=getchar();
	while(c<'0'||c>'9'){
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		s*=10;
		s+=c-'0';
		c=getchar();
	}
}
void add(int u,int v){
	k++;
	a[k].v=v;
	a[k].next=head[u];
	head[u]=k;
}
int main(){
	int n,m,u,v,i,s=0,last=0;
	read(n);
	read(m);
	while(m--){
		read(u);
		read(v);
		add(u,v);
		add(v,u);
	}
	for(i=1;i<=n;i++){
		q[0].push(i);
	}
	while(last>=0){
		while(!q[last].empty()){
			u=q[last].front();
			q[last].pop();
			if(!vis[u]){
				bool ok=0;
				vis[u]=1;
				for(i=head[u];i;i=a[i].next){
					v=a[i].v;
					if(!vis[v]){
						ans[v]++;
						q[ans[v]].push(v);
						ok=1;
					}
				}
				last+=ok;
			}
		}
		while(q[last].empty()){
			last--;
		}
	}
	for(i=1;i<=n;i++){
		s=max(s,ans[i]);
	}
	printf("%d\n",s+1);
	return 0;
}

优先队列：

#include<stdio.h>
#include<queue>
#include<iostream>
using namespace std;
const int N=10005,M=1000005;
struct edge{
	int v,next;
}a[2*M];
struct node{
	int x,s;
	bool operator <(node p)const{
		return s==p.s?x>p.x:s<p.s;
	}
};
priority_queue<node>q;
int k,head[N],ans[N];
bool vis[N];
void read(int &s){
	s=0;
	char c=getchar();
	while(c<'0'||c>'9'){
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		s*=10;
		s+=c-'0';
		c=getchar();
	}
}
void add(int u,int v){
	k++;
	a[k].v=v;
	a[k].next=head[u];
	head[u]=k;
}
int main(){
	int n,m,u,v,i;
	read(n);
	read(m);
	while(m--){
		read(u);
		read(v);
		add(u,v);
		add(v,u);
	}
	for(i=1;i<=n;i++){
		node tp;
		tp.s=0;
		tp.x=i;
		q.push(tp);
	}
	while(!q.empty()){
		node tp=q.top();
		q.pop();
		u=tp.x;
		if(!vis[u]){
			vis[u]=1;
			for(i=head[u];i;i=a[i].next){
				v=a[i].v;
				if(!vis[v]){
					ans[v]++;
					tp.s=ans[v];
					tp.x=v;
					q.push(tp);
				}
			}
		}
	}
	for(i=1,u=0;i<=n;i++){
		u=max(u,ans[i]);
	}
	printf("%d\n",u+1);
	return 0;
}

暴力(O(n^2)竟然也能过。。。)：

#include<stdio.h>
#include<iostream>
using namespace std;
const int N=10005,M=1000005;
struct edge{
	int v,next;
}a[2*M];
int head[N],k,ans[N];
bool vis[N];
void read(int &s){
	s=0;
	char c=getchar();
	while(c<'0'||c>'9'){
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		s*=10;
		s+=c-'0';
		c=getchar();
	}
}
void add(int u,int v){
	k++;
	a[k].v=v;
	a[k].next=head[u];
	head[u]=k;
}
int main(){
	int n,m,u,v,i,j;
	read(n);
	read(m);
	while(m--){
		read(u);
		read(v);
		add(u,v);
		add(v,u);
	}
	for(i=1;i<=n;i++){
		for(j=1,v=-1;j<=n;j++){
			if(!vis[j]&&ans[j]>v){
				v=ans[j];
				u=j;
			}
		}
		vis[u]=1;
		for(j=head[u];j;j=a[j].next){
			v=a[j].v;
			if(!vis[v]){
				ans[v]++;
			}
		}
	}
	for(i=1,u=0;i<=n;i++){
		u=max(u,ans[i]);
	}
	printf("%d\n",u+1);
	return 0;
}