首页 > 其他 > 详细

2_sat

时间:2019-05-02 14:01:54      阅读:115      评论:0      收藏:0      [点我收藏+]

要求字典序的情况的话,爆搜

不要求的话

1:建图,有向边A--->B的意义为选择A则必须选择B,一般一个点的两种取值情况会拆点。

2:缩点。

3:建反向图,跑拓扑排序(有说不用建再跑,但我不懂为什么)。

4:根据实际情况输出。

例题:https://www.luogu.org/problemnew/show/P4782

代码

#include<cstdio>
#include<ctype.h>
#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
#define ll long long
inline ll rd()
{
    ll x=0,f=1;char c=getchar();
    while(!isdigit(c)){if(c==-) f=-f;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return x*f;
}
const int N=2e6+13;
vector<int>g[N];
struct zx{int nx,to;}e[N<<1];
int dfn[N],low[N],top,num,nm,cnt,a[N],v[N],bj[N],du[N],h[N],opp[N];
void TJ(int x)
{
    dfn[x]=low[x]=++num;a[++top]=x;
    for(int i=0,y;i<g[x].size();i++) if(!dfn[y=g[x][i]]) TJ(y),low[x]=min(low[x],low[y]);else if(!v[y]) low[x]=min(low[x],dfn[y]);
    if(low[x]!=dfn[x]) return ;v[x]=++cnt;
    while(a[top]!=x) v[a[top--]]=cnt;top--;
}
void add(int x,int y){e[++nm]={h[x],y};h[ax]=nm;du[y]++;}
int main()
{
    int n=rd(),m=rd(),i,c,b,j,l=0,r=0,x,y;
    for(int k=1;k<=m;k++)  i=rd(),c=rd(),j=rd(),b=rd(),g[i<<1|c].push_back(j<<1|(!b)),g[j<<1|b].push_back(i<<1|(!c));
    for(i=2;i<=(n<<1|1);i++) if(!dfn[i]) TJ(i);
    
    for(x=2;x<=(n<<1|1);x++) for(i=0;i<g[x].size();i++) if(v[g[x][i]]!=v[x]) add(v[g[x][i]],v[x]);
    for(i=1;i<=n;i++)if(v[i<<1]==v[i<<1|1]){printf("IMPOSSIBLE");return 0;}
    for(i=2;i<=(n<<1|1);i++) opp[v[i]]=v[i^1];
    for(int i=1;i<=cnt;i++) if(!du[i]) a[r++]=i;
    while(l<r)
    {
        x=a[l++];bj[x]=1;bj[opp[x]]=-1;
        for(i=h[x];i;i=e[i].nx) if(!--du[y=e[i].to]&&!bj[y]) a[r++]=y;
    }
    printf("POSSIBLE\n");
    for(int i=1;i<=n;i++) if(bj[v[i<<1]]==1)printf("1 ");else printf("0 ");
    return 0;
}

 

2_sat

原文:https://www.cnblogs.com/LWL--Figthing/p/10802186.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!