Maze （HDU

  题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4035

 

  题解：

　　一道经典的循环期望，难点在于公式的推导。

　　推导太复杂，看着题解都推了好久，不想再用公式编辑器写一遍了，想看的点这里。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
#define RI register int
#define eps 1e-10
using namespace std;
const int INF = 0x7ffffff ;
const int N = 10000 + 10 ;
const int M = 1e7 + 10 ;

inline int read() {
    int k = 0 , f = 1 ; char c = getchar() ;
    for( ; !isdigit(c) ; c = getchar())
      if(c == '-') f = -1 ;
    for( ; isdigit(c) ; c = getchar())
      k = k*10 + c-'0' ;
    return k*f ;
}
struct Edge {
    int to, next ;
}ee[M] ;
int n, cnt = 0 ; int head[N] ; double k[N], e[N], A[N], B[N], C[N], in[N] ;
inline void add_edge(int x,int y) {
    ee[++cnt].to = y, ee[cnt].next = head[x], head[x] = cnt ;
}

void dfs(int x,int f) {
    A[x] = k[x]*in[x], B[x] = 1.0-k[x]-e[x], C[x] = (1-k[x]-e[x])*in[x] ; double fm = in[x] ;
    for(int i=head[x];i;i=ee[i].next) {
        int y = ee[i].to ; if(y == f) continue ;
        dfs(y,x) ;
        A[x] += B[x]*A[y] ; C[x] += B[x]*C[y] ; fm -= B[x]*B[y] ;
    }
    A[x] /= fm, B[x] /= fm, C[x] /= fm ;
}
inline void solve(int w) {
    n = read() ; memset(head,0,sizeof(head)) ; memset(in,0,sizeof(in)) ; cnt = 0 ;
    for(int i=1;i<n;i++) {
        int x = read(), y = read() ; in[x]++, in[y]++ ;
        add_edge(x,y) ; add_edge(y,x) ;
    }
    for(int i=1;i<=n;i++) k[i] = (double)read()/100.0, e[i] = (double)read()/100.0 ;
    dfs(1,0) ;
    double ans = C[1]/(1.0-A[1]) ;
    printf("Case %d: ",w) ; 
    if(abs(A[1]-1.0) < eps) printf("impossible
") ;
    else printf("%.6lf
",ans) ;
}

int main() {
    int t = read() ;
    for(int w=1;w<=t;w++) solve(w) ;
    return 0 ;
}

另：

•  这题精度一定要开够，一开始开了1e-5 WA了，看别人博客发现也有人遇到了同样的问题=-=，开到1e-10即可。
•  库中有的函数不要自己写（尤其是与库中函数同名）。一开始没加cmath库，自己写了个abs（）函数，疯狂CE =-=