#include<bits/stdc++.h>
using namespace std;
inline int read()
{
	int n=0,f=1,ch=getchar();
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		n=n*10+ch-'0';
		ch=getchar();
	}
	return n*f;
}
int p[200005];
int sl[200005];
bool cz[200005],vis[200005];
vector<int>v[200005];
vector<int>yy[200005],qs[200005];
int fa[200005];
int findf(int n)
{
	if(fa[n]==n)return n;
	return fa[n]=findf(fa[n]);
}
int xl1[200005],cnt1;
int xl2[200005],cnt2;
bool visd[200005];
int n,m;
vector<int>wz[200005];
int now[1000005];
int pp[1000005];
int sta[1000005],ttop;
void dfs1(int x)
{
	for(int i=now[x];i<v[x].size();i=now[x])
	{
		now[x]=i+1;
		dfs1(v[x][i]);
	}
	sta[++ttop]=x;
}
int tsl[405][405];
int sc1[1000005],sc2[1000005];
int xl[1000005],cnt;
int main()
{
	n=read();
	m=read();
	for(int i=1;i<=n*m;i++)
	{
		p[i]=read();
		pp[i]=p[i];
	}
	for(int i=1;i<=m;i++)fa[i]=i;
	int ans=0;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)sl[j]=0;
		for(int j=1;j<=m;j++)
		{
			sl[p[(i-1)*m+j]]++;
			if(sl[p[(i-1)*m+j]]>=2)
			{
				cz[p[(i-1)*m+j]]=true;
				qs[i].push_back(p[(i-1)*m+j]);
				wz[p[(i-1)*m+j]].push_back((i-1)*m+j);
			}
		} 
		for(int j=1;j<=m;j++)if(sl[j]==0)yy[i].push_back(j);
		for(int j=0;j<yy[i].size();j++)
		{
			v[yy[i][j]].push_back(i+m);
			v[i+m].push_back(qs[i][j]);
		}
		if(!yy[i].empty())
		{
			for(int j=1;j<yy[i].size();j++)fa[findf(yy[i][0])]=findf(yy[i][j]);
			for(int j=0;j<qs[i].size();j++)fa[findf(yy[i][0])]=findf(qs[i][j]);
		}
	}
	for(int i=1;i<=m;i++)
	{
		if(!cz[i])continue;
		if(findf(i)==i)
		{
			v[m+n+1].push_back(i);
			v[i].push_back(m+n+1);
		}
	}
	dfs1(m+n+1);
	reverse(sta+1,sta+ttop+1);
	cnt1=0;
	for(int i=ttop;i>=1;i--)
	{
		if(sta[i]==n+m+1)
		{
			xl[++cnt]=n*m+1;
		}
		else if(sta[i]>m)
		{
			continue;
		}
		else
		{
			if(sta[i-1]==n+m+1)continue;
			int bh=sta[i-1]-m;
			for(int j=1;j<=m;j++)
			{
				if(p[(bh-1)*m+j]==sta[i])
				{
					xl[++cnt]=(bh-1)*m+j;
					p[(bh-1)*m+j]=-1;
					break;
				}
			}
		}
	}
	printf("%d\n",cnt-1);
	for(int i=1;i<=cnt-1;i++)printf("%d %d\n",xl[i+1],xl[i]);
	return 0;
	for(int i=1;i<=n*m;i++)p[i]=pp[i];
	for(int i=1;i<=cnt-1;i++)
	{
		int now1=xl[i+1],now2=xl[i];
		assert(p[now2]==0);
		p[now2]=p[now1];
		p[now1]=0;
	}
	for(int i=1;i<=n*m;i++)printf("%d ",p[i]);
	printf("\n");
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)sl[j]=0;
		for(int j=1;j<=m;j++)sl[p[(i-1)*m+j]]++;
		for(int j=1;j<=m;j++)assert(sl[j]==1);
	}
	return 0;
}