[bzoj] 1242: Zju1015 Fishing Net弦图判定

依然是cdq的论文，用最大势求出完美消除序列，然后判断是否合法。

用学长的话来概括：

一句话概括 MCS 算法：每次寻找没有标号的点中连接已标号的点最多的，标号该点后加入序列。

一句话概括判断序列是否合法：,判断序列中每个元素相邻的所有点 j1,j2,j3,...,jn 中序号最小的点 j1 和其他点 j2,j3,...,jn 是否相邻。

复杂度为： O(n+m)

#include <cstdio>
#define C (c=getchar())
using namespace std;
const int N=1005;
int n,m,tot,qaq;
int head[N],sa[N];
int sh[N],num[N];
bool vis[N];
int t[N];
int map[N][N],hi;
inline void read(int &x){  
    static char c;
	int b=1;x=0;C;
    while (c<'0'||c>'9') {if (c=='-') b=-1;C;}
	while (c>='0'&&c<='9') x=x*10+c-'0',C;
	x*=b;
}  
struct edge{
	int to,nxt;
}e[5000005];
inline void add(int u,int v)
{
	e[++tot].to=v;
	e[tot].nxt=head[u];
	head[u]=tot;
	return;
}

int main()
{
	register int i,j;
	read(n);read(m);
	for (i=1;i<=m;i++)
	{
		int u,v;
		read(u);read(v);
		add(u,v);add(v,u);
		map[u][v]=map[v][u]=1;
	}
	for (i=1;i<=n;i++) add(n+1,i);
	for (i=n;i;i--)
	{
		int v;
		while (1)
		{
			v=0;
			for (j=head[hi+1+n];j&&!v;j=e[j].nxt)
			{
				if (!vis[e[j].to]) v=e[j].to;
				else head[hi+1+n]=e[j].nxt;
			}
			if (v)
			{
				vis[v]=1;sa[i]=v;
				for (j=head[v];j;j=e[j].nxt)
				if (!vis[e[j].to])
				{
					sh[e[j].to]++;
					add(sh[e[j].to]+1+n,e[j].to);
					if(sh[e[j].to]>hi) hi=sh[e[j].to];
					}
				break;
			}
			else hi--;
		}
	}
	for (i=1;i<=n;i++) num[sa[i]]=i;
	for (i=n;i;i--)
	{
		int x=sa[i],mi=N+1; qaq=0;
		for (j=head[x];j;j=e[j].nxt)
		{
			if (num[e[j].to]>num[x])
			{
				t[++qaq]=num[e[j].to];
				if (num[e[j].to]<mi) mi=num[e[j].to];
			}
		}
		for (j=1;j<=qaq;j++)
		{
			if (mi!=t[j]&&!map[sa[mi]][sa[t[j]]])
			{
				printf("Imperfect\n");
				return 0;
			}
		}
	}
	printf("Perfect\n");  
    return 0;  
}