#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return true; } return false; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return true; } return false; }
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define si(x) int(x.size())
#define double long double
const int mod=1000000007,MAX=105,INF=100000;
const double eps=1e-10;
typedef vector<double> vec;
typedef vector<vec> mat;

int N;
int dis[MAX][MAX];
int edge[MAX][MAX];

void floyd(){
    for(int k=0;k<N;k++){
        for(int i=0;i<N;i++){
            for(int j=0;j<N;j++){
                chmin(dis[i][j],dis[i][k]+dis[k][j]);
            }
        }
    }
}

vec gauss_jordan(const mat& A,const vec &b){
    int n=A.size();
    mat B(n,vec(n+1));
    
    for(int i=0;i<n;i++) for(int j=0;j<n;j++) B[i][j]=A[i][j];
    
    for(int i=0;i<n;i++) B[i][n]=b[i];
    
    for(int i=0;i<n;i++){
        int pivot=i;
        for(int j=i;j<n;j++){
            if(abs(B[j][i])>abs(B[pivot][i])) pivot=j;
        }
        swap(B[i],B[pivot]);
        
        if(abs(B[i][i])<eps) return vec();
        
        for(int j=i+1;j<=n;j++) B[i][j]/=B[i][i];
        for(int j=0;j<n;j++){
            if(i!=j){
                for(int k=i+1;k<=n;k++) B[j][k]-=B[j][i]*B[i][k];
            }
        }
    }
    
    vec x(n);
    
    for(int i=0;i<n;i++) x[i]=B[i][n];
    return x;
}

int main(){
    
    std::ifstream in("text.txt");
    std::cin.rdbuf(in.rdbuf());
    cin.tie(0);
    ios::sync_with_stdio(false);
    
    while(1){
        int s,t;cin>>N>>s>>t;
        memset(edge,0,sizeof(edge));
        if(N+s+t==0) break;
        s--;t--;
        vector<int> X(N);
        for(int i=0;i<N;i++) cin>>X[i];
        for(int i=0;i<N;i++){
            for(int j=0;j<N;j++){
                cin>>dis[i][j];
                if(dis[i][j]==0){
                    if(i!=j) dis[i][j]=INF;
                }else{
                    edge[i][j]=dis[i][j];
                }
            }
        }
        floyd();
        
        mat A(N,vec(N));
        vec B(N);
        
        for(int i=0;i<N;i++){
            if(i==t){
                A[i][i]=1;
            }else{
                vector<int> use;
                int sum=0;
                if(X[i]){
                    for(int j=0;j<N;j++){
                        if(edge[i][j]&&edge[i][j]+dis[j][t]==dis[i][t]){
                            use.push_back(j);
                            sum+=edge[i][j];
                        }
                    }
                    
                    for(int u:use) A[i][u]=-1.0;
                    A[i][i]=si(use);
                    
                    B[i]=sum;
                    
                }else{
                    for(int j=0;j<N;j++){
                        if(edge[i][j]){
                            use.push_back(j);
                            sum+=edge[i][j];
                        }
                    }
                    
                    if(edge[i][i]==0){
                        for(int u:use) A[i][u]=-1.0;
                        A[i][i]=si(use);
                        
                        B[i]=sum;
                    }else{
                        for(int u:use){
                            if(u==i) continue;
                            A[i][u]=-1.0;
                        }
                        
                        A[i][i]=si(use)-1.0;
                        
                        B[i]=sum;
                    }
                }
            }
        }
        
        vec c=gauss_jordan(A,B);
        
        cout<<setprecision(25)<<c[s]<<endl;
        
    }
}





Measured memory_bytes: 3572