poj 3683 Priest John's Busiest Day【2-SAT+tarjan+拓扑】

转换成2-SAT模型，建边是如果时间(i,j)冲突就连边(i,j')，其他同理
tarjan缩点，判可行性
返图拓扑，输出方案

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int N=5005;
int n,a[N],b[N],h[N],cnt,dfn[N],low[N],tot,s[N],top,bl[N],col,du[N],id[N],t[N];
bool v[N];
vector<pair<int,int> >p;
struct qwe
{
    int ne,no,to;
}e[N*N];
int read()
{
    int r=0,f=1;
    char p=getchar();
    while(p>'9'||p<'0')
    {
        if(p=='-')
            f=-1;
        p=getchar();
    }
    while(p>='0'&&p<='9')
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r*f;
}
void add(int u,int v)
{
    cnt++;
    e[cnt].ne=h[u];
    e[cnt].no=u;
    e[cnt].to=v;
    h[u]=cnt;
}
void tarjan(int u)
{
    dfn[u]=low[u]=++tot;
    v[s[++top]=u]=1;
    for(int i=h[u];i;i=e[i].ne)
    {
        if(!dfn[e[i].to])
        {
            tarjan(e[i].to);
            low[u]=min(low[u],low[e[i].to]);
        }
        else if(v[e[i].to])
            low[u]=min(low[u],dfn[e[i].to]);
    }
    if(dfn[u]==low[u])
    {
        col++;
        while(s[top]!=u)
        {
            bl[s[top]]=col;
            v[s[top--]]=0;
        }
        bl[s[top]]=col;
        v[s[top--]]=0;
    }
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        int x=read(),xx=read(),y=read(),yy=read(),z=read();
        x=x*60+xx,y=y*60+yy;
        a[i]=x,b[i]=x+z,a[i+n]=y-z,b[i+n]=y;//cerr<<a[i]<<" "<<b[i]<<"   "<<c[i]<<" "<<d[i]<<endl;
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(i!=j)
            {
                if(!(a[i]>=b[j]||a[j]>=b[i]))
                    add(i,j+n);
                if(!(a[i]>=b[j+n]||a[j+n]>=b[i]))
                    add(i,j);
                if(!(a[i+n]>=b[j]||a[j]>=b[i+n]))
                    add(i+n,j+n);
                if(!(a[i+n]>=b[j+n]||a[j+n]>=b[i+n]))
                    add(i+n,j);
            }
    for(int i=1;i<=n+n;i++)
        if(!dfn[i])
            tarjan(i);
    for(int i=1;i<=n;i++)
    {
        id[bl[i]]=bl[i+n];
        id[bl[i+n]]=bl[i];
        if(bl[i]==bl[i+n])
        {
            puts("NO");
            return 0;
        }
    }
    puts("YES");
    for(int i=1;i<=cnt;i++)
        if(bl[e[i].no]!=bl[e[i].to])
            p.push_back(make_pair(bl[e[i].to],bl[e[i].no]));
    memset(h,0,sizeof(h));
    cnt=0;
    for(int i=0,len=p.size();i<len;i++)
        add(p[i].first,p[i].second),du[p[i].second]++;
    queue<int>q;
    for(int i=1;i<=col;i++)
        if(!du[i])
            q.push(i);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(!t[u])
            t[u]=1,t[id[u]]=2;
        for(int i=h[u];i;i=e[i].ne)
            if(!(--du[e[i].to]))
                q.push(e[i].to);
    }
    for(int i=1;i<=n;i++)
    {
        if(t[bl[i]]==1)
            printf("%02d:%02d %02d:%02d\n",a[i]/60,a[i]%60,b[i]/60,b[i]%60);
        else
            printf("%02d:%02d %02d:%02d\n",a[i+n]/60,a[i+n]%60,b[i+n]/60,b[i+n]%60);
    }
    return 0;
}

转载于:https://www.cnblogs.com/lokiii/p/9978811.html