QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#578387#5668. Cell Nuclei DetectionCelestialCoder#TL 4266ms30124kbC++203.4kb2024-09-20 18:47:082024-09-20 18:47:08

Judging History

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

  • [2024-09-20 18:47:08]
  • 评测
  • 测评结果:TL
  • 用时:4266ms
  • 内存:30124kb
  • [2024-09-20 18:47:08]
  • 提交

answer

#include <bits/stdc++.h>
#ifdef SHARAELONG
#include "../../cpp-header/debug.hpp"
#endif
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;

struct HopcroftKarp {
    int n, m;
    vector<vector<int>> g;
    vector<int> dst, le, ri;
    vector<char> visit, track;
    HopcroftKarp(int n, int m): n(n), m(m), g(n), dst(n), le(n, -1), ri(m, -1), visit(n), track(n+m) {}

    void add_edge(int s, int e) {
        g[s].push_back(e);
    }

    bool bfs() {
        bool res = false;
        queue<int> que;
        fill(dst.begin(), dst.end(), 0);
        for (int i=0; i<n; ++i) {
            if (le[i] == -1) que.push(i), dst[i] = 1;
        }
        while (!que.empty()) {
            int v = que.front();
            que.pop();
            for (auto i: g[v]) {
                if (ri[i] == -1) res = true;
                else if (!dst[ri[i]]) dst[ri[i]] = dst[v]+1, que.push(ri[i]);
            }
        }
        return res;
    }

    bool dfs(int v) {
        if (visit[v]) return false;
        visit[v] = 1;
        for (auto i: g[v]) {
            if (ri[i] == -1 || !visit[ri[i]] && dst[ri[i]] == dst[v]+1 && dfs(ri[i])) {
                le[v] = i;
                ri[i] = v;
                return true;
            }
        }
        return false;
    }

    int maximum_matching() {
        int res = 0;
        fill(le.begin(), le.end(), -1);
        fill(ri.begin(), ri.end(), -1);
        while (bfs()) {
            fill(visit.begin(), visit.end(), 0);
            for (int i=0; i<n; ++i) {
                if (le[i] == -1) res += dfs(i);
            }
        }
        return res;
    }
};

void solve() {
    int m, n;
    cin >> m >> n;
    vector<array<int, 4>> gr(m), de(n);
    for (int i=0; i<m; ++i) {
        for (int j=0; j<4; ++j) cin >> gr[i][j];
    }
    for (int i=0; i<n; ++i) {
        for (int j=0; j<4; ++j) cin >> de[i][j];
    }

    map<pii, vector<int>> mp;
    for (int i=0; i<n; ++i) {
        auto&[x1, y1, x2, y2] = de[i];
        mp[{ x1,y1 }].push_back(i);
        mp[{ x1,y2 }].push_back(i);
        mp[{ x2,y1 }].push_back(i);
        mp[{ x2,y2 }].push_back(i);
    }
    
    auto inter = [](const array<int, 4>& g, const array<int, 4>& d) {
        if (min(g[2], d[2]) < max(g[0], d[0])) return false;
        if (min(g[3], d[3]) < max(g[1], d[1])) return false;
        int sz = (min(g[2], d[2]) - max(g[0], d[0])) * (min(g[3], d[3]) - max(g[1], d[1]));
        return (2*sz >= (g[2]-g[0])*(g[3]-g[1]));
    };

    auto search = [&](int i, pii c, set<int>& s) {
        for (int x=c.first-4; x<=c.first+4; ++x) {
            for (int y=c.second-4; y<=c.second+4; ++y) {
                if (mp.find({ x,y }) != mp.end()) {
                    auto& v = mp[{ x,y }];
                    for (int j: v) {
                        if (inter(gr[i], de[j])) s.insert(j);
                    }
                }
            }
        }
    };

    HopcroftKarp g(m, n);
    for (int i=0; i<m; ++i) {
        auto[x1, y1, x2, y2] = gr[i];
        set<int> s;
        search(i, { x1,y1 }, s);
        search(i, { x1,y2 }, s);
        search(i, { x2,y1 }, s);
        search(i, { x2,y2 }, s);
        for (int j: s) g.add_edge(i, j);
    }
    cout << g.maximum_matching() << '\n';
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int cases;
    cin >> cases;
    while (cases--) solve();
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 3564kb

input:

3
2 2
1 1 3 3
3 3 5 5
2 2 4 4
4 4 6 6
2 3
1 1 3 3
3 3 5 5
1 3 3 5
2 1 4 5
3 1 5 3
3 3
1 1 2 2
2 2 3 3
3 3 4 4
1 1 3 3
2 2 4 4
3 3 5 5

output:

0
1
3

result:

ok 3 lines

Test #2:

score: 0
Accepted
time: 0ms
memory: 3488kb

input:

3
2 2
1 1 3 3
3 3 5 5
2 2 4 4
4 4 6 6
2 3
1 1 3 3
3 3 5 5
1 3 3 5
2 1 4 5
3 1 5 3
3 3
1 1 2 2
2 2 3 3
3 3 4 4
1 1 3 3
2 2 4 4
3 3 5 5

output:

0
1
3

result:

ok 3 lines

Test #3:

score: 0
Accepted
time: 3703ms
memory: 16676kb

input:

5
50000 50000
0 0 4 4
4 0 8 4
8 0 12 4
12 0 16 4
16 0 20 4
20 0 24 4
24 0 28 4
28 0 32 4
32 0 36 4
36 0 40 4
40 0 44 4
44 0 48 4
48 0 52 4
52 0 56 4
56 0 60 4
60 0 64 4
64 0 68 4
68 0 72 4
72 0 76 4
76 0 80 4
80 0 84 4
84 0 88 4
88 0 92 4
92 0 96 4
96 0 100 4
100 0 104 4
104 0 108 4
108 0 112 4
112 ...

output:

50000
50000
0
50000
3150

result:

ok 5 lines

Test #4:

score: 0
Accepted
time: 3672ms
memory: 30124kb

input:

5
50000 50000
0 0 1 1
1 0 2 1
2 0 3 1
3 0 4 1
4 0 5 1
5 0 6 1
6 0 7 1
7 0 8 1
8 0 9 1
9 0 10 1
10 0 11 1
11 0 12 1
12 0 13 1
13 0 14 1
14 0 15 1
15 0 16 1
16 0 17 1
17 0 18 1
18 0 19 1
19 0 20 1
20 0 21 1
21 0 22 1
22 0 23 1
23 0 24 1
24 0 25 1
25 0 26 1
26 0 27 1
27 0 28 1
28 0 29 1
29 0 30 1
30 0 ...

output:

50000
25050
12500
16000
8000

result:

ok 5 lines

Test #5:

score: 0
Accepted
time: 4266ms
memory: 29516kb

input:

5
50000 50000
0 0 2 4
4 0 7 1
8 0 10 1
12 0 15 3
16 0 19 1
20 0 22 2
24 0 26 4
28 0 30 4
32 0 36 3
36 0 40 1
40 0 44 1
44 0 47 2
48 0 49 3
52 0 54 1
56 0 59 4
60 0 64 3
64 0 68 3
68 0 70 1
72 0 76 4
76 0 80 3
80 0 84 4
84 0 87 2
88 0 90 1
92 0 94 4
96 0 98 1
100 0 104 1
104 0 107 2
108 0 110 4
112 0...

output:

10594
10779
10618
10381
10779

result:

ok 5 lines

Test #6:

score: -100
Time Limit Exceeded

input:

5
50000 50000
0 0 4 4
1 0 5 4
2 0 6 4
3 0 7 4
4 0 8 4
5 0 9 4
6 0 10 4
7 0 11 4
8 0 12 4
9 0 13 4
10 0 14 4
11 0 15 4
12 0 16 4
13 0 17 4
14 0 18 4
15 0 19 4
16 0 20 4
17 0 21 4
18 0 22 4
19 0 23 4
20 0 24 4
21 0 25 4
22 0 26 4
23 0 27 4
24 0 28 4
25 0 29 4
26 0 30 4
27 0 31 4
28 0 32 4
29 0 33 4
30...

output:


result: