#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define N 810
ll n,m;
ll a[N];
vector<ll>t[N];
ll fir[N],v[N*N],nxt[N*N],w[N*N];ll cnt=0; 
inline void add(ll x,ll y,ll z){
	v[++cnt]=y;w[cnt]=z;nxt[cnt]=fir[x];fir[x]=cnt;
	return ;
}
bool vis[N];
ll xu[N*N],cn=0;
inline void dfs(ll x){
	vis[x]=1;
	for(;fir[x];){
		ll o=v[fir[x]];ll p=fir[x];
		fir[x]=nxt[fir[x]];dfs(o);
		xu[++cn]=p;
	}return ;
} 
ll x[N*N],y[N*N],cnn=0;
ll tt[N*N];ll ji=0;
ll jit[N];
int main()
{
//	freopen("test1.in","r",stdin);
	//freopen(".in","r",stdin);
//	freopen("test1.out","w",stdout);
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	cin>>n>>m;ll g=n*m+1;for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++)t[j].clear();
		for(int j=1;j<=m;j++){
			cin>>a[j];tt[++ji]=a[j];t[a[j]].push_back(j);
		}for(int j=1;j<=m;j++){
			if(t[j].empty()){
				add(j+n,i,0);continue;
			}for(int k=1;k<t[j].size();k++)add(i,j+n,t[j][k]+(i-1)*m);
		}
	}for(int i=1;i<=n;i++){
		if(vis[i])continue;
		cn=0;dfs(i);if(!cn)continue;
		reverse(xu+1,xu+cn+1);
		assert(v[xu[cn]]==i);
		reverse(xu+1,xu+cn);
		x[++cnn]=w[xu[1]];y[cnn]=g;
		for(int i=3;i<=cn;i+=2){
			x[++cnn]=w[xu[i-2]];y[cnn]=w[xu[i]];
		}x[++cnn]=g;y[cnn]=w[xu[cn-1]];
	}cout<<cnn<<'\n';
	for(int i=1;i<=cnn;i++){
		cout<<x[i]<<' '<<y[i]<<'\n';swap(tt[x[i]],tt[y[i]]);
	}for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			jit[tt[(i-1)*m+j]]++;
		}for(int j=1;j<=m;j++){
			if(jit[j]!=i){
				cout<<1;
			}
			assert(jit[j]==i);
		}
	}
	return 0;
}