[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;  }