#include <bits/stdc++.h>
#ifdef CKISEKI
#define safe cerr<<__PRETTY_FUNCTION__<<" line "<<__LINE__<<" safe\n"
#define debug(a...) debug_(#a, a)
#define orange(a...) orange(#a, a)
using std::cerr;
template <typename ...T>
void debug_(const char *s, T ...a) {
    cerr << "\e[1;32m(" << s << ") = (";
    int cnt = sizeof...(T);
    (..., (cerr << a << (--cnt ? ", " : ")\e[0m\n")));
}
template <typename I>
void orange_(const char *s, I L, I R) {
    cerr << "\e[1;32m[ " << s << " ] = [ ";
    for (int f = 0; L != R; ++L)
        cerr << (f++ ? ", " : "") << *L;
    cerr << " ]\e[0m\n";
}
#else
#define safe ((void)0)
#define debug(...) ((void)0)
#define orange(...) ((void)0)
#endif

#define all(v) begin(v),end(v)

using namespace std;

constexpr int maxc = 1'000'000 + 5;
constexpr int maxn = 500 + 5;

int a[maxn][maxn];

bool exi[maxc];
int x_min[maxc], x_max[maxc];
int y_min[maxc], y_max[maxc];

int v[maxn][maxn];
int vis[maxc], visc;

int main() {
    cin.tie(nullptr) -> sync_with_stdio(false);
    int n, m; cin >> n >> m;
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < m; ++j)
            cin >> a[i][j];

    for (int i = 0; i < maxc; ++i) {
        y_min[i] = m;
        x_min[i] = n;
    }

    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < m; ++j) {
            const int c = a[i][j];
            x_min[c] = min(x_min[c], i);
            x_max[c] = max(x_max[c], i);
            y_min[c] = min(y_min[c], j);
            y_max[c] = max(y_max[c], j);
            exi[c] = true;
        }
    }

    for (int c = 0; c < maxc; ++c) {
        if (not exi[c]) continue;
        v[x_min[c]][y_min[c]] += 1;
    }
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < m; ++j) {
            if (i > 0) v[i][j] += v[i - 1][j];
            if (j > 0) v[i][j] += v[i][j - 1];
            if (i > 0 and j > 0) v[i][j] -= v[i - 1][j - 1];
        }
    }

    int64_t ans = 0;
    for (int c = 0; c < maxc; ++c) {
        if (not exi[c]) continue;
        int s = v[x_max[c]][y_max[c]];
        if (x_min[c] > 0) s -= v[x_min[c] - 1][y_max[c]];
        if (y_min[c] > 0) s -= v[x_max[c]][y_min[c] - 1];
        if (x_min[c] > 0 and y_min[c] > 0)
            s += v[x_min[c] - 1][y_min[c] - 1];
        
        visc += 1;
        auto upd = [&](int x, int y) {
            int o = a[x][y];
            if (vis[o] == visc) return;
            vis[o] = visc;
            if (x_min[c] > x_min[o] or x_min[o] > x_max[c])
                return;
            if (y_min[c] > y_min[o] or y_min[o] > y_max[c])
                return;
            if (x_min[c] > x_min[o] or x_min[o] > x_max[c] or y_min[c] > y_min[o] or y_min[o] > y_max[c])
                s -= 1;
        };

        for (int j = y_min[c]; j <= y_max[c]; ++j) {
            upd(x_min[c], j);
            upd(x_max[c], j);
        }
        for (int i = x_min[c]; i <= x_max[c]; ++i) {
            upd(i, y_min[c]);
            upd(i, y_max[c]);
        }
        ans += s - 1;
    }
    cout << ans << '\n';
    return 0;
}