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