#include <bits/stdc++.h>
#pragma GCC optimize ("O2,unroll-loops")
//#pragma GCC optimize("no-stack-protector,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> pii;
typedef pair<pii, int> piii;
typedef pair<ll, ll> pll;
#define debug(x) cerr<<#x<<'='<<(x)<<endl;
#define debugp(x) cerr<<#x<<"= {"<<(x.first)<<", "<<(x.second)<<"}"<<endl;
#define debug2(x, y) cerr<<"{"<<#x<<", "<<#y<<"} = {"<<(x)<<", "<<(y)<<"}"<<endl;
#define debugv(v) {cerr<<#v<<" : ";for (auto x:v) cerr<<x<<' ';cerr<<endl;}
#define all(x) x.begin(), x.end()
#define pb push_back
#define kill(x) return cout<<x<<'\n', 0;

const ld eps=1e-7;
const int inf=1000000010;
const ll INF=10000000000000010LL;
const int mod=1000000007;
const int N=105;

int n, m, k, u, v, x, y, t, a, b, ans;
int A[N][N];
int B[N][N];
int C[N][N];
int match[N], tmp[N];
bool mark[N], ok[N];
bool G[N][N];

bool dfs(int node){
	mark[node]=1;
	for (int v=1; v<=n; v++) if (G[node][v]){
		if (!match[v] || (!mark[match[v]] && dfs(match[v]))){
			match[v]=node;
			return 1;
		}
	}
	return 0;
}

int main(){
	ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	//freopen("input.txt", "r", stdin);
	//freopen("output.txt", "w", stdout);
	cin>>n>>m;
	for (int i=1; i<=n; i++) for (int j=1; j<=m; j++) cin>>A[i][j];
	for (int j=1; j<=m; j++){
		memset(G, 0, sizeof(G));
		for (int x=1; x<=n; x++) for (int y=1; y<=m; y++) if (A[x][y]) G[x][(A[x][y]+m-1)/m]=1;
		memset(match, 0, sizeof(match));
		memset(ok, 0, sizeof(ok));
//		debug("start")
		int flag=1;
		while (flag--){
			memset(mark, 0, sizeof(mark));
			for (int i=1; i<=n; i++) if (!mark[i] && !ok[i] && dfs(i)) flag=ok[i]=1;
		}
		for (int i=1; i<=n; i++){
			assert(match[i]);
//			debug2(i, match[i])
			int x=match[i];
			for (int y=1; y<=m; y++) if ((A[x][y]+m-1)/m==i){
				B[x][j]=A[x][y];
				A[x][y]=0;
				break ;
			}
		}
	}
	for (int j=1; j<=m; j++){
		for (int i=1; i<=n; i++) tmp[i]=B[i][j];
		sort(tmp+1, tmp+n+1);
		for (int i=1; i<=n; i++) C[i][j]=tmp[i];
	}
	for (int i=1; i<=n; i++) for (int j=1; j<=m; j++) cout<<B[i][j]<<" \n"[j==m];
	for (int i=1; i<=n; i++) for (int j=1; j<=m; j++) cout<<C[i][j]<<" \n"[j==m];
	
	return 0;
}
