#include <bits/stdc++.h>
#define MOD 998244353
using namespace std;

int n,m,p[410],id[410],num[410],a[410][410],b[410][410],f[410][410],c[410][410],dp[410],pw1[100010],pw2[100010];
vector<int> vec[410];

inline int find_set(int x)
{
	return p[x]==x?x:p[x]=find_set(p[x]);
}

inline void add1(int &x,int y)
{
	x+=y;
	if(x>=MOD) x-=MOD;
}

inline void add2(int &x,int y)
{
	x+=y;
	if(x<0) x+=MOD; 
}

inline int my_pow(int x,int y)
{
	if(y==0) return 1;
	if(y==1) return x;
	int res=my_pow(x,y/2);
	if(y%2==0) return 1ll*res*res%MOD;
	else return 1ll*res*res%MOD*x%MOD; 
}

signed main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int t;cin>>t;
	for(int i=0;i<=400;i++) c[i][0]=1;
	for(int i=1;i<=400;i++){
		for(int j=1;j<=400;j++){
			c[i][j]=c[i-1][j];
			add1(c[i][j],c[i-1][j-1]);
		}
	}
	while(t--){
		cin>>n>>m;
		for(int i=1;i<=n;i++) p[i]=i,vec[i].clear();
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				cin>>a[i][j];
				if(a[i][j]==0){
					int rootx=find_set(i),rooty=find_set(j);
					if(rootx!=rooty) p[rooty]=rootx;
				}
			}
		}
		for(int i=1;i<=n;i++) vec[find_set(i)].push_back(i);
		bool flag=true;
		for(int i=1;i<=n;i++){
			for(int j=i+1;j<=n;j++) if(a[i][j]!=a[j][i]) flag=false;
		}
		int tot=0;
		for(int i=1;i<=n;i++) if(find_set(i)==i){
			tot++;num[tot]=0;
			for(int j=0;j<(int)vec[i].size();j++) id[vec[i][j]]=tot,num[tot]++;
			for(int j=0;j<(int)vec[i].size();j++){
				for(int k=j+1;k<(int)vec[i].size();k++){
					if(a[vec[i][j]][vec[i][k]]!=0) flag=false;
				}
			}
		}
		for(int i=1;i<=tot;i++) for(int j=1;j<=tot;j++) b[i][j]=0; 
		for(int i=1;i<=n;i++){
			for(int j=i+1;j<=n;j++){
				if(b[id[i]][id[j]]==0) b[id[i]][id[j]]=b[id[j]][id[i]]=a[i][j];
				else{
					if(a[i][j]!=b[id[i]][id[j]]) flag=false;
				}
			}
		}
		if(!flag){
			cout<<0<<'\n';
			continue;
		}
		for(int i=1;i<=tot;i++){
			for(int j=1;j<=tot;j++){
				f[i][j]=b[i][j];
			}
		}
		for(int k=1;k<=tot;k++){
			for(int i=1;i<=tot;i++){
				for(int j=1;j<=tot;j++){
					f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
				}
			}
		}
		int ans=0,res=1;
		dp[1]=1;
		for(int i=2;i<=n;i++){
			dp[i]=my_pow(m+1,i*(i-1)/2);
			for(int j=0;j<i-1;j++){
				add2(dp[i],-1ll*dp[j+1]*c[i-1][j]%MOD*my_pow(m,(j+1)*(i-j-1))%MOD*my_pow(m+1,(i-j-1)*(i-j-2)/2)%MOD);
			}
		}
		for(int i=1;i<=tot;i++) res=1ll*res*dp[num[i]]%MOD;
		for(int i=1;i<=tot;i++){
			for(int j=i+1;j<=tot;j++){
				if(f[i][j]<b[i][j]){
					res=0;continue;
				}
				bool ok=false;
				for(int k=1;k<=tot;k++) if(k!=i&&k!=j){
					if(f[i][k]+f[k][j]==b[i][j]){
						ok=true;break;
					}
				}
				if(m<b[i][j]){
					res=0;continue;
				}
				if(ok){
					res=1ll*res*my_pow(m-b[i][j]+1,num[i]*num[j])%MOD;
				}
				else{
					int cur=my_pow(m-b[i][j]+1,num[i]*num[j]);
					add2(cur,-my_pow(m-b[i][j],num[i]*num[j]));
					res=1ll*res*cur%MOD;
				}
			}
		}
		add1(ans,res);
		if(m>0){
			m--;res=1;
			dp[1]=1;
			for(int i=2;i<=n;i++){
				dp[i]=my_pow(m+1,i*(i-1)/2);
				for(int j=0;j<i-1;j++){
					add2(dp[i],-1ll*dp[j+1]*c[i-1][j]%MOD*my_pow(m,(j+1)*(i-j-1))%MOD*my_pow(m+1,(i-j-1)*(i-j-2)/2)%MOD);
				}
			}
			for(int i=1;i<=tot;i++) res=1ll*res*dp[num[i]]%MOD;
			for(int i=1;i<=tot;i++){
				for(int j=i+1;j<=tot;j++){
					if(f[i][j]<b[i][j]){
						res=0;continue;
					}
					bool ok=false;
					for(int k=1;k<=tot;k++) if(k!=i&&k!=j){
						if(f[i][k]+f[k][j]==b[i][j]){
							ok=true;break;
						}
					}
					if(m<b[i][j]){
						res=0;continue;
					}
					if(ok){
						res=1ll*res*my_pow(m-b[i][j]+1,num[i]*num[j])%MOD;
					}
					else{
						int cur=my_pow(m-b[i][j]+1,num[i]*num[j]);
						add2(cur,-my_pow(m-b[i][j],num[i]*num[j]));
						res=1ll*res*cur%MOD;
					}
				}
			}
			add2(ans,-res);
		}
		cout<<ans<<'\n';
	}
	return 0;
}