#include <cstdio>
#include <cstring>
#include <cmath>
const int INF=0x3fffffff;
typedef long long ll;
const int mod=5201314;
const int maxn=505;
int n,dp[2][maxn][maxn];
char val[maxn][maxn];
void add(int &x, int y){
x=(x+y)%mod;
}
void input(){
int i;
scanf("%d",&n);
for (i=1;i<=n;i++)
scanf("%s",val[i]+1);
}
void solve(){
int i,j,k,x1,x2,y1,y2;
dp[0][1][n]=(val[1][1]==val[n][n]);
k=1;
for(i=2;i<=n;i++){
memset(dp[k],0,sizeof(dp[k]));
for(x1=1;x1<=i;x1++)
for(x2=n;x2>=n-i+1;x2--){
y1=i-x1+1;
y2=2*n-x2-i+1;
if(val[x1][y1]!=val[x2][y2])
continue;
add(dp[k][x1][x2],dp[!k][x1][x2]);
add(dp[k][x1][x2],dp[!k][x1-1][x2]);
add(dp[k][x1][x2],dp[!k][x1][x2+1]);
add(dp[k][x1][x2],dp[!k][x1-1][x2+1]);
}
k=!k;
}
int ans=0;
for(i=1;i<=n;i++)
add(ans,dp[!k][i][i]);
printf("%d\n",ans);
}
int main(){
int t;
scanf("%d", &t);
while (t--){
input();
solve();
}
return 0;
}