| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 982003 | vjudge1 | Aliens (IOI16_aliens) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
#define int long long
struct Line {
int a, b;
int idx;
double x(const auto& rhs) { return 1.0 * (rhs.b - b) / (a - rhs.a); };
int cal(const auto& x) { return a * x + b; }
};
struct CHT {
Line s[100'000 + 10];
int st = 1, ed = 0;
void add(Line nline) {
while (ed - st + 1 >= 2 && s[ed - 1].x(s[ed]) > s[ed - 1].x(nline)) ed -= 1;
s[++ed] = nline;
}
Line get(int x) {
while (ed - st + 1 >= 2 && s[st].cal(x) > s[st + 1].cal(x)) st += 1;
return s[st];
}
};
int take_photos(int n, int m, int k, vector<int> row, vector<int> column) {
using pii = pair<int, int>;
vector<pii> p(n);
for (int i = 0; i < n; ++i) {
if (row[i] > column[i]) swap(row[i], column[i]);
p[i] = {row[i], column[i]};
}
auto dis = [&](pii a) { return abs(a.first - a.second) + 1; };
{ //sort_the_points
sort(p.begin(), p.end(), [&](const auto& a, const auto& b) {
return make_pair(a.first, -dis(a)) < make_pair(b.first, -dis(b));
});
vector<bool> mk(n, true);
pii point = {-1, -1};
for (int i = 0; i < n; ++i) {
if (point.first <= p[i].first && p[i].second <= point.second) mk[i] = false;
else point = p[i];
}
int cnt = 0;
for (int i = 0; i < n; ++i) {
cnt += mk[i];
if (mk[i]) p.push_back(p[i]);
}
p.erase(p.begin(), p.begin() + n);
p.push_back({-1, -1});
k = min(k, n = cnt);
sort(p.begin(), p.end(), [&](const auto& a, const auto& b) { return a.first < b.first; });
}
auto cost = [&](const auto& C) {
vector<int> f(n + 1, 1'000'000'000'000'000), g(n + 1);
CHT hull;
f[0] = 0;
for (int i = 1; i <= n; ++i) {
int pl = -p[i].first + 1;
int pre = max(0ll, p[i - 1].second + pl); pre *= pre;
Line nline = {2 * pl, f[i - 1] + pl * pl - pre, i - 1};
hull.add(nline);
Line opt = hull.get(p[i].second);
g[i] = g[opt.idx] + 1;
f[i] = opt.cal(p[i].second) + 1ll * p[i].second * p[i].second + C;
}
return make_pair(f[n], g[n]);
};
int l = 0, r = 1'000'000'000'000'000'000;
int answer = 1'000'000'000'000'000'000;
while (l <= r) {
int mid = l + r >> 1;
auto [val, cnt] = cost(mid);
if (cnt <= k) {
r = mid - 1;
answer = val - k * mid;
} else l = mid + 1;
}
return answer;
}
signed main() {
cin.tie(0)->sync_with_stdio(0);
int n, m, k; cin >> n >> m >> k;
vector<int> r(n), c(n);
for (int i = 0; i < n; ++i) cin >> r[i] >> c[i];
cout << take_photos(n, m, k, r, c) << "\n";
}