#【NOI2017】游戏(搜索,2-sat) ##题面 BZOJ的SPJ是假的 兹磁洛谷ios
##题解 若是没有$x$地图的影响 这就是一个裸的$2-sat$问题优化
可是如今有不超过$8$个$x$地图的影响 咱们不难想到枚举$x$地图的状态再来$2-sat$判断剩余是否可行。 这样的复杂度是$O(3^dn)$,稍微算一下发现这个复杂度有点假ui
考虑如何优化,咱们枚举$x$地图不连什么 表面上看起来仍是$O(3^dn)$ 可是,当他等价于$a,b$两种地图时,$2^d$种方案中必定会包含着最后的解,也就是说,$a,b$两种地图替换$x$的全部状态,等价于$a,b,c$三种地图替换$x$ 这样的复杂度是$O(2^dn)$spa
如今考虑$2-sat$如何连边 对于限制$i->j$ 显然直接链接$i->j$ 而且,若是$j'$被选择,那么必须选择$i'$ 因此,同时链接$j'->i'$ 一些没有意义的连边能够直接省略。code
至于输出方案,$Tarjan$求出来的强连通份量已经和拓扑序相关,所以没有必要从新建图拓扑排序排序
#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> #include<set> #include<map> #include<vector> #include<queue> using namespace std; #define ll long long #define RG register #define MAX 333333 inline int read() { RG int x=0,t=1;RG char ch=getchar(); while((ch<'0'||ch>'9')&&ch!='-')ch=getchar(); if(ch=='-')t=-1,ch=getchar(); while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar(); return x*t; } char g[MAX]; struct Limit{int i,j;char a,b;}q[MAX]; int n,D,m,pos[20]; struct Line{int v,next;}e[MAX<<3]; int h[MAX],cnt; int p[MAX]; inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;} bool fl; void Build() { memset(h,0,sizeof(h));cnt=1;fl=true; for(int i=1;i<=n;++i) { p[i]=p[i+n]=p[i+n+n]=0; if(g[i]=='a')p[i+n]=i+n+n,p[i+n+n]=i+n; else if(g[i]=='b')p[i]=i+n+n,p[i+n+n]=i; else p[i]=i+n,p[i+n]=i; } for(int i=1;i<=m;++i) { int u=q[i].i,v=q[i].j; int a=q[i].a-65,b=q[i].b-65; if(u==v&&a==b)continue; if(g[u]-97==a)continue; if(g[v]-97==b){Add(u+a*n,p[u+a*n]);continue;} Add(u+a*n,v+b*n); Add(p[v+b*n],p[u+a*n]); } } int dfn[MAX],low[MAX],S[MAX],top; int gr,G[MAX],tim,ans[MAX]; bool Ins[MAX]; void Tarjan(int u) { dfn[u]=low[u]=++tim; S[++top]=u;Ins[u]=true; int v; for(int i=h[u];i;i=e[i].next) { v=e[i].v; if(!dfn[v])Tarjan(v),low[u]=min(low[u],low[v]); else if(Ins[v])low[u]=min(low[u],dfn[v]); } if(dfn[u]==low[u]) { ++gr; do{v=S[top--];G[v]=gr;Ins[v]=false;}while(u!=v); } } void Calc() { Build();gr=tim=0;if(!fl)return; memset(dfn,0,sizeof(dfn)); memset(low,0,sizeof(low)); memset(G,0,sizeof(G)); for(int i=1;i<=n+n+n;++i)if(!dfn[i])Tarjan(i); for(int i=1;i<=n+n+n;++i)if(G[i]==G[p[i]])return; for(int i=1;i<=n;++i) { if(g[i]=='a')(G[i+n]<G[i+n+n])?ans[i]=1:ans[i]=2; if(g[i]=='b')(G[i]<G[i+n+n])?ans[i]=0:ans[i]=2; if(g[i]=='c')(G[i]<G[i+n])?ans[i]=0:ans[i]=1; } for(int i=1;i<=n;++i)putchar(ans[i]+65);puts(""); exit(0); } void dfs(int x) { if(x==D+1){Calc();return;} g[pos[x]]='a';dfs(x+1); g[pos[x]]='b';dfs(x+1); } int main() { n=read();D=read();D=0; scanf("%s",g+1); m=read(); for(int i=1;i<=m;++i) q[i].i=read(),q[i].a=getchar(),q[i].j=read(),q[i].b=getchar(); for(int i=1;i<=n;++i) if(g[i]=='x')pos[++D]=i; dfs(1); puts("-1"); return 0; }