[2-SAT][POJ2683][HDU1814]2-SAT两种模板

2-SAT有两种模板,一种 O(M) 可以求任意可行解,另一种 O(MN) 十分暴力,所以可以求字典序最小或最大解

POJ2683

给定限制求任意可行解

这题就是 O(M) 的模板。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <queue>
#define N 4010
#define mp(a,b) make_pair(a,b)

using namespace std;
typedef pair<int,int> prs;
typedef pair<prs,prs> prs4;

prs4 p[N];
int n,cnt,tp,ti,g;
int G[N],sta[N],dfn[N],low[N],vis[N],bel[N],nG[N],v[N],d[N],opp[N];
struct edge{
  int t,nx;
}E[N*N<<1];
queue<int> Q;

inline bool cover(prs a,prs b){
  if(a.first>b.first) swap(a,b);
  return a.second>b.first;
}

inline void InserT(int x,int y){
  E[++cnt].t=x;E[cnt].nx=G[y];G[y]=cnt;
}

inline void Insert(int x,int y){
  E[++cnt].t=y;E[cnt].nx=nG[x];nG[x]=cnt;d[y]++;
}

void tarjan(int x){
  vis[x]=1;
  dfn[x]=low[x]=++ti;
  sta[++tp]=x;
  for(int i=G[x];i;i=E[i].nx){
    if(!vis[E[i].t]) tarjan(E[i].t);
    if(vis[E[i].t]==1) low[x]=min(low[x],low[E[i].t]);
  }
  if(low[x]==dfn[x]){
    ++g; int k;
    do{
      k=sta[tp--];
      bel[k]=g;
      vis[k]=2;
    }while(tp&&k!=x);
  }
}

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&n);
  for(int i=1;i<=n;i++){
    int s,s1,t,t1,l;
    scanf("%d:%d %d:%d %d",&s,&s1,&t,&t1,&l);
    s=s*60+s1; t=t*60+t1;
    p[i]=mp(mp(s,s+l),mp(t-l,t));
  }
  for(int i=1;i<=n;i++)
    for(int j=i+1;j<=n;j++){
      if(cover(p[i].first,p[j].first)) InserT(i<<1,j<<1|1),InserT(j<<1,i<<1|1);
      if(cover(p[i].first,p[j].second)) InserT(i<<1,j<<1),InserT(j<<1|1,i<<1|1);
      if(cover(p[i].second,p[j].first)) InserT(i<<1|1,j<<1|1),InserT(j<<1,i<<1);
      if(cover(p[i].second,p[j].second)) InserT(i<<1|1,j<<1),InserT(j<<1|1,i<<1);
    }
  for(int i=2;i<=(n<<1|1);i++)
    if(!dfn[i]) tarjan(i);
  for(int i=1;i<=n;i++) if(bel[i<<1]==bel[i<<1|1]) {puts("NO");return 0;}
  for(int k=2;k<=(n<<1|1);k++){
    for(int i=G[k];i;i=E[i].nx)
      if(bel[E[i].t]!=bel[k])
    Insert(bel[E[i].t],bel[k]);
    opp[bel[k]]=bel[k^1];
  }
  memset(vis,0,sizeof(vis));
  for(int i=1;i<=g;i++)
    if(!d[i]) Q.push(i);
  while(!Q.empty()){
    int x=Q.front(); Q.pop();
    if(!vis[x]) vis[x]=1,vis[opp[x]]=2;
    for(int i=nG[x];i;i=E[i].nx)
      if(!(--d[E[i].t])) Q.push(E[i].t);
  }
  puts("YES");
  for(int i=1;i<=n;i++){
    int s,s1,t,t1;
    if(vis[bel[i<<1|1]]==1){
      s=p[i].first.first/60,s1=p[i].first.first%60;
      t=p[i].first.second/60,t1=p[i].first.second%60;
      printf("%02d:%02d %02d:%02d\n",s,s1,t,t1);
    }
    else{
      s=p[i].second.first/60,s1=p[i].second.first%60;
      t=p[i].second.second/60,t1=p[i].second.second%60;
      printf("%02d:%02d %02d:%02d\n",s,s1,t,t1);
    }
  }
}

HDU 1814

给定限制,求字典序最小解

O(NM) 模板

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#define N 20010
#define M 40010

using namespace std;

int n,m,x,y,cnt,tp;
int G[N],sta[N],vis[N];
struct edge{
  int t,nx;
}E[M<<1];

inline void InserT(int x,int y){
  E[++cnt].t=y;E[cnt].nx=G[x];G[x]=cnt;
}

bool dfs(int x){
  if(vis[x]) return true;
  if(vis[x^1]) return false;
  sta[++tp]=x; vis[x]=1;
  for(int i=G[x];i;i=E[i].nx)
    if(!dfs(E[i].t)) return false;
  return true;             
}

inline bool solve(){
  memset(vis,0,sizeof(vis));
  for(int i=1;i<=n;i++)
    if(!vis[i<<1]&&!vis[i<<1|1]){
      tp=0;
      if(dfs(i<<1)) continue;
      for(;tp;tp--)vis[sta[tp]]=0;
      if(dfs(i<<1|1)) continue;
      return false;
    }
  for(int i=2;i<=(n<<1|1);i++) if(vis[i]) printf("%d\n",i-1);
  return true;
}

int main(){
  while(scanf("%d%d",&n,&m)==2){
    cnt=0;
    for(int i=1;i<=m;i++){
      scanf("%d%d",&x,&y);
      x++;y++;
      InserT(x,y^1);InserT(y,x^1);
    }
    if(!solve()) puts("NIE");
    memset(G,0,sizeof(G));
  }
  return 0;
}

表示本蒟蒻只会打模板,之前还把2-SAT搞成2-SET……
之后还得学各种2-SAT乱搞

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值