HDU 4115 Eliminate the Conflict

2-SAT，拆成六个点。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<stack>
#include<queue>
#include<algorithm>
using namespace std;

const int maxn=+;

int T,N,M;
stack<int>S;
vector<int>G[maxn];
vector<int>FG[maxn];
int Belong[maxn];
int flag[maxn];
int Block;

int X[maxn];

void init()
{
    for(int i=; i<maxn; i++) G[i].clear();
    for(int i=; i<maxn; i++) FG[i].clear();
    memset(Belong,,sizeof Belong);
    memset(flag,,sizeof flag);
    while(!S.empty()) S.pop();
    Block=;
}

void addEdge(int x,int y)
{
    G[x].push_back(y);
    FG[y].push_back(x);
}

void dfs1(int now)
{
    flag[now]=;
    for(int i=; i<G[now].size(); i++)
        if(!flag[G[now][i]])
            dfs1(G[now][i]);
    S.push(now);
}

void dfs2(int now)
{
    Belong[now]=Block;
    for(int i=; i<FG[now].size(); i++)
        if(!Belong[FG[now][i]])
            dfs2(FG[now][i]);
}

bool judge()
{
    for(int i=; i<*N; i++) if(!flag[i]) dfs1(i);
    while(!S.empty())
    {
        int Top=S.top();
        S.pop();
        if(!Belong[Top])
        {
            Block++;
            dfs2(Top);
        }
    }
    for(int i=; i<*N; i++)
        if(Belong[i]==Belong[i+*N])
            return ;
    return ;
}

void read()
{
    int A,B,C;
    for(int i=;i<N;i++)
    {
        addEdge(*i+,*i++*N);
        addEdge(*i+,*i++*N);

        addEdge(*i+,*i++*N);
        addEdge(*i+,*i++*N);

        addEdge(*i+,*i++*N);
        addEdge(*i+,*i++*N);
    }
    for(int i=;i<N;i++)
    {
        scanf("%d",&X[i]); X[i]--;
        if(X[i]==)
        {
            addEdge(*i++*N,*i+);
            addEdge(*i++*N,*i+);
        }
        else if(X[i]==)
        {
            addEdge(*i++*N,*i+);
            addEdge(*i++*N,*i+);
        }
        else if(X[i]==)
        {
            addEdge(*i++*N,*i+);
            addEdge(*i++*N,*i+);
        }
    }
    for(int i=;i<=M;i++)
    {
        scanf("%d%d%d",&A,&B,&C);
        A--,B--;
        if(C==)
        {
            addEdge(*A+,*B+);
            addEdge(*A+,*B+);
            addEdge(*A+,*B+);

            addEdge(*B+,*A+);
            addEdge(*B+,*A+);
            addEdge(*B+,*A+);
        }

        else if(C==)
        {
            addEdge(*A+,*B++*N);
            addEdge(*A+,*B++*N);
            addEdge(*A+,*B++*N);

            addEdge(*B+,*A++*N);
            addEdge(*B+,*A++*N);
            addEdge(*B+,*A++*N);
        }
    }
}

int main()
{
    scanf("%d",&T);
    for(int Case=;Case<=T;Case++)
    {
        scanf("%d%d",&N,&M);
        init();
        read();
        printf("Case #%d: ",Case);
        if(judge()) printf("yes\n");
        else printf("no\n");
    }
    return ;
}