//Author: Kevin
#include<bits/stdc++.h>
//#pragma GCC optimize("O2")
using namespace std;
#define ll long long
#define ull unsigned ll
#define pb emplace_back
#define mp make_pair
#define ALL(x) (x).begin(),(x).end()
#define rALL(x) (x).rbegin(),(x).rend()
#define srt(x) sort(ALL(x))
#define rev(x) reverse(ALL(x))
#define rsrt(x) sort(rALL(x))
#define sz(x) (int)(x.size())
#define inf 0x3f3f3f3f
#define pii pair<int,int>
#define lb(v,x) (int)(lower_bound(ALL(v),x)-v.begin())
#define ub(v,x) (int)(upper_bound(ALL(v),x)-v.begin())
#define uni(v) v.resize(unique(ALL(v))-v.begin())
#define longer __int128_t
void die(string S){puts(S.c_str());exit(0);}
int isp[5005];
int val[55][55];
map<int,vector<int>> Mp;
int ksm(int a,int b,int mod)
{
	int ans=1;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		b>>=1;
		a=a*a%mod;
	}
	return ans;
}
int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	for(int i=2;i<5005;i++)
		isp[i]=1;
	for(int i=2;i<5005;i++)
		for(int j=i+i;j<5005;j+=i)
			isp[j]=0;
	int t;
	cin>>t;
	while(t--)
	{
		int n,r;
		cin>>n>>r;
		if(n==3)
		{
			if(r)
			{
				cout<<"No\n";
				continue;
			}
			cout<<"Yes\n";
			cout<<"1 2 3\n";
			cout<<"4 5 6\n";
			cout<<"7 8 9\n";
			continue;
		}
		else if(isp[n*n+1]&&r)
		{
			int g=2;
			while(true)
			{
				int ok=1;
				vector<int> vec;
				int tmp=1;
				for(int i=1;i<=n*n;i++)
				{
					tmp=tmp*g%(n*n+1);
					vec.pb(tmp);
				}
				srt(vec);
				uni(vec);
				if(sz(vec)==n*n) break;
				g++;
			}
			int tmp=1;
			int v=0;
			while(true)
			{
				if(tmp==r)
				{
					r=v;
					break;
				}
				tmp=tmp*g%(n*n+1);
				v++;
			}
			if(v%4!=2)
			{
				cout<<"No\n";
				continue;
			}
			for(int i=0;i<n;i++)
				val[i][0]=i*2;
			for(int i=0;i<n;i++)
				val[0][i]=(i/2)*(n*2)+(i&1);
			int sum=6;
			for(int i=1;i<n;i++)
				for(int j=1;j<n;j++)
					val[i][j]=((sum-val[i-1][j]-val[i][j-1]-val[i-1][j-1])%(n*n)+n*n)%(n*n);
			int x=((v-sum)/4+n*n)%(n*n);
			for(int i=0;i<n;i++)
				for(int j=0;j<n;j++)
					val[i][j]+=x;
			cout<<"Yes\n";
			for(int i=0;i<n;i++)
			{
				for(int j=0;j<n;j++)
				{
					cout<<ksm(g,val[i][j],n*n+1)<<" ";
				}
				cout<<'\n';
			}
		}
		else
			cout<<"No\n";
	}
	return 0;
}