QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#635833 | #9412. Stop the Castle | FangYifan | AC ✓ | 1ms | 4052kb | C++17 | 6.3kb | 2024-10-12 21:06:56 | 2024-10-12 21:06:57 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using db = double;
constexpr int mod = 998244353;
constexpr int N = 1e6 + 10;
template <class T>
struct Flow {
T inf;
int n;
struct Edge {
int v;
T cap;
Edge(int v, T cap) : v(v), cap(cap) {}
};
vector<Edge> e; // e[i] 存编号为 i 的边所指向的节点 v / 流量 cap
vector<vector<int>> g; // g[u] 存储了所有以 u 为起点的边 u -> v 的编号 i
vector<int> cur, deep;
Flow(int n) : inf(numeric_limits<T>::max()), n(n), g(n + 1) {}
bool bfs(int s, int t) {
deep.assign(n + 1, -1);
queue<int> q;
deep[s] = 0;
q.push(s);
while (!q.empty()) {
int u = q.front();
q.pop();
for (int i : g[u]) {
auto [v, c] = e[i];
if (c > 0 && deep[v] == -1) {
deep[v] = deep[u] + 1;
if (v == t) return true;
q.push(v);
}
}
}
return false;
}
T dfs(int u, int t, T flow) {
if (u == t) return flow;
T remain = flow;
for (int &i = cur[u]; i < (int)g[u].size(); i++) {
int j = g[u][i];
auto [v, c] = e[j];
if (c > 0 && deep[v] == deep[u] + 1) {
T cost = dfs(v, t, min(remain, c));
e[j].cap -= cost;
e[j ^ 1].cap += cost;
remain -= cost;
if (remain == 0) return flow;
}
}
return flow - remain;
}
void add(int u, int v, T cap) {
g[u].push_back(e.size());
e.emplace_back(v, cap);
g[v].push_back(e.size());
e.emplace_back(u, 0);
}
// 最大流: 从源点 s 到汇点 t 的最大流量
T max_flow(int s, int t) {
T ans = 0;
while (bfs(s, t)) {
cur.assign(n + 1, 0);
ans += dfs(s, t, inf);
}
return ans;
}
// 最小割: 分割源点 s 到汇点 t 容量和最小的割边集合,最小割 = 最大流
vector<int> min_cut_ST() { // 标记最小割 (S, T)
vector<int> c(n + 1);
for (int i = 1; i <= n; i++) {
c[i] = (deep[i] != -1); // 源点 s 可到达节点 i
}
return c;
}
struct _Edge {
int u, v;
T cap;
T flow;
};
vector<_Edge> edges() {
vector<_Edge> t;
for (int i = 0; i < e.size(); i += 2) {
_Edge x;
x.u = e[i + 1].v;
x.v = e[i].v;
x.cap = e[i].cap + e[i + 1].cap;
x.flow = e[i + 1].cap;
t.push_back(x);
}
return t;
}
};
void solve() {
int n; cin >> n;
map<int, vector<int>> rp, cp;
vector<int> r(n + 1), c(n + 1);
// 每一行 / 每一列 维护坐标
for (int i = 1; i <= n; i++) {
cin >> r[i] >> c[i];
rp[r[i]].push_back(c[i]);
cp[c[i]].push_back(r[i]);
}
int m; cin >> m;
for (int i = 1; i <= m; i++) {
int x, y; cin >> x >> y;
rp[x].push_back(y);
cp[y].push_back(x);
}
for (auto &[x, y] : rp) sort(y.begin(), y.end());
for (auto &[x, y] : cp) sort(y.begin(), y.end());
// 枚举两车位置关系,如果相邻 则无解,如果碰撞 则创建新的碰撞线段
vector<array<int, 3>> rline, cline;
for (int i = 1; i <= n; i++) {
for (int j = i + 1; j <= n; j++) {
if (r[i] == r[j]) {
if (abs(c[i] - c[j]) == 1) {
cout << "-1\n";
return;
}
auto &v = rp[r[i]];
auto it = lower_bound(v.begin(), v.end(), c[i]);
if ((it != v.begin() && *prev(it) == c[j]) || (it != v.end() && *next(it) == c[j])) {
rline.push_back({r[i], min(c[i], c[j]), max(c[i], c[j])});
}
}
if (c[i] == c[j]) {
if (abs(r[i] - r[j]) == 1) {
cout << "-1\n";
return;
}
auto &v = cp[c[i]];
auto it = lower_bound(v.begin(), v.end(), r[i]);
if ((it != v.begin() && *prev(it) == r[j]) || (it != v.end() && *next(it) == r[j])) {
cline.push_back({c[i], min(r[i], r[j]), max(r[i], r[j])});
}
}
}
}
// 对于行列碰撞的交点 (x, y),建立 行碰撞线段 -> 列碰撞线段 的有向边,容量为 1,跑最大流即可
int numr = rline.size(), numc = cline.size();
Flow<int> f(numr + numc + 2);
int S = numr + numc + 1, T = numr + numc + 2;
for (int i = 1; i <= numr; i++) f.add(S, i, 1);
for (int i = 1; i <= numc; i++) f.add(numr + i, T, 1);
for (int i = 1; i <= numr; i++) {
for (int j = 1; j <= numc; j++) {
auto [x, y1, y2] = rline[i - 1];
auto [y, x1, x2] = cline[j - 1];
if (x1 < x && x < x2 && y1 < y && y < y2) {
// cerr << i << ' ' << j << "\n";
// cerr << x << ' ' << y1 << ' ' << y2 << ' ' << y << ' ' << x1 << ' ' << x2 << "\n";
f.add(i, numr + j, 1);
}
}
}
// 处理行列碰撞线段 and 对应的放置方案
f.max_flow(S, T);
auto e = f.edges();
set<pair<int, int>> vis;
vector<pair<int, int>> ans;
for (auto [u, v, cap, flow] : e) {
if (u == S || v == T) continue;
if (flow == 1) {
auto [x, y1, y2] = rline[u - 1];
auto [y, x1, x2] = cline[v - numr - 1];
ans.emplace_back(x, y);
vis.insert({x, y1 + 1});
vis.insert({x1 + 1, y});
}
}
for (auto [x, y1, y2] : rline) {
if (!vis.count({x, y1 + 1})) {
ans.emplace_back(x, y1 + 1);
}
}
for (auto [y, x1, x2] : cline) {
if (!vis.count({x1 + 1, y})) {
ans.emplace_back(x1 + 1, y);
}
}
cout << ans.size() << "\n";
for (auto [x, y] : ans) cout << x << ' ' << y << "\n";
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
int tt = 1;
cin >> tt;
while (tt--) {
solve();
}
return 0;
}
这程序好像有点Bug,我给组数据试试?
详细
Test #1:
score: 100
Accepted
time: 1ms
memory: 3584kb
input:
4 7 1 3 6 6 4 7 2 1 6 3 4 1 2 6 3 3 4 6 4 3 1 2 1 1 2 2 0 3 1 1 1 3 3 3 1 1 2 3 1 1 1 3 2 3 0
output:
2 4 6 2 3 0 1 2 3 -1
result:
ok ok 4 cases (4 test cases)
Test #2:
score: 0
Accepted
time: 1ms
memory: 3588kb
input:
21 11 3 5 18 6 19 12 27 48 28 38 30 12 33 18 42 18 43 38 45 46 48 34 3 2 6 24 4 41 45 15 11 6 15 27 16 9 16 14 16 48 19 26 25 12 27 26 28 4 31 40 32 6 33 50 37 50 46 11 50 29 17 1 49 4 9 9 15 9 22 11 31 11 46 14 28 17 5 17 49 18 43 20 31 30 46 33 11 33 33 34 15 36 6 47 10 2 16 46 36 30 11 7 34 7 41 ...
output:
3 20 12 29 38 34 18 5 16 10 16 15 12 6 20 26 34 50 0 1 16 10 0 1 43 22 5 1 3 1 13 33 10 44 45 42 44 0 5 27 41 29 26 44 4 8 1 21 15 1 32 9 0 0 0 0 6 23 10 35 34 12 11 23 44 29 21 24 46 0 3 20 30 43 25 31 17 0 -1 3 16 36 25 7 17 39 6 5 5 8 9 6 4 7 3 3 10 2 11
result:
ok ok 21 cases (21 test cases)
Test #3:
score: 0
Accepted
time: 1ms
memory: 3636kb
input:
2 200 7 52 9 160 10 4 12 61 18 436 19 430 20 499 24 139 25 416 29 53 33 153 35 275 35 310 37 25 38 477 39 349 42 120 44 158 45 330 49 486 50 204 51 67 52 138 52 305 56 139 63 132 66 4 67 327 70 484 71 67 72 308 74 218 76 479 81 139 82 90 86 201 87 335 91 35 91 220 92 496 94 29 94 436 96 359 97 299 1...
output:
46 52 139 91 61 94 153 126 138 126 275 132 67 154 160 154 496 185 147 187 467 199 432 224 35 270 356 277 305 306 374 311 177 334 186 349 112 352 61 393 307 35 276 187 433 224 393 260 223 261 468 274 67 277 189 311 455 390 308 440 179 453 251 11 4 19 436 21 499 38 25 39 477 52 67 57 139 138 429 240 3...
result:
ok ok 2 cases (2 test cases)
Test #4:
score: 0
Accepted
time: 1ms
memory: 3888kb
input:
20 13 89287395 441913071 89287395 590314987 142917372 582720391 142917372 598813561 232930851 582720391 263974835 468188689 263974835 490702144 543529670 152471388 543529670 219371706 647691663 598813561 772865038 598813561 773363571 482933265 773363571 561453301 8 128947560 120560639 264980960 4742...
output:
8 89287395 441913072 142917372 582720392 263974835 468188690 543529670 152471389 773363571 482933266 142917373 582720391 142917373 598813561 647691664 598813561 7 808969359 818037531 808969359 386523246 808969359 891229782 59832704 871951216 92745430 358274214 469117951 762966373 808969360 354711322...
result:
ok ok 20 cases (20 test cases)
Test #5:
score: 0
Accepted
time: 1ms
memory: 3600kb
input:
2 184 35763790 435860426 152681166 443483111 261932524 745277126 261932524 787982418 314305867 342102981 314305867 522593781 314305867 747023891 314305867 811404535 314305867 816914009 314305867 857896659 319901983 797458033 321347896 342102981 332550928 902179526 343207116 914408792 350635050 15515...
output:
160 314305867 435860426 314305867 787982418 350635050 443483111 350635050 747023891 350635050 811404535 407380694 342102981 469496529 159320601 469496529 438907614 469496529 526518374 469496529 797458033 469496529 816914009 469496529 837963933 495537855 732085440 496813081 342102981 496813081 787982...
result:
ok ok 2 cases (2 test cases)
Test #6:
score: 0
Accepted
time: 1ms
memory: 3660kb
input:
2 193 78368109 329329995 87932514 657474465 87932514 903492853 94608574 845872655 103602204 88649154 124431127 405860533 145262472 108198774 145262472 907127202 154400439 136367504 165270984 773451933 184157521 409828915 197728547 291882089 197728547 609623645 211750768 843867422 246178069 108396139...
output:
145 145262472 903492853 197728547 405860533 246178069 291882089 291242118 564718673 365329221 136367504 365329221 211655411 365329221 329329995 365329221 409828915 365329221 657489855 365329221 773451933 365329221 845872655 365329221 909607734 382558552 843867422 385064292 480642952 391226954 502303...
result:
ok ok 2 cases (2 test cases)
Test #7:
score: 0
Accepted
time: 0ms
memory: 3660kb
input:
2 192 1075648 579875717 1075648 937648715 1393060 191616432 3957616 541362608 3957616 971879257 5415801 541362608 5887448 541362608 5887448 624842533 25767120 610618164 26193206 214214028 26193206 231445627 26727662 374314271 29550509 755417826 29550509 917592925 30033126 610618164 31134588 58188305...
output:
131 31134588 191616432 31134588 541362608 33207011 231445627 46425938 755417826 223207847 374314271 223207847 624842533 223207847 762642201 235369723 288390629 250896470 264208611 261816954 741119775 295641139 541362608 340743563 381919591 341366249 231445627 360971009 755417826 417168695 231445627 ...
result:
ok ok 2 cases (2 test cases)
Test #8:
score: 0
Accepted
time: 1ms
memory: 3600kb
input:
20 14 378192988 147499872 378192988 837331700 449924898 147499872 449924898 837331700 592684636 147499872 592684636 837331700 783332532 147499872 783332532 837331700 378192988 197547597 783332532 197547597 378192988 343857831 783332532 343857831 378192988 576579017 783332532 576579017 2 449924898 19...
output:
15 378192988 147499873 378192988 576579018 783332532 147499873 783332532 576579018 378192988 197547598 783332532 197547598 378192988 343857832 783332532 343857832 378192989 147499872 378192989 837331700 449924899 147499872 449924899 837331700 592684637 147499872 592684637 837331700 378192989 5765790...
result:
ok ok 20 cases (20 test cases)
Test #9:
score: 0
Accepted
time: 1ms
memory: 4052kb
input:
2 200 8623893 3649468 8623893 989134714 24820918 3649468 24820918 989134714 43705451 3649468 43705451 989134714 45500146 3649468 45500146 989134714 58265855 3649468 58265855 989134714 112263159 3649468 112263159 989134714 126496650 3649468 126496650 989134714 137574156 3649468 137574156 989134714 14...
output:
249 24820918 6255307 43705451 11610914 45500146 23366504 58265855 28938721 112263159 31914492 126496650 36029445 137574156 66226935 140089445 84152868 156036789 132710904 170794245 149083019 176576709 174305881 186774359 180413073 194114590 195080121 222368048 231772824 237249112 234663336 250839457...
result:
ok ok 2 cases (2 test cases)
Test #10:
score: 0
Accepted
time: 1ms
memory: 3712kb
input:
2 200 2691939 27674656 2691939 984614319 30225807 27674656 30225807 984614319 39388076 27674656 39388076 984614319 55883389 27674656 55883389 984614319 77258752 27674656 77258752 984614319 106874917 27674656 106874917 984614319 144383292 27674656 144383292 984614319 145360776 27674656 145360776 9846...
output:
216 279200011 206375958 347252805 243902652 446662965 257900116 496057717 293316179 514655989 363531956 523731402 376231503 538628510 383623456 562810007 390555088 630132600 441004197 631858791 490065400 643299496 540189713 754008138 545422987 782710197 565468975 2691939 27674657 2691939 983896590 9...
result:
ok ok 2 cases (2 test cases)
Test #11:
score: 0
Accepted
time: 1ms
memory: 3684kb
input:
2 200 1005937 7 3412691 7 5764908 7 7897785 7 8061992 7 12801501 7 18529696 7 22574605 7 25203351 7 34574958 7 34998216 7 40928451 7 49792193 7 55962320 7 59632864 7 62353373 7 62611196 7 65346591 7 71459320 7 71550121 7 72215989 7 82584263 7 86119549 7 96204862 7 96940713 7 97693112 7 102067050 7 1...
output:
199 1005938 7 3412692 7 5764909 7 7897786 7 8061993 7 12801502 7 18529697 7 22574606 7 25203352 7 34574959 7 34998217 7 40928452 7 49792194 7 55962321 7 59632865 7 62353374 7 62611197 7 65346592 7 71459321 7 71550122 7 72215990 7 82584264 7 86119550 7 96204863 7 96940714 7 97693113 7 102067051 7 105...
result:
ok ok 2 cases (2 test cases)
Extra Test:
score: 0
Extra Test Passed