# | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
---|---|---|---|---|---|---|---|
640492 | piOOE | 도시들 (IOI15_towns) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
#include "towns.h"
using namespace std;
int hubDistance(int n, int sub) {
map<pair<int, int>, int> mp;
auto mygetDistance = [&](int x, int y) {
if (mp.count({x, y})) {
return mp[{x, y}];
} else {#include <bits/stdc++.h>
#include "towns.h"
using namespace std;
int hubDistance(int n, int sub) {
map<pair<int, int>, int> mp;
auto mygetDistance = [&](int x, int y) {
if (mp.count({x, y})) {
return mp[{x, y}];
} else {
return mp[{x, y}] = mp[{y, x}] = getDistance(x, y);
}
};
if (sub <= 2) {
vector<int> dist(n);
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, 0);
}
int v = max_element(dist.begin(), dist.end()) - dist.begin();
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, v);
}
int u = max_element(dist.begin(), dist.end()) - dist.begin();
vector<int> dist2(n);
for (int i = 0; i < n; ++i) {
dist2[i] = mygetDistance(i, u);
}
int ans = numeric_limits<int>::max();
for (int i = 0; i < n; ++i) {
int e = dist[i] + dist2[i] - dist[u] >> 1;
ans = min(ans, max(dist[i] - e, dist2[i] - e));
}
return ans;
} else {
vector<int> dist(n);
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, 0);
}
int a = max_element(dist.begin(), dist.end()) - dist.begin();
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, a);
}
int b = max_element(dist.begin(), dist.end()) - dist.begin();
vector<int> dist2(n);
for (int i = 0; i < n; ++i) {
dist2[i] = mygetDistance(i, b);
}
int ans = numeric_limits<int>::max();
for (int i = 0; i < n; ++i) {
int e = dist[i] + dist2[i] - dist[b] >> 1;
ans = min(ans, max(dist[i] - e, dist2[i] - e));
}
vector<int> mega;
for (int i = 0; i < n; ++i) {
if (i != b && i != a) {
int e = dist[i] + dist2[i] - dist[b] >> 1;
mega.push_back(dist[i] - e);
}
}
sort(mega.begin(), mega.end());
mega.resize(unique(mega.begin(), mega.end()) - mega.begin());
const int m = mega.size();
vector<int> adj(n), cnt(m);
vector<vector<int>> kids(m);
for (int i = 0; i < n; ++i) {
if (i != b && i != a) {
int e = dist[i] + dist2[i] - dist[b] >> 1;
adj[i] = lower_bound(mega.begin(), mega.end(), dist[i] - e) - mega.begin();
kids[adj[i]].push_back(i);
++cnt[adj[i]];
}
}
for (int i = 0; i < m; ++i) {
if (max(mega[i], dist[b] - mega[i]) == ans) {
int L = 1, R = 1;
for (int j = 0; j < i; ++j) {
L += cnt[j];
}
for (int j = i + 1; j < m; ++j) {
R += cnt[j];
}
if (max({L, R, cnt[i]}) <= n / 2 && 0) {
assert(false);
return -ans;
}
}
}
return ans;
}
}
return mp[{x, y}] = mp[{y, x}] = getDistance(x, y);
}
};
if (sub <= 2) {
vector<int> dist(n);
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, 0);
}
int v = max_element(dist.begin(), dist.end()) - dist.begin();
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, v);
}
int u = max_element(dist.begin(), dist.end()) - dist.begin();
vector<int> dist2(n);
for (int i = 0; i < n; ++i) {
dist2[i] = mygetDistance(i, u);
}
int ans = numeric_limits<int>::max();
for (int i = 0; i < n; ++i) {
int e = dist[i] + dist2[i] - dist[u] >> 1;
ans = min(ans, max(dist[i] - e, dist2[i] - e));
}
return ans;
} else {
vector<int> dist(n);
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, 0);
}
int a = max_element(dist.begin(), dist.end()) - dist.begin();
for (int i = 0; i < n; ++i) {
dist[i] = mygetDistance(i, a);
}
int b = max_element(dist.begin(), dist.end()) - dist.begin();
vector<int> dist2(n);
for (int i = 0; i < n; ++i) {
dist2[i] = mygetDistance(i, b);
}
int ans = numeric_limits<int>::max();
for (int i = 0; i < n; ++i) {
int e = dist[i] + dist2[i] - dist[b] >> 1;
ans = min(ans, max(dist[i] - e, dist2[i] - e));
}
vector<int> mega;
for (int i = 0; i < n; ++i) {
if (i != b && i != a) {
int e = dist[i] + dist2[i] - dist[b] >> 1;
mega.push_back(dist[i] - e);
}
}
sort(mega.begin(), mega.end());
mega.resize(unique(mega.begin(), mega.end()) - mega.begin());
const int m = mega.size();
vector<int> adj(n), cnt(m);
vector<vector<int>> kids(m);
for (int i = 0; i < n; ++i) {
if (i != b && i != a) {
int e = dist[i] + dist2[i] - dist[b] >> 1;
adj[i] = lower_bound(mega.begin(), mega.end(), dist[i] - e) - mega.begin();
kids[adj[i]].push_back(i);
++cnt[adj[i]];
}
}
for (int i = 0; i < m; ++i) {
if (max(mega[i], dist[b] - mega[i]) == ans) {
int L = 1, R = 1;
for (int j = 0; j < i; ++j) {
L += cnt[j];
}
for (int j = i + 1; j < m; ++j) {
R += cnt[j];
}
if (max({L, R, cnt[i]}) <= n / 2) {
assert(false);
return -ans;
}
}
}
return ans;
}
}