hdu 4421 BitMagic

这是一道区域赛的题目，解法有许多，这边是2-sat的做法

题目大意:自己看题

分析：对于A[i]的每一位做2-SAT，判断是否可行。

主要是建图：

对于a&b=0  有 a-> ┐b, b-> ┐a

a&b=1            ┐a->a , ┐b->b

a|b=0            a-> ┐a,b-> ┐b

a|b=1　　  　   ┐a->b, ┐b->a

a^b=0　　　　a->b,b->a, ┐a-> ┐b, ┐b-> ┐a

a^b=1　　　　a-> ┐b,b-> ┐a, ┐a->b, ┐b->a

用Kosaraju算法会T(也许我写渣了)用Tarjan算法就ok╮(╯▽╰)╭

#include <iostream>
#include <cstdio>
#include <cstring>
#include <set>
#include <algorithm>
#include <map>
#include<vector>
using namespace std;
const int MAXN = 2010;
const int MAXM = 1010*501*2;
struct Edge
{
    int to,next;
} edge[MAXM];
int head[MAXN],tot;
void addedge(int u,int v)
{
    edge[tot].to = v;
    edge[tot].next = head[u];
    head[u] = tot++;
}
int Low[MAXN],DFN[MAXN],Stack[MAXN],Belong[MAXN];//Belong数组的值1~scc
int Index,top;
int scc;
bool Instack[MAXN];
int num[MAXN];
void init()
{
    tot = 0;
    memset(head,-1,sizeof(head));
}
void Tarjan(int u)
{
    int v;
    Low[u] = DFN[u] = ++Index;
    Stack[top++] = u;
    Instack[u] = true;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        v = edge[i].to;
        if( !DFN[v] )
        {
            Tarjan(v);
            if(Low[u] > Low[v])Low[u] = Low[v];
        }
        else if(Instack[v] && Low[u] > DFN[v])
            Low[u] = DFN[v];
    }
    if(Low[u] == DFN[u])
    {
        scc++;
        do
        {
            v = Stack[--top];
            Instack[v] = false;
            Belong[v] = scc;
            num[scc]++;
        }
        while(v != u);
    }
}
int maps[1010][1010];
bool solve(int m){
    for(int i=0;i<m;++i){
        for(int j=0;j<m;++j){
            if(i==j&&maps[i][j])return false;
            if(maps[i][j]!=maps[j][i])return false;
        }
    }
    for(int k=0;k<31;++k){
        init();
        for(int i=0;i<m;++i){
            for(int j=0;j<m;++j){
                if(i==j)continue;
                int is = maps[i][j]&(1<<k);
                if((i&1)&&(j&1)){//|
                    if(is){//a|b=1
                        addedge(i,j+m);
                        addedge(j,i+m);
                    }
                    else{//a|b=0
                        addedge(i+m,i);
                        addedge(j+m,j);
                    }
                }
                else if((i%2==0)&&(j%2==0)){//&
                    if(is){//a&b=1
                        addedge(i,i+m);
                        addedge(j,j+m);
                    }
                    else{//a&b=0
                        addedge(i+m,j);
                        addedge(j+m,i);
                    }
                }
                else {//^
                    if(is){//a^b=1
                        addedge(i+m,j);
                        addedge(j+m,i);
                        addedge(i,j+m);
                        addedge(j,i+m);
                    }
                    else{//a^b=0
                        addedge(i+m,j+m);
                        addedge(j+m,i+m);
                        addedge(i,j);
                        addedge(j,i);

                    }
                }
            }
        }
        memset(DFN,0,sizeof(DFN));
        for(int i=0;i<2*m;++i){
            if(!DFN[i])Tarjan(i);
        }
        for(int i=0;i<m;++i){
            if(Belong[i]==Belong[m+i]){
                return false;
            }
        }
    }
    return true;
}
int main (){
    int m;
    while(scanf("%d",&m)!=EOF){
        for(int i=0;i<m;++i){
            for(int j=0;j<m;++j){
                scanf("%d",&maps[i][j]);
            }
        }
        if(solve(m))printf("YES
");
        else printf("NO
");
    }
}