BZOJ 3470: Freda’s Walk

没人写题解的样子。。 有生之年第一次刷了个rank1 望dalao们轻虐

正着走一遍 然后求出从1到每个点的概率g
反着走一遍 求出每个点对1的贡献f（过程跟求期望没什么两样）

最后考虑割一条边
算一个点x以及其以下的点对1的贡献是 ∑(1+f[y])∗w[k]∑w[y] （有x到y的边k 边权为w[k]）
如果割掉一条边
我们对于每个x算出分子和分母 那对y就等于让分子和分母都减去一个数再算一下就好了
答案就是 原 ans+max( (max_x-g[x])*g[x] ) ) max_x就是上面那个分数的最大值。

#include<bits/stdc++.h>using namespace std;const int N=10002,M=100002;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<<1)+(x<<3)+ch-'0'; ch=getchar();}    return x*f;}struct E{int x,y,c,next;}a[M],e[M]; int first[N],fir[N];int in[N];queue<int>q;double f[N],g[N],mx,s[N];void solve(int x){    double t=0,d=f[x];    for(int k=first[x];k;k=a[k].next){        int y=a[k].y;        t+=(1+f[y])*a[k].c;    }    for(int k=first[x];k;k=a[k].next){        int y=a[k].y;        double tt=t-(1+f[y])*a[k].c;        if(s[x]>a[k].c)d=max(d,tt/(s[x]-a[k].c));    }    mx=max(mx,(d-f[x])*g[x]);}int main(){    int n=read(),m=read(),i;    for(i=1;i<=m;i++){        int x=read()+1,y=read()+1,c=read(); s[x]+=c;        a[i]=(E){x,y,c,first[x]},first[x]=i,in[y]++;    }    for(i=1;i<=n;i++)if(!in[i])q.push(i);    g[1]=1;    while(!q.empty()){        int x=q.front(); q.pop();        for(int k=first[x];k;k=a[k].next){            int y=a[k].y;            g[y]+=g[x]*a[k].c/s[x],in[y]--;            if(!in[y])q.push(y);        }    }    for(i=1;i<=m;i++)        e[i]=(E){a[i].y,a[i].x,a[i].c,fir[a[i].y]},fir[a[i].y]=i,in[a[i].x]++;    for(i=1;i<=n;i++)if(!in[i])q.push(i);    while(!q.empty()){        int x=q.front(); q.pop();        for(int k=fir[x];k;k=e[k].next){            int y=e[k].y;            f[y]+=(f[x]+1)*e[k].c/s[y],in[y]--;            if(!in[y])q.push(y);        }    }    double ans=f[1]; mx=0;    for(i=1;i<=n;i++)solve(i);    printf("%.6lf\n",ans+mx);    return 0;}