[2-SAT 任意解 Tarjan 模板题] POJ 3683 Priest John's Busiest Day

得好好说说2-SAT这个东西啦
我们研究下这篇论文 由对称性解2-SAT问题–伍昱
论文中提到的一种方法是 对建完的图跑一边Tarjan 跑出强连通分量
如果两个相关变量也就是 x x在同一个SCC里 那么就是无解
否则提出了一种基于缩点后反图拓扑排序和染色的一种输出方案的方法
这样我觉得麻烦了点
首先 反图的拓扑序 就是我们生成SCC的顺序 这很显然 而且反图不反图也无所谓 反正也就是拓扑序反一下的事
然后 完全不必要染色
我们考虑 x 所在的分量s x 所在的分量 s
我们选择原图拓扑序在后的分量包含的那个条件成立 这样就不会出现由这个条件连向相反条件的一条路径

代码短短的很兹瓷啊

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef pair<int,int> abcd;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
  return *p1++;
}
inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=2005;
struct edge{
  int u,v,next;
}G[N*N];
int head[N],inum;
inline void add(int u,int v,int p){
  G[p].u=u; G[p].v=v; G[p].next=head[u]; head[u]=p;
}

int pre[N],low[N],clk; int cnt,scc[N]; int Stack[N],pnt;
#define V G[p].v
inline void dfs(int u){
  pre[u]=low[u]=++clk; Stack[++pnt]=u;
  for (int p=head[u];p;p=G[p].next)
    if (!pre[V]){
      dfs(V);
      low[u]=min(low[u],low[V]);
    }else if (!scc[V])
      low[u]=min(low[u],pre[V]);
  if (low[u]==pre[u]){
    ++cnt;
    while (Stack[pnt]!=u) scc[Stack[pnt--]]=cnt; scc[Stack[pnt--]]=cnt;
  }
}

abcd eve[N];

inline int jud(abcd a,abcd b){  
  return !(a.second<=b.first || b.second<=a.first); 
}

int n,ans[N];

int main(){
  int a1,b1,a2,b2,d;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(n);
  for (int i=0;i<n;i++){
    read(a1),read(b1),read(a2),read(b2),read(d);
    eve[i<<1]=abcd(a1*60+b1,a1*60+b1+d),eve[i<<1|1]=abcd(a2*60+b2-d,a2*60+b2);
  }
  for (int i=0;i<2*n;i++)  
    for (int j=i+1;j<2*n;j++)
      if (((i>>1)!=(j>>1)) && jud(eve[i],eve[j]))
    add(i,j^1,++inum),add(j,i^1,++inum);
  for (int i=0;i<2*n;i++)
    if (!pre[i])
      dfs(i);
  for (int i=0;i<2*n;i+=2)
    if (scc[i]==scc[i^1]){
      printf("NO\n"); return 0;
    }else
      ans[i>>1]=scc[i]<scc[i+1]?i:i+1;
  printf("YES\n");
  for (int i=0;i<n;i++){  
    a1=eve[ans[i]].first/60,b1=eve[ans[i]].first%60;  
    a2=eve[ans[i]].second/60,b2=eve[ans[i]].second%60;  
    printf("%02d:%02d %02d:%02d\n",a1,b1,a2,b2);  
  }
  return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值