#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=210,mod=1e9+7;
int n,a[N],f[N][N],g[N][N];
char s[N][N];
signed main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	int T;
	cin>>T;
	while(T--)
	{
		memset(f,0,sizeof(f));
		cin>>n;
		for(int i=1;i<=n;i++)cin>>a[i];
		for(int i=1;i<=n;i++)cin>>(s[i]+1);
		for(int i=1;i<=n;i++)
		{
			memset(g,0,sizeof(g));
			for(int j=1;j<=n;j++)
				for(int k=1;k<=n;k++)
					g[j][k]=((g[j-1][k]+g[j][k-1]-g[j-1][k-1]+f[j-1][k-1]*(s[a[j-1]][a[i]]^48))%mod+mod)%mod;
			for(int j=1;j<=n;j++)
				for(int k=1;k<=n;k++)
					f[j][k]=(f[j][k]+(g[j][k]+1)*(a[i]==a[j])*(a[i]==a[k]))%mod;
		}
		int ans=0;
		for(int j=1;j<=n;j++)
			for(int k=1;k<=n;k++)
				ans=(ans+f[j][k])%mod;
		cout<<ans<<"\n";
	}
	return 0;
}