图论之2-sat模板

#include<cstdio>
#include<cstring>
#include<iostream>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
int n;
const int maxn=2011;
struct Node{
int s,e;
}Tim[2*maxn];

vector<int>e[2*maxn];
vector<int>e2[2*maxn];
vector<int>ans;
stack<int>st;

int id,col;
int dfn[2*maxn],low[2*maxn],flag[2*maxn], in[2*maxn];
bool inst[2*maxn];
int cf[2*maxn];//cf[x]保存的是和编号为 x 的连通分量矛盾的连通分量的编号
int Col[2*maxn];
//利用矛盾关系建图、缩点、反向建图、利用拓扑关系重新着色。
void addedge(int u,int v){
e[u].push_back(v);
}
bool Judge(int x,int y){
if(Tim[x].e<=Tim[y].s||Tim[x].s>=Tim[y].e)return 1;
return 0;
}
void buildmap(){
for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
if(i!=j){
if(!Judge(2*i,2*j))
addedge(2*i,2*j+1);
if(!Judge(2*i,2*j+1))
addedge(2*i,2*j);
if(!Judge(2*i+1,2*j))
addedge(2*i+1,2*j+1);
if(!Judge(2*i+1,2*j+1))
addedge(2*i+1,2*j);
}
}
}

}
void Tarjan(int u){
dfn[u]=low[u]=id++;
st.push(u);
inst[u]=1;
for(int i=0;i<e[u].size();i++){
int k=e[u][i];
if(dfn[k]==0){
Tarjan(k);
low[u]=std::min(low[u],low[k]);
}
else if(inst[k]){
low[u]=std::min(low[u],dfn[k]);
}
}
if(dfn[u]==low[u]){
while(1){
int x=st.top();
st.pop();
flag[x]=col;
inst[x]=0;
if(x==u)break;
}
col++;
}
}
int solve(){
for(int i=0;i<2*n;i++){
if(!dfn[i]){
Tarjan(i);
}
}
for(int i=0;i<2*n;i+=2){
if(flag[i]==flag[i^1])return 1;
cf[flag[i]]=flag[i^1];
cf[flag[i^1]]=flag[i];
}
return 0;
}
void Topsort(){
queue<int>q;
for(int i=1;i<col;i++){
if(in[i]==0)q.push(i);
}
while(!q.empty()){
int u=q.front();
q.pop();
if(Col[u]==0){
Col[u]=1;
Col[cf[u]]=-1;
}
for(int i=0;i<e2[u].size();i++){
int v=e2[u][i];
in[v]--;
if(in[v]==0)q.push(v);
}
}
}
void Getans(){
for(int i=0;i<2*n;i++){//建立反向图
for(int j=0;j<e[i].size();j++){
if(flag[i]!=flag[e[i][j]]){
e2[flag[e[i][j]]].push_back(flag[i]);
in[flag[i]]++;
}
}
}
Topsort();//求结果
}
void init(){
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
memset(flag,0,sizeof(flag));
memset(inst,0,sizeof(inst));
}
int main(){
init();
buildmap();
solve();
return 0;
}