bzoj3470 Freda’s Walk（期望dp+概率dp+拓扑序+数学）

先正着拓扑序一遍，求出p[i],表示从1出发到达i点的概率。
然后倒着拓扑序一遍，求出f[i]，表示从i开始走的期望步数。
这时候的f[1]就是不删边时的期望步数了。
我们考虑如果删掉边x->y,w，对答案有什么影响。
显然对x之前的点是没有任何影响的。
原本的$f[x]=\frac{(f[y]+1)*w(x,y)+\sum_v (f[v]+1)*w(x,v)}{sum[x]}$
而现在的$f[x]=\frac{\sum_v (f[v]+1)*w(x,v)}{sum[x]-w(x,y)}$
他们的差值*p[x]就是对答案的影响。因此我们循环m条边，取最大的即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 10010
#define M 100010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,h[N],num=0,in[N],sum[N];
double p[N],f[N];
struct Edge{
    int x,y,val;
}e[M];
struct edge{
    int to,next,val;
}data[M];
inline void add(int x,int y,int val){
    data[++num].to=y;data[num].next=h[x];h[x]=num;
    data[num].val=val;in[y]++;
}
queue<int>q;
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=m;++i){
        e[i].x=read()+1,e[i].y=read()+1,e[i].val=read();
        add(e[i].x,e[i].y,e[i].val);sum[e[i].x]+=e[i].val;
    }p[1]=1;
    for(int i=1;i<=n;++i) if(!in[i]) q.push(i);
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;p[y]+=p[x]*data[i].val/sum[x];
            if(--in[y]==0) q.push(y);
        }
    }memset(h,0,sizeof(h));num=0;memset(in,0,sizeof(in));
    for(int i=1;i<=m;++i) add(e[i].y,e[i].x,e[i].val);
    for(int i=1;i<=n;++i) if(!in[i]) q.push(i);
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;f[y]+=(f[x]+1)*data[i].val/sum[y];
            if(--in[y]==0) q.push(y);
        }
    }double ans=f[1];
    for(int i=1;i<=m;++i){
        int x=e[i].x,y=e[i].y,val=e[i].val;
        double del=(f[x]-(f[y]+1)*val/sum[x])*sum[x]/(sum[x]-val)-f[x];
        ans=max(ans,f[1]+p[x]*del);
    }printf("%.6lf\n",ans);
    return 0;
}