#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int MAXN = 800 + 5;

vector<pii> g[MAXN];

vector<int> sta;
int cur[MAXN];

void dfs(int u)
{
	while(cur[u] < (int)g[u].size())
	{
		int i = cur[u]; ++cur[u];
		
		dfs(g[u][i].first);
		
		if(g[u][i].second > 0)
			sta.emplace_back(g[u][i].second);
	}
}

int main(void)
{
	int n,m;
	scanf("%d%d",&n,&m);
	
	for(int i=1; i<=n; ++i)
	{
		static vector<int> pos[MAXN];
		for(int j=1; j<=m; ++j)
			pos[j].clear();
		for(int j=1; j<=m; ++j)
		{
			int x;
			scanf("%d",&x);
			pos[x].emplace_back((i-1) * m + j);
		}
		
		for(int j=1; j<=m; ++j)
		{
			if((int)pos[j].size() == 1) continue;
			if((int)pos[j].size() > 1)
			{
				for(int k=1; k<(int)pos[j].size(); ++k)
					g[i].emplace_back(j + n, pos[j][k]);
			}
			else
			{
				g[j + n].emplace_back(i, -1);
			}
		}
	}
	
	vector<pii> ans;
	
	int lstpos = n * m + 1;
	auto OP = [&] (int i)
	{
		ans.emplace_back(i, lstpos);
		lstpos = i;
	};
	
	for(int rt=1; rt<=n; ++rt)
	{
		sta.clear();
		dfs(rt);
		if((int)sta.size() == 0) continue;
		
		for(auto t: sta)
			OP(t);
		OP(n * m + 1);
	}
	
	printf("%d\n",(int)ans.size());
	for(auto t: ans)
		printf("%d %d\n",t.first,t.second);
	return 0;
}