#include<bits/stdc++.h>
using namespace std;
const int N=221,inf=0x3f3f3f3f;
struct edge{
	int y,inv,nxt,w;
	edge(int y1,int i1,int n1,int w1):y(y1),inv(i1),nxt(n1),w(w1){}
	edge(){}
}e[N*N];
int cnt,st,ed;
int n,m;
int head[3*N],cur[3*N],d[3*N];
int s[N][N];
int sz[N][3];
int vis[N];
void add(int x,int y,int i,int w){
	e[cnt]=edge(y,i,head[x],w);
	head[x]=cnt;
}
void build(int x,int y,int w){
	cnt++;add(x,y,cnt+1,w);
	cnt++;add(y,x,cnt-1,0);
}
int dinic(int x,int flow){
	if(!flow||x==ed) return flow;
	int res=0;
	for(int i=cur[x];~i;i=e[i].nxt,cur[x]=i){
		int y=e[i].y,w=e[i].w;
		if(!w||d[y]!=d[x]+1) continue;
		int tfl=dinic(y,min(flow,w));
		e[i].w-=tfl,e[e[i].inv].w+=tfl;
		flow-=tfl,res+=tfl;
		if(!flow) break;
	}
	return res;
}
bool bfs(){
	queue<int> q;
	memset(d,0x3f,sizeof(d));
	memset(vis,0,sizeof(vis));
	for(int i=st;i<=ed;i++) cur[i]=head[i];
	d[st]=1,q.push(st);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=head[x];~i;i=e[i].nxt){
			int y=e[i].y,w=e[i].w;
			if(!w||d[y]<inf) continue;
			d[y]=d[x]+1;
			q.push(y);
		}
	}
	return d[ed]!=inf;
}
signed main(){
	cin>>n>>m;
	for(int i=1;i<=n;i++) for(int j=1;j<=n-i;j++){
		scanf("%d",&s[i][i+j]);
		s[i+j][i]=s[i][i+j];
	}
	if((m+1)&1){
		printf("No");
		return 0;
	}
	for(int i=1;i<=n;i++){
		memset(vis,0,sizeof(vis));
		for(int j=1;j<=n;j++){
			if(!s[i][j]) continue;
			if(vis[s[i][j]]){
				printf("No");
				return 0;
			}
			vis[s[i][j]]++;
		}
		for(int j=1;j<=m;j++){
			if(vis[j]) sz[j][2]++;
			else sz[j][1]++;
		}
	}
	for(int i=1;i<=m;i++){
		if(sz[i][2]/2+sz[i][1]>(m+1)/2){
			printf("No");
			return 0;
		}
	}
	printf("Yes\n");
	for(int x=n+1;x<=m+1;x++){
		memset(head,-1,sizeof(head));
		memset(sz,0,sizeof(sz));
		cnt=0,st=0,ed=n+m+2;
		for(int i=1;i<=n;i++){
			memset(vis,0,sizeof(vis));
			for(int j=1;j<=n;j++){
				if(!s[i][j]) continue;
				if(vis[s[i][j]]){
					printf("No");
					return 0;
				}
				vis[s[i][j]]++;
			}
			for(int j=1;j<=m;j++){
				if(vis[j]) sz[j][2]++;
				else sz[j][1]++;
			}
		}
		for(int i=1;i<=m;i++){
			build(0,i,1);
			for(int j=1;j<=2*((m+1)/2-sz[i][1]-sz[i][2]/2);j++){
				build(i,m+n+1,1);
			}
		}
		for(int i=1;i<=n;i++) build(m+i,ed,1);
		build(m+n+1,ed,m-n);
		for(int i=1;i<=n;i++){
			memset(vis,0,sizeof(vis));
			for(int j=1;j<=n;j++){
				if(!s[i][j]) continue;
				vis[s[i][j]]++;
			}
			for(int j=1;j<=m;j++){
				if(!vis[j]) build(j,m+i,1);
			}
		}
		int ans=0;
		while(bfs()){
			ans+=dinic(st,inf);
		}
		//printf("%lld\n",ans);
		for(int t=1;t<=m;t++){
			for(int i=head[t];~i;i=e[i].nxt){
				if(e[i].y!=st&&e[i].w==0){
					if(e[i].y<=n+m) s[e[i].y-m][x]=s[x][e[i].y-m]=t;
					break;
				}
			}
		}
		n++;
	}
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			printf("%d ",s[i][j]);
		}
		printf("\n");
	}
	
}