#include "bits/stdc++.h"
#ifdef DEBUG
#include "PrettyDebug.hpp"
#else
#define debug(...) [](auto...){}(__VA_ARGS__)
#define debuga(...) [](auto...){}(__VA_ARGS__)
#endif
#define For(i,a,b) for(int i=a;i<=b;i++)
#define Rev(i,a,b) for(int i=a;i>=b;i--)
#define Fin(file) freopen(file,"r",stdin)
#define Fout(file) freopen(file,"w",stdout)
using namespace std; typedef long long ll;
// N is the number of vertices, M is the number of edges you're about to add
// cst_type is what you think it is
template<int N,int M,typename cst_type=int>
class Dinic{
    int n,S,T,head[N],nxt[M*2],to[M*2],tot,q[N],h,t,dep[N],now[N]; cst_type cap[M*2],maxflow;
    bool bfs(){
        h=t=0; q[t++]=S; For(i,1,n) dep[i]=0; ; dep[S]=1; now[S]=head[S];
        while(h<t){
            int u=q[h++]; for(int e=head[u];e;e=nxt[e]){
                if(cap[e]==0) continue; ; int v=to[e];
                if(!dep[v]) { dep[v]=dep[u]+1; now[v]=head[v]; q[t++]=v; if(v==T) return true; }
            }
        } return false;
    }
    cst_type dfs(int u,cst_type flow){
        if(u==T) return flow; ; cst_type fw=0;
        for(int& e=now[u];e;e=nxt[e]){
            if(cap[e]==0) continue; ; int v=to[e]; if(dep[v]!=dep[u]+1) continue;
            cst_type res=dfs(v,min(flow,cap[e])); if(res==0) { dep[v]=0; continue; }
            cap[e]-=res; cap[e^1]+=res; fw+=res; flow-=res; if(flow==0) break;
        } return fw;
    }
public:
    void init(int _n) { n=_n; tot=1; } 
    void clear(int _n) { For(i,1,n) head[i]=0; ; maxflow=0; n=_n; tot=1; }
    int add_edge(int x,int y,cst_type z) { nxt[++tot]=head[x]; head[x]=tot; to[tot]=y; cap[tot]=z; nxt[++tot]=head[y]; head[y]=tot; to[tot]=x; cap[tot]=0; return tot-1; }
    cst_type solve(int _S,int _T,cst_type inf=numeric_limits<cst_type>::max()){
        if(!n||tot%2==0) { cerr<<"INIT YOUR DINIC!!!!!!!!!!!!"<<endl; exit(-1); }
        S=_S; T=_T; while(bfs()) { cst_type res=-1; while((res=dfs(S,inf))!=0) maxflow+=res; } ; return maxflow;
    }
    cst_type flowing(int e) { return cap[e^1]; }
    vector<int> get_cut() { bfs(); vector<int> res; for(int e=2;e<=tot;e+=2) if(dep[to[e^1]]&&!dep[to[e]]) res.push_back(e); ; return res; }
};
constexpr int N=210;
int n,m,a[N][N],vis[N][N],id[N][N]; Dinic<N*2,N*N*2> G;
signed main(){
    atexit([](){cerr<<"Time = "<<clock()<<" ms"<<endl;});
    cin>>n>>m; if(~m&1) return cout<<"No\n",0;
    For(i,1,n) For(j,i+1,n){
        cin>>a[i][j],a[j][i]=a[i][j];
        if(vis[i][a[i][j]]) return cout<<"No\n",0; else vis[i][a[i][j]]=1;
        if(vis[j][a[i][j]]) return cout<<"No\n",0; else vis[j][a[i][j]]=1;
    }
    while(n<m+1){
        debug(n);
        int S=m+n+2,T=m+n+3; G.clear(T);
        For(i,1,m) G.add_edge(S,i,1);
        For(i,1,n) G.add_edge(m+i,T,1);
        G.add_edge(m+n+1,T,m-n);
        For(i,1,m){
            int cc=0; For(j,1,n) if(!vis[j][i]) id[i][j]=G.add_edge(i,m+j,1); else cc++,id[i][j]=0;
            if(n-cc/2!=(m+1)/2) G.add_edge(i,m+n+1,1);
        }
        if(G.solve(S,T)!=m) return cout<<"No\n",0;
        For(i,1,m) For(j,1,n) if(!vis[j][i]&&G.flowing(id[i][j])) a[n+1][j]=a[j][n+1]=i,vis[n+1][i]=vis[j][i]=1;
        n++; For(i,1,n) debuga(a[i],1,n);
    }
    cout<<"Yes\n";
    For(i,1,m+1) For(j,i+1,m+1) cout<<a[i][j]<<" \n"[j==m+1];
    return 0;
}

// START TYPING IF YOU DON'T KNOW WHAT TO DO
// STOP TYPING IF YOU DON'T KNOW WHAT YOU'RE DOING
// CONTINUE, NON-STOPPING, FOR CHARLIEVINNIE

// Started Coding On: January 26 Fri, 16 : 07 : 13