hdu 4421——Bit Magic

http://www.cnblogs.com/wuminye/p/3274609.html

第一道2-sat，模仿着写。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define maxn 1100
#define maxm (maxn*maxn)
int n;

int head[maxn],v[maxm],next[maxm],cnt;

int low[maxn],dfn[maxn],step;
int ID[maxn],IdNum;
int sta[maxn],vis[maxn],top;

int b[510][510];
void add(int x,int y)
{
	v[cnt]=y;
	next[cnt]=head[x];
	head[x]=cnt++;
}
void init()
{
	memset(head,-1,sizeof(head));
	memset(low,0,sizeof(low));
	memset(dfn,0,sizeof(dfn));
	memset(vis,0,sizeof(vis));
	step=top=cnt=IdNum=0;
}
void addedge(int x)
{
	for(int i=0;i<n;i++)// i代表0，i+n代表1 
		for(int j=i+1;j<n;j++)
		{
			int m=b[i][j]&(1<<x);
			if((i&1)&&(j&1))// |
			{
				if(m)
				{
					add(i,j+n);
					add(j,i+n);
				}
				else
				{
					add(i+n,i);//   i必须是0 
					add(j+n,j);//   j必须是0 
				}
			} else
			{
				if((i%2==0)&&(j%2==0))// &
				{
					if(m)
					{
						add(i,i+n);
						add(j,j+n);
					}
					else
					{
						add(i+n,j);
						add(j+n,i);
					}
				}
				else// ^
				{
					if(m)
	                {
	                    add(i,j+n);
	                    add(j,i+n);
	                    add(j+n,i);
	                    add(i+n,j);
	                }else  //==
	                {
	                    add(i,j);
	                    add(j,i);
	                    add(i+n,j+n);
	                    add(j+n,i+n);
	                }
				}
			}
		}
}
void tarjan(int u)
{
	dfn[u]=low[u]=++step;
	sta[++top]=u;
	vis[u]=1;
	for(int i=head[u];~i;i=next[i])
	{
		int to=v[i];
		if(!dfn[to])//没进过栈 
		{
			tarjan(to);
			low[u]=min(low[u],low[to]);
		} 
		else if(vis[to])
				low[u]=min(low[u],dfn[to]);	
	}
	if(dfn[u]==low[u])
	{
		IdNum++;
		int x;
		do
		{
			x=sta[top--];
			vis[x]=0;
			ID[x]=IdNum;
		}while(u!=x);
	}
}
int judge(int x)
{
	init();
	addedge(x);//建边

	for(int i=0;i<2*n;i++)
		if(!dfn[i])
			tarjan(i);//求强连通分量
			
	for(int i=0;i<n;i++)
		if(ID[i]==ID[i+n])
			return 0;
	return 1; 
}
int solve()
{
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		{
			if(i==j&&b[i][j])
				return 0;
			else if(b[i][j]!=b[j][i])
				return 0;
		}
	for(int i=0;i<31;i++)//从第1位到第31位
		if(judge(i)==0)
			return 0;
	return 1;
}
int main()
{
	while(cin>>n)
	{
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				scanf("%d",&b[i][j]);
		if(solve())
			printf("YES\n");
		else
			printf("NO\n");
	}
	return 0;
}