結果

問題 No.1028 闇討ち
コンテスト
ユーザー Jumbo_kpr
提出日時 2020-04-17 21:58:23
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 155 ms / 2,000 ms
コード長 4,153 bytes
コンパイル時間 1,327 ms
コンパイル使用メモリ 119,472 KB
実行使用メモリ 20,736 KB
最終ジャッジ日時 2024-10-03 12:51:44
合計ジャッジ時間 3,503 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

//#pragma GCC optimize("Ofast")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
//#pragma GCC optimize("unroll-loops")
//#pragma warning(disable : 4996)
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<math.h>
#include<iomanip>
#include<set>
#include<numeric>
#include<cstring>
#include<cstdio>
#include<functional>
#include<bitset>
#include<limits.h>
#include<cassert>
#include<iterator>
#include<complex>
#include<stack>
#include<unordered_map>
#include<unordered_set>
#include<time.h>
#include<random>
#include<array>


using namespace std;

#define REP(i, n) for(int i = 0;i < n;++i)
#define REPR(i, n) for(int i = n-1;i >= 0;--i)
#define FOR(i, m, n) for(int i = m;i < n;++i)
#define FORR(i, m, n) for(int i = m-1;i >= n;--i)
#define SORT(v, n) sort(v, v+n);
#define VSORT(v) sort(v.begin(), v.end());
#define REVERSE(v,n) reverse(v,v+n);
#define VREVERSE(v) reverse(v.begin(), v.end())
#define ll long long
#define print(x) cout<<(x)<<'\n'
#define pe(x) cout<<(x)<<" "
#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define lb(v,n) lower_bound(v.begin(), v.end(), (n))
#define ub(v,n) upper_bound(v.begin(), v.end(), (n))
#define int long long
//#define double long double
#define all(x) (x).begin(), (x).end()
#define print_space(v) REP(i,v.size())cout << v[i] << ((i == v.size() - 1) ? "\n" : " ")
template<typename T1, typename T2> inline void chmin(T1 & a, T2 b) { if (a > b) a = b; }
template<typename T1, typename T2> inline void chmax(T1& a, T2 b) { if (a < b) a = b; }
typedef pair<int, int>pii;
typedef array<int, 3> arr3;
std::random_device rd;
std::mt19937 mt(rd());
constexpr ll MOD = 1e9 + 7; constexpr int MAX = 2000020;
const double pi = acos(-1); constexpr double EPS = 1e-8;
constexpr ll INF = 1e18;
void yes(bool c) { if (c)print("Yes"); else print("No"); };

//const int mod = 1000000007;
const int mod = 998244353;
struct mint {
	ll x; // typedef long long ll;
	mint(ll x = 0) :x((x%mod + mod) % mod) {}
	mint operator-() const { return mint(-x); }
	mint& operator+=(const mint a) {
		if ((x += a.x) >= mod) x -= mod;
		return *this;
	}
	mint& operator-=(const mint a) {
		if ((x += mod - a.x) >= mod) x -= mod;
		return *this;
	}
	mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this; }
	mint operator+(const mint a) const { return mint(*this) += a; }
	mint operator-(const mint a) const { return mint(*this) -= a; }
	mint operator*(const mint a) const { return mint(*this) *= a; }
	mint pow(ll t) const {
		if (!t) return 1;
		mint a = pow(t >> 1);
		a *= a;
		if (t & 1) a *= *this;
		return a;
	}

	// for prime mod
	mint inv() const { return pow(mod - 2); }
	mint& operator/=(const mint a) { return *this *= a.inv(); }
	mint operator/(const mint a) const { return mint(*this) /= a; }
};
istream& operator>>(istream& is, const mint& a) { return is >> a.x; }
ostream& operator<<(ostream& os, const mint& a) { return os << a.x; }


int A[1010][1010];
void solve() {
	int N; cin >> N;
	//vector<int>ans(N);
	vector<vector<int>>X(N),Y(N);
	REP(i, N) {
		REP(j, N) {
			int a; cin >> a;
			a--;
			X[a].push_back(i);
			Y[a].push_back(j);
		}
	}
	ll ans = 0;
	REP(i, N) {
		int l = 0, r = N;
		int ret = 1e9;
		while (r - l > 10) {
			int m1 = (l * 2 + r) / 3, m2 = (l + r * 2)/3;
			int res1 = 0, res2 = 0;
			REP(j, N) {
				res1 += Y[i][j];
				int x1 = m1 + Y[i][j],x2= m1 - Y[i][j];
				if (X[i][j] > x1)res1 += X[i][j] - x1;
				else if (X[i][j] < x2)res1 += x2 - X[i][j];
			}
			REP(j, N) {
				res2 += Y[i][j];
				int x1 = m2 + Y[i][j], x2 = m2 - Y[i][j];
				if (X[i][j] > x1)res2+= X[i][j] - x1;
				else if (X[i][j] < x2)res2 += x2 - X[i][j];
			}
			if (res1 < res2)r = m2;
			else l = m1;
			ret = min({ ret,res1,res2 });
		}
		FOR(m, l, r + 1) {
			int res1 = 0;
			REP(j, N) {
				res1 += Y[i][j];
				int x1 = m + Y[i][j], x2 = m - Y[i][j];
				if (X[i][j] > x1)res1 += X[i][j] - x1;
				else if (X[i][j] < x2)res1 += x2 - X[i][j];
			}
			ret = min(ret, res1);
		}
		ans += ret;
	}
	print(ans);


}

signed main() {
	cin.tie(0);
	ios::sync_with_stdio(false);
	//int q; cin >> q;
	//while (q--)
	solve();
}
0