HDU4035 Maze 【树形DP】【期望DP】

题目分析:

以前一直不会这个方法, 我好菜啊。

转移分为三个部分,一个是直接成功,一个是转移到E1,还有一个是转移到自己周围的一圈儿点。 如果是叶子那么只能转移到父亲,如果不是叶子可以把非叶子的转移代换,这样也只转移到父亲,判一下无解就行了。

代码:

 #include<bits/stdc++.h>
using namespace std; const int maxn = ;
const double eps = 1e-; int n,Tmp,cas;
double a[maxn],b[maxn];
vector<int> g[maxn]; double A[maxn],B[maxn],C[maxn]; void dfs(int now,int f){
int cnt = ;
for(int i=;i<g[now].size();i++){
if(g[now][i] == f) continue;
cnt++; dfs(g[now][i],now);
}
if(cnt == ){A[now] = a[now]; B[now] = -a[now]-b[now]; C[now] = B[now];}
else{
double alpha = ,beta = ,gamma = ;
for(int i=;i<g[now].size();i++){
if(g[now][i] == f) continue;
alpha += A[g[now][i]];beta += B[g[now][i]];gamma += C[g[now][i]];
}
if(f != ){
gamma += cnt+;
alpha *= (-a[now]-b[now])/(cnt+);
beta *= (-a[now]-b[now])/(cnt+);
gamma *= (-a[now]-b[now])/(cnt+);
B[now] = (-a[now]-b[now])/(cnt+);
beta = -beta; C[now] = gamma; A[now] = alpha+a[now];
C[now] /= beta; A[now] /= beta; B[now] /= beta;
}else{
gamma += cnt;
alpha *= (-a[now]-b[now])/cnt;
beta *= (-a[now]-b[now])/cnt;
gamma *= (-a[now]-b[now])/cnt;
B[now] = ; C[now] = gamma; beta = -(alpha+beta+a[now]);
if(beta <= eps){C[now] = -;}
else C[now] /= beta;
}
}
} void read(){
scanf("%d",&n);
memset(A,,sizeof(A));
memset(B,,sizeof(B));
memset(C,,sizeof(C));
memset(a,,sizeof(a));
memset(b,,sizeof(b));
for(int i=;i<=n;i++) g[i].clear();
for(int i=;i<n;i++){
int u,v; scanf("%d%d",&u,&v);
g[u].push_back(v); g[v].push_back(u);
}
for(int i=;i<=n;i++) {
int u,v; scanf("%d%d",&u,&v);
a[i] = u/100.0,b[i] = v/100.0;
}
} void work(){
dfs(,);
printf("Case %d: ",cas);
if(C[] == -){puts("impossible");}
else{printf("%.6lf\n",C[]);}
} int main(){
scanf("%d",&Tmp);
while(cas++,Tmp--){
read();
work();
}
return ;
}
上一篇:[CQOI2018]破解D-H协议


下一篇:String的intern方法的用处