POJ 3683 2-SAT +输出方案

输出方案好烦啊 代码超级多~第一次写这种代码,照着别人的代码打的。。

 

建图:

对于婚礼i和婚礼j,i表示开始主持,i`表示结束主持,j也是一样

枚举每一对不同的i和j,

如果i和j冲突,连接ij`

如果i和j`冲突,连接ij

如果i`和j冲突,连接i`j`

如果i`和j冲突,连接i`j

 

输出方案:

我们可以把处在同一个强连通分量中的点和边缩成一个点,得到新的有向图G1。然后,我们把G中的所有弧反向,得到图G2。
我们把G中所有顶点置为“未着色”。按照拓扑顺序重复下面的操作:
1、选择第一个未着色的顶点x。把x染成红色。
2、把所有与x矛盾的顶点(如果存在bi,~bi∈B ,且bi属于x代表的强连通分量,~bi属于代表y的强连通分量,那么x和y就是互相矛盾的顶点)及其子孙全部全部染成蓝色。
3、重复操作1和2,直到不存在未着色的点为止。此时,G2中被染成红色的点在图G中对应的顶点集合,就对应着该2-SAT的一组解。
此实现可以递归。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<stack>
#define maxn 2010
#define maxe  maxn*maxn*3
using namespace std;

struct node
{
    int v,next;
}edge[maxe];

int n;
int s[maxn],t[maxn],seq[maxn];
int cnt,head1[maxn],head2[maxn],head3[maxn];
int scc,idex,dfn[maxn],low[maxn],belong[maxn];
int top,sstack[maxn];
bool instack[maxn],vis[maxn];
int d[maxn],col[maxn];

void addedge(int u,int v,int *head)
{
    edge[cnt].v=v;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}

void tarjan(int u)
{
    dfn[u]=low[u]=++idex;
    sstack[++top]=u;
    instack[u]=1;
    for(int i=head1[u];i!=-1;i=edge[i].next)
    {
        int v=edge[i].v;
        if(dfn[v]==0)
        {
            tarjan(v);
            low[u]=min(low[v],low[u]);
        }
        else if(instack[v])
        {
            low[u]=min(low[u],dfn[v]);
        }
    }
    if(dfn[u]==low[u])
    {
        scc++;
        while(1)
        {
            int tmp=sstack[top--];
            addedge(scc,tmp,head3);
            instack[tmp]=0;
            belong[tmp]=scc;
            if(tmp==u)
            break;
        }
    }
}

void build()
{
    int i,j;
    for(i=1;i<=2*n;i++)
    {
        for(j=head1[i];j!=-1;j=edge[j].next)
        {
            int v=edge[j].v;
            if(belong[i]!=belong[v])
            {
                addedge(belong[v],belong[i],head2);
            }
        }
    }
}

bool ok(int i,int j)
{
    if(s[i]>=t[j]||t[i]<=s[j])
    return false;
    return true;
}

bool check()
{
    int i;
    for(i=1;i<=n;i++)
    {
        if(belong[i<<1]==belong[i*2-1])
        return false;
    }
    return true;
}

void topo()
{
    int i,j;
    top=0;
    stack<int> sta;
    for(i=1;i<=scc;i++)
    {
        for(j=head2[i];j!=-1;j=edge[j].next)
        {
            d[edge[j].v]++;
        }
    }
    for(i=1;i<=scc;i++)
    {
        if(!d[i])
        sta.push(i);
    }
    while(!sta.empty())
    {
        int u=sta.top();
        sta.pop();
        seq[++top]=u;
        for(i=head2[u];i!=-1;i=edge[i].next)
        {
            int v=edge[i].v;
            d[v]--;
            if(!d[v])
            sta.push(v);
        }
    }
}

void paint(int u)
{
    int k;
    for(k=head2[u];k!=-1;k=edge[k].next)
    {
        if(col[edge[k].v]==0)
        {
            col[edge[k].v]=-1;
            paint(edge[k].v);
        }
    }
}

void color()
{
    int i,j,k;
    for(j=1;j<=top;j++)
    {
        int u=seq[j];
        if(col[u]==0)
        {
            col[u]=1;
            for(i=head3[u];i!=-1;i=edge[i].next)
            {
                int v=edge[i].v;
                if(col[belong[((v-1)^1)+1]]==0)
                {
                    col[belong[((v-1)^1)+1]]=-1;
                    paint(belong[((v-1)^1)+1]);
                }
            }
        }
    }
}

void print()
{
    int i,j;
    for(i=1;i<=n;i++)
    {
        if(col[belong[i*2]]==1)
        {
            printf("%02d:%02d %02d:%02d\n",s[i*2]/60,s[i*2]%60,t[i*2]/60,t[i*2]%60);
        }
        else
        {
            printf("%02d:%02d %02d:%02d\n",s[i*2-1]/60,s[i*2-1]%60,t[i*2-1]/60,t[i*2-1]%60);
        }
    }
}

int main()
{
    freopen("in.txt","r",stdin);
    int i,j;
    scanf("%d",&n);
    memset(head1,-1,sizeof(head1));
    memset(head2,-1,sizeof(head2));
    memset(head3,-1,sizeof(head3));
    cnt=0;
    for(i=1;i<=n;i++)
    {
        int u,v,tmp;
        scanf("%d:%d",&u,&v);
        s[2*i-1]=u*60+v;
        scanf("%d:%d",&u,&v);
        t[2*i]=u*60+v;
        scanf("%d",&tmp);
        t[2*i-1]=s[2*i-1]+tmp;
        s[2*i]=t[2*i]-tmp;
    }
    for(i=1;i<=n;i++)
    {
        for(j=i+1;j<=n;j++)
        {
            if(ok(2*i,2*j))
            {
                addedge(2*i,2*j-1,head1);
                addedge(2*j,2*i-1,head1);
            }
            if(ok(2*i-1,2*j)){
				addedge(2*i-1,2*j-1,head1);
				addedge(2*j,2*i,head1);
			}
			if(ok(2*i,2*j-1)){
				addedge(2*i,2*j,head1);
				addedge(2*j-1,2*i-1,head1);
			}
			if(ok(2*i-1,2*j-1)){
				addedge(2*i-1,2*j,head1);
				addedge(2*j-1,2*i,head1);
			}
        }
    }
    scc=0;
    for(i=1;i<=2*n;i++)
    {
        if(!dfn[i])
        tarjan(i);
    }
    if(!check())
    printf("NO\n");
    else
    {
        printf("YES\n");
        build();
        topo();
        color();
        print();
    }
    return 0;
}







 

 

 

 

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值