#include<iostream>
#include<cstdio>
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
using namespace std;
#include<vector>
#include<cmath>
#define int ll
typedef long long ll;
int T;
int find_closest_k(ll n){
	ll k=sqrt(2*n)-1;
	while(1){
		if((ll)k*(k-1)/2<=n&&(ll)k*(k+1)/2>n){
			return k;
		}
		k++;
	}
}

void solve(ll n,ll m){
	if(((m*(m+1)/2)%2==1)&&((n*(n+1)/2)%2==1)){
		cout<<"No"<<endl;
	}
	else{
        if((m*m-3*m)%4==0){
            cout<<"Yes"<<endl;
		    vector<int> v;
			int res=m;
			m=(m*m-3*m)/4;
			while(m){
				int k=find_closest_k(m);
                //cout<<k<<endl;
				v.push_back(k);
				res-=k;
				m-=k*(k-1)/2;
			}
            int len=v.size();
            // for(int i=0;i<len;i++){
            //     cout<<v[i]<<endl;
            // }
            // cout<<res<<endl;
		for(int i=1;i<=n;i++){
			int len=v.size();
			bool flag=1;
			for(int j=0;j<len;j++){
				for(int k=1;k<=v[j];k++){
					cout<<(int)flag<<" ";
				}
				flag=!flag;
			}
            int temp=res;
			while(temp--){
				cout<<(int)flag<<" ";
				flag=!flag;
			}
            cout<<endl;
		}
        }
		else if((n*n-3*n)%4==0){
            cout<<"Yes"<<endl;
            vector<int> v;
            int res=n;
            n=(n*n-3*n)/4;
            while(n){
                int k=find_closest_k(n);
                v.push_back(k);
                res-=k;
                n-=k*(k-1)/2;
            }
            int len=v.size();
            // for(int i=0;i<len;i++){
            //     cout<<v[i]<<endl;
            // }
            bool flag=1;
            for(int i=0;i<len;i++){
                for(int k=1;k<=v[i];k++){
                    for(int j=1;j<=m;j++){
                        cout<<(int)flag<<" ";
                    }
                    cout<<endl;
                }
                flag=!flag;
            }
            while(res--){
                for(int i=1;i<=m;i++){
                    cout<<(int)flag<<" ";
                }
                cout<<endl;
                flag=!flag;
            }
        }
	}
}
signed main(){
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	cin>>T;
	while(T--){
		// ll n,m;
		// cin>>m;
        // vector<int> v;
		// 	ll res=m;
		// 	m=(m*m-3*m)/4;
		// 	while(m){
		// 		ll k=find_closest_k(m);
        //         //cout<<k<<endl;
		// 		v.push_back(k);
		// 		res-=k;
		// 		m-=k*(k-1)/2;
		// 	}
        //     int len=v.size();
        //     ll ans=0;
        //     for(int i=0;i<len;i++){
        //         //cout<<v[i]<<" ";
        //         ans+=v[i];
        //     }
        //     ans+=res;
        //     cout<<endl;
        //     cout<<ans<<endl;
        ll n,m;
        cin>>n>>m;
        solve(n,m);
	}
	return 0;
}