#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 210,mod = 1e9+7;
int addv(int x,int y){return (x+y>=mod)?(x+y-mod):(x+y);}
int subv(int x,int y){return (x>=y)?(x-y):(x-y+mod);}
void add(int& x,int y){x = addv(x,y);}
void sub(int& x,int y){x = subv(x,y);}
int qpow(int a,int n){
    int res = 1; a%=mod;
    while(n){
        if(n&1)res = 1ll*res*a%mod;
        n>>=1; a = 1ll*a*a%mod;
    }
    return res;
}
int qinv(int a){return qpow(a,mod-2);}
//
int T,n,a[N],e[N][N],f[N][N][N],g[N][N],h[N][N],ans;

void solve(){
    cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    memset(e,0,sizeof e);
    for(int i=1;i<=n;i++){
        string s; cin>>s; for(int j=1;j<=n;j++)e[i][j] = s[j-1]-'0';
    }
    memset(f,0,sizeof f); 
    memset(g,0,sizeof g);
    ans = 0;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)for(int k=1;k<=n;k++)if(a[i] == a[j] && a[j] == a[k]){
            f[i][j][k] = addv(g[j][k],1);
            add(ans,f[i][j][k]);
        }
        memcpy(h,f[i],sizeof h);
        for(int j=1;j<=n;j++)for(int k=1;k<=n;k++)add(h[j][k],h[j][k-1]);
        for(int j=1;j<=n;j++)for(int k=1;k<=n;k++)add(h[j][k],h[j-1][k]);
        for(int j=1;j<=n;j++)for(int k=1;k<=n;k++)if(a[j] == a[k] && e[a[i]][a[j]]){
            add(g[j][k],h[j-1][k-1]);
        }
     }
    cout<<ans<<endl;
}

int main(){
    cin>>T;
    while(T--)solve();

    return 0;
}