dfs会T，只好正反两遍拓扑了……
#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
using namespace std;
const int N=100005;
int n,m,h[N],cnt,d1[N],d2[N];
double f[N],p[N],s[N],ans;
struct qwe
{
int ne,no,to;
double va;
}e[N];
vector<pair<int,double> >a[N];
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
void add(int u,int v,int w)
{
cnt++;
e[cnt].ne=h[u];
e[cnt].no=u;
e[cnt].to=v;
e[cnt].va=w;
h[u]=cnt;
}
int main()
{
n=read(),m=read();
for(int i=1;i<=m;i++)
{
int x=read()+1,y=read()+1,z=read();
s[x]+=z;
a[y].push_back(make_pair(x,z)),d1[x]++;
add(x,y,z),d2[y]++;
}
queue<int>q;
for(int i=1;i<=n;i++)
if(!d1[i])
q.push(i);
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=0;i<a[u].size();i++)
{
f[a[u][i].first]+=(f[u]+1.0)*a[u][i].second/s[a[u][i].first];
if(!(--d1[a[u][i].first]))
q.push(a[u][i].first);
}
}
p[1]=1;
for(int i=1;i<=n;i++)
if(!d2[i])
q.push(i);
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=h[u];i;i=e[i].ne)
{
p[e[i].to]+=p[u]*e[i].va/s[u];
if(!(--d2[e[i].to]))
q.push(e[i].to);
}
}
ans=f[1];
for(int i=1;i<=cnt;i++)
ans=max(ans,f[1]+((f[e[i].no]-(f[e[i].to]+1)*e[i].va/s[e[i].no])*s[e[i].no]/(s[e[i].no]-e[i].va)-f[e[i].no])*p[e[i].no]);
printf("%.6f
",ans);
return 0;
}