#include<bits/stdc++.h>

#define pb push_back
#define fi first
#define se second
#define mp make_pair

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef long double ld;

template <typename T> bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template <typename T> bool chkmax(T &x,T y){return x<y?x=y,1:0;}

int readint(){
	int x=0,f=1; char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

int n,m,tot,cnt;
int v[160005],c[160005],nxt[160005],h[405],ans[160005],a[160005],num[405][405],f[405],fr[405][405];
bool vis[160005];
vector<int> hv[405][405];
vector<pii> res;

void addedge(int x,int y,int z){
	v[++tot]=y; c[tot]=z; nxt[tot]=h[x]; h[x]=tot;
}

int getf(int x){return x==f[x]?x:f[x]=getf(f[x]);}

void merge(int x,int y){
	int fx=getf(x),fy=getf(y);
	if(fx==fy) return;
	f[fx]=fy;
}

void dfs(int u){
	for(int &p=h[u];p;p=nxt[p]){
		if(vis[p]) continue;
		vis[p]=1;
		int edge=p;
		dfs(v[p]);
		ans[++cnt]=edge;
	}
}

int main(){
	n=readint(); m=readint();
	for(int i=1;i<=n*m;i++) a[i]=readint();
	for(int i=1;i<=n;i++){
		for(int j=(i-1)*m+1;j<=i*m;j++) num[i][a[j]]++;
		for(int j=(i-1)*m+1;j<=i*m;j++) hv[i][a[j]].pb(j);
	}
	for(int i=1;i<=m;i++){
		vector<int> vec;
		for(int j=1;j<=n;j++) if(!num[j][i]) vec.pb(j);
		for(int j=1;j<=n;j++){
			for(int k=2;k<=num[j][i];k++){
				addedge(j,vec.back(),i);
				fr[i][vec.back()]=tot;
				vec.pop_back();
			}
		}
	}
	for(int i=1;i<=n;i++) f[i]=i;
	for(int i=1;i<=n;i++) for(int p=h[i];p;p=nxt[p]) merge(i,v[p]);
	for(int i=1;i<=m;i++){
		vector<int> vec;
		for(int j=1;j<=n;j++) if(!num[j][i]) vec.pb(j);
		for(int j=1;j<vec.size();j++){
			if(getf(vec[j])!=getf(vec[0])){
				int t1=fr[i][vec[0]],t2=fr[i][vec[j]];
				swap(v[t1],v[t2]);
				merge(vec[0],vec[j]);
			}
		}
	}
	for(int i=1;i<=n;i++){
		cnt=0;
		dfs(i);
		if(!cnt) continue;
		// cout<<"##### "<<endl;
		// for(int j=cnt;j>=1;j--) cout<<v[ans[j]]<<' ';
		// cout<<endl;
		int lst=hv[i][c[ans[cnt]]].back();
		hv[i][c[ans[cnt]]].pop_back();
		res.pb(mp(lst,n*m+1));
		for(int i=1;i<cnt;i++){
			int tmp=hv[v[ans[i+1]]][c[ans[i]]].back();
			hv[v[ans[i+1]]][c[ans[i]]].pop_back();
			res.pb(mp(tmp,lst));
			lst=tmp;
		}
		res.pb(mp(n*m+1,lst));
	}
	printf("%d\n",res.size());
	for(auto r:res) printf("%d %d\n",r.fi,r.se);
	return 0;
}