【UVA11825】Hackers' Crackdown

题意

  Miracle Corporations has a number of system services running in a distributed computer system which is a prime target for hackers. The system is basically a set of N computer nodes with each of them running a set of N services. Note that, the set of services running on every node is same everywhere in the network. A hacker can destroy a service by running a specialized exploit for that service in all the nodes.
  One day, a smart hacker collects necessary exploits for all these N services and launches an attack on the system. He finds a security hole that gives him just enough time to run a single exploit in each computer. These exploits have the characteristic that, its successfully infects the computer where it was originally run and all the neighbor computers of that node.
  Given a network description, find the maximum number of services that the hacker can damage.

题面

  有$n$个节点，每一个节点都有$n$种标记，现在可以对每一个节点分别进行一次操作，每一次操作可以选择一种标记，使这个节点及与其相邻的节点覆盖这种标记（一个节点可以被覆盖多种标记），问经过$n$操作后，最多有多少种标记覆盖了所有节点
  $n≤16$，节点编号从0开始

解法

子集$DP$：
  第一次做子集$DP$的题……
  这个题目的意思就是说：全集为$U$，$|U|=n$，有$n$个集合，为{$P_1，P_2……P_n$}，现在要对这$n$个集合进行分组，要求$P_1∪P_2∪……∪P_k=U$，$P$为该组内的集合，求最多能分多少组
  设$sta_i$为一个二进制数，表示第$i$个节点相邻的点及自身，例如$i=10$_（2），$sta_i=10110$_（2），表示与1号节点相邻的有4，2两个节点。由此可知，$sta_i$表示对$i$进行操作能够标记的点的集合
  设$p_i$表示选择$i$表示的集合能够标记的点集，$i$为一个二进制数，例如$i=10010$_（2）表示选择{0，4}这个集合能够标记的点集

  设$f_i$表示选择$i$集合能够标记的全集个数，即答案，$i$也是一个二进制数。状态已经设出来了，下面就是如何转移：
  对于一种状态$i$，$i$表示的就是一个集合，考虑如何枚举出$i$的所有子集，可以参见刘汝佳的紫书

复杂度

O（$2^n$*平均非空子集个数），当有$n$个物品时，总非空子集数为：$\sum_{i=1}^n C_n^i$

代码

#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#define Lint long long int
using namespace std;
const int N=20;
const int L=(1<<17);
int f[L+5],p[L+5];
int sta[N];
int n,A;
int main()
{
    int len,x,C=0;
    while( scanf("%d",&n) && n )
    {
        memset( sta,0x0,sizeof sta );
        memset( f,0x0,sizeof f ),memset( p,0x0,sizeof p );
        for(int i=0;i<=n-1;i++)
        {
            scanf("%d",&len);
            for(int j=1;j<=len;j++)
            {
                scanf("%d",&x);
                sta[i]|=(1<<x);
            }
            sta[i]|=(1<<i);
        }
        A=(1<<n)-1;
        for(int i=0;i<=A;i++)
        {
            p[i]=0;
            for(int j=0;j<=n-1;j++)
                if( i&(1<<j) )
                    p[i]|=sta[j];
        }
        for(int i=0;i<=A;i++)
        {
            if( p[i]!=A )   f[i]=0;
            else   f[i]=1;
            for(int j=i; j ;j=(j-1)&i)
                f[i]=max( f[i],f[j]+f[i^j] );
        }
        printf("Case %d: %d\n",++C,f[A]);
    }
    return 0;
}