uva 11600(期望dp)

题意：有n个城市，编号为1～n，城市两两直接都有一条双向道路（一共n(n-1)/2条），有m条路是安全的，其他的路上都有妖怪。为了让城市两两可达，有一个捉妖者第一天晚上住在城市1，然后每个白天随机选择一个新的城市，顺着能到达新城市的道路上把妖怪收服，晚上住在新城市。问需要多少个白天才能让任意两个城市之间存在一条路径上没有妖怪而相互可达。
题解：让任意两两城市可达，那就是让图成为连通图的问题。城市最多有30个，那么就可以从状态压缩来考虑，那么刚开始最多有30个连通分量，那么状态最大2^32-1，数组肯定放不下，考虑用map映射，每种状态对应一个期望。所以要先用dfs统计初始连通分量，把每个连通分量的城市的数量存起来，然后dp中通过把和当前状态s中未连通的连通分量进行连通计算期望值。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <map>
#include <vector>
using namespace std;
const int N = 35;
int n, m, vis[N], cnt, num[N];
map<int, double> f;
vector<int> g[N];

void init() {
    memset(vis, 0, sizeof(vis));
    memset(num, 0, sizeof(num));
    f.clear();
    for (int i = 0; i <= n; i++)
        g[i].clear();
}

void input() {
    int u, v;
    for (int i = 0; i < m; i++) {
        scanf("%d%d", &u, &v);
        g[u].push_back(v);      
        g[v].push_back(u);
    }
}

int dfs(int u) {
    int count = 1;
    vis[u] = 1;
    for (int i = 0; i < g[u].size(); i++) {
        int v = g[u][i];
        if (!vis[v])
            count += dfs(v);
    }
    return count;
}

double dp(int s) {
    if (f[s] > 1e-9)
        return f[s];
    int con = 0;
    for (int i = 0; i < cnt; i++)
        if (s & (1 << i))
            con += num[i];
    if (con == n)
        return f[s] = 0;
    f[s] = (n - 1) * 1.0 / (n - con);
    for (int i = 0; i < cnt; i++)
        if (!(s & (1 << i)))
            f[s] += dp(s | (1 << i)) * num[i] * 1.0 / (n - con);
    return f[s];
}

void solve() {
    cnt = 0;
    for (int i = 1; i <= n; i++)
        if (!vis[i])
            num[cnt++] = dfs(i);
}

int main() {
    int t, cas = 1;
    scanf("%d", &t);
    while (t--) {
        scanf("%d%d", &n, &m);
        init();
        input();
        solve();
        printf("Case %d: %lf\n", cas++, dp(1));
    }
    return 0;
}