#include <iostream>
#include <algorithm>
#include <iomanip>
#include <vector>
#include <queue>
#include <set>
#include <map>

using namespace std;
typedef long long ll;

typedef pair<int, int> P;
int A[100][100];
int N, M;

const int INF =  100000000;

struct edge {int to, cap, rev;};


class FordFolkerson{
    public:
        FordFolkerson(int n){
            N = n;
            G = vector<vector<edge>>(n, vector<edge>());
            used = vector<bool>(n, false);
        }

        void add_edge(int from, int to, int cap){
            G[from].push_back((edge{to, cap, (int)G[to].size()}));
            G[to].push_back((edge{from, 0, (int)G[from].size()-1}));
        }

        int max_flow(int s, int t){
            int flow = 0;
            while(true){
                clear_used();
                int f = dfs(s, t, INF);
                if(f == 0){
                    break;
                }
                flow += f;
            }
            return flow;
        }
        int match(int v){
            for(edge e : G[v]){
                if(e.to == 0) continue;
                if(e.cap == 0) return e.to;
            }
        }

    private:
        vector<vector<edge>> G;
        vector<bool> used;
        int N;
        void clear_used(){
            for(int i = 0; i < N; i++) used[i] = false;
        }

        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 to_row(int a){
    return (a-1)/M;
}

int B[100][100], C[100][100];
bool used[100][100];

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << setprecision(10) << fixed;
    cin >> N >> M;
    for(int i = 0; i < N; i++){
        for(int j = 0; j < M; j++){
            cin >> A[i][j];
        }
    }
    for(int t = 0; t < M-1; t++){
        FordFolkerson din(N*2+2);
        for(int i = 1; i <= N; i++) din.add_edge(0, i, 1);
        for(int i = 1; i <= N; i++) din.add_edge(N+i, 2*N+1, 1);
        for(int i = 0; i < N; i++){
            for(int j = 0; j < M; j++){
                if(!used[i][j]) {
                    int r = to_row(A[i][j]);
                    din.add_edge(r+1, N+i+1, 1);
                }
            }
        }

        int f = din.max_flow(0, 2*N+1);
        for(int i = 0; i < N; i++){
            int m = din.match(i+1)-N-1;
            for(int j = 0; j < M; j++){
                if(!used[m][j] &&  to_row(A[m][j]) == i) {
                    B[m][t] = A[m][j];
                    used[m][j] = true;
                    break;
                }
            }
        }
    }

    for(int i = 0; i < N; i++){
        for(int j = 0; j < M; j++){
            if(!used[i][j]) {
                B[i][M-1] = A[i][j];
                used[i][j] = true;
                break;
            }
        }
    }
    for(int i = 0; i < N; i++){
        for(int j = 0; j < M; j++){
            cout << B[i][j] << ' ';
        }
        cout << endl;
    }
    for(int i = 0; i < M; i++){
        vector<int> v;
        for(int j = 0; j < N; j++) v.push_back(B[j][i]);
        sort(v.begin(), v.end());
        for(int j = 0; j < N; j++) C[j][i] = v[j];
    }
    for(int i = 0; i < N; i++){
        for(int j = 0; j < M; j++){
            cout << C[i][j] << ' ';
        }
        cout << endl;
    }
}   

Measured memory_bytes: 820