poj 3740 Easy Finding(Dancing Links)

3 3
0 1 0
0 0 1
1 0 0
4 4
0 0 0 1
1 0 0 0
1 1 0 1
0 1 0 0

Yes, I found it
It is impossible

题意:

给你一个n*m(n<=16,m<=300)的矩阵。矩阵的每一个元素仅仅能是0或1.如今问你能不能从里面选一些列出来使的没一列有且仅有一个1.

思路:

開始逗比了。把题看成每行仅仅有16个。一看16就乐了。这是典型的壮压嘛。

把每行压成长度为16的二进制数。

然后就是可行性判定了。写完一交re了。

这还用说。

m,n都看错了不re才怪。因为每行值个数太多。不能壮压了。无赖。仅仅能百度。得知是Dancing Links.于是百度其资料。

链接一枚。

具体见代码：

#include<iostream>
#include<stdio.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=6010;
int U[maxn],D[maxn],L[maxn],R[maxn],C[maxn];//上下左右指针。
c[i]结点i相应的列。
int S[310],H[310];//S[i]为i列1的个数。H[i]为i行的尾指针。

int n,m,cnt;
void init()
{
    int i;
    for(i=1;i<=n;i++)
        H[i]=-1;
    for(i=0;i<=m;i++)
    {
         S[i]=0;
         L[i+1]=i;
         R[i]=i+1;
         U[i]=D[i]=i;
    }
    R[m]=0;
    cnt=m+1;
}
void Insert(int r,int c)
{                        //头插法建链表
    U[cnt]=c,D[cnt]=D[c];//确定新增结点上下指针信息
    U[D[c]]=cnt,D[c]=cnt;//恢复链表信息
    if(H[r]==-1)         //确定左右指针信息
        H[r]=L[cnt]=R[cnt]=cnt;//增加头
    else
    {
        L[cnt]=H[r],R[cnt]=R[H[r]];//头插法
        L[R[H[r]]]=cnt,R[H[r]]=cnt;
    }
    S[c]++;//更新附加信息
    C[cnt++]=c;
}
void Remove(int c)//移除c列。

{
    int i,j;
    R[L[c]]=R[c],L[R[c]]=L[c];
    for(i=D[c];i!=c;i=D[i])
        for(j=R[i];j!=i;j=R[j])
            D[U[j]]=D[j],U[D[j]]=U[j],S[C[j]]--;
}
void Resume(int c)//还原c列。

{
    int i,j;
    R[L[c]]=L[R[c]]=c;
    for(i=D[c];i!=c;i=D[i])
        for(j=R[i];j!=i;j=R[j])
            D[U[j]]=U[D[j]]=j,S[C[j]]++;
}
bool dfs()
{
    if(R[0]==0)
        return true;
    int i,j,c,miv=INF;
    for(i=R[0];i;i=R[i])
        if(S[i]<miv)
            miv=S[i],c=i;
    Remove(c);//处理第c列
    for(i=D[c];i!=c;i=D[i])
    {
        for(j=R[i];j!=i;j=R[j])
            Remove(C[j]);
        if(dfs())
            return true;
        for(j=L[i];j!=i;j=L[j])
            Resume(C[j]);
    }
    Resume(c);
    return false;
}
int main()
{
    int op,i,j;
    while(~scanf("%d%d",&n,&m))
    {
        init();
        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++)
            {
                scanf("%d",&op);
                if(op)
                    Insert(i,j);
            }
        if(dfs())
            printf("Yes, I found it\n");
        else
            printf("It is impossible\n");
    }
    return 0;
}