結果
| 問題 |
No.1028 闇討ち
|
| コンテスト | |
| ユーザー |
wk
|
| 提出日時 | 2020-04-19 17:22:42 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 220 ms / 2,000 ms |
| コード長 | 4,082 bytes |
| コンパイル時間 | 1,711 ms |
| コンパイル使用メモリ | 170,168 KB |
| 実行使用メモリ | 11,328 KB |
| 最終ジャッジ日時 | 2024-10-05 19:32:08 |
| 合計ジャッジ時間 | 5,648 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 20 |
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:168:17: warning: 'c' may be used uninitialized [-Wmaybe-uninitialized]
168 | c[A[0][0]][1]++;
| ~~~~~~~~~~~~^
main.cpp:158:10: note: 'c' declared here
158 | long c[1005][2];
| ^
main.cpp:168:13: warning: 'A' may be used uninitialized [-Wmaybe-uninitialized]
168 | c[A[0][0]][1]++;
| ~~~~~~^
main.cpp:155:10: note: 'A' declared here
155 | long A[1005][1005];
| ^
main.cpp:168:17: warning: 'c' may be used uninitialized [-Wmaybe-uninitialized]
168 | c[A[0][0]][1]++;
| ~~~~~~~~~~~~^
main.cpp:158:10: note: 'c' declared here
158 | long c[1005][2];
| ^
ソースコード
#include <bits/stdc++.h>
#define REP(i, n) for(int i = 0; (i) < (n); (i)++)
#define INF 5000000000000000
using namespace std;
struct Edge{
int to;
long weight;
Edge(int t, long w) : to(t), weight(w) {}
};
long modpow(long a, long n, long mod) {
long res = 1;
while (n > 0) {
if (n & 1) res = res * a % mod;
a = a * a % mod;
n >>= 1;
}
return res;
}
// a^{-1} mod を計算する
long modinv(long a, long mod) {
return modpow(a, mod - 2, mod);
}
struct compare1 {
bool operator()(const pair<long, long>& value,
const long& key)
{
return (value.first < key);
}
bool operator()(const long& key,
const pair<long, long>& value)
{
return (key < value.first);
}
};
struct UnionFind {
vector<int> par;
vector<int> rank;
UnionFind(int n = 1){
init(n);
}
void init(int n = 1){
par.resize(n);
rank.resize(n);
REP(i, n) par[i] = i, rank[i] = 0;
}
int root(int x){
if(par[x] == x) return x;
else return par[x] = root(par[x]);
}
bool issame(int x, int y){
return root(x) == root(y);
}
bool merge(int x, int y){
x = root(x); y = root(y);
if(x == y) return false;
if(rank[x] < rank[y]) swap(x, y);
if(rank[x] == rank[y]) rank[x]++;
par[y] = x;
return true;
}
};
template<class Abel> struct weightedUnionFind{
vector<int> par;
vector<int> rank;
vector<Abel> diff_weight;
weightedUnionFind(int n = 1, Abel SUM_UNITY = 0){
init(n, SUM_UNITY);
}
void init(int n = 1, Abel SUM_UNITY = 0){
par.resize(n); rank.resize(n); diff_weight.resize(n);
REP(i, n) par[i] = i, rank[i] = 0, diff_weight[i] = SUM_UNITY;
}
int root(int x){
if(par[x] == x) return x;
else{
int r = root(par[x]);
diff_weight[x] += diff_weight[par[x]];
return par[x] = r;
}
}
Abel weight(int x){
root(x);
return diff_weight[x];
}
bool issame(int x, int y){
return root(x) == root(y);
}
bool merge(int x, int y, Abel w){
w += weight(x); w -= weight(y);
x = root(x); y = root(y);
if(x == y) return false;
if(rank[x] < rank[y]) swap(x, y), w = -w;
if(rank[x] == rank[y]) rank[x]++;
par[y] = x;
diff_weight[y] = w;
return true;
}
Abel diff(int x, int y){
return weight(y) - weight(x);
}
};
using Graph = vector<vector<int>>;
using P = pair<int, int>;
/*
void dijkstra(int s, int V, Graph &G, long* d){
priority_queue<P, vector<P>, greater<P>> pque;
fill(d, d + V, INF);
d[s] = 0;
pque.push(P(0, s));
while(!pque.empty()){
P p = pque.top(); pque.pop();
int now = p.second;
if(d[now] < p.first) continue;
REP(i, G[now].size()){
Edge e = G[now][i];
if(d[e.to] > d[now] + e.weight){
d[e.to] = d[now] + e.weight;
pque.push(P(d[e.to], e.to));
}
}
}
}
*/
long GCD(long a, long b){
if(b == 0l) return a;
if(a < b) return GCD(b, a);
else return GCD(b, a%b);
}
int main()
{
int N; cin >> N;
long A[1005][1005];
long b[1005];
long ans[1005];
long c[1005][2];
REP(i, N) REP(j, N){
int a; cin >> a; a--;
A[i][j] = a;
}
REP(i, N) b[i] = 0l;
REP(i, N) REP(j, N) b[A[i][j]] += max(i, j);
REP(i, N) ans[i] = b[i];
REP(i, N) c[i][0] = 0l, c[i][1] = 0l;
REP(i, N) REP(j, N) if(i > j) c[A[i][j]][0]++;
c[A[0][0]][1]++;
REP(i, N-1){
REP(j, N) b[j] -= c[j][0];
REP(j, N) b[j] += c[j][1];
REP(j, N) ans[j] = min(ans[j], b[j]);
for(int j = i+1; j >= 0; j--) c[A[j][i+1-j]][1]++;
for(int j = i+1; j < N; j++) c[A[j][j-i-1]][0]--;
}
long su = 0l;
REP(i, N) su += ans[i];
cout << su << endl;
return 0;
}
wk