// #pragma GCC optimize("O3,unroll-loops,Ofast")
// #pragma GCC target("avx2")
#include<bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace __gnu_pbds;
#define int long long
#define endl '\n'

using namespace std;
using pii = pair<int, int>;
using tup = tuple<int, int, int>;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
template <class T> using ordered_set = tree<T, null_type,
                         less<T>, rb_tree_tag, tree_order_statistics_node_update>;

const int inf = 1e18;
const int mod = 1000000007;
const double eps = 1e-9;
const int N = 2005;

void preprocess() {}

int n, k;
vector<tup> edges;

int p[N][5], sz[N], par[N];

int dist(int pi, int pj) {
    int ret = 0;
    for(int i=0; i<k; i++) ret += abs(p[pi][i] - p[pj][i]);
    return ret;
}

int find(int x) {
    return par[x] = (par[x] == x ? x : find(par[x]));
}

void join(int u, int v) {
    u = find(u);
    v = find(v);
    if(u == v) return;
    sz[u] += sz[v];
    par[v] = u;
}

int kruskal() {
    sort(edges.rbegin(), edges.rend());
    for(int i=1; i<=n; i++) par[i] = i, sz[i] = 1;

    int ans = 0;
    for(auto [d, u, v] : edges) {
        if(find(u) == find(v)) continue;
        ans += d;
        join(u, v);
    } 
    return ans;
}

void solve(int tc) {
    cin >> n >> k;

    for(int i=1; i<=n; i++) {
        for(int j=0; j<k; j++) cin >> p[i][j];

        for(int j=1; j<i; j++) {
            edges.push_back({dist(i, j), i, j});
        }
    }

    cout << kruskal() << endl;
}
     
int32_t main() {
    cin.tie(NULL)->sync_with_stdio(false);
    cout.precision(10);

    preprocess(); 

    int T = 1;
    // cin >> T;

    for (int i = 1; i <= T; i++) {
        solve(i);
    }

    return 0;
}