2-sat——operator_

1.什么是 2 − s a t 2-sat 2sat

个人理解 : : :

2 − s a t 2-sat 2sat问题指的就是有若干个集合,每个集合内部有 2 2 2个点,这两个点只能且必须取一个(其实就是一个 b o o l bool bool变量的 2 2 2种取值),在集合外,还有若干个关系,形如 : : : x 1 x_1 x1 y 2 y_2 y2 中至少取一个,求是否有满足所有关系的可行解(可能还会要求构造)。

形式化一点,就是求构造 b o o l bool bool数列 { x } \{x\} {x}使得 ( x 1 ∨ x 2 ) ∧ ( x 1 ∨ ¬ x 4 ) ∧ ( ¬ x 3 ∨ x 4 ) ∧ ( x 2 ∨ x 5 ) ∧ ⋯ = = t r u e (x_1 \vee x_2)\land(x_1 \vee \neg x_4)\land(\neg x_3 \vee x_4)\land(x_2 \vee x_5)\land\cdots==true (x1x2)(x1¬x4)(¬x3x4)(x2x5)==true

2.怎么解决?

考虑图论建模:一条边 u → v u\rightarrow v uv就代表选了 u u u就必须选 v v v
那么一个 b o o l bool bool变量就需要分为 2 2 2个点,由此可得:

( x 1 ∨ x 2 ) : ¬ x 1 → x 2 , ¬ x 2 → x 1 (x_1 \vee x_2):\neg x_1\rightarrow x_2,\neg x_2\rightarrow x_1 (x1x2):¬x1x2,¬x2x1

( ¬ x 1 ∨ x 2 ) : x 1 → x 2 , ¬ x 2 → ¬ x 1 (\neg x_1 \vee x_2):x_1\rightarrow x_2,\neg x_2\rightarrow \neg x_1 (¬x1x2):x1x2,¬x2¬x1

( ¬ x 1 ∨ ¬ x 2 ) : x 1 → ¬ x 2 , x 2 → ¬ x 1 (\neg x_1 \vee \neg x_2):x_1\rightarrow \neg x_2,x_2\rightarrow \neg x_1 (¬x1¬x2):x1¬x2,x2¬x1

至此建模完成。下面考虑求解。

显然有朴素的dfs,基于该操作:选一个点就把它的后继全选上,把它的对面(即 ¬ x \neg x ¬x)的前驱全删除。但这样不仅麻烦,复杂度还特高。

所以想优化,注意到其中会形成环,每个环中的点必然状态相同,要么全选,要么全不选。可以想到强连通分量的缩点。于是用 t a r j a n tarjan tarjan跑一遍缩点模板,如果一个变量的 2 2 2个状态的对应点都在同一个强连通分量中,就不可能有解,反之,则有解。

可能还需要构造,那么观察到缩点后就可以形成有向无环图,又想到拓扑序,但拓扑序是无前继,所以反着处理,就可以达到无后继的效果,按拓扑序选择一遍就行了(其实若 t a r j a n tarjan tarjan写得优雅都不用再来一遍拓扑)。

这样,就可以完美解决 2 − s a t 2-sat 2sat问题了。

3.满汉全席

模板,代码放一下吧。

#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read() {
	int s=0,m=0;char ch=getchar();
	while(!isdigit(ch)) {if(ch=='-')m=1;ch=getchar();}
	while( isdigit(ch)) s=(s<<3)+(s<<1)+(ch^48),ch=getchar();
	return m?-s:s;
}
int h[205],cnt;
struct QWQ{int v,nxt;} e[2005];
void add(int u,int v) {e[++cnt]={v,h[u]},h[u]=cnt;}
int t,n,m,x1,x2;
char ch1,ch2;//m为2*x,h为2*x-1 
stack<int> st;
int f[205],dfn[205],low[205],c,k,scc[205],flag;
void tarjan(int u,int fa) {
	st.push(u);
	f[u]=1,dfn[u]=low[u]=++c;
	for(int i=h[u];i;i=e[i].nxt) {
		int v=e[i].v;
		if(!dfn[v])
			tarjan(v,u),low[u]=min(low[u],low[v]);
		else if(f[v])
			low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u]) {
		k++;
		while(1) {
			int top=st.top();st.pop();
			scc[top]=k,f[top]=0;//top在强连通分量k中 
			if(top==u) break;
		}
	}
} 
signed main() {
	cin>>t;
	while(t--) {
		cnt=k=c=flag=0;
		memset(h,0,sizeof(h));
		memset(f,0,sizeof(f));
		memset(dfn,0,sizeof(dfn));
		memset(low,0,sizeof(low));
		n=read(),m=read();
		for(int i=1;i<=m;i++) {
			cin>>ch1,x1=read(),cin>>ch2,x2=read();
			add(2*x1-(ch1=='m'),2*x2-(ch2=='h'));
			add(2*x2-(ch2=='m'),2*x1-(ch1=='h'));
		}
		for(int i=1;i<=2*n;i++) 
			if(!dfn[i]) tarjan(i,i);
		for(int i=1;i<=2*n;i+=2)
			if(scc[i]==scc[i+1]) flag=1;
		if(!flag) printf("GOOD\n");
		else printf("BAD\n");
	}
	return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值