原文链接http://www.cnblogs.com/zhouzhendong/p/8146041.html
题目传送门 - 洛谷3825
题解
我们考虑到地图中x的个数很少,最多只有8个。
所以我们可以考虑穷举。
我们只需要把x变成a和b,这样就涵盖了选择A,B,C的三种情况。
所以我们状压枚举每一个x可以变成什么情况。
然后对于每一种情况,几乎就是2-sat裸题了。
然后我们考虑特殊情况:
如果选了A就得选择B。
如果A的状态不合法,那么显然这条边是不用建立的。
否则:
如果B的状态不合法,那么显然不可以到达A,于是我们将状态A连向他的对立点。
如果B的状态合法,那么我们要连接2条边。一条是A到B的,一条是opp(B)到opp(A)的。(其中opp(x)表示状态x的对立状态)
同学们注意了,Tarjan缩点千万别写错,我已经连续3个程序死在Tarjan缩点上面了。
代码
#include <cstring> #include <cstdio> #include <algorithm> #include <cmath> #include <cstdlib> using namespace std; const int N=500005*2,M=1000005*2; struct Gragh{ int cnt,x[M],y[M],nxt[M],fst[N]; void clear(){ cnt=0; memset(fst,0,sizeof fst); } void add(int a,int b){ x[++cnt]=a,y[cnt]=b,nxt[cnt]=fst[a],fst[a]=cnt; } }g,g2; int n,m,d,pos[10],w[N][3]; int dfn[N],low[N],vis[N],inst[N],st[N],bh[N],top,time,cnt; int o[N],in[N],q[N],res[N],head,tail; char place[N]; struct Limit{ int i,j,hi,hj; }l[M]; void Get_Pos(){ int cnt=0; for (int i=1;i<=n;i++) if (place[i]=='x') pos[++cnt]=i; } int opp(int x){ return x+n*(x<=n?1:-1); } void Tarjan_Prepare(){ top=time=cnt=0; memset(st,0,sizeof st); memset(bh,0,sizeof bh); memset(dfn,0,sizeof dfn); memset(low,0,sizeof low); memset(vis,0,sizeof vis); memset(inst,0,sizeof inst); } void Tarjan(int x){ dfn[x]=low[x]=++time; st[++top]=x; inst[x]=vis[x]=1; for (int i=g.fst[x];i;i=g.nxt[i]) if (!vis[g.y[i]]){ Tarjan(g.y[i]); low[x]=min(low[x],low[g.y[i]]); } else if (inst[g.y[i]]) low[x]=min(low[x],low[g.y[i]]); if (dfn[x]==low[x]){ cnt++; bh[st[top]]=cnt; inst[st[top]]=0; while (st[top--]!=x){ bh[st[top]]=cnt; inst[st[top]]=0; } } } bool check(){ for (int i=1;i<=n;i++) if (bh[i]==bh[i+n]) return 0; else o[bh[i]]=bh[i+n],o[bh[i+n]]=bh[i]; return 1; } void solve(int s){ for (int i=1;i<=d;i++) place[pos[i]]='a'+((s>>(i-1))&1); for (int i=1;i<=n;i++){ if (place[i]=='a')w[i][0]=0,w[i][1]=i,w[i][2]=i+n; if (place[i]=='b')w[i][0]=i,w[i][1]=0,w[i][2]=i+n; if (place[i]=='c')w[i][0]=i,w[i][1]=i+n,w[i][2]=0; } g.clear(); for (int I=1;I<=m;I++){ int i=l[I].i,j=l[I].j,hi=l[I].hi,hj=l[I].hj; if (!w[i][hi]) continue; if (w[j][hj]){ g.add(w[i][hi],w[j][hj]); g.add(opp(w[j][hj]),opp(w[i][hi])); } else g.add(w[i][hi],opp(w[i][hi])); } Tarjan_Prepare(); for (int i=1;i<=n*2;i++) if (!vis[i]) Tarjan(i); if (!check()) return; memset(in,0,sizeof in); g2.clear(); for (int i=1;i<=g.cnt;i++) if (bh[g.x[i]]!=bh[g.y[i]]) g2.add(bh[g.y[i]],bh[g.x[i]]),in[bh[g.x[i]]]++; head=tail=0; for (int i=1;i<=cnt;i++) if (!in[i]) q[++tail]=i; memset(res,0,sizeof res); while (head<tail){ int x=q[++head]; if (!res[x]) res[x]=1,res[o[x]]=-1; for (int i=g2.fst[x];i;i=g2.nxt[i]){ in[g2.y[i]]--; if (!in[g2.y[i]]) q[++tail]=g2.y[i]; } } for (int i=1;i<=n;i++){ int v=res[bh[i]]==1?i:(i+n),c; for (int j=0;j<3;j++) if (w[i][j]==v) c=j; printf("%c",c+'A'); } exit(0); } int main(){ scanf("%d%d%s%d",&n,&d,place+1,&m); for (int i=1;i<=m;i++){ char ch1[3],ch2[3]; scanf("%d%s%d%s",&l[i].i,ch1,&l[i].j,ch2); l[i].hi=ch1[0]-'A',l[i].hj=ch2[0]-'A'; } Get_Pos(); for (int i=0;i<(1<<d);i++) solve(i); printf("-1"); return 0; }