首先求出SCC缩点,E[T]=0,按拓扑序计算
对于无边连出的块,如果不是T所在块,则称该块是死路块
对于一个块,如果其中的点连出的边是死路块,则它也是死路块
否则对于每块进行高斯消元求出期望
如果S点所在块为死路块,则答案为INF
#include<cstdio>
#include<cmath>
const int N=10010,M=1000010;
int n,m,x,y,i,j,S,T;
int g[3][N],nxt[3][M],v[3][M],ed,G[N],NXT[N],V[N],d[N],q[N],f[N],h,t,cnt,id[N],tot;bool vis[N],inf[N];
double e[N],a[110][110],ans[110],tmp;
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
inline void add(int x,int y){
v[0][++ed]=y;nxt[0][ed]=g[0][x];g[0][x]=ed;
v[1][ed]=x;nxt[1][ed]=g[1][y];g[1][y]=ed;
}
inline void ADD(int x,int y){d[y]++,v[2][++ed]=y;nxt[2][ed]=g[2][x];g[2][x]=ed;}
inline void addnode(int x,int y){V[++ed]=y;NXT[ed]=G[x];G[x]=ed;}
void dfs1(int x){
vis[x]=1;
for(int i=g[0][x];i;i=nxt[0][i])if(!vis[v[0][i]])dfs1(v[0][i]);
q[++t]=x;
}
void dfs2(int x,int y){
vis[x]=0,addnode(f[x]=y,x);
for(int i=g[1][x];i;i=nxt[1][i])if(vis[v[1][i]])dfs2(v[1][i],y);
}
inline void cal(int x){
int i,j,k,u,y,cnt;double t;
for(tot=0,i=G[x];i;i=NXT[i])id[V[i]]=++tot;
for(cnt=0,i=G[x];i;i=NXT[i]){
for(u=V[i],++cnt,j=1;j<=tot+1;j++)a[cnt][j]=0;
if(u==T){a[cnt][cnt]=1;continue;}
for(j=g[0][u];j;j=nxt[0][j]){
if(inf[f[y=v[0][j]]]){inf[x]=1;return;}
a[cnt][cnt]+=1,a[cnt][tot+1]+=1;
if(f[y]==x)a[cnt][id[y]]-=1;else a[cnt][tot+1]+=e[y];
}
}
for(i=1;i<=tot;i++){
for(k=i,j=i+1;j<=tot;j++)if(std::fabs(a[j][i])>std::fabs(a[k][i]))k=j;
if(k!=i)for(j=i;j<=tot+1;j++)tmp=a[i][j],a[i][j]=a[k][j],a[k][j]=tmp;
for(j=i+1;j<=tot;j++)for(t=a[j][i]/a[i][i],k=i;k<=tot+1;k++)a[j][k]-=a[i][k]*t;
}
for(ans[tot]=a[tot][tot+1]/a[tot][tot],i=tot-1;i;i--){
for(ans[i]=a[i][tot+1],j=tot;j>i;j--)ans[i]-=ans[j]*a[i][j];
ans[i]/=a[i][i];
}
for(i=G[x];i;i=NXT[i])e[V[i]]=ans[id[V[i]]];
}
int main(){
read(n),read(m),read(S),read(T);
while(m--)read(x),read(y),add(x,y);
for(i=1;i<=n;i++)if(!vis[i])dfs1(i);
for(ed=0,i=n;i;i--)if(vis[q[i]])dfs2(q[i],++cnt);
for(ed=0,i=1;i<=n;i++)for(j=g[0][i];j;j=nxt[0][j])if(f[i]!=f[v[0][j]])ADD(f[v[0][j]],f[i]);
for(i=1;i<=cnt;i++)if(!d[i])if(f[T]!=i)inf[i]=1;
for(h=i=1,t=0;i<=cnt;i++)if(!d[i]){
q[++t]=i;
if(f[T]==i)cal(i);
}
while(h<=t)for(i=g[2][x=q[h++]];i;i=nxt[2][i])if(!(--d[v[2][i]]))cal(q[++t]=v[2][i]);
if(inf[f[S]])puts("INF");else printf("%.3f",e[S]);
return 0;
}