QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#635833#9412. Stop the CastleFangYifanAC ✓1ms4052kbC++176.3kb2024-10-12 21:06:562024-10-12 21:06:57

Judging History

你现在查看的是最新测评结果

  • [2024-10-12 21:06:57]
  • 评测
  • 测评结果:AC
  • 用时:1ms
  • 内存:4052kb
  • [2024-10-12 21:06:56]
  • 提交

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,我给组数据试试?

Details

Tip: Click on the bar to expand more detailed information

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