#include <bits/stdc++.h>
using namespace std;
#define a first
#define b second
#define pb push_back
typedef long long llo;
//#define endl '\n'
int n,m;
int it[401][401];
vector<int> pre[401][401];
int co[401][401];
int tt[401*401];
vector<pair<int,int>> adj[801];
int vis[801];
vector<pair<int,int>> xx[801];
vector<pair<int,int>> cot2;
void dfs(int no){
	vector<pair<int,int>> zz=xx[no];
	xx[no].clear();
	for(auto j:zz){
		if(xx[j.a].size()){
			dfs(j.a);
		}
		cot2.pb(j);
	}
}
int main(){
	ios_base::sync_with_stdio(false);
	cin.tie(NULL);
	cin>>n>>m;
	for(int i=0;i<n;i++){
		for(int j=0;j<m;j++){
			cin>>it[i][j];
			tt[i*m+j]=it[i][j];
			it[i][j]--;
			pre[i][it[i][j]].pb(j);
			co[i][it[i][j]]++;
		}
		for(int j=0;j<m;j++){
			co[i][j]-=1;
			if(co[i][j]==-1){
				adj[j+n].pb({i,-1});
			}
			else if(co[i][j]>0){
				for(int k=0;k+1<pre[i][j].size();k++){
					adj[i].pb({j+n,pre[i][j][k]+i*m});
				}
			}
		}
	}
	vector<pair<int,int>> ans;
	for(int i=0;i<n;i++){
		if(adj[i].size()){
			cot2.clear();
			int cur=i;
			vector<pair<int,int>> cot;
			while(true){
				if(cur==i and adj[cur].size()==0){
					break;
				}
				cot.pb({cur,adj[cur].back().b});
				int ne=adj[cur].back().a;
				adj[cur].pop_back();
				cot.pb({ne,adj[ne].back().b});
				cur=adj[ne].back().a;
				adj[ne].pop_back();
			}
			vector<int> yy;
			for(int i=0;i<cot.size();i++){
				yy.pb(cot[i].a);
				vis[cot[i].a]=1;
			}
			while(yy.size()){
				vector<int> zz;
				for(int i=0;i<yy.size();i++){
					if(adj[yy[i]].size()){
						cur=yy[i];
						while(true){
							if(cur==yy[i] and adj[cur].size()==0){
								break;
							}
							if(vis[cur]==0){
								vis[cur]=1;
								zz.pb(cur);
							}
							xx[yy[i]].pb({cur,adj[cur].back().b});
							
					//		cot.pb({cur,adj[cur].back().b})
							int ne=adj[cur].back().a;
							adj[cur].pop_back();
							if(vis[ne]==0){
								vis[ne]=1;
								zz.pb(ne);

							}
							xx[yy[i]].pb({ne,adj[ne].back().b});
							cur=adj[ne].back().a;
							adj[ne].pop_back();
						}
					}
				}
				for(auto j:zz){
					vis[j]=1;
				}
				yy=zz;
			}
			
			for(auto j:cot){
				dfs(j.a);
				cot2.pb(j);
			}
			cot.clear();
			for(int j=0;j<cot2.size();j+=2){
				cot.pb(cot2[j]);
			}
			for(auto j:cot){
				assert(adj[j.a].size()==0);
			}

			ans.pb({cot[0].b,n*m});
			for(int i=cot.size()-1;i>=1;i--){
				ans.pb({cot[i].b,cot[(i+1)%((int)cot.size())].b});
			}
			ans.pb({n*m,cot[1].b});
		}
	}
	cout<<ans.size()<<endl;
	for(auto j:ans){
		tt[j.b]=tt[j.a];
		cout<<j.a+1<<" "<<j.b+1<<endl;
	}
/*	for(int j=0;j<n*m;j++){
		cout<<tt[j]<<",";
	}
	cout<<endl;
*/








	return 0;
}