#include<bits/stdc++.h>
using namespace std;
const int N=4e5+10;
int T,n,m,k;
long long mx[N];
vector<int>b[N];
struct Layn {
    int s,cnt,id;
    long long sum;
    vector<int>p;
}a[N];
bool cmp1(Layn x,Layn y){return x.s<y.s;}
bool cmp2(Layn x,Layn y){return x.id<y.id;}
int main() {
    scanf("%d",&T);
    while(T--) {
        scanf("%d%d%d",&n,&m,&k);
        for(int i=1;i<=n;i++)a[i].p.clear(),a[i].p.push_back(0),b[i].clear();
        for(int i=1;i<=n;i++) {
            scanf("%d",&a[i].s),a[i].cnt=0,a[i].sum=0;
            for(int j=1,x;j<=m;j++) {
                scanf("%d",&x),a[i].p.push_back(x);
                if(x<0)a[i].cnt++;
                else a[i].sum+=x;
                a[i].id=i;
            }
        }
        sort(a+1,a+n+1,cmp1);
        int siz=0;
        for(int i=1;i<=n;i++) {
            if(a[i].s>a[i-1].s)++siz;
            b[siz].push_back(i);
        }bool flag=0;
        for(int i=1;i<=siz;i++) {
            mx[i]=0;
            for(int j:b[i]) {
                //printf("%d %d %d %d----\n",a[j].s,a[j].cnt,a[j].sum,mx[i-1]);
                if(a[j].sum+1ll*k*a[j].cnt<=mx[i-1]){flag=1;break;}
                for(int l=1;l<=m;l++)if(a[j].p[l]<0) {
                    if(a[j].sum<=mx[i-1])a[j].p[l]=min(mx[i-1]+1-a[j].sum,1ll*k);
                    else a[j].p[l]=0;
                    a[j].sum+=a[j].p[l];
                }
                mx[i]=max(mx[i],a[j].sum);
            }
            if(flag)break;
        }
        sort(a+1,a+n+1,cmp2);
        if(flag)puts("No");
        else {
            puts("Yes");
            for(int i=1;i<=n;i++,puts(""))for(int j=1;j<=m;j++)printf("%d ",a[i].p[j]);
        }
    }
}
// 5
// 3 4 5
// 5 1 3 -1 -1
// 2 -1 5 -1 5
// 3 3 -1 -1 4
// 2 3 10
// 10000 5 0 -1
// 1 10 10 10
// 2 3 10
// 10 1 2 3
// 100 4 5 6
// 2 3 10
// 100 1 2 3
// 10 4 5 6
// 2 3 10000
// 100 -1 -1 -1
// 1 -1 -1 -1