題目鏈接：點擊查看

題目大意：給出一個 n 個點的樹，初始時位于點 1 ，每個點都會給出一個概率 k[ i ] 和一個概率 e[ i?] ，其代表的意義分別是：

k[ i ] ：在點 i 有 k[ i ] 的幾率死亡，復活到點 1

e[ i ] ：在點 i 有 e[ i ] 的幾率逃生

1 - k[ i ] - e[ i ] ：等概率經過周圍的邊到達相鄰的點

問從點 1 開始逃生的期望步數

題目分析：

樹形dp維護一下概率dp所需要的變量即可，記得特判無解的情況，有個小坑就是會卡精度，eps調成 1e-9 或 1e-10 才能過

代碼：
?

#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cmath> #include<cstring> #include<algorithm> #include<stack> #include<climits> #include<queue> #include<map> #include<set> #include<sstream> #include<cassert> #include<bitset> using namespace std;typedef long long LL;typedef unsigned long long ull;const int inf=0x3f3f3f3f;const int N=1e4+100;const double eps=1e-10;int sgn(double x) {if(fabs(x)<eps)return 0;if(x<0)return -1;else return 1; }vector<int>node[N];double k[N],e[N],t[N],A[N],B[N],C[N];void init(int n) {for(int i=1;i<=n;i++)node[i].clear(); }bool dfs(int u,int fa) {int m=node[u].size();A[u]=k[u];B[u]=t[u]/m;C[u]=t[u];double temp=0;for(auto v:node[u]){if(v==fa)continue;if(!dfs(v,u))return false;A[u]+=t[u]/m*A[v];C[u]+=t[u]/m*C[v];temp+=t[u]/m*B[v];}if(sgn(1-temp)==0)return false;A[u]/=(1-temp);B[u]/=(1-temp);C[u]/=(1-temp);return true; }int main() { #ifndef ONLINE_JUDGE // freopen("data.in.txt","r",stdin); // freopen("data.out.txt","w",stdout); #endif // ios::sync_with_stdio(false);int w;cin>>w;int kase=0;while(w--){int n;scanf("%d",&n);init(n);for(int i=1;i<n;i++){int u,v;scanf("%d%d",&u,&v);node[u].push_back(v);node[v].push_back(u);}for(int i=1;i<=n;i++){scanf("%lf%lf",k+i,e+i);k[i]/=100;e[i]/=100;t[i]=1-k[i]-e[i];}if(dfs(1,-1)&&sgn(1-A[1])!=0)printf("Case %d: %.10f\n",++kase,C[1]/(1-A[1]));elseprintf("Case %d: impossible\n",++kase);}return 0; }

?

?

總結

以上是生活随笔為你收集整理的HDU - 4035 Maze(概率dp)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。