LightOJ 1287 Where to Run (期望+状压dp+记忆化搜索)

题意：给出无环图，小偷一开始在起点0，走过的点不能再走，当还有路可以走时，可以选择停留5min，相当于加一个自环，求最后无路可走的时间期望。

题解：期望+状压dp+记忆化搜索
因为n只有15，用状压dp。
我们用$dp[u][sta]$表示在点u，状态为sta（即表示走过了那些路口）的时间期望，知道起点，然后逆推。
用cnt表示可走的路数量，cnt为0时$dp[u][sta]=0$，否则
$dp[u][sta]=(dp[u][sta]+5+\Sigma (dp[v][sta2]+w))/(1+cnt)$
化简后得$dp[u][sta]=(\Sigma (dp[v][sta2]+w)+5)/cnt$，然后记忆化搜索。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<set>
#include<map>
#include<sstream>
#include<iomanip>
#define ll long long
using namespace std;
int t, n, m, u, v, w, k, head[222], vis[22][1<<15];
double dp[22][1 << 15];
struct node {
	int v, w, next;
}edge[222];
void add(int u, int v, int w) { //链式前向星存图
	edge[++k].next = head[u];
	edge[k].v = v;
	edge[k].w = w;
	head[u] = k;
}
int dfs(int u, int sta) {
	if (sta == (1 << n) - 1) {  //当前状态为末状态
		dp[u][sta] = 0;
		return 1;
	}
	if (vis[u][sta]) {
		if (dp[u][sta]) return 1;  //走过但有路可走
		else return 0;  //走过但无路可走
	}
	vis[u][sta] = 1;
	dp[u][sta] = 5;
	int cnt = 0; //可走数量
	for (int i = head[u]; i; i = edge[i].next) {
		int v = edge[i].v, w = edge[i].w;
		if ((sta & (1 << v)) == 0 && dfs(v, sta | (1 << v))) {  //当前状态不包含v且可从v走
			cnt++;
			dp[u][sta] += dp[v][sta | (1 << v)] + w;
		}
	}
	if (!cnt) {  //无路可走 置0
		dp[u][sta] = 0;
		return 0;
	}
	else {
		dp[u][sta] /= cnt;
		return 1;
	}
}
int main() {
	scanf("%d", &t);
	int cas = 0;
	while (t--) {
		k = 0;
		memset(head, 0, sizeof(head));
		memset(dp, 0, sizeof(dp));
		memset(vis, 0, sizeof(vis));
		scanf("%d%d", &n, &m);
		for (int i = 1; i <= m; i++) {
			scanf("%d%d%d", &u, &v, &w);
			add(u, v, w);
			add(v, u, w);
		}
		dfs(0, 1);
		printf("Case %d: %f\n", ++cas, dp[0][1]);
	}
	return 0;
}