HDU 6149 Valley Numer II（状压DP）

Description

众所周知，度度熊非常喜欢图。

它最近发现了图中也是可以出现 $valley$ —— 山谷的，像下面这张图。

为了形成山谷，首先要将一个图的顶点标记为高点或者低点。标记完成后如果一个顶点三元组 <X,Y,Z> < X , Y , Z > <script type="math/tex" id="MathJax-Element-2"> </script>中， $X$和$Y$之间有边， $Y$与$Z$之间也有边，同时 $X$和$Z$是高点， $Y$是低点，那么它们就构成一个$valley$。

度度熊想知道一个无向图中最多可以构成多少个$valley$，一个顶点最多只能出现在一个$valley$中。

Input

第一行为$T$，表示输入数据组数。

每组数据的第一行包含三个整数$N，M，K$，分别表示顶点个数，边的个数，标记为高点的顶点个数。

接着的M行，每行包含两个两个整数$Xi，Yi$，表示一条无向边。

最后一行包含$K$个整数$Vi$，表示这些点被标记为高点，其他点则都为低点。

● $1\le T\le 20$

● $1\le N\le 30$

● $1\le M\le N*(N-1)/2$

● $0\le K\le min(N,15)$

● $1\le Xi, Yi\le N, Xi\neq Yi$

● $1\le Vi\le N$

Output

对每组数据输出最多能构成的$valley$数目。

Sample Input

3
3 2 2
1 2
1 3
2 3
3 2 2
1 2
1 3
1 2
7 6 5
1 2
1 3
1 4
2 3
2 6
2 7
3 4 5 6 7

Sample Output

1
0
2

Solution

状压$DP$，用$k$位$01$表示$k$个高点的使用情况，$1$表示该高点已经被用过了，$0$表示没有被用过，设前$l$个低点使用的高点状态为$s$可以组成的最多$valley$数目，对于第$l+1$个低点，要么不用这个低点，要么枚举第$l+1$个低点可以配对的高点组$(x,y)$，注意$s$中第$x,y$位需为$0$，进而有转移

$dp[l+1][s+2^x+2^y]=max(dp[l+1][s+2^x+2^y],dp[l][s]+1)$

$max(dp[res][s],0\le s<2^k)$即为答案，其中$res$为低点数量

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
typedef pair<int,int>P;
int T,n,m,k,low[31],high[16],dp[31][(1<<15)+5];
bool a[31][31],vis[31];
vector<P>g[31];
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&m,&k);
        memset(a,0,sizeof(a));
        while(m--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            a[x][y]=a[y][x]=1;
        }
        memset(vis,0,sizeof(vis));
        for(int i=0;i<k;i++)
            scanf("%d",&high[i]),vis[high[i]]=1;
        int res=0;
        for(int i=1;i<=n;i++)
            if(!vis[i])low[res++]=i;
        for(int l=0;l<res;l++)
        {
            g[l].clear();
            for(int i=0;i<k;i++)
                if(a[low[l]][high[i]])
                    for(int j=i+1;j<k;j++)
                        if(a[low[l]][high[j]])
                            g[l].push_back(P(i,j));
        }
        memset(dp,0,sizeof(dp));
        int K=1<<k;
        for(int l=0;l<res;l++)
            for(int s=0;s<K;s++)
            {
                dp[l+1][s]=max(dp[l+1][s],dp[l][s]);
                for(int i=0;i<g[l].size();i++)
                {
                    int x=g[l][i].first,y=g[l][i].second;
                    if((s&(1<<x))||(s&(1<<y)))continue;
                    dp[l+1][s^(1<<x)^(1<<y)]=max(dp[l+1][s^(1<<x)^(1<<y)],dp[l][s]+1);
                } 
            }
        int ans=0;
        for(int s=0;s<K;s++)ans=max(ans,dp[res][s]);
        printf("%d\n",ans);
    }
    return 0;
}