【BZOJ4945】【NOI2017】游戏(搜索,2-sat)
【NOI2017】游戏(搜索,2-sat)
题面
BZOJ的SPJ是假的
兹磁洛谷
题解
如果没有\(x\)地图的影响
这就是一个裸的\(2-sat\)问题
但是现在有不超过\(8\)个\(x\)地图的影响
我们不难想到枚举\(x\)地图的状态再来\(2-sat\)判断剩余是否可行。
这样的复杂度是\(O(3^dn)\),稍微算一下发现这个复杂度有点假
考虑如何优化,我们枚举\(x\)地图不连什么
表面上看起来还是\(O(3^dn)\)
但是,当他等价于\(a,b\)两种地图时,\(2^d\)种方案中一定会包含着最后的解,也就是说,\(a,b\)两种地图替换\(x\)的所有状态,等价于\(a,b,c\)三种地图替换\(x\)
这样的复杂度是\(O(2^dn)\)
现在考虑\(2-sat\)如何连边
对于限制\(i->j\)
显然直接连接\(i->j\)
并且,如果\(j'\)被选择,那么必须选择\(i'\)
所以,同时连接\(j'->i'\)
一些没有意义的连边可以直接省略。
至于输出方案,\(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;
}