poj3683

Posted mikuknight

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了poj3683相关的知识,希望对你有一定的参考价值。

2-SAT算法流程

1.建立选择边

2.tarjan缩点(同一块内必然同时选择),判断是否存在可行解(若存在一对对立点{2*i,2*i-1}在同一块内,则不可行)

3.把缩点后的块之间建立反向选择边,并传递不选择命令(code中用op[ ]数组记录)

4.按照拓扑序找到一为选择的点(块),标为选择,传递不选择标记

//poj3683
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#define inf 2000000000
#define ll long long 
#define mod 1000000007
using namespace std;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<0||ch>9){if(ch==-)f=-1;ch=getchar();}
    while(ch>=0&&ch<=9){x=x*10+ch-0;ch=getchar();}
    return x*f;
}
int n,cnt,top,scc,ind;
int a[2005],b[2005],belong[2005],op[2005];
int dfn[2005],low[2005],q[2005];
bool inq[2005];
struct edge{int to,next;}e[2000005],ed[2000005];
int last[2005],last2[2005],d[2005];
int color[2005];
bool jud(int x,int y){
    if(b[x]<=a[y]||a[x]>=b[y])return 0;
    return 1;
}
void insert(int u,int v){e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;}
void insert2(int u,int v){d[v]++;ed[++cnt].to=v;ed[cnt].next=last2[u];last2[u]=cnt;}
void build(){
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++){
            if(jud(2*i,2*j)){insert(2*i,2*j-1);insert(2*j,2*i-1);}
            if(jud(2*i,2*j-1)){insert(2*i,2*j);insert(2*j-1,2*i-1);}
            if(jud(2*i-1,2*j)){insert(2*i-1,2*j-1);insert(2*j,2*i);}
            if(jud(2*i-1,2*j-1)){insert(2*i-1,2*j);insert(2*j-1,2*i);}
        }
}
void tarjan(int x){
    dfn[x]=low[x]=++ind;
    q[++top]=x;inq[x]=1;
    for(int i=last[x];i;i=e[i].next)
        if(!dfn[e[i].to]){tarjan(e[i].to);low[x]=min(low[e[i].to],low[x]);}
        else if(inq[e[i].to])low[x]=min(dfn[e[i].to],low[x]);
    if(low[x]==dfn[x]){
        int now=0;scc++;
        while(now!=x){
            now=q[top--];
            inq[now]=0;
            belong[now]=scc;
        }
    }
}
void rebuild(){
    cnt=0;
    for(int x=1;x<=2*n;x++)
        for(int i=last[x];i;i=e[i].next)
            if(belong[x]!=belong[e[i].to])
                insert2(belong[e[i].to],belong[x]);
}
void dfs(int x){
    if(color[x])return;
    color[x]=-1;
    for(int i=last2[x];i;i=ed[i].next)dfs(ed[i].to);
}
void topsort(){
    for(int i=1;i<=scc;i++)if(!d[i])q[++top]=i;
    while(top){
        int now=q[top--];
        if(color[now])continue;
        color[now]=1;dfs(op[now]);
        for(int i=last2[now];i;i=ed[i].next){
            d[ed[i].to]--;
            if(!d[ed[i].to])q[++top]=ed[i].to;
        }
    }
}
void print(int x){printf("%.2d:",x/60);printf("%.2d ",x%60);}
int main(){
    n=read();
    int x;
    for(int i=1;i<=n;i++){
        a[2*i]=read();a[2*i]=a[2*i]*60+read();
        b[2*i-1]=read();b[2*i-1]=b[2*i-1]*60+read();
        x=read();b[2*i]=a[2*i]+x;a[2*i-1]=b[2*i-1]-x;
    }
    build();
    for(int i=1;i<=2*n;i++)if(!dfn[i])tarjan(i);
    for(int i=1;i<=n;i++)
        if(belong[2*i]==belong[2*i-1]){puts("NO");return 0;}
    puts("YES");
    rebuild();
    for(int i=1;i<=n;i++){
        op[belong[2*i]]=belong[2*i-1];
        op[belong[2*i-1]]=belong[2*i];
    }
    topsort();
    for(int i=1;i<=n;i++)
        if(color[belong[2*i]]==1)
            print(a[2*i]),print(b[2*i]),puts("");
        else print(a[2*i-1]),print(b[2*i-1]),puts("");
    return 0;
}

 

以上是关于poj3683的主要内容,如果未能解决你的问题,请参考以下文章

POJ3683 Falsita

poj3683

poj3683 2-SAT 同上一道

poj 3683(2-SAT+SCC)

POJ 3683 Priest John's Busiest Day(2-SAT)

poj 3683 Priest John's Busiest Day - 2-sat