loj 1251(2-sat + 输出一组可行解）

题目链接：http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=26961

思路：u表示留下，~u表示离开，同理v,对于+u,-v，我们可以这样来定义：若u离开，则v必须留下，如v离开，则u必须留下,于是我们可以连边u+n->v,v+n->u,后面的同理。

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<algorithm>
  5 #include<stack>
  6 #include<vector>
  7 #include<queue>
  8 using namespace std;
  9 #define MAXN 20000
 10 #define MAXM 444444
 11 
 12 struct Edge{
 13     int v,next;
 14 }edge[MAXM];
 15 
 16 int n,m,NE;
 17 int head[MAXN];
 18 
 19 void Insert(int u,int v)
 20 {
 21     edge[NE].v=v;
 22     edge[NE].next=head[u];
 23     head[u]=NE++;
 24 }
 25 
 26 int cnt,bcc_count;
 27 int low[MAXN],dfn[MAXN],color[MAXN];
 28 bool mark[MAXN];
 29 stack<int>S;
 30 
 31 void Tarjan(int u)
 32 {
 33     low[u]=dfn[u]=++cnt;
 34     mark[u]=true;
 35     S.push(u);
 36     for(int i=head[u];i!=-1;i=edge[i].next){
 37         int v=edge[i].v;
 38         if(dfn[v]==0){
 39             Tarjan(v);
 40             low[u]=min(low[u],low[v]);
 41         }else if(mark[v]){
 42             low[u]=min(low[u],dfn[v]);
 43         }
 44     }
 45     if(low[u]==dfn[u]){
 46         int v;
 47         bcc_count++;
 48         do{
 49             v=S.top();
 50             S.pop();
 51             mark[v]=false;
 52             color[v]=bcc_count;
 53         }while(u!=v);
 54     }
 55 }
 56 
 57 int opp[MAXN];
 58 bool Judge()
 59 {
 60     for(int i=1;i<=n;i++){
 61         if(color[i]==color[i+n])return false;
 62         //opp[x]保存的是和编号为x的连通分量矛盾的连通分量的编号
 63         opp[color[i]]=color[i+n];
 64         opp[color[i+n]]=color[i];
 65     }
 66     return true;
 67 }
 68 
 69 vector<vector<int> >g;
 70 vector<int>ans;
 71 int degree[MAXN];
 72 int vis[MAXN];
 73 
 74 void TopSort()
 75 {
 76     queue<int>que;
 77     for(int i=1;i<=bcc_count;i++){
 78         if(degree[i]==0)que.push(i);
 79     }
 80     while(!que.empty()){
 81         int u=que.front();
 82         que.pop();
 83         if(vis[u]==0){
 84             //染色
 85             vis[u]=1;
 86             vis[opp[u]]=-1;
 87         }
 88         for(int i=0;i<g[u].size();i++){
 89             degree[g[u][i]]--;
 90             if(degree[g[u][i]]==0)que.push(g[u][i]);
 91         }
 92     }
 93 }
 94 
 95 
 96 int main()
 97 {
 98     int _case,t=1,a,b,u,v;
 99     scanf("%d",&_case);
100     while(_case--){
101         scanf("%d %d",&m,&n);
102         NE=0;
103         memset(head,-1,sizeof(head));
104         while(m--){
105             scanf("%d%d",&a,&b);
106             u=abs(a),v=abs(b);
107             if(a>0&&b>0)Insert(u+n,v),Insert(v+n,u);
108             else if(a>0&&b<0)Insert(u+n,v+n),Insert(v,u);
109             else if(a<0&&b>0)Insert(u,v),Insert(v+n,u+n);
110             else Insert(u,v+n),Insert(v,u+n);
111         }
112         cnt=bcc_count=0;
113         memset(mark,false,sizeof(mark));
114         memset(dfn,0,sizeof(dfn));
115         memset(color,0,sizeof(color));
116         for(int i=1;i<=2*n;i++)if(dfn[i]==0)Tarjan(i);
117         printf("Case %d: ",t++);
118         if(!Judge()){
119             puts("No");
120             continue;
121         }
122         puts("Yes");
123         //反向建图
124         memset(degree,0,sizeof(degree));
125         g.clear();
126         g.resize(2*n+2);
127         for(int u=1;u<=2*n;u++){
128             for(int i=head[u];i!=-1;i=edge[i].next){
129                 int v=edge[i].v;
130                 if(color[u]!=color[v]){
131                     g[color[v]].push_back(color[u]);
132                     degree[color[u]]++;
133                 }
134             }
135         }
136         memset(vis,0,sizeof(vis));
137         TopSort();
138         ans.clear();
139         for(int i=1;i<=n;i++){
140             if(vis[color[i]]==1)ans.push_back(i);
141         }
142         printf("%d",(int)ans.size());
143         sort(ans.begin(),ans.end());
144         for(int i=0;i<(int)ans.size();i++){
145             printf(" %d",ans[i]);
146         }
147         puts("");
148      }
149     return 0;
150 }

View Code