#include <bits/stdc++.h>
using namespace std;
struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_;
#define FOR(i, begin, end) for(int i=(begin);i<(end);i++)
#define REP(i, n) FOR(i,0,n)
#define IFOR(i, begin, end) for(int i=(end)-1;i>=(begin);i--)
#define IREP(i, n) IFOR(i,0,n)
#define Sort(v) sort(v.begin(), v.end())
#define Reverse(v) reverse(v.begin(), v.end())
#define all(v) v.begin(),v.end()
#define SZ(v) ((int)v.size())
#define Lower_bound(v, x) distance(v.begin(), lower_bound(v.begin(), v.end(), x))
#define Upper_bound(v, x) distance(v.begin(), upper_bound(v.begin(), v.end(), x))
#define Max(a, b) a = max(a, b)
#define Min(a, b) a = min(a, b)
#define bit(n) (1LL<<(n))
#define bit_exist(x, n) ((x >> n) & 1)
#define debug(x) cout << #x << "=" << x << endl;
#define vdebug(v) cout << #v << "=" << endl; REP(i, v.size()){ cout << v[i] << ","; } cout << endl;
#define mdebug(m) cout << #m << "=" << endl; REP(i, m.size()){ REP(j, m[i].size()){ cout << m[i][j] << ","; } cout << endl;}
#define pb push_back
#define f first
#define s second
//#define int long long
#define INF 1000000000
template<typename T> istream &operator>>(istream &is, vector<T> &v){ for (auto &x : v) is >> x; return is; }
template<typename T> ostream &operator<<(ostream &os, vector<T> &v){ for(int i = 0; i < v.size(); i++) { cout << v[i]; if(i != v.size() - 1) cout << endl; }; return os; }
template<typename T> void Out(T x) { cout << x << endl; }
template<typename T1, typename T2> void Ans(bool f, T1 y, T2 n) { if(f) Out(y); else Out(n); }

using vec = vector<int>;
using mat = vector<vec>;
using Pii = pair<int, int>;
using PiP = pair<int, Pii>;
using PPi = pair<Pii, int>;
using bools = vector<bool>;
using pairs = vector<Pii>;

//int dx[4] = {1,0,-1,0};
//int dy[4] = {0,1,0,-1};
//char d[4] = {'D','R','U','L'};

const int mod = 1000000007;
//const int mod = 998244353;
#define Add(x, y) x = (x + (y)) % mod
#define Mult(x, y) x = (x * (y)) % mod

struct SparseTable
{
    using T = int;
    vector<vector<T>> table;

    T F(T a, T b){
        return max(a, b);
    }
    
    SparseTable(int N){
        int n = 32 - __builtin_clz((signed)N);
        table = vector<vector<T>>(n, vector<T>(N));
    }

    SparseTable(int N, vector<T> &A){
        int n = 32 - __builtin_clz((signed)N);
        table = vector<vector<T>>(n, vector<T>(N));
        initialize(A);
    }

    void initialize(vector<T> &A){
        for(int i = 0; i < A.size(); i++) table[0][i] = A[i];
        for(int n = 0; n < table.size() - 1; n++){
            for(int i = 1 << n; i < A.size(); i++){
                table[n + 1][i - (1 << n)] = F(table[n][i - (1 << n)], table[n][i]);
            }
        }
    }

    // the range is [l, r]
    // 0 <= l <= r < N must be satisfied
    T get(int l, int r){
        int n = 31 - __builtin_clz((signed)(r - l + 1));
        return F(table[n][l], table[n][r - (1 << n) + 1]);
    }

};

int calcsum(int x1, int x2, int y1, int y2, mat &S){
    return S[x2 + 1][y2 + 1] - S[x2 + 1][y1] - S[x1][y2 + 1] + S[x1][y1];
}

signed main(){

    int H, W; cin >> H >> W;
    vector<string> A(H); cin >> A;

    mat S(H + 1, vec(W + 1, 0));
    REP(i, H) REP(j, W) if(A[i][j] == '#') S[i + 1][j + 1]++;
    REP(i, H + 1) REP(j, W) S[i][j + 1] += S[i][j];
    REP(j, W + 1) REP(i, H) S[i + 1][j] += S[i][j];
    //mdebug(S);

    int M = 16;
    vector<mat> dp1(H, mat(H, vec(W, -INF))), dp2(W, mat(W, vec(H, -INF)));
    REP(i1, H) FOR(i2, i1, H) REP(j, W){
        int v = 0;
        if(A[i1][j] == '#') v = 1;
        if(calcsum(i1, i2, j, j, S) == v * (i2 - i1 + 1)){
            int x1 = j, x2 = W;
            while(x2 - x1 > 1){
                int x = (x1 + x2) / 2;
                if(calcsum(i1, i2, j, x, S) == v * (i2 - i1 + 1) * (x - j + 1)) x1 = x;
                else x2 = x;
            }
            dp1[i1][i2][j] = x1;
        }
    }
    REP(j1, W) FOR(j2, j1, W) REP(i, H){
        int v = 0;
        if(A[i][j1] == '#') v = 1;
        if(calcsum(i, i, j1, j2, S) == v * (j2 - j1 + 1)){
            int x1 = i, x2 = H;
            while(x2 - x1 > 1){
                int x = (x1 + x2) / 2;
                if(calcsum(i, x, j1, j2, S) == v * (j2 - j1 + 1) * (x - i + 1)) x1 = x;
                else x2 = x;
            }
            dp2[j1][j2][i] = x1;
        }
    }
    if(dp1[0][H - 1][0] == W - 1 || dp2[0][W - 1][0] == H - 1){
        Out(0);
        return 0;
    }
    
    SparseTable ST1(W), ST2(H);
    FOR(k, 1, M + 1){
        vector<mat> dq1(H, mat(H, vec(W, -INF))), dq2(W, mat(W, vec(H, -INF)));
        REP(i1, H) FOR(i2, i1, H){
            
            ST1.initialize(dp1[i1][i2]);

            REP(j, W) if(dp1[i1][i2][j] >= 0){
                int tmp = ST1.get(j, min(dp1[i1][i2][j] + 1, W - 1));
                Max(dq1[i1][i2][j], tmp);
            }
        }

        REP(j1, W) FOR(j2, j1, W){
            
            ST2.initialize(dp2[j1][j2]);

            REP(i, H) if(dp2[j1][j2][i] >= 0){
                int tmp = ST2.get(i, min(dp2[j1][j2][i] + 1, H - 1));
                Max(dq2[j1][j2][i], tmp);
            }
        }
        dp1.swap(dq1); dp2.swap(dq2);
        REP(i1, H) FOR(i2, i1, H){
            int jmax = 0;
            REP(j, W){
                if(dp2[j][j][i1] < i2) continue;
                if(jmax < j) jmax = j;
                while(jmax < W - 1 && dp2[j][jmax + 1][i1] >= i2) jmax++;
                Max(dp1[i1][i2][j], jmax);
            }
        }
            
        REP(j1, W) FOR(j2, j1, W){
            int imax = 0;
            REP(i, H){
                if(dp1[i][i][j1] < j2) continue;
                if(imax < i) imax = i;
                while(imax < H - 1 && dp1[i][imax + 1][j1] >= j2) imax++;
                Max(dp2[j1][j2][i], imax);
            }
        }

        if(dp1[0][H - 1][0] == W - 1 || dp2[0][W - 1][0] == H - 1){
            Out(k);
            return 0;
        }
    }

    return 0;
}


Measured memory_bytes: 106452