bzoj 1415 无环期望

 #include <cstdio>
 #include <vector>
 #include <queue>
 #include <algorithm>
 #define maxn 1010
 using namespace std;

 int n, m, a, b;
 vector<int> g[maxn];
 int s[maxn][maxn];
 double f[maxn][maxn];

 void bfs( int src ) {
     queue<int> qu;
     s[src][src] = src;
     for( int t=; t<g[src].size(); t++ ) {
         int v = g[src][t];
         s[src][v] = v;
         qu.push( v );
     }
     while( !qu.empty() ) {
         int u = qu.front();
         qu.pop();
         for( int t=; t<g[u].size(); t++ ) {
             int v = g[u][t];
             if( s[src][v] ) continue;
             s[src][v] = s[src][u];
             qu.push(v);
         }
     }
 }
 void dodp( int i, int j ) {
     double &dv = f[i][j];
     if( dv>-0.5 ) return;
     if( i==j ) {
         dv = 0.0;
         return;
     }
     if( s[i][j]==j || s[s[i][j]][j]==j ) {
         dv = 1.0;
         return;
     }
     dv = 0.0;
     for( int t=; t<g[j].size(); t++ ) {
         int v = g[j][t];
         dodp( s[s[i][j]][j], v );
         dv += f[s[s[i][j]][j]][v];
     }
     dv /= g[j].size();
     dv += 1.0;
 }

 int main() {
     scanf( "%d%d%d%d", &n, &m, &a, &b );
     for( int i=,u,v; i<=m; i++ ) {
         scanf( "%d%d", &u, &v );
         g[u].push_back(v);
         g[v].push_back(u);
     }
     for( int i=; i<=n; i++ ) sort( g[i].begin(), g[i].end() );
     for( int i=; i<=n; i++ ) bfs(i);
     for( int i=; i<=n; i++ )
         g[i].push_back(i);
     for( int i=; i<=n; i++ )
         for( int j=; j<=n; j++ ) f[i][j] = -1.0;
     dodp( a, b );
     printf( "%.3lf\n", f[a][b] );
 }