hdu4778 状态压缩

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <vector>
#include <string.h>
using namespace std;
bool use[1<<22];
int dp[1<<23];
int G,B,S;
int P[22][10];
int nu[10];
int work3()
{
    int ans=0;
    for(int i=1; i<=G; i++)
    {
        int ge=0;
        for(int j=0;j <B; j++)
            ge+=P[j][i];
         ans+=ge/S;
    }
    return ans;
}
int jia(int j)
{
    int ans=0;
    for(int i=1; i<=G; i++)
        {
            nu[i]+=P[j][i];
            ans+=nu[i]/S;
           nu[i]%=S;
        }
    return ans;
}
void jian(int j)
{
    for(int i=1; i<=G; i++)
    {
        nu[i]=((nu[i]-P[j][i])%S+S)%S;
    }
}
void dfs(int D,int cntge)
{
    if(use[D])return ;
    use[D]=true;
    dp[D]=0;
         for(int j=0;j<B; j++)
         {
             if((D&(1<<j))==0)continue;
            int num=jia(j);
             dfs(D^(1<<j),cntge-num);
             if(num>0)
              dp[D]=max(dp[D^(1<<j)]+num,dp[D]);
             else{
               dp[D]=max(cntge-dp[D^(1<<j)],dp[D]);
             }
             jian(j);
         }
}
int main()
{

    while(scanf("%d%d%d",&G,&B,&S)&&G+B+S!=0)
    {
        if(B==0){
            printf("0\n");continue;
        }
        for(int i=0; i<B; i++)
        {
            int d;
            scanf("%d",&d);
            memset(P[i],0,sizeof(P[i]));
            for(int j=0;j<d; j++){
                int k;scanf("%d",&k);
                 P[i][k]++;
            }
        }
        int K=1<<B;
        memset(use,false,sizeof(use));
        use[0]=true;
        int ge1=work3();
        dfs(K-1,ge1);
        int ge=dp[K-1];
        printf("%d\n",ge*2-ge1);
    }

    return 0;
}