#include<bits/stdc++.h>
using namespace std;
typedef pair<int,int> pii;
int main(){
  ios::sync_with_stdio(false);
  cin.tie(0); cout.tie(0);
  int n,m; cin>>n>>m;
  vector<vector<int> > a(n+1);
  for(int i=0;i<n;i++){
    a[i].resize(m);
    for(auto &j:a[i])cin>>j,j--;
  }
  vector<pii> r;
  auto op=[&](int x,int y){
    r.emplace_back(x,y);
    a[y].emplace_back(a[x].back()),a[x].pop_back();
  };
  function<void(int,int)> dfs=[&](int l,int r){
    if(l==r)return;
    int d=l+r>>1;
    for(int i=0,p=-1;i<n;i++){
      int mn=*min_element(a[i].begin(),a[i].end()),
        mx=*max_element(a[i].begin(),a[i].end());
      if(l<=mn&&mx<=r){
        if(~p){
          int s=0,c=-1;
          for(int j:a[p])s+=j>d;
          for(int j:a[i])s+=j>d;
          c=s>m,s=0;
          auto pd=[&](int x){return c^(x<=d);};
          for(int j:a[p])s+=pd(j);
          while(a[i].size()>m-s)op(i,n);
          while(!a[p].empty())pd(a[p].back())?op(p,i):op(p,n);
          while(a[i].size()>m-s)op(i,p);
          while(a[n].size()>s)op(n,p);
          while(!a[i].empty())op(i,n);
          while(a[p].size()>s)op(p,i);
          while(!a[n].empty())pd(a[n].back())&&a[p].size()<m?op(n,p):op(n,i);
        }
        p=i;
      }
    }
    dfs(l,d),dfs(d+1,r);
  };
  dfs(0,n-1);
  cout<<r.size()<<endl;
  for(auto [x,y]:r)
    cout<<x+1<<' '<<y+1<<'\n';
  return 0;
}