DLX

// hdu 2828
// poj 3740
int mx[N][M];
int U[nm],D[nm],L[nm],R[nm],X[nm],Y[nm];
int S[M],H[N];
int sz,ans;
void ver(int x,int y){///上 下
    D[x]=y;U[y]=x;
}
void hor(int x,int y){////左 右
    R[x]=y;L[y]=x;
}

namespace AA{/////////可重复覆盖
    int vis[M];
    void Remove(int c){
        for(int i = D[c]; i != c; i = D[i])
            hor(L[i],R[i]);
    }

    void Resume(int c){
        for(int i = D[c]; i != c; i = D[i])
            hor(L[i],i);hor(i,R[i]);
    }

    int h(){
        memset(vis,0,sizeof(vis));
        int ans = 0;
        for(int c = R[0] ; c ; c = R[c]){
            if(!vis[c]){
                ans++;
                for(int i = D[c]; i != c; i = D[i])
                    for(int j = R[i]; j != i; j = R[j])
                        vis[Y[j]] = 1;
            }
        }
        return ans;
    }

    bool Dance(int dep){
        int c=R[0];
        if(c == 0)
            return 1;
        if(dep + h() > ans)       //A_Star剪枝
            return 0;


        for(int i=R[0]; i; i=R[i])
            if(S[i]<S[c]) c=i;

        for(i=D[c]; i!=c; i=D[i]){
            Remove(i);
            for(j=R[i]; j!=i; j=R[j])Remove(j);
            if(Dance(dep+1))return 1;
            for(j=L[i]; j!=i; j=L[j])Resume(j);//在这里我遵循的原则是先删除的后还原，后删除的先还原。
            Resume(i);
        }

        return 0;
    }
}

namespace BB{/////////不可重复覆盖
    void Remove(int c){
        hor(L[c],R[c]);
        for(int i=D[c];i!=c;i=D[i]){
            for(int j=R[i];j!=i;j=R[j]){
                ver(U[j],D[j]);
                S[Y[j]]--;
            }
        }
    }
    void Resume(int c){
        hor(L[c],c);
        hor(c,R[c]);
        for(int i=U[c];i!=c;i=U[i]){
            for(int j=L[i];j!=i;j=L[j]){
                ver(U[j],j);
                ver(j,D[j]);
                S[Y[j]]++;
            }
        }
    }
    int Dance(int k){
        int c=R[0];
        if(c==0){
            return 1;
        }
        for(int i=R[0];i!=0;i=R[i]){
            if(S[i]<S[c])c=i;
        }
        Remove(c);
        for(int i=D[c];i!=c;i=D[i]){
            for(int j=R[i];j!=i;j=R[j])Remove(Y[j]);
            if(Dance(k+1))return 1;
            for(int j=L[i];j!=i;j=L[j])Resume(Y[j]);
        }
        Resume(c);
        return 0;
    }
}

void Link(int r,int c){
    ver(sz,D[c]);
    ver(c,sz);
    if(H[r]<0)H[r]=L[sz]=R[sz]=sz;
    else {
        hor(sz,R[H[r]]);
        hor(H[r],sz);
    }
    S[c]++;
    Y[sz]=c;
    X[sz++]=r;
}

void Build(int n,int m){
    for(int i=0;i<=m;i++){
        hor(i,i+1);
        D[i]=U[i]=i;
        S[i]=0;
    }
    hor(m,0);
    sz=m+1;
    for(int i=1;i<=n;i++)H[i]=-1;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(mx[i][j])Link(i,j);
        }
    }
}

void solve(int n,int m){
    Build(n,m);
    if(AA::Dance(0))printf("yes\n");
    if(BB::Dance(0))printf("yes\n");
}
int main(){
    return 0;
}