hust 1017

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

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

 #include <cstdio>

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

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

 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] == -)
     {
         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[] == )
     {
         len = k;
         ans[k] = -;
         return true;
     }
     else
     {
         int m = 0xffffff,num;
         for(int i = R[];i != ;i = R[i])
         {
             if(S[i] == )
             {
                 return false;
             }
             if(m > S[i])
             {
                 m = S[i];
                 num = i;
                 if(m == )
                 {
                     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 + ))
             {
                 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 = ;i <= N;i++)
         {
             int k;
             scanf("%d",&k);
             H[i] = -;
             for(int j = ;j <= k;j++)
             {
                 int c;
                 scanf("%d",&c);
                 ins(i,c);
             }
         }
         if(dfs())
         {
             printf("%d",len);
             for(int i = ;ans[i] != -;i++)
             {
                 printf(" %d",ans[i]);
             }
             printf("\n");
         }
         else
         {
             printf("NO\n");
         }
     }
     return ;
 }