// Hydro submission #66ffeb8555b5e82bc017aaa7@1728048006090
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=205,p=1e9+7;
char c[N][N];
int t,n,a[N],f[N][N],sum[N][N];
signed main() {
	cin>>t;
	while(t--) {
		cin>>n;
		for(int i=1; i<=n; i++) cin>>a[i];
		for(int i=1; i<=n; i++) for(int j=1; j<=n; j++) cin>>c[i][j];
		memset(f,0,sizeof(f));
		int ans=0;
		for(int i=1; i<=n; i++) {
			for(int j=1; j<=n; j++) {
				for(int k=1; k<=n; k++) {
					if(c[a[j]][a[i]]=='1') sum[j+1][k+1]=f[j][k];
					else sum[j+1][k+1]=0;
					sum[j+1][k+1]=(sum[j+1][k+1]+sum[j][k+1]+sum[j+1][k]-sum[j][k]+p)%p;
				}
			}
			for(int j=1; j<=n; j++) {
				for(int k=1; k<=n; k++) {
					if(a[i]==a[j]&&a[i]==a[k]) {
						ans=(ans+sum[j][k]+1)%p;
						f[j][k]=(f[j][k]+sum[j][k]+1)%p;
					}
				}
			}
		}
		cout<<ans<<endl;
	}
	return 0;
}