#include <bits/stdc++.h>
using namespace std;
vector<vector<long long int>>szl;
vector<bool>jart;
vector<int>be;

void dfs(int c){
    jart[c]=true;
    for(int x : szl[c]){
        if(!jart[x]) dfs(x);
    }

}
int main() {
	int n, m;
    cin>>n>>m;
    vector<vector<long long int>>t(n+1, vector<long long int>(m+1, 0));
    szl.resize(n*m+1);
    jart.resize(n*m+1);
    be.resize(n*m+1);
    for(int i=1; i<=n; i++){
        for(int j=1; j<=m; j++){
            cin>>t[i][j];
            int k=(i-1)*m+j%m;
            if(j%m==0) k+=m;
            if(t[i][j]>0 && t[i][j-1]>0){
                if(t[i][j]>=t[i][j-1]){
                    szl[k-1].push_back(k);
                    be[k]++;
                }
                if(t[i][j]<=t[i][j-1]){
                    szl[k].push_back(k-1);
                    be[k-1]++;
                }
            }
            if(t[i][j]>0 && t[i-1][j]>0){
                if(t[i][j]>=t[i-1][j]){
                    szl[k-m].push_back(k);
                    be[k]++;
                }
                if(t[i][j]<=t[i-1][j]){
                    szl[k].push_back(k-m);
                    be[k-m]++;
                }
            }
        }
    }
    /*for(int i=1; i<=n*m; i++){
        for(int x : szl[i]) cout<<x<<" ";
    cout<<endl;
    }*/
   int db=0;
   vector<pair<int,int>>mo;
    for(int i=1; i<=n; i++){
        for(int j=1; j<=m; j++){
            if(t[i][j]>0){
                int k=(i-1)*m+j%m;
                if(j%m==0) k+=m;
                if(!jart[k] && be[k]==0){
                    mo.push_back({i, j});
                    dfs(k);
                    db++;
                }
            }
    }
    }
    cout<<db<<endl;
    for(auto x : mo) cout<<x.first<<" "<<x.second<<endl;

}