#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define llint long long #define inf 1e18 #define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++) #define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++) #define chmin(x, y) (x) = min((x), (y)) #define chmax(x, y) (x) = max((x), (y)) using namespace std; typedef pair P; llint n; llint a[1005][1005]; vector vec[1005][1005]; llint d[1005][1005]; llint s[1005][1005]; int main(void) { ios::sync_with_stdio(0); cin.tie(0); cin >> n; for(int y = 1; y <= n; y++){ for(int x = 1; x <= n; x++){ cin >> a[x][y]; } } for(int y = 1; y <= n; y++){ for(int x = 1; x <= n; x++){ if(y+x-1 <= n) vec[a[x][y]][y+x-1].push_back(x-1); } } for(int i = 1; i <= n; i++){ llint sum = 0, cnt = 0; for(int y = 1; y <= n; y++){ s[i][y] += sum; for(int j = 0; j < vec[i][y].size(); j++){ llint x = vec[i][y][j]; sum += x; cnt++; } sum += cnt; } } for(int y = 1; y <= n; y++){ for(int x = 1; x <= n; x++){ vec[x][y].clear(); } } for(int y = 1; y <= n; y++){ for(int x = 1; x <= n; x++){ if(y-(x-1) >= 1) vec[a[x][y]][y-(x-1)].push_back(x-1); } } for(int i = 1; i <= n; i++){ llint sum = 0, cnt = 0; for(int y = n; y >= 1; y--){ s[i][y] += sum; for(int j = 0; j < vec[i][y].size(); j++){ llint x = vec[i][y][j]; sum += x; cnt++; } sum += cnt; } } for(llint y = 1; y <= n; y++){ for(llint x = 1; x <= n; x++){ llint id = a[x][y]; llint l = max(1LL, y-(x-1)), r = min(n, y+(x-1)); if(l > r) continue; d[id][l] += x-1, d[id][r+1] -= x-1; } } for(int i = 1; i <= n; i++){ for(int y = 1; y <= n; y++) d[i][y] += d[i][y-1]; } for(int i = 1; i <= n; i++){ for(int y = 1; y <= n; y++) s[i][y] += d[i][y]; } llint ans = 0; for(int i = 1; i <= n; i++){ llint mn = inf; for(int y = 1; y <= n; y++) mn = min(mn, s[i][y]); ans += mn; } cout << ans << endl; return 0; }