#include <bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
void add(int &x,int y){
    if((x+=y)>=mod)x-=mod;
}
int solve()
{
    int n;cin>>n;
    vector<int>a(n);
    for(int i=0;i<n;i++)cin>>a[i],a[i]--;
    vector<vector<pair<int,int>>> p(n);
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++)if(a[i]==a[j])p[a[i]].emplace_back(i,j);
    }   
    vector<string> s(n);
    for(int i=0;i<n;i++)cin>>s[i];
    vector<vector<int>>val(n,vector<int>(n));
    for(int i=0;i<n;i++){
        vector<vector<int>>now(n,vector<int>(n));
        for(int j=0;j<n;j++)if(s[j][a[i]]=='1'){
            for(auto [x,y]:p[j])now[x][y]=val[x][y];
        }
        for(int j=0;j<n;j++){
            for(int k=1;k<n;k++){
                add(now[j][k],now[j][k-1]);
            }
        }
        for(int j=1;j<n;j++){
            for(int k=0;k<n;k++)add(now[j][k],now[j-1][k]);
        }
        for(auto [x,y]:p[a[i]]){
            add(val[x][y],1);
            if(x>0&&y>0)add(val[x][y],now[x-1][y-1]);
        }
    }
    int ans=0;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)add(ans,val[i][j]);
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    int T;cin>>T;while(T--)cout<<solve()<<"\n";
    return 0;
}