// competitive-verifier: PROBLEM https://yukicoder.me/problems/no/1028 #pragma GCC optimize("Ofast,fast-math,unroll-all-loops") #include #if !defined(ATCODER) && !defined(EVAL) #pragma GCC target("sse4.2,avx2,bmi2") #endif template constexpr bool chmax(T &a, const U &b) { return a < (T)b ? a = (T)b, true : false; } template constexpr bool chmin(T &a, const U &b) { return (T)b < a ? a = (T)b, true : false; } constexpr std::int64_t INF = 1000000000000000003; constexpr int Inf = 1000000003; constexpr double EPS = 1e-7; constexpr double PI = 3.14159265358979323846; #define FOR(i, m, n) for (int i = (m); i < int(n); ++i) #define FORR(i, m, n) for (int i = (m) - 1; i >= int(n); --i) #define FORL(i, m, n) for (std::int64_t i = (m); i < std::int64_t(n); ++i) #define rep(i, n) FOR (i, 0, n) #define repn(i, n) FOR (i, 1, n + 1) #define repr(i, n) FORR (i, n, 0) #define repnr(i, n) FORR (i, n + 1, 1) #define all(s) (s).begin(), (s).end() struct Sonic { Sonic() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); std::cout << std::fixed << std::setprecision(20); } constexpr void operator()() const {} } sonic; struct increment_impl { template const increment_impl &operator>>(std::vector &v) const { for (auto &x : v) ++x; return *this; } } Inc; struct decrement_impl { template const decrement_impl &operator>>(std::vector &v) const { for (auto &x : v) --x; return *this; } } Dec; struct sort_impl { template const sort_impl &operator>>(std::vector &v) const { std::sort(v.begin(), v.end()); return *this; } } Sort; struct unique_impl { template const unique_impl &operator>>(std::vector &v) const { std::sort(v.begin(), v.end()); v.erase(std::unique(v.begin(), v.end()), v.end()); return *this; } } Uniq; using namespace std; using ll = std::int64_t; using ld = long double; template std::istream &operator>>(std::istream &is, std::pair &p) { return is >> p.first >> p.second; } template std::istream &operator>>(std::istream &is, std::vector &v) { for (T &i : v) is >> i; return is; } template std::ostream &operator<<(std::ostream &os, const std::pair &p) { return os << '(' << p.first << ',' << p.second << ')'; } template std::ostream &operator<<(std::ostream &os, const std::vector &v) { for (auto it = v.begin(); it != v.end(); ++it) os << (it == v.begin() ? "" : " ") << *it; return os; } template void co(Head &&head, Tail &&...tail) { if constexpr (sizeof...(tail) == 0) std::cout << head << '\n'; else std::cout << head << ' ', co(std::forward(tail)...); } template void ce(Head &&head, Tail &&...tail) { if constexpr (sizeof...(tail) == 0) std::cerr << head << '\n'; else std::cerr << head << ' ', ce(std::forward(tail)...); } void Yes(bool is_correct = true) { std::cout << (is_correct ? "Yes\n" : "No\n"); } void No(bool is_not_correct = true) { Yes(!is_not_correct); } void YES(bool is_correct = true) { std::cout << (is_correct ? "YES\n" : "NO\n"); } void NO(bool is_not_correct = true) { YES(!is_not_correct); } void Takahashi(bool is_correct = true) { std::cout << (is_correct ? "Takahashi" : "Aoki") << '\n'; } void Aoki(bool is_not_correct = true) { Takahashi(!is_not_correct); } #include /// @brief slope trick template struct slope_trick { slope_trick() = default; [[nodiscard]] T min_x() const { return l.top(); } [[nodiscard]] T min_value() const { return min_f; } /// @brief Add f(x) = a void add_const(T a) { min_f += a; } /// @brief Add f(x) = max(0, x - a) void add_x_minus_a(T a) { if (!l.empty() && l.top() > a) { min_f += l.top() - a; } l.push(a); r.push(l.top()); l.pop(); } /// @brief Add f(x) = max(0, a - x) void add_a_minus_x(T a) { if (!r.empty() && a > r.top()) { min_f += a - r.top(); } r.push(a); l.push(r.top()); r.pop(); } /// @brief Add f(x) = abs(x - a) = max(0, x - a) + max(0, a - x) void add_abs(T a) { add_x_minus_a(a); add_a_minus_x(a); } /// @brief f(x) <- min_{y <= x} f(y) void clear_right() { r = decltype(r)(); } /// @brief f(x) <- min_{y >= x} f(y) void clear_left() { l = decltype(l)(); } private: T min_f{}; std::priority_queue l; std::priority_queue, std::greater> r; }; int main(void) { int n; std::cin >> n; std::vector a(n, std::vector(n, 0)); for (auto& v : a) { for (auto& e : v) std::cin >> e; } std::vector> v(n); for (int i = 0; i < n; ++i) { for (int j = 0; j < n; ++j) { v[a[i][j] - 1].add_const(j); v[a[i][j] - 1].add_x_minus_a(i + j); v[a[i][j] - 1].add_a_minus_x(i - j); } } int ans = 0; for (int i = 0; i < n; ++i) ans += v[i].min_value(); std::cout << ans << '\n'; return 0; }