TYVJ P1933 「Poetize3」绿豆蛙的归宿 数学期望DP

传送门

---

很明显，如果我们正着推是错的，所以应该倒着推（额(⊙o⊙)…）
f[i]代表从i到n的路径期望值，f[i]+=(f[root]+w(i,root))/k[i]（注意：边是反着存的Notice！！！）,然后按照拓扑序更新即可

---

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
const int maxn=100000+5,maxm=2*maxn;
int n,m,hd[maxn],to[maxm],nxt[maxm],w[maxm],out[maxn],k[maxn],cnt;
double f[maxn],s[maxn];//f[i]代表从1到i的路径期望长度，s 
void add(int x,int y,int z){
    to[cnt]=y;
    w[cnt]=z;
    nxt[cnt]=hd[x];
    hd[x]=cnt++;
}
queue<int> q;
void bfs(int root){
    q.push(root);
    while(!q.empty()){
        int v=q.front();
        q.pop();
        for(int i=hd[v];i!=-1;i=nxt[i]){
            f[to[i]]+=(f[v]+w[i])/(double)(k[to[i]]);
            out[to[i]]--;
            if(out[to[i]]==0)
                q.push(to[i]); 
        }
    } 
}
signed main(void){
    scanf("%d%d",&n,&m),cnt=0;
    memset(out,0,sizeof(out));
    memset(hd,-1,sizeof(hd));
    memset(k,0,sizeof(k));
    memset(f,0,sizeof(f));
    for(int i=1,a,b,c;i<=m;i++)
        scanf("%d%d%d",&a,&b,&c),add(b,a,c),out[a]++,k[a]++;
    bfs(n);
    printf("%.2f\n",f[1]);
    return 0;
}

---

by >o< neighthorn