#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define cmin(a,b) a=min(a,b)
#define cmax(a,b) a=max(a,b)
#define modd(a,b) a=a%b
#define pb push_back
#define db double
#define lowbit(x) x&(-x)
#define cerr(x) cout<<#x<<"="<<x<<endl
const ll mod=1e9+7;
ll ksm(ll a,ll b) {ll ans=1,bs=a;while(b) {if(b&1) ans=ans*bs%mod;bs=bs*bs%mod;b>>=1;}return ans;}
#define fi first
#define se second
#define N 1010
int n,m;
int c[N][N],a[N];
vector<pair<int,int>>g[N];
void solve() {
    cin>>n>>m;
    int cnt=0;
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cin>>c[i][j],a[++cnt]=c[i][j];
    sort(a+1,a+cnt+1);
    int tot=unique(a+1,a+cnt+1)-a-1;
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) {
        c[i][j]=lower_bound(a+1,a+tot+1,c[i][j])-a;
    }
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) {
        g[c[i][j]].pb(make_pair(i,j));
    }
    // for(int i=1;i<=n;i++,cout<<endl) for(int j=1;j<=m;j++) cout<<c[i][j]<<" ";
    ll ans=0;
    for(int i=1;i<=tot;i++) {
        if(g[i].size()<=1) continue;
        sort(g[i].begin(),g[i].end(),[&](pair<int,int>a,pair<int,int>b){
            return a.first<b.first;
        });
        for(int j=1;j<g[i].size();j++) {
            ans+=1ll*j*((int)g[i].size()-1-j+1)*(g[i][j].first-g[i][j-1].first);
        }
        sort(g[i].begin(),g[i].end(),[&](pair<int,int>a,pair<int,int>b) {
            return a.second<b.second;
        });
        for(int j=1;j<g[i].size();j++) {
            ans+=1ll*j*((int)g[i].size()-1-j+1)*(g[i][j].second-g[i][j-1].second);
        }
    }
    cout<<ans*2<<endl;
}
int main() {
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    // int T;
    // cin>>T;
    // while(T--) 
    solve();
    return 0;
}