hdu4421 Bit Magic 2-sat

點擊打開鏈接～～

主要就是下圖中，給出b[n][n],求是否有a[n]符合，有則“YES”，否則“NO”

至於2-sat，就是處理兩個關係以下的算法，比如a&b==0,那麼如果a==1，那麼 b一定是0，如果b==1，那麼a一定等於0，具體看代碼

可以將i看作i==0時，i+n看作i==1時

#include 
#include 
#include 
#include 
using namespace std;
const int N=400010;
const int M=400010;
int b[505][505];
struct twosat
{
    int num[N],ad[M],v[M],pre[M];   
    int dfn[N],low[N],s[N];
    int cnt,top,tol,n,scnt;
    void init()
    {
        tol=0;
        memset(ad,0,sizeof ad);
    }
    void add(int x,int y)   //建圖
    {
        pre[++tol]=ad[x];
        ad[x]=tol;
        v[tol]=y;
    }
/*    void add_c(int x,int xi,int y,int yi)
    {
        add(x+xi*n,y+(!yi)*n);
        add(y+yi*n,x+(!xi)*n);
    }*/
    void dfs(int u)
    {
        dfn[u]=low[u]=cnt++; //排序     
        s[++top]=u;   //入棧
        for(int i=ad[u];i;i=pre[i])
        {
            if(!dfn[v[i]])   //如果沒有訪問過，則再繼續找
            {
                dfs(v[i]);
                low[u]=min(low[u],low[v[i]]); //最小公共祖先
            }
            else if(!num[v[i]])
                low[u]=min(low[u],dfn[v[i]]);  //有可能low[v[i]]是已經被刪除 ？
        }
        if(low[u]==dfn[u])
        {
            ++scnt;
            while(1)
            {
                int x=s[top--];
                num[x]=scnt;
                if(x==u) break;
            }
        }
    }
    bool solve()
    {
        memset(num,0,sizeof num);
        memset(dfn,0,sizeof dfn);
        memset(low,0,sizeof low);
        cnt=top=scnt=0;
        for(int i=0;i<2*n;i++)
        {
            if(!dfn[i])
                dfs(i);
        }
        for(int i=0;i 選了i，就一定要選j
                        sat.add(j,i+n);
                    }
                    else
                    {
                        sat.add(i+n,i);
                        sat.add(j+n,j);
                    }
                }
                else if(i%2==0 && j%2==0) //&
                {
                    if(b[i][j]&(1<