#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll,ll> pl;
typedef pair<int,int> pii;
typedef tuple<int,int,int> tt;
#define all(a) a.begin(), a.end()
#define filter(a) a.erase(unique(all(a)), a.end())
vector<int> cmp;
int getID (int targ) {
return lower_bound(all(cmp), targ) - cmp.begin();
}
const int mn = 4e6 + 6;
int lo[2002][2002], hi[2002][2002], loCol[2002][2002], hiCol[2002][2002], hiRow[2002][2002];
int calcLo[mn], calcHi[mn];
vector<pii> contain[mn];
bool check (int delta) {
int pos = -1, curLo = INT_MAX, curHi = INT_MIN;
for (int k = 0; k < cmp.size(); k++) {
while (pos + 1 < cmp.size() && cmp[pos + 1] < cmp[k] - delta) {
pos++;
for (auto [row, col] : contain[pos]) {
curLo = min(curLo, lo[row][col]);
curHi = max(curHi, hi[row][col]);
}
}
int cHi = max(calcHi[k], curHi);
int cLo = min(calcLo[k], curLo);
if (cHi == INT_MIN || cHi - cLo <= delta) return 1;
}
return 0;
}
int solve (vector<vector<int>> &vec, int n, int m) {
for (int j = 1; j <= m; j++) {
int curLo = INT_MAX, curHi = INT_MIN;
for (int i = n; i >= 1; i--) {
loCol[i][j] = curLo = min(curLo, vec[i][j]);
hiCol[i][j] = curHi = max(curHi, vec[i][j]);
}
}
for (int i = 1; i <= n; i++) {
int curLo = INT_MAX, curHi = INT_MIN;
for (int j = 1; j <= m; j++) {
lo[i][j] = curLo = min(curLo, loCol[i][j]);
hi[i][j] = curHi = max(curHi, hiCol[i][j]);
}
curHi = INT_MIN;
for (int j = m; j >= 1; j--)
hiRow[i][j] = curHi = max(curHi, vec[i][j]);
}
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++)
contain[getID(vec[i][j])].push_back(make_pair(i, j));
int curLo = INT_MAX, curHi = INT_MIN;
for (int k = cmp.size() - 1; k >= 0; k--) {
calcLo[k] = curLo, calcHi[k] = curHi;
for (auto [row, col] : contain[k]) {
curLo = min(curLo, lo[row][col]);
curHi = max(curHi, hi[row][col]);
}
}
int ans = 0;
if (!check(0)) {
for (int i = 29; i >= 0; i--) {
int mask = ans | (1 << i);
if (!check(mask)) ans |= (1 << i);
}
ans++;
}
for (int k = 0; k < cmp.size(); k++)
contain[k].clear();
return ans;
}
int main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int n, m; cin >> n >> m;
vector<vector<int>> vec(n + 1, vector<int>(m + 1));
vector<vector<int>> rot(m + 1, vector<int>(n + 1));
//cmp.push_back(0);
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
cin >> vec[i][j];
rot[j][n - i + 1] = vec[i][j];
cmp.push_back(vec[i][j]);
}
}
sort(all(cmp)); filter(cmp);
cout << min(solve(vec, n, m), solve(rot, m, n));
return 0;
}
