没人写题解的样子。。 有生之年第一次刷了个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;
}