#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx")
#pragma GCC optimize("unroll-loops")
//#pragma warning(disable : 4996)

#ifdef _MSC_VER
#include <intrin.h>
#define __builtin_popcount __popcnt
#define __builtin_popcountll __popcnt64
#endif

#include <limits.h>
#include <math.h>
#include <time.h>

#include <algorithm>
#include <array>
#include <bitset>
#include <cassert>
#include <complex>
#include <cstdio>
#include <cstring>
#include <functional>
#include <iomanip>
#include <iostream>
#include <iterator>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <stack>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <vector>

using namespace std;

#define REP(i, n) for (int i = 0; i < (n); ++i)
#define REPR(i, n) for (int i = n - 1; i >= 0; --i)
#define FOR(i, m, n) for (int i = m; i < n; ++i)
#define FORR(i, m, n) for (int i = m - 1; i >= n; --i)
#define SORT(v, n) sort(v, v + n);
#define VSORT(v) sort(v.begin(), v.end());
#define REVERSE(v, n) reverse(v, v + n);
#define VREVERSE(v) reverse(v.begin(), v.end())
#define ll long long
#define print(x) cout << (x) << '\n'
#define pe(x) cout << (x) << " "
#define DEBUG(x) cout << #x << ": " << x << endl
#define lb(v, n) lower_bound(v.begin(), v.end(), (n))
#define ub(v, n) upper_bound(v.begin(), v.end(), (n))
#define int long long
//#define double long double
#define all(x) (x).begin(), (x).end()
#define print_space(v) REP(i, v.size()) cout << v[i] << ((i == v.size() - 1) ? "\n" : " ")
template <typename T1, typename T2> inline void chmin(T1& a, T2 b) { if (a > b) a = b; }
template <typename T1, typename T2> inline void chmax(T1& a, T2 b) { if (a < b) a = b; }
typedef pair<int, int> pii;
typedef pair<long long, long long> pll;
std::random_device rd;
std::mt19937 mt(rd());
constexpr ll MOD = 1e9 + 7;
constexpr int MAX = 500050;
const double pi = acos(-1);
constexpr double EPS = 1e-8;
constexpr ll LINF = 1e18 + 1;
constexpr int INF = 1e18 + 1;

//Ford-Fulkerson法で最大流を求める O(F|E|)
const int MAX_V = 500;
struct edge { int to, cap, rev,num; };
vector<edge>G[MAX_V];
bool used[MAX_V];

void add_edge(int from, int to, int cap,int num) {
	G[from].push_back({ to, cap, G[to].size(),num });
	G[to].push_back({ from,0,G[from].size() - 1,num });
}
int dfs(int v, int t, int f) {
	if (v == t)return f;
	used[v] = true;
	for (int i = 0; i < G[v].size(); i++) {
		edge &e = G[v][i];
		if (!used[e.to] && e.cap > 0) {
			int d = dfs(e.to, t, min(f, e.cap));
			if (d > 0) {
				e.cap -= d;
				G[e.to][e.rev].cap += d;
				return d;
			}
		}
	}
	return 0;
}
int max_flow(int s, int t) {
	int flow = 0;
	for (;;) {
		memset(used, 0, sizeof(used));
		int f = dfs(s, t, INF);
		if (f == 0) {
			return flow;
		}
		flow += f;
	}
}

int N, M;
int A[101][101];

bool vis[101 * 101];

int B[101][101], C[101][101];

void solve() {
	cin >> N >> M;
	REP(i, N) {
		REP(j, M) {
			cin >> A[i][j];
			A[i][j]--;
		}
	}
	int V = N * 2 + 2;
	int s = N * 2, t = N * 2 + 1;
	REP(k, M) {
		REP(i, V) {
			used[i] = false;
			G[i].clear();
		}
		REP(i, N) {
			add_edge(s, i, 1, -1);
			add_edge(N + i, t, 1, -1);
		}
		REP(i, N) {
			REP(j, M) {
				if (!vis[A[i][j]]) {
					int id = A[i][j] / M;
					int v = N + i;
					add_edge(id, v, 1, A[i][j]);
				}
			}
		}
		int f = max_flow(s, t);
		REP(i, N) {
			for (auto e : G[i]) {
				if (e.cap == 0) {
					int num = e.num;
					int row = e.to - N;
					B[row][k] = num;
					vis[num] = true;
				}
			}
		}
	}
	REP(i, N) {
		REP(j, M) {
			cout << B[i][j] + 1 << " \n"[j == M - 1];
		}
	}
	REP(j, M) {
		vector<int>v;
		REP(i, N)v.push_back(B[i][j] + 1);
		VSORT(v);
		REP(i, N)C[i][j] = v[i];
	}
	REP(i, N) {
		REP(j, M) {
			cout << C[i][j] << " \n"[j == M - 1];
		}
	}
}


signed main() {
	cin.tie(0);
	ios::sync_with_stdio(false);
	//while (test());
	//int q;
	//cin >> q;
	//while (q--)
	solve();
}


Measured memory_bytes: 4604