poj 3678 Katu Puzzle

题目链接:http://poj.org/problem?id=3678

题目大意:判断是否有一种分配(给结点一个值0或1)方案可以满足输入的要求(某两个结点的与,或,异或操作为c).

题目思路:
经典2-SAT问题(表示自己比较弱,对没接触过的算法,想破脑袋想不到...看了一下午才明白点).
既然知道是2-SAT,重点就在构图了,我们以结点a,b(~a,~b为其非)为例.
a AND b == 1 (a,b同时为1)       a<-->b,  ~a-->a,  ~b-->b;
a AND b == 0 (a,b不同时为1)   a-->~b,  b-->~a;
a OR  b == 1  (a,b不同时为0)   ~a-->b,  ~b-->a;
a OR  b == 0  (a,b同时为0)      
~a<-->~b,  a-->~a,  b-->~b;
a XOR b == 1 (a,b不相同)      
a<-->~b,  b<-->~a;
a XOR b == 0 (a,b相同)           a<-->b,~a<-->~b;
对于这题,不需要求方案,只要求方案是否存在,所以我们不需要缩点,而且还可以把环
(黄字部分)从邻接表中去掉.


代码:

#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <time.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;

#define ull unsigned __int64
//#define ll __int64
//#define ull unsigned long long
#define ll long long
#define son1 New(p.xl,xm,p.yl,ym),(rt<<2)-2
#define son2 New(p.xl,xm,min(ym+1,p.yr),p.yr),(rt<<2)-1
#define son3 New(min(xm+1,p.xr),p.xr,p.yl,ym),rt<<2
#define son4 New(min(xm+1,p.xr),p.xr,min(ym+1,p.yr),p.yr),rt<<2|1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define middle (l+r)>>1
#define MOD 1000000007
#define esp (1e-8)
const int INF=0x3F3F3F3F;
const double DINF=10000.00;
//const double pi=acos(-1.0);
const int M=2010;
int min(int x,int y){return x<y? x:y;}
int max(int x,int y){return x>y? x:y;}
void swap(int& x,int& y){int t=x;x=y;y=t;}
int T,cas;

vector<int>init[M];
int low[M],pre[M],idx[M],ss[M];
int n,nn,Index,top,m,N;

void dfs(int u){
	low[u]=pre[u]=++Index;
	ss[++top]=u;
	for(int i=0;i<init[u].size();i++){
		int v=init[u][i];
		if(pre[v]==-1) dfs(v),low[u]=min(low[u],low[v]);
		else if(idx[v]==-1) low[u]=min(low[u],pre[v]);
	}
	if(low[u]==pre[u]){
		int v=-1;
		++nn;
		while(u!=v) idx[v=ss[top--]]=nn;
	}
}

void Tarjan(){
	memset(idx,-1,sizeof(int)*(n+1));
	memset(pre,-1,sizeof(int)*(n+1));
	Index=top=nn=0,ss[0]=-1;
	for(int u=1;u<=n;u++) if(pre[u]==-1) dfs(u);
}

bool Judge(){
	for(int u=1;u<=N;u++) if(idx[u]==idx[u+N]) return false;
	return true;
}

void _sof(){
	Tarjan();
	printf("%s\n",Judge()? "YES":"NO");
}

void _init(){
	int a,b,c;
	char op[4];
	n=N<<1|1;
	for(int i=1;i<=n;i++) init[i].clear();
	while(m--){
		scanf("%d%d%d%s",&a,&b,&c,op);
		a++,b++;
		if(c){
			if(op[0]=='A'){
				init[a+N].push_back(a);
				init[b+N].push_back(b);
			}else if(op[0]=='O'){
				init[a+N].push_back(b);
				init[b+N].push_back(a);
			}
		}else{
			if(op[0]=='A'){
				init[a].push_back(b+N);
				init[b].push_back(a+N);
			}else if(op[0]=='O'){
				init[a].push_back(a+N);
				init[b].push_back(b+N);
			}
		}
	}
	_sof();
}

int main(){
	//freopen("1.in","r",stdin);
	//freopen("1.out","w",stdout);
	//_init()
	//scanf("%d",&T);for(cas=1;cas<=T;cas++) _init();
	while(~scanf("%d%d",&N,&m))  _init();
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值