[BZOJ1242][Zju1015][弦图判定]Fishing Net弦图判定

题意

弦图判定

弦图判定裸题……然而小号WA了一页

用MCS算出完美消除序列

设{vi+1,…,vn}中所有与vi相邻的点依次为vj1, …, vjk。

只需判断vj1是否与vj2, …, vjk相邻即可。

——CDQ

这里用一点要注意，vj1是在序列中标号最小的与vi相邻的点，不是任意点……我太菜没有注意到这点

#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 1010

using namespace std;

int n,m,u,v,cnt;
int G
,p
,vis
,c
,tmp
,rank
;
int Map

;
struct edge{
int t,nx;
}E[5000005];
struct stp{
int g;
stp *l,*r;
}A
,*B
;

inline void reaD(int &x){
char c=getchar(); x=0;
for(;c>57||c<48;c=getchar());for(;c>=48&&c<=57;x=x*10+c-48,c=getchar());
}

inline void Insert(int x,int y){
Map[x][y]=Map[y][x]=1;
E[++cnt].t=y; E[cnt].nx=G[x]; G[x]=cnt;
E[++cnt].t=x; E[cnt].nx=G[y]; G[y]=cnt;
}

inline void MCS(){
for(int i=1;i<=n;i++){
if(i>1) A[i].l=&A[i-1];
if(i<n) A[i].r=&A[i+1];
A[i].g=i;
}
B[0]=&A[1]; int maxx=0;
for(int i=n;i;i--){
while(!B[maxx]) maxx--;
int x=B[maxx]->g;
p[i]=x; vis[x]=1;
B[maxx]=B[maxx]->r;
if(B[maxx]) B[maxx]->l=0;
for(int j=G[x];j;j=E[j].nx)
if(!vis[E[j].t]){
int v=E[j].t;
if(A[v].l) A[v].l->r=A[v].r;
else B[c[v]]=B[c[v]]->r;
if(A[v].r) A[v].r->l=A[v].l;
c[v]++; A[v].l=0;
if(B[c[v]]) B[c[v]]->l=&A[v];
A[v].r=B[c[v]];
B[c[v]]=&A[v];
if(c[v]>maxx) maxx=c[v];
}
}
}

int main(){
freopen("1242.in","r",stdin);
freopen("1242.out","w",stdout);
reaD(n); reaD(m);
for(int i=1;i<=m;i++)
reaD(u),reaD(v),Insert(u,v);
MCS();
for(int i=1;i<=n;i++) rank[p[i]]=i;
for(int i=1;i<=n;i++){
int g=0,Min=n+1,v=p[i],x;
for(int j=G[v];j;j=E[j].nx)
if(rank[E[j].t]>rank[v]){
tmp[++g]=E[j].t;
if(rank[tmp[g]]<Min) Min=rank[tmp[g]],x=tmp[g];
}
for(int j=1;j<=g;j++)
if(tmp[j]!=x&&!Map[x][tmp[j]]) return puts("Imperfect"),0;
}
return puts("Perfect"),0;
}