HDU 6149 Valley Numer II (背包+状压+滚动数组)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6149

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
const int  maxn =30+5;
const int mod=1e9+7;
const int ub=10000;
void gmax(int& x,int y) {if(y>x) x=y;}
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}

/*
题目大意：给定一个图，
事先标出一些高点，然后稳有多少个三元组满足，
一低二高，低和高之间有边相连。
每个点只能属于这一个三元组。

对每个低点进行背包，选择的集合是当前低点的边集合，
因为k只有15所以对k进行状压（再加一层下标映射），
背包里面一层是普通的状压，不过要考虑当前更新过的不能重复使用，
所以牺牲点空间进行滚动数组。
（刚开始还不以为需要滚动数组，注意对于一个单一点有多种二元组选择。）

*/

int n,m,k,x,y;
int dp[2][1<<16],vis[40];
int a[20],mp[100];
vector<int> g[maxn];

int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d",&n,&m,&k);
        for(int i=0;i<n;i++) g[i].clear();
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&x,&y);x--,y--;
            g[x].push_back(y),g[y].push_back(x);///
        }

        mst(vis,0);mst(mp,0);
        int s=0;///总状态
        for(int i=0;i<k;i++)
        {
            scanf("%d",&a[i]);
            a[i]--;vis[a[i]]=1;
            mp[a[i]]=i;///映射
        }
        s=(1<<k)-1;

        int now=0;
        mst(dp,-1);dp[now][s]=0;///当前状态为零
        for(int i=0;i<n;i++) if(vis[i]==0){///挑选那些低点
            now^=1;memcpy(dp[now],dp[now^1],sizeof(dp[now]));
            for(int j=0;j<g[i].size();j++)
            {
                if(vis[g[i][j]]==0) continue;
                for(int q=j+1;q<g[i].size();q++)///选择点对
                {
                    if(vis[g[i][q]]==0) continue;
                    int tmp=(1<<mp[g[i][j]])+(1<<mp[g[i][q]]);
                    for(int p=0;p<(1<<k);p++)
                        if((p&tmp)==tmp && dp[now^1][p]!=-1)
                            dp[now][p^tmp]=max(dp[now][p^tmp],dp[now^1][p]+1);
                }
            }
        }

        int ans=0;for(int i=0;i<(1<<k);i++) gmax(ans,dp[now][i]);
        printf("%d\n",ans);
    }
    return 0;
}