poj3678 Katu Puzzle 2-sat

http://poj.org/problem?id=3678

变量名写错了检查好久。。。

2018年寒假的时候过于摸鱼，没有看懂2-sat，现在看了算法竞赛进阶指南上，发现讲得很好很容易懂，而且跟并查集的虚点进行了对比，当时并查集的虚点还是TB教我的，非常巧妙，印象比较深刻。

这题我们把几种位运算分一下情况。

a and b =0 ，此时我们知道2个都是0时可行的，2-sat的连边是表示有了条件1，必须有条件2，这样的意思

所以当a=1时，b必须为0，b=1时，a必须为0，连边就是add(a+n,b),add(b+n,a);

a and b =1，时，我们知道a，b都必须为1，那么反过来想，a=0，b=0就是不可能的情况，于是add(a,a+n),add(b,b+n)

在2-sat中我们要注意一个命题的逆否命题也需要加入图中，而a=0,a=1的逆否命题就是本身，所以不用加额外的条件。

接下来就是一样的讨论

a | b = 1 ,add(a,b+n),add(b,a+n)

a | b = 0 ,add(a+n,a),add(b+n,b);

a ^ b = 0,add(a,b),add(b,a),add(a+n,b+n) ,add(b+n,a+n)

a ^ b = 1,add(a,b+n),add(b+n,a),add(a+n,b),add(a,b+n)

可以发现异或的时候是可以用并查集的虚点的方法解决的，即只决定了a,b相同还是a,b不同的关系，而并查集的虚点实际上就是一条双向边，在并查集中a^b=0就是a，b相同，即union(a,b),union(a+n,b+n)代表其命题和逆否命题，在2-sat中就是4条边。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int maxl=2e3+10;
const int maxm=4e6+10;

int n,m,cnt,top,ind,ff,ans;
int ehead[maxl],s[maxl],dfn[maxl],low[maxl],f[maxl];
struct ed
{
	int to,nxt;
}e[maxm];
bool in[maxl];
char ch[5];

inline void add(int u,int v)
{
	e[++cnt].to=v;e[cnt].nxt=ehead[u];ehead[u]=cnt;
}

inline void prework()
{
	ind=0;top=0;ff=0;cnt=0;
	for(int i=1;i<=2*n;i++)
		ehead[i]=0,dfn[i]=0,low[i]=0,in[i]=false; 
	int a,b,c;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d%s",&a,&b,&c,ch);
		a++;b++;
		if(ch[0]=='A')
		{
			if(c==0)
				add(a+n,b),add(b+n,a);
			else
				add(a,a+n),add(b,b+n);
		}
		else if(ch[0]=='O')
		{
			if(c==1)
				add(a,b+n),add(b,a+n);
			else
				add(a+n,a),add(b+n,b);
		} 
		else
		{
			if(c==0)
				add(a+n,b+n),add(a,b),add(b+n,a+n),add(b,a);
			else
				add(a,b+n),add(a+n,b),add(b,a+n),add(b+n,a);
		}
	}
}

inline void tarjan(int u)
{
	dfn[u]=low[u]=++ind;s[++top]=u;in[u]=true;
	int v;
	for(int i=ehead[u];i;i=e[i].nxt)
	{
		v=e[i].to;
		if(!dfn[v])
		{
			tarjan(v);
			low[u]=min(low[v],low[u]);
		}else if(in[v])
			low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u])
	{
		ff++;
		do
		{
			v=s[top--];
			f[v]=ff;in[v]=false;
		}while(v!=u);
	}
}

inline void mainwork()
{
	ans=1;
	for(int i=1;i<=2*n;i++)
	if(!dfn[i])
		tarjan(i);
	for(int i=1;i<=n;i++)
	if(f[i]==f[i+n])
		ans=0;
}

inline void print()
{
	if(ans)
		puts("YES");
	else
		puts("NO");
}

int main()
{
	//freopen("4.in","r",stdin);
	while(~scanf("%d%d",&n,&m))
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}