// Hydro submission #66ffc8ea55b5e82bc01754fe@1728039147013
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1e9+7;
int n;
int a[205];
bool e[205][205];
int f[205][205],g[205][205];
void test()
{
	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++)
		{
			char x;
			cin>>x;
			e[i][j]=x-'0';
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			f[i][j]=g[i][j]=0;
	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(a[i]==a[j]&&a[j]==a[k])
					f[j][k]=g[j][k]+1;
				else
					f[j][k]=0;
		for(int j=1;j<=n;j++)
			for(int k=1;k<=n;k++)
				f[j][k]=(f[j][k]+f[j-1][k]+f[j][k-1]-f[j-1][k-1]+mod)%mod;
		(ans+=f[n][n])%=mod;
		for(int j=1;j<=n;j++)
			for(int k=1;k<=n;k++)
				if(a[j]==a[k]&&e[a[i]][a[j]])
					(g[j][k]+=f[j-1][k-1])%=mod;
	}
	cout<<ans<<"\n";
}
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	int t;
	cin>>t;
	for(int i=1;i<=t;i++)
		test();
}
/*
1
4
1 2 1 2
1111
1111
1111
1111
*/