HDU 3062 Party

Party

Problem Description

有n对夫妻被邀请参加一个聚会，因为场地的问题，每对夫妻中只有1人可以列席。在2n 个人中，某些人之间有着很大的矛盾（当然夫妻之间是没有矛盾的），有矛盾的2个人是不会同时出现在聚会上的。有没有可能会有n 个人同时列席？

 

Input

n： 表示有n对夫妻被邀请 (n<= 1000)
m： 表示有m 对矛盾关系 ( m < (n - 1) * (n -1))

在接下来的m行中，每行会有4个数字，分别是 A1,A2,C1,C2 
A1,A2分别表示是夫妻的编号 
C1,C2 表示是妻子还是丈夫 ，0表示妻子 ，1是丈夫
夫妻编号从 0 到 n -1 

 

Output

如果存在一种情况 则输出YES 
否则输出 NO 

 

Sample Input

  
  

   
   2 
1
0 1 1 1 
  
  

 

Sample Output

  
  

   
   YES
  
  

 

Source

2009 Multi-University Training Contest 16 - Host by NIT

 

分析：简单的2-SAT。很容易想到2-SAT模型。用ai表示第i对的妻子，用ai+n表示第i对的丈夫。然后就是根据输入信息建图。由于不需要输出解，可以用直接的DFS算法或者Tarjan判断一个强连通的合法性就可以了。下面给出2-SAT的两种算法。

DFS:

#include<bits/stdc++.h>
using namespace std;

int n,m;
vector<int> vec[2005];


void input(){
    int a1,a2,c1,c2;
    for(int i=0;i<=2*n;i++){
        vec[i].clear();
    }
    for(int i=1;i<=m;i++){
        scanf("%d%d%d%d",&a1,&a2,&c1,&c2);
        if(c1==0&&c2==0){//a1和a2有矛盾
            vec[a1].push_back(a2+n);
            vec[a2].push_back(a1+n);
        }
        else if(c1==1&&c2==0){//a1+n和a2有矛盾
            vec[a1+n].push_back(a2+n);
            vec[a2].push_back(a1);
        }
        else if(c1==0&&c2==1){//a1和a2+n有矛盾
            vec[a1].push_back(a2);
            vec[a2+n].push_back(a1+n);
        }
        else if(c1==1&&c2==1){//a1+n和a2+n有矛盾
            vec[a1+n].push_back(a2);
            vec[a2+n].push_back(a1);
        }
    }
}

int id[2005];
int top;
int Stack[2005];

bool dfs(int u){
    if(u>=n){
        if(id[u-n]){
            return false;
        }
    }
    else{
        if(id[u+n]){
            return false;
        }
    }
    if(id[u]){
        return true;
    }
    id[u]=1;
    Stack[top++]=u;
    for(int i=0;i<vec[u].size();i++){
        if(!dfs(vec[u][i])){
            return false;
        }
    }
    return true;
}

bool solve(){
    memset(id,0,sizeof(id));
    for(int i=0;i<n;i++){
        if(id[i]==0&&id[i+n]==0){
            top=0;
            if(!dfs(i)){
                while(top){
                    id[Stack[--top]]=0;
                }
                if(!dfs(i+n)){
                    return false;
                }
            }
        }
    }
    return true;
}

int main(){
    while(~scanf("%d%d",&n,&m)){
        input();
        if(solve()){
            printf("YES\n");
        }
        else{
            printf("NO\n");
        }
    }
    return 0;
}

Tarjan+缩点+拓扑+染色：

#include<bits/stdc++.h>
using namespace std;

int n,m;
vector<int> vec[2005];
int top;
int Stack[2005];

int dfn[2005];
int low[2005];
int clk;
int scc;
bool instack[2005];
int scc_num[2005];
int func[2005];
vector<int> Vec[2005];
int indegree[2005];
int color[2005];

inline void init(){
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(instack,false,sizeof(instack));
    memset(scc_num,0,sizeof(scc_num));
    memset(func,0,sizeof(func));
    memset(indegree,0,sizeof(indegree));
    memset(color,0,sizeof(color));
    for(int i=0;i<=2*n;i++){
        Vec[i].clear();
        vec[i].clear();
    }
}

inline void input(){
    int a1,a2,c1,c2;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d%d",&a1,&a2,&c1,&c2);
        if(c1==0&&c2==0){//a1和a2有矛盾
            vec[a1].push_back(a2+n);
            vec[a2].push_back(a1+n);
        }
        else if(c1==1&&c2==0){//a1+n和a2有矛盾
            vec[a1+n].push_back(a2+n);
            vec[a2].push_back(a1);
        }
        else if(c1==0&&c2==1){//a1和a2+n有矛盾
            vec[a1].push_back(a2);
            vec[a2+n].push_back(a1+n);
        }
        else if(c1==1&&c2==1){//a1+n和a2+n有矛盾
            vec[a1+n].push_back(a2);
            vec[a2+n].push_back(a1);
        }
    }
}

inline void tarjan(int u,int pre){
    dfn[u]=low[u]=clk++;
    Stack[top++]=u;
    instack[u]=true;
    int len=vec[u].size();
    for(int i=0;i<len;i++){
        int v=vec[u][i];
        //if(v!=pre){
            if(dfn[v]==0){
                tarjan(v,u);
                low[u]=min(low[v],low[u]);
            }
            else if(instack[v]){
                low[u]=min(low[u],dfn[v]);
            }
        //}
    }
    if(dfn[u]==low[u]){
        scc++;
        do{
            int v=Stack[--top];
            instack[v]=false;
            scc_num[v]=scc;
        }while(Stack[top]!=u);
    }
}

inline void suodian(){
    for(int i=0;i<2*n;i++){
        int len=vec[i].size();
        for(int j=0;j<len;j++){
            int k=vec[i][j];
            int u=scc_num[i];
            int v=scc_num[k];
            if(u!=v){
                Vec[v].push_back(u);//反向建图
                indegree[u]++;
            }
        }
    }
}

inline void toposort(){
    queue<int> q;
    for(int i=1;i<=scc;i++){
        if(indegree[i]==0){
            q.push(i);
        }
    }
    while(q.empty()==0){
        int u=q.front();
        q.pop();
        int len=vec[u].size();
        if(color[u]==0){
            color[u]=1;
            color[func[u]]=2;
        }
        for(int i=0;i<len;i++){
            int v=vec[u][i];
            indegree[v]--;
            if(indegree[v]==0){
                q.push(v);
            }
        }
    }
}

inline bool solve(){
    clk=0;
    scc=0;
    for(int i=0;i<2*n;i++){
        if(dfn[i]==0){
            tarjan(i,-1);
        }
    }
    for(int i=0;i<n;i++){
        if(scc_num[i]==scc_num[i+n]){
            return false;
        }
        else{
            func[scc_num[i]]=scc_num[i+n];
            func[scc_num[i+n]]=scc_num[i];
        }
    }
    suodian();
    toposort();
    return true;
}

int main(){
    while(~scanf("%d%d",&n,&m)){
        init();
        input();
        if(solve()){
            printf("YES\n");
        }
        else{
            printf("NO\n");
        }
    }
    return 0;
}