#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
inline int read()
{
char c=getchar();int x=0;bool f=0;
for(;!isdigit(c);c=getchar())f^=!(c^45);
for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+(c^48);
return f?-x:x;
}
int mod;
struct modint{
unsigned int x;
modint(int o=0){x=o;}
modint &operator = (int o){return x=o,*this;}
modint &operator +=(modint o){return x=x+o.x>=mod?x+o.x-mod:x+o.x,*this;}
modint &operator -=(modint o){return x=x<o.x?x-o.x+mod:x-o.x,*this;}
modint &operator *=(modint o){return x=1ll*x*o.x%mod,*this;}
modint &operator ^=(int b){
modint a=*this,c=1;
for(;b;b>>=1,a*=a)if(b&1)c*=a;
return x=c.x,*this;
}
modint &operator /=(modint o){return *this *=o^=mod-2;}
friend modint operator +(modint a,modint b){return a+=b;}
friend modint operator -(modint a,modint b){return a-=b;}
friend modint operator *(modint a,modint b){return a*=b;}
friend modint operator /(modint a,modint b){return a/=b;}
friend modint operator ^(modint a,int b){return a^=b;}
friend bool operator ==(modint a,modint b){return a.x==b.x;}
friend bool operator !=(modint a,modint b){return a.x!=b.x;}
bool operator ! () {return !x;}
modint operator - () {return x?mod-x:0;}
bool operator <(const modint&b)const{return x<b.x;}
};
inline modint qpow(modint x,int y){return x^y;}
vector<modint> fac,ifac,iv;
inline void initC(int n)
{
if(iv.empty())fac=ifac=iv=vector<modint>(2,1);
int m=iv.size(); ++n;
if(m>=n)return;
iv.resize(n),fac.resize(n),ifac.resize(n);
For(i,m,n-1){
iv[i]=iv[mod%i]*(mod-mod/i);
fac[i]=fac[i-1]*i,ifac[i]=ifac[i-1]*iv[i];
}
}
inline modint C(int n,int m){
if(m<0||n<m)return 0;
return initC(n),fac[n]*ifac[m]*ifac[n-m];
}
inline modint sign(int n){return (n&1)?(mod-1):(1);}
#define fi first
#define se second
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pii;
typedef vector<int>vi;
#define maxn 1000005
#define inf 0x3f3f3f3f
/*
*/
int n,r;
vi o={2,8,9,15};
int g=3;
int lg[maxn],pw[maxn];
bool isp(int p){
For(i,2,p-1) if(p%i==0) return 0;
return 1;
}
void sols(){
if(n==2){
if(r==4) {
puts("Yes");
puts("1 2");
puts("3 4");
return;
}
puts("No");
return;
}
if(n==3){
if(r==0){
puts("Yes");
int idx=0;
For(i,1,3)For(j,1,3)cout<<(++idx)<<" \n"[j==3];
return;
}
}
puts("No");
return;
}
int mat[405][405]={
{16 ,14 ,1 ,11},
{5,15,4,2},
{8,6,9,3},
{13,7,12,10}
};
int a[405][405];
void work()
{
n=read();
mod=n*n+1;
r=read();
int p=n*n+1;
if(!isp(p) || n==2){
sols();
return;
}
bool qwq=0;
for(g=2;g<n*n;++g){
int t=1;
For(i,0,n*n) lg[i]=0;
lg[1]=0;
For(j,1,n*n-1){
t=t*g%mod;
lg[t]=j;
}
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
bool ok=1;
For(i,2,n*n) if(!lg[i]) ok=0;
if(ok){
pw[0]=1;
For(i,1,n*n) pw[i]=pw[i-1]*g%mod;
qwq=1; break;
}
}
cerr<<"qwq "<<qwq<<"\n";
cerr<<"g: "<<g<<"\n";
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
// For(i,0,n*n) cout<<pw[i]<<" "; cout<<"\n";
assert(qwq);
{
int idx=-1;
For(i,1,n)For(j,1,n){
if(i%2==1 && j%2==1){
a[i][j]=++idx;
a[i][j+1]=n*n-idx-1;
}
if(i%2==0 && j%2==0){
a[i][j]=++idx;
a[i][j-1]=n*n-idx-1;
}
}
r=lg[r];
if(r%4!=2){
puts("No");
return;
}
int add=(r-(n*n-2)+n*n*2)%(n*n);
add/=4;
For(i,1,n)For(j,1,n){
a[i][j]+=add;
a[i][j]%=(n*n);
}
puts("Yes");
For(i,1,n)For(j,1,n) cout<<pw[a[i][j]]<<" \n"[j==n];
}
}
signed main()
{
int T=read();
while(T--)work();
return 0;
}
/*#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
inline int read()
{
char c=getchar();int x=0;bool f=0;
for(;!isdigit(c);c=getchar())f^=!(c^45);
for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+(c^48);
return f?-x:x;
}
int mod;
struct modint{
unsigned int x;
modint(int o=0){x=o;}
modint &operator = (int o){return x=o,*this;}
modint &operator +=(modint o){return x=x+o.x>=mod?x+o.x-mod:x+o.x,*this;}
modint &operator -=(modint o){return x=x<o.x?x-o.x+mod:x-o.x,*this;}
modint &operator *=(modint o){return x=1ll*x*o.x%mod,*this;}
modint &operator ^=(int b){
modint a=*this,c=1;
for(;b;b>>=1,a*=a)if(b&1)c*=a;
return x=c.x,*this;
}
modint &operator /=(modint o){return *this *=o^=mod-2;}
friend modint operator +(modint a,modint b){return a+=b;}
friend modint operator -(modint a,modint b){return a-=b;}
friend modint operator *(modint a,modint b){return a*=b;}
friend modint operator /(modint a,modint b){return a/=b;}
friend modint operator ^(modint a,int b){return a^=b;}
friend bool operator ==(modint a,modint b){return a.x==b.x;}
friend bool operator !=(modint a,modint b){return a.x!=b.x;}
bool operator ! () {return !x;}
modint operator - () {return x?mod-x:0;}
bool operator <(const modint&b)const{return x<b.x;}
};
inline modint qpow(modint x,int y){return x^y;}
vector<modint> fac,ifac,iv;
inline void initC(int n)
{
if(iv.empty())fac=ifac=iv=vector<modint>(2,1);
int m=iv.size(); ++n;
if(m>=n)return;
iv.resize(n),fac.resize(n),ifac.resize(n);
For(i,m,n-1){
iv[i]=iv[mod%i]*(mod-mod/i);
fac[i]=fac[i-1]*i,ifac[i]=ifac[i-1]*iv[i];
}
}
inline modint C(int n,int m){
if(m<0||n<m)return 0;
return initC(n),fac[n]*ifac[m]*ifac[n-m];
}
inline modint sign(int n){return (n&1)?(mod-1):(1);}
#define fi first
#define se second
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pii;
typedef vector<int>vi;
#define maxn 1000005
#define inf 0x3f3f3f3f
/*
*/
int n,r;
vi o={2,8,9,15};
int g=3;
int lg[maxn],pw[maxn];
bool isp(int p){
For(i,2,p-1) if(p%i==0) return 0;
return 1;
}
void sols(){
if(n==2){
if(r==4) {
puts("Yes");
puts("1 2");
puts("3 4");
return;
}
puts("No");
return;
}
if(n==3){
if(r==0){
puts("Yes");
int idx=0;
For(i,1,3)For(j,1,3)cout<<(++idx)<<" \n"[j==3];
return;
}
}
puts("No");
return;
}
int mat[405][405]={
{16 ,14 ,1 ,11},
{5,15,4,2},
{8,6,9,3},
{13,7,12,10}
};
int a[405][405];
void work()
{
n=read();
mod=n*n+1;
r=read();
int p=n*n+1;
if(!isp(p) || n==2){
sols();
return;
}
bool qwq=0;
for(g=2;g<n*n;++g){
int t=1;
For(i,0,n*n) lg[i]=0;
lg[1]=0;
For(j,1,n*n-1){
t=t*g%mod;
lg[t]=j;
}
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
bool ok=1;
For(i,2,n*n) if(!lg[i]) ok=0;
if(ok){
pw[0]=1;
For(i,1,n*n) pw[i]=pw[i-1]*g%mod;
qwq=1; break;
}
}
cerr<<"qwq "<<qwq<<"\n";
cerr<<"g: "<<g<<"\n";
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
// For(i,0,n*n) cout<<pw[i]<<" "; cout<<"\n";
assert(qwq);
{
int idx=-1;
For(i,1,n)For(j,1,n){
if(i%2==1 && j%2==1){
a[i][j]=++idx;
a[i][j+1]=n*n-idx-1;
}
if(i%2==0 && j%2==0){
a[i][j]=++idx;
a[i][j-1]=n*n-idx-1;
}
}
r=lg[r];
if(r%4!=2){
puts("No");
return;
}
int add=(r-(n*n-2)+n*n*2)%(n*n);
add/=4;
For(i,1,n)For(j,1,n){
a[i][j]+=add;
a[i][j]%=(n*n);
}
puts("Yes");
For(i,1,n)For(j,1,n) cout<<pw[a[i][j]]<<" \n"[j==n];
}
}
signed main()
{
int T=read();
while(T--)work();
return 0;
}#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
inline int read()
{
char c=getchar();int x=0;bool f=0;
for(;!isdigit(c);c=getchar())f^=!(c^45);
for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+(c^48);
return f?-x:x;
}
int mod;
struct modint{
unsigned int x;
modint(int o=0){x=o;}
modint &operator = (int o){return x=o,*this;}
modint &operator +=(modint o){return x=x+o.x>=mod?x+o.x-mod:x+o.x,*this;}
modint &operator -=(modint o){return x=x<o.x?x-o.x+mod:x-o.x,*this;}
modint &operator *=(modint o){return x=1ll*x*o.x%mod,*this;}
modint &operator ^=(int b){
modint a=*this,c=1;
for(;b;b>>=1,a*=a)if(b&1)c*=a;
return x=c.x,*this;
}
modint &operator /=(modint o){return *this *=o^=mod-2;}
friend modint operator +(modint a,modint b){return a+=b;}
friend modint operator -(modint a,modint b){return a-=b;}
friend modint operator *(modint a,modint b){return a*=b;}
friend modint operator /(modint a,modint b){return a/=b;}
friend modint operator ^(modint a,int b){return a^=b;}
friend bool operator ==(modint a,modint b){return a.x==b.x;}
friend bool operator !=(modint a,modint b){return a.x!=b.x;}
bool operator ! () {return !x;}
modint operator - () {return x?mod-x:0;}
bool operator <(const modint&b)const{return x<b.x;}
};
inline modint qpow(modint x,int y){return x^y;}
vector<modint> fac,ifac,iv;
inline void initC(int n)
{
if(iv.empty())fac=ifac=iv=vector<modint>(2,1);
int m=iv.size(); ++n;
if(m>=n)return;
iv.resize(n),fac.resize(n),ifac.resize(n);
For(i,m,n-1){
iv[i]=iv[mod%i]*(mod-mod/i);
fac[i]=fac[i-1]*i,ifac[i]=ifac[i-1]*iv[i];
}
}
inline modint C(int n,int m){
if(m<0||n<m)return 0;
return initC(n),fac[n]*ifac[m]*ifac[n-m];
}
inline modint sign(int n){return (n&1)?(mod-1):(1);}
#define fi first
#define se second
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pii;
typedef vector<int>vi;
#define maxn 1000005
#define inf 0x3f3f3f3f
/*
*/
int n,r;
vi o={2,8,9,15};
int g=3;
int lg[maxn],pw[maxn];
bool isp(int p){
For(i,2,p-1) if(p%i==0) return 0;
return 1;
}
void sols(){
if(n==2){
if(r==4) {
puts("Yes");
puts("1 2");
puts("3 4");
return;
}
puts("No");
return;
}
if(n==3){
if(r==0){
puts("Yes");
int idx=0;
For(i,1,3)For(j,1,3)cout<<(++idx)<<" \n"[j==3];
return;
}
}
puts("No");
return;
}
int mat[405][405]={
{16 ,14 ,1 ,11},
{5,15,4,2},
{8,6,9,3},
{13,7,12,10}
};
int a[405][405];
void work()
{
n=read();
mod=n*n+1;
r=read();
int p=n*n+1;
if(!isp(p) || n==2){
sols();
return;
}
bool qwq=0;
for(g=2;g<n*n;++g){
int t=1;
For(i,0,n*n) lg[i]=0;
lg[1]=0;
For(j,1,n*n-1){
t=t*g%mod;
lg[t]=j;
}
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
bool ok=1;
For(i,2,n*n) if(!lg[i]) ok=0;
if(ok){
pw[0]=1;
For(i,1,n*n) pw[i]=pw[i-1]*g%mod;
qwq=1; break;
}
}
cerr<<"qwq "<<qwq<<"\n";
cerr<<"g: "<<g<<"\n";
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
// For(i,0,n*n) cout<<pw[i]<<" "; cout<<"\n";
assert(qwq);
{
int idx=-1;
For(i,1,n)For(j,1,n){
if(i%2==1 && j%2==1){
a[i][j]=++idx;
a[i][j+1]=n*n-idx-1;
}
if(i%2==0 && j%2==0){
a[i][j]=++idx;
a[i][j-1]=n*n-idx-1;
}
}
r=lg[r];
if(r%4!=2){
puts("No");
return;
}
int add=(r-(n*n-2)+n*n*2)%(n*n);
add/=4;
For(i,1,n)For(j,1,n){
a[i][j]+=add;
a[i][j]%=(n*n);
}
puts("Yes");
For(i,1,n)For(j,1,n) cout<<pw[a[i][j]]<<" \n"[j==n];
}
}
signed main()
{
int T=read();
while(T--)work();
return 0;
}#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
inline int read()
{
char c=getchar();int x=0;bool f=0;
for(;!isdigit(c);c=getchar())f^=!(c^45);
for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+(c^48);
return f?-x:x;
}
int mod;
struct modint{
unsigned int x;
modint(int o=0){x=o;}
modint &operator = (int o){return x=o,*this;}
modint &operator +=(modint o){return x=x+o.x>=mod?x+o.x-mod:x+o.x,*this;}
modint &operator -=(modint o){return x=x<o.x?x-o.x+mod:x-o.x,*this;}
modint &operator *=(modint o){return x=1ll*x*o.x%mod,*this;}
modint &operator ^=(int b){
modint a=*this,c=1;
for(;b;b>>=1,a*=a)if(b&1)c*=a;
return x=c.x,*this;
}
modint &operator /=(modint o){return *this *=o^=mod-2;}
friend modint operator +(modint a,modint b){return a+=b;}
friend modint operator -(modint a,modint b){return a-=b;}
friend modint operator *(modint a,modint b){return a*=b;}
friend modint operator /(modint a,modint b){return a/=b;}
friend modint operator ^(modint a,int b){return a^=b;}
friend bool operator ==(modint a,modint b){return a.x==b.x;}
friend bool operator !=(modint a,modint b){return a.x!=b.x;}
bool operator ! () {return !x;}
modint operator - () {return x?mod-x:0;}
bool operator <(const modint&b)const{return x<b.x;}
};
inline modint qpow(modint x,int y){return x^y;}
vector<modint> fac,ifac,iv;
inline void initC(int n)
{
if(iv.empty())fac=ifac=iv=vector<modint>(2,1);
int m=iv.size(); ++n;
if(m>=n)return;
iv.resize(n),fac.resize(n),ifac.resize(n);
For(i,m,n-1){
iv[i]=iv[mod%i]*(mod-mod/i);
fac[i]=fac[i-1]*i,ifac[i]=ifac[i-1]*iv[i];
}
}
inline modint C(int n,int m){
if(m<0||n<m)return 0;
return initC(n),fac[n]*ifac[m]*ifac[n-m];
}
inline modint sign(int n){return (n&1)?(mod-1):(1);}
#define fi first
#define se second
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pii;
typedef vector<int>vi;
#define maxn 1000005
#define inf 0x3f3f3f3f
/*
*/
int n,r;
vi o={2,8,9,15};
int g=3;
int lg[maxn],pw[maxn];
bool isp(int p){
For(i,2,p-1) if(p%i==0) return 0;
return 1;
}
void sols(){
if(n==2){
if(r==4) {
puts("Yes");
puts("1 2");
puts("3 4");
return;
}
puts("No");
return;
}
if(n==3){
if(r==0){
puts("Yes");
int idx=0;
For(i,1,3)For(j,1,3)cout<<(++idx)<<" \n"[j==3];
return;
}
}
puts("No");
return;
}
int mat[405][405]={
{16 ,14 ,1 ,11},
{5,15,4,2},
{8,6,9,3},
{13,7,12,10}
};
int a[405][405];
void work()
{
n=read();
mod=n*n+1;
r=read();
int p=n*n+1;
if(!isp(p) || n==2){
sols();
return;
}
bool qwq=0;
for(g=2;g<n*n;++g){
int t=1;
For(i,0,n*n) lg[i]=0;
lg[1]=0;
For(j,1,n*n-1){
t=t*g%mod;
lg[t]=j;
}
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
bool ok=1;
For(i,2,n*n) if(!lg[i]) ok=0;
if(ok){
pw[0]=1;
For(i,1,n*n) pw[i]=pw[i-1]*g%mod;
qwq=1; break;
}
}
cerr<<"qwq "<<qwq<<"\n";
cerr<<"g: "<<g<<"\n";
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
// For(i,0,n*n) cout<<pw[i]<<" "; cout<<"\n";
assert(qwq);
{
int idx=-1;
For(i,1,n)For(j,1,n){
if(i%2==1 && j%2==1){
a[i][j]=++idx;
a[i][j+1]=n*n-idx-1;
}
if(i%2==0 && j%2==0){
a[i][j]=++idx;
a[i][j-1]=n*n-idx-1;
}
}
r=lg[r];
if(r%4!=2){
puts("No");
return;
}
int add=(r-(n*n-2)+n*n*2)%(n*n);
add/=4;
For(i,1,n)For(j,1,n){
a[i][j]+=add;
a[i][j]%=(n*n);
}
puts("Yes");
For(i,1,n)For(j,1,n) cout<<pw[a[i][j]]<<" \n"[j==n];
}
}
signed main()
{
int T=read();
while(T--)work();
return 0;
}#include<bits/stdc++.h>
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
inline int read()
{
char c=getchar();int x=0;bool f=0;
for(;!isdigit(c);c=getchar())f^=!(c^45);
for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+(c^48);
return f?-x:x;
}
int mod;
struct modint{
unsigned int x;
modint(int o=0){x=o;}
modint &operator = (int o){return x=o,*this;}
modint &operator +=(modint o){return x=x+o.x>=mod?x+o.x-mod:x+o.x,*this;}
modint &operator -=(modint o){return x=x<o.x?x-o.x+mod:x-o.x,*this;}
modint &operator *=(modint o){return x=1ll*x*o.x%mod,*this;}
modint &operator ^=(int b){
modint a=*this,c=1;
for(;b;b>>=1,a*=a)if(b&1)c*=a;
return x=c.x,*this;
}
modint &operator /=(modint o){return *this *=o^=mod-2;}
friend modint operator +(modint a,modint b){return a+=b;}
friend modint operator -(modint a,modint b){return a-=b;}
friend modint operator *(modint a,modint b){return a*=b;}
friend modint operator /(modint a,modint b){return a/=b;}
friend modint operator ^(modint a,int b){return a^=b;}
friend bool operator ==(modint a,modint b){return a.x==b.x;}
friend bool operator !=(modint a,modint b){return a.x!=b.x;}
bool operator ! () {return !x;}
modint operator - () {return x?mod-x:0;}
bool operator <(const modint&b)const{return x<b.x;}
};
inline modint qpow(modint x,int y){return x^y;}
vector<modint> fac,ifac,iv;
inline void initC(int n)
{
if(iv.empty())fac=ifac=iv=vector<modint>(2,1);
int m=iv.size(); ++n;
if(m>=n)return;
iv.resize(n),fac.resize(n),ifac.resize(n);
For(i,m,n-1){
iv[i]=iv[mod%i]*(mod-mod/i);
fac[i]=fac[i-1]*i,ifac[i]=ifac[i-1]*iv[i];
}
}
inline modint C(int n,int m){
if(m<0||n<m)return 0;
return initC(n),fac[n]*ifac[m]*ifac[n-m];
}
inline modint sign(int n){return (n&1)?(mod-1):(1);}
#define fi first
#define se second
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pii;
typedef vector<int>vi;
#define maxn 1000005
#define inf 0x3f3f3f3f
/*
*/
int n,r;
vi o={2,8,9,15};
int g=3;
int lg[maxn],pw[maxn];
bool isp(int p){
For(i,2,p-1) if(p%i==0) return 0;
return 1;
}
void sols(){
if(n==2){
if(r==4) {
puts("Yes");
puts("1 2");
puts("3 4");
return;
}
puts("No");
return;
}
if(n==3){
if(r==0){
puts("Yes");
int idx=0;
For(i,1,3)For(j,1,3)cout<<(++idx)<<" \n"[j==3];
return;
}
}
puts("No");
return;
}
int mat[405][405]={
{16 ,14 ,1 ,11},
{5,15,4,2},
{8,6,9,3},
{13,7,12,10}
};
int a[405][405];
void work()
{
n=read();
mod=n*n+1;
r=read();
int p=n*n+1;
if(!isp(p) || n==2){
sols();
return;
}
bool qwq=0;
for(g=2;g<n*n;++g){
int t=1;
For(i,0,n*n) lg[i]=0;
lg[1]=0;
For(j,1,n*n-1){
t=t*g%mod;
lg[t]=j;
}
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
bool ok=1;
For(i,2,n*n) if(!lg[i]) ok=0;
if(ok){
pw[0]=1;
For(i,1,n*n) pw[i]=pw[i-1]*g%mod;
qwq=1; break;
}
}
cerr<<"qwq "<<qwq<<"\n";
cerr<<"g: "<<g<<"\n";
// For(i,1,n*n) cout<<lg[i]<<" "; cout<<"\n";
// For(i,0,n*n) cout<<pw[i]<<" "; cout<<"\n";
assert(qwq);
{
int idx=-1;
For(i,1,n)For(j,1,n){
if(i%2==1 && j%2==1){
a[i][j]=++idx;
a[i][j+1]=n*n-idx-1;
}
if(i%2==0 && j%2==0){
a[i][j]=++idx;
a[i][j-1]=n*n-idx-1;
}
}
r=lg[r];
if(r%4!=2){
puts("No");
return;
}
int add=(r-(n*n-2)+n*n*2)%(n*n);
add/=4;
For(i,1,n)For(j,1,n){
a[i][j]+=add;
a[i][j]%=(n*n);
}
puts("Yes");
For(i,1,n)For(j,1,n) cout<<pw[a[i][j]]<<" \n"[j==n];
}
}
signed main()
{
int T=read();
while(T--)work();
return 0;
}*/