#include<bits/stdc++.h>
using namespace std;

#define rep(i,a,b) for(int i=(a),i##end=(b);i<=i##end;++i)
#define per(i,a,b) for(int i=(a),i##end=(b);i>=i##end;--i)

#define pb push_back

typedef vector<int>vi;
typedef pair<int,int>pii;

int n,m,A[405*405];
vi pos[405][405];

int ecnt,h[805];
struct edges{
  int nxt,to;
}E[405*405*2];
void addline(int u,int v){
  E[++ecnt]={h[u],v},h[u]=ecnt;
}
vi path;
bool vis[805];
void dfs(int u){
  vis[u]=1;
  for(int&i=h[u];i;){
    int v=E[i].to;
    i=E[i].nxt;
    dfs(v);
  }
  path.pb(u);
}

int main(){
  cin.tie(0);
  ios::sync_with_stdio(0);
  cin>>n>>m;
  rep(i,1,n*m)cin>>A[i];
  rep(i,1,n){
    vi vis(m+1,0);
    rep(j,1,m){
      int x=A[(i-1)*m+j];
      if(!vis[x]){
        vis[x]=1;
      }else{
        pos[i][x].pb((i-1)*m+j);
        addline(i,n+x);
      }
    }
    rep(x,1,m){
      if(!vis[x]){
        addline(n+x,i);
      }
    }
  }
  vector<pii>as;
  rep(s,1,n)if(!vis[s]){
    path.clear();
    dfs(s);
    reverse(path.begin(),path.end());
    vi cir;
    for(int i=1;i<(int)path.size();i+=2){
      int u=path[i-1];
      int x=path[i]-n;
//      printf("#buk %d, val %d\n",u,x);
      cir.pb(pos[u][x].back());
      pos[u][x].pop_back();
    }
    as.pb({n*m+1,cir.back()});
    per(i,cir.size()-1,1){
      as.pb({cir[i],cir[i-1]});
    }
    as.pb({n*m+1,cir[0]});
  }
  cout<<as.size()<<endl;
  for(auto[x,y]:as)cout<<x<<' '<<y<<endl;
  return 0;
}