hust 1017

题意：求01矩阵的精确覆盖。

分析：本来想学习dancing links来解决数独问题，发现dancing links最初解决的问题是精确覆盖，于是就找到这道题来做了。这种NPC问题只能用DFS暴搜的情况下，很适合的一种优化方式就是用dancing links加速状态的改变，利用双向循环十字链表使元素的删除与恢复操作非常简便快捷。

#include <cstdio>

int U[102005],D[102005],L[102005],R[102005];
int X[100005],Y[100005],H[1005],S[1005],ans[1005],len,M,N,sz;

void init(int n,int m)
{
    for(int i = 0;i <= m;i++)
    {
        U[i] = D[i] = i;
        L[i + 1] = i;
        R[i] = i + 1;
        S[i] = 0;
    }
    R[m] = 0;
    L[0] = m;
    sz = m + 1;
}

void remove(int c)
{
    //删除一整列
    R[L[c]] = R[c];
    L[R[c]] = L[c];
    //删除行
    for(int i = D[c];i != c;i = D[i])
    {
        for(int j = R[i];j != i;j = R[j])
        {
            D[U[j]] = D[j];
            U[D[j]] = U[j];
            S[X[j]]--;
        }
    }
}

void resume(int c)
{
    //恢复一整列
    L[R[c]] = c;
    R[L[c]] = c;
    //恢复行
    for(int i = D[c];i != c;i = D[i])
    {
        for(int j = R[i];j != i;j = R[j])
        {
            D[U[j]] = j;
            U[D[j]] = j;
            S[X[j]]++;
        }
    }
}

void ins(int r,int c)
{
    S[c]++;
    //纵向插入
    D[U[c]] = sz;
    U[sz] = U[c];
    D[sz] = c;
    U[c] = sz;
    X[sz] = c;
    Y[sz] = r;
    //横向插入
    if(H[r] == -1)
    {
        H[r] = L[sz] = R[sz] = sz;
    }
    else
    {
        R[L[H[r]]] = sz;
        L[sz] = L[H[r]];
        R[sz] = H[r];
        L[H[r]] = sz;
    }
    sz++;
}

bool dfs(int k)
{
    if(R[0] == 0)
    {
        len = k;
        ans[k] = -1;
        return true;
    }
    else
    {
        int m = 0xffffff,num;
        for(int i = R[0];i != 0;i = R[i])
        {
            if(S[i] == 0)
            {
                return false;
            }
            if(m > S[i])
            {
                m = S[i];
                num = i;
                if(m == 1)
                {
                    break;
                }
            }
        }
        remove(num);
        for(int i = D[num];i != num;i = D[i])
        {
            ans[k] = Y[i];
            for(int j = R[i];j != i;j = R[j])
            {
                remove(X[j]);
            }
            if(dfs(k + 1))
            {
                return true;
            }
            for(int j = R[i];j != i;j = R[j])
            {
                resume(X[j]);
            }
        }
        resume(num);
    }
    return false;
}

int main()
{
    while(~scanf("%d%d",&N,&M))
    {
        init(N,M);
        for(int i = 1;i <= N;i++)
        {
            int k;
            scanf("%d",&k);
            H[i] = -1;
            for(int j = 1;j <= k;j++)
            {
                int c;
                scanf("%d",&c);
                ins(i,c);
            }
        }
        if(dfs(0))
        {
            printf("%d",len);
            for(int i = 0;ans[i] != -1;i++)
            {
                printf(" %d",ans[i]);
            }
            printf("\n");
        }
        else
        {
            printf("NO\n");
        }
    }
    return 0;
}