#include<bits/stdc++.h>
using namespace std;
constexpr int M=404;
int n,m,tp,stk[M*M];
vector<int>vec[M];
vector<pair<int,int>>adj[M<<1],ans;
int read(){
    int x=0;char ch=getchar();
    while (!isdigit(ch)) ch=getchar();
    while (isdigit(ch)) x=x*10+ch-48,ch=getchar();
    return x;
}
void dfs(int x){
    while (!adj[x].empty()){
        auto [y,z]=adj[x].back();adj[x].pop_back();
        dfs(y); if (z) stk[++tp]=z;
    }
}
int main(){
    n=read();m=read();
    for (int i=1;i<=n;i++){
        for (int j=1;j<=m;j++) vec[j].clear();
        for (int j=1;j<=m;j++) vec[read()].emplace_back((i-1)*m+j);
        for (int j=1;j<=m;j++)
            if (vec[j].empty())
                adj[i].emplace_back(n+j,0);
            else for (int k=1;k<vec[j].size();k++)
                adj[n+j].emplace_back(i,vec[j][k]);
    }
    for (int i=1;i<=n;i++)
        if (!adj[i].empty()){
            tp=0;dfs(i);
            reverse(stk+1,stk+tp+1);
            ans.emplace_back(n*m+1,stk[1]);
            for (int i=1;i<tp;i++)
                ans.emplace_back(stk[i],stk[i+1]);
            ans.emplace_back(stk[tp],n*m+1);
        }
    printf("%d\n",ans.size());
    for (auto [x,y]:ans) printf("%d %d\n",x,y);
    return 0;
}