pku3678 2_SAT现有N个数,每个数不是0就是1。现给出一些有关这些数的与,或,异或的关系,问可能不?

        题意:现有N个数,每个数不是0就是1。现给出一些有关这些数的与,或,异或的关系,问可能不?

        分析:对于每个数而言,不是0就是1,对于每个数给定两个状态,一个表示0,一个表示1,只要每两个数最后不在同一个强连通分量里面就行。不过这题给出与,或,异或的关系,如果两数与为1的话,明显就可以确定那两个数都是1了,所以我把已经能确定的数标记出来。先对已经确定的数深搜一遍,看是不是存在同一个数的两个状态在同一个强连通分量中。。。而后再用Tarjan算法划分强连通分量,只要每个状态属于不同的强连通分量即输出YES.]

        优化:其实对于每个已经确定的数,比如第i个数确定0,那么只要加一条从第i个数为1的状态到第i个数为0的状态的有向边即可。即如果x必须选,那么就加一条~x到x的有向边。

         心得:用链接表存储图时,用vector数组比自己写方便多了,以后要改过来;tarjan与kosaraju会的话,时间复杂度上没多大区别,因为kosaraju需要建正反图,建图简单的可以用kosaraju,建图复杂的用tarjan;最近状态还不错,一些题目通过自己深入思考,可以解决,不仅锻炼了能力,而且心情也蛮爽的。深入思考才是王道啊。。。不像以前,想不了多久就看解题报告,或看discuss(也许这是做usaco锻炼出来的吧)。。。

//注释了的为第一种方法(没优化的)
#include<stdio.h>
#include<vector>
#include<iostream>
using namespace std;

int low[2100],DFN[2100],cnt[2100],q[2100],top,tnum,num,n,m;
bool flag[2100],vis[2100],inq[2100];
vector<int> g[2100];

//将确定的状态不断扩大
/*void DFS1(int t)
{
	flag[t]=true;
	int i;
	for(i=0;i<g[t].size();i++)
		if(!flag[g[t][i]])
			DFS1(g[t][i]);
}*/
//tarjan算法求强连通分量
void DFS(int t)
{
	int i,j;
	vis[t]=true;
	low[t]=DFN[t]=++tnum;
	q[top++]=t;
	inq[t]=true;
	for(i=0;i<g[t].size();i++)
	{
		j=g[t][i];
		if(!vis[j])
		{
			DFS(j);
			if(low[t]>low[j])
				low[t]=low[j];
		}
		else if(inq[j]&&low[t]>DFN[j])
			low[t]=DFN[j];
	}

	if(low[t]==DFN[t])
	{
		++num;
		do
		{
			cnt[q[--top]]=num;
			inq[q[top]]=false;
		}while(q[top]!=t);
	}
}
int main()
{
	int i,j,a,b,c;
	char ch[10];
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		memset(flag,false,sizeof(flag));
		for(i=0;i<m;i++)
		{
			scanf("%d %d %d %s",&a,&b,&c,ch);
			if(ch[0]=='A')
			{
				if(c==1)
				{
				//	flag[2*a+1]=flag[2*b+1]=true;//相与为1,肯定都为1
					g[2*a].push_back(2*a+1);
					g[2*b].push_back(2*b+1);
				}
				else
				{
					g[2*a+1].push_back(2*b);//相与为0,所以一个为1时,另一个必须为0
					g[2*b+1].push_back(2*a);
				}
			}
			else if(ch[0]=='O')
			{
				if(c==1)
				{
					g[2*a].push_back(2*b+1);//相或为1时,一个为0,另一个必须为1
					g[2*b].push_back(2*a+1);
				}
				else
				{
				//	flag[2*a]=flag[2*b]=true;
					g[2*a+1].push_back(2*a);//相或为0时,两个都为0
					g[2*b+1].push_back(2*b);
				}
			}
			else
			{
				if(c==0)
				{
					g[2*a].push_back(2*b);//异或为0时,两个数状态一定相同
					g[2*b].push_back(2*a);
					g[2*a+1].push_back(2*b+1);
					g[2*b+1].push_back(2*a+1);
				}
				else
				{
					g[2*a].push_back(2*b+1);//异或为1时,两个数状态一定相异
					g[2*b].push_back(2*a+1);
					g[2*a+1].push_back(2*b);
					g[2*b+1].push_back(2*a);
				}
			}
		}

	    //已经确定取0还是1的接着确定
	/*	for(i=0;i<2*n;i++)
			if(flag[i])
				DFS1(i);

		for(i=0;i<2*n;i+=2)
			if(flag[i]&&flag[i+1])
			{
				printf("NO\n");
				break;
			}
		if(i<2*n) continue;*/
		
		memset(vis,false,sizeof(vis));
		memset(inq,false,sizeof(inq));
		memset(cnt,-1,sizeof(cnt));
		for(top=0,num=0,tnum=0,i=0;i<2*n;i++)
			if(!vis[i])
				DFS(i);
		//会不会在同一个强连通分量
		for(i=0;i<2*n;i+=2)
			if(cnt[i]!=-1&&cnt[i]==cnt[i+1])
			{
				printf("NO\n");
				break;
			}
		if(i>=2*n) printf("YES\n");
	}
	return 0;
}


 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值