bzoj3836

状压dp

图上怎么跑dp？我们跑三进制状压dp，0表示选了，1表示既没选也没覆盖，2表示没选但是被覆盖了。

状态是dp[dep][S]表示当前走到了深度为dep的节点，状态为S,按照dfs序转移

每次转移就是计算这个点选了没选，然后像树形dp一样更新节点

返祖边也要处理

#include<bits/stdc++.h>
using namespace std;
const int N = 5e4 + ;
int n, m, ans;
vector<int> G[N];
int dp[][N], c[N], vis[N], bin[], st[N], d[N];
inline int rd()
{
    int x = , f = ; char c = getchar();
    while(c < '' || c > '') { if(c == '-') f = -; c = getchar(); }
    while(c >= '' && c <= '') { x = x *  + c - ''; c = getchar(); }
    return x * f;
}
int bit(int S, int t)
{
    return S / bin[t] % ;
}
void dfs(int u, int dep)
{
    vis[u] = ;
    d[u] = dep;
    if(!dep)
    {
        dp[][] = c[u];
        dp[][] = ;
        dp[][] = 1e9;
    }
    else
    {
        int top = ;
        for(int i = ; i < G[u].size(); ++i)
        {
            int v = G[u][i];
            if(d[v] < d[u] && vis[v]) st[++top] = d[v];
        }
        for(int i = ; i < bin[dep + ]; ++i) dp[dep][i] = 1e9;
        for(int i = ; i < bin[dep]; ++i)
        {
            int U = , V = i;
            for(int j = ; j <= top; ++j) if(bit(i, st[j]) == ) U = ; else if(bit(i, st[j]) == ) V += bin[st[j]];
            dp[dep][i + U * bin[dep]] = min(dp[dep][i + U * bin[dep]], dp[dep - ][i]);
            dp[dep][V] = min(dp[dep][V], dp[dep - ][i] + c[u]);
        }
    }
    for(int i = ; i < G[u].size(); ++i)
    {
        int v = G[u][i];
        if(vis[v]) continue;
        dfs(v, dep + );
        for(int j = ; j < bin[dep + ]; ++j) dp[dep][j] = min(dp[dep + ][j], dp[dep + ][j +  * bin[dep + ]]);
    }
}
int main()
{
    n = rd();
    m = rd();
    bin[] = ;
    for(int i = ; i <= ; ++i) bin[i] = bin[i - ] * ;
    for(int i = ; i <= n; ++i) c[i] = rd();
    for(int i = ; i <= m; ++i)
    {
        int u = rd(), v = rd();
        G[u].push_back(v);
        G[v].push_back(u);
    }
    for(int i = ; i <= n; ++i) if(!vis[i])
    {
        dfs(i, );
        ans += min(dp[][], dp[][]);
    }
    printf("%d\n", ans);
    return ;
}