#include <bits/stdc++.h>
using namespace std;
using pii = pair<int, int>;

int n,m;
vector<pii> xy={{0,-1},{-1,0},{0,1},{1,0}};

void dfs(vector<vector<int>>& nm,int x,int y){
    int cur=nm[x][y];
    nm[x][y]=0;
    for(auto [a,b]:xy){
        int i=x+a,j=y+b;
        if(i>n||i<1||j>m||j<1||!nm[i][j])continue;
        if(nm[i][j]>=cur)dfs(nm,i,j);
    }
}

int main(){
    cin>>n>>m;
    vector<vector<int>> nm(n+1,vector<int>(m+1));
    vector<array<int,3>> q;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            cin>>nm[i][j];
            if(nm[i][j])q.push_back({nm[i][j],i,j});
        }
    sort(q.begin(),q.end());
    vector<pii> ans;
    for(auto [a,x,y]:q){
        if(nm[x][y]){
            dfs(nm,x,y);
            ans.push_back({x,y});
        }
    }
    cout<<ans.size()<<'\n';
    for(auto [a,b]:ans)cout<<a<<' '<<b<<'\n';
}