light1287Where to Run：概率DP

lightoj1287

题意：如果能通过下一个点走完所有的点，那么这个点就叫EJ。如果EJ个数为0,不行动。否则,可以选择往任意一个EJ走或者在原地停留5分钟。求警察抓到你的期望时间。

题解

• 状态压缩+概率DP。
• 假设在u点有k种选择，包括停留5分钟。那么那么状态转移方程为：
• 则。注意当往其它点走时，要保证能够走完所有的点。
• dp[s][u]表示已经走完的状态为s，目前处于i点，走完所有还需要的期望时间。

代码

#include <bits/stdc++.h>
using namespace std;
int const N = 16;
int n,m,mp[N][N];
double dp[1<<N][N];
bool vis[1<<N][N],canreach[1<<N][N];
bool DP(int s,int u){   //已经遍历的点的状态为s,目前处于点u;
	if(s == (1<<n)-1)	return dp[s][u] = 0,true;  //可以到达
	if(vis[s][u])	return canreach[s][u];
	vis[s][u] = true;
	int cnt = 0;
	double res = 0;
	for(int i=0;i<n;i++){
		if(!mp[u][i] || s & (1<<i))	continue;  //没有通路，或者已经走过
		if(DP(s|(1<<i),i)){   //判断能否到达，即使有u通往i的路，但是不一定能够到达所有的点
			cnt++;
			canreach[s][u] = true;
			res += (mp[u][i] + dp[s|1<<i][i]);
		}
	}
	if(!canreach[s][u])	return dp[s][u] = 0,false;   //如果不能到达，那么即不能等待也不能继续往前走
	else return dp[s][u] = (res + 5.0) / cnt;
}
int main(){
	int T,caser = 0;
	scanf("%d",&T);
	while(T--){
		memset(mp,0,sizeof(mp));
		memset(dp,0,sizeof(dp));
		memset(vis,false,sizeof(vis));
		memset(canreach,false,sizeof(canreach));
		scanf("%d%d",&n,&m);
		for(int i=1;i<=m;i++){
			int u,v,c;
			scanf("%d%d%d",&u,&v,&c);
			mp[u][v] = mp[v][u] = c;
		}
		DP(1,0);
		printf("Case %d: %.6f\n",++caser,dp[1][0]);
	}
	return 0;
}