2 - SAT 建边

建边

1.我们利用一条有向边<i,j>,来表示选i的情况下,一定要选j;

2.用i表示某个点是true,那么i’表示某个点是false

3.建边:

A. x 和 y 必选(都为1):^x -> x, ^y->y
B. x 和 y 都不能选(都为0):x -> ^x, y -> ^y
C.x 和 y 至少有一个为 1 ( x 为 1 时 y 任取或者 y 为 1 时 x 任取 / x和y不能都选0):^x->y, ^y->x
D. x和y至少有一个为0(x为0时y任取或者y为0时x任取 / x和y不能都选1):y->^x, x -> ^y
E. x和y不同(一个为1/0另一个为0/1):x -> ^y, y -> ^x , ^y -> x, ^x -> y
F. x和y相同(同为1/0) :x -> y ,^x -> ^y ,y -> x , ^y -> ^x
G. x选1时,y不能选0 :x->y, ^y-> ^x
H. x选0时,y不能选1 :^x -> ^y, y -> x

解决方案

1.求字典序最小的解的方法:

暴力dfs求解(复杂度O(N*M))

2.判断当前的2-sa问题t是否有解

tarjan强连通缩点,加判断(复杂度O(N+M))

3.求出当前的2-sat问题的任意一组解

tarjan强连通缩点+拓扑排序+构建一组解(复杂度O(N+M))

至于想要知道当前元素归到了哪一边,看bl,哪个小就是哪个。

1.求字典序最小的解的方法:

题目:HDU1814

题意

有n队人(2i和2i+1一队),有m种限制关系,第x个人出现第y个人不能出现,输出字典序最小的一种解

#include<iostream>
#include<vector>
using namespace std;
const int maxn=1e5+5;
int n,d,h1,h2,m,dfn[maxn],low[maxn],stk[maxn],vis[maxn],color[maxn],cnt,num,top,r,t;
vector<int>gap[maxn];
bool dfs(int x){
    if(color[x^1]) return 0;//当前节点的反点已经走过,则不合法
    if(color[x]) return 1;//
    stk[top++]=x;
    color[x]=1;
	for(int i=0;i<gap[x].size();i++){
		int y=gap[x][i];
		if(!dfs(y)) return 0;
	}
	return 1;
}
bool check(){
	for(int i=0;i<2*n;i+=2){
		if(!color[i]&&!color[i^1]){//都没被染过色
			top=0;
		    if(!dfs(i)){//从i开始dfs,不成立的话,初始化dfs过的点的颜色
			    for(int j=0;j<top;j++) 
			        color[stk[j]]=color[stk[j]^1]=0;
			    if(!dfs(i^1))//从i^1(i+1)开始dfs,若还不成功,则无解,返回 0
			    return 0;
	   	    }
		}
	} 
    return 1;
}
int main(){
	while(scanf("%d%d",&n,&m)!=EOF){//n从0开始到n-1
	    int a,b;
	    for(int i=0;i<=2*n;i++) gap[i].clear(),color[i]=stk[i]=0;
	    num=top=cnt=0;
	    for(int i=1;i<=m;i++){
	        scanf("%d %d",&a,&b);
	        a--,b--;
		    gap[a].push_back(b^1),gap[b].push_back(a^1); 
	    }
	    if(check()){
		    for(int i=0;i<2*n;i++)
		        if(color[i])
		    printf("%d\n",i+1);
	    }
        else
        printf("NIE\n");
    }
	return 0;    
}

2.3 判断是否有解,输出任意一解

洛谷UVA1391 Astronauts
思路
选择 i i i 和 ^ i i i 中拓扑序大的,因为拓扑序大对结果的影响小 。tarjan后的color数组存的是逆拓扑序,所有可以直接选择 c o l o r [ i ] color[i] color[i] c o l o r [ color[ color[^ i ] i] i] 中拓扑序“小”的

#include<iostream>
#include<vector>
using namespace std;
const int maxn=4e5+5;
int n,d,h1,h2,m,dfn[maxn],low[maxn],stk[maxn],vis[maxn],color[maxn],cnt,num,top,r,t,sum,age[maxn];
vector<int>gap[maxn];
void tarjan(int x){
	dfn[x]=low[x]= ++cnt;
	stk[++top]=x;
	vis[x]=1;
	for(int i=0;i<gap[x].size();i++){
		int y=gap[x][i];
		if(!dfn[y]){
			tarjan(y);
			low[x]=min(low[x],low[y]);
		}
		else if(vis[y])
		low[x]=min(low[x],dfn[y]);
	}
	if(low[x]==dfn[x]){
		int tp=0;
		++num;
		do{
			tp=stk[top--];
			vis[tp]=0;
			color[tp]=num;
		}while(x!=tp);
	}
} 
int main(){
	while(scanf("%d%d",&n,&m)!=EOF&&n!=0&&m!=0){
	    int a,b;
	    for(int i=1;i<=2*n;i++) gap[i].clear(),color[i]=stk[i]=low[i]=dfn[i]=vis[i]=0;
	    num=top=cnt=sum=0;
	    for(int i=1;i<=n;i++){
	    	scanf("%d",&age[i]);
	    	sum+=age[i];
		}
		for(int i=1;i<=n;i++)
		    if(age[i]*n>=sum)age[i]=1;
		    else age[i]=0;
	    for(int i=1;i<=m;i++){
	        scanf("%d %d",&a,&b);
		    if(age[a]==age[b])gap[a+n].push_back(b),gap[b+n].push_back(a);
		    gap[a].push_back(b+n),gap[b].push_back(a+n);
	    }
	    for(int i=1;i<=2*n;i++)
	        if(!dfn[i])
	        tarjan(i);
	    int s=1;
		for( s=1;s<=n;s++){
		    if(color[s]==color[s+n]){ 
		    	printf("No solution.\n"); 
				break;
			} 
		}
		if(s==n+1){ 
		    for(int i=1;i<=n;i++) 
			    if(color[i]>color[i+n])//反拓扑序,i+n的拓扑序小  
			    	if(age[i])printf("A\n"); 
					else printf("B\n"); 	
				else { 
					printf("C\n");
				}
		}
    }
	return 0;    
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值