QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#180617#5668. Cell Nuclei Detectionarseny_yAC ✓1471ms313716kbC++234.6kb2023-09-16 02:02:372023-09-16 02:02:37

Judging History

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

  • [2023-09-16 02:02:37]
  • 评测
  • 测评结果:AC
  • 用时:1471ms
  • 内存:313716kb
  • [2023-09-16 02:02:37]
  • 提交

answer

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef long double ld;
#define all(a) (a).begin(), (a).end()
#define X first
#define Y second

const ll INF = LLONG_MAX;
const int MOD = 1e9 + 7, MAXN = 2e5 + 1337;
const ld EPS = 1e-6;

void solve();

int main() {
#ifdef LOCAL
    freopen("../main_cpp/input.txt", "r", stdin);
#else
    cin.tie(nullptr)->sync_with_stdio(false);
#endif
    int _ = 1;
    cin >> _;
    while (_--) {
        solve();
    }
}

vector<pair<int, int>> g[MAXN];

struct Rect {
    int x1, y1, x2, y2;

    int sq() {
        if (min(x2 - x1, y2 - y1) <= 0) return 0;
        return (x2 - x1) * (y2 - y1);
    }

    int intersec(const Rect oth) const {
        Rect ans;
        ans.x1 = max(x1, oth.x1);
        ans.y1 = max(y1, oth.y1);
        ans.x2 = min(x2, oth.x2);
        ans.y2 = min(y2, oth.y2);
        return ans.sq();
    }
};

const int MX = 2000 + 10;

vector<int> have[MX][MX];

struct Dinic {
    struct edge {
        int v, to;
        ll f, c;

        edge(int a, int b, ll C, ll d) : v(a), to(b), f(C), c(d) {}

        edge() = default;
    };

    vector<vector<int>> g;
    vector<edge> e;
    vector<int> dist, pnt;
    ll low = 1;
    int s, t, n;

    Dinic() = default;

    Dinic(int n, int ss = -1, int tt = -1) : n(n + 2), s(ss == -1 ? n : ss), t(tt == -1 ? n + 1 : tt), dist(n + 2),
                                             pnt(n + 2, 0) {
        g.assign(n + 2, {});
    }

    void init(int nn, int ss, int tt) {
        n = nn + 2, s = (ss == -1 ? nn : ss), t = (tt == -1 ? nn + 1 : tt);
        dist.resize(n), pnt.assign(n, 0), g.assign(n, {});
    }

    void create(int v, int to, ll c) {
        g[v].push_back(e.size());
        e.emplace_back(v, to, 0, c);
        g[to].push_back(e.size());
        e.emplace_back(to, v, 0, 0);
    }

    bool bfs() {
        fill(dist.begin(), dist.end(), INT_MAX);
        dist[s] = 0;
        queue<int> q;
        q.push(s);
        while (!q.empty()) {
            int v = q.front();
            q.pop();
            for (int i: g[v])
                if (e[i].c - e[i].f >= low && dist[e[i].to] == INT_MAX) {
                    dist[e[i].to] = dist[v] + 1;
                    q.push(e[i].to);
                }
        }
        return dist[t] != INT_MAX;
    }

    ll dfs(int v, ll minC = INF) {
        if (v == t) return minC;
        while (pnt[v] < g[v].size()) {
            int i = g[v][pnt[v]];
            if (dist[e[i].to] == dist[v] + 1 && e[i].c - e[i].f >= low) {
                ll delta = dfs(e[i].to, min(minC, e[i].c - e[i].f));
                if (delta) {
                    e[i].f += delta;
                    e[i ^ 1].f -= delta;
                    return delta;
                }
            }
            ++pnt[v];
        }
        return 0;
    }

    ll flow(int S = -1, int Tt = -1) {
        low = 1;
        s = (S == -1 ? s : S), t = (S == -1 ? t : Tt);
        ll ans = 0;
        while (bfs()) {
            fill(pnt.begin(), pnt.end(), 0);
            while (ll y = dfs(s))
                ans += y;
        }
        return ans;
    }
};

void solve() {
    for (int i = 0; i < MX; ++i) for (int j = 0; j < MX; ++j) have[i][j].clear();
    int n, m;
    cin >> n >> m;
    int sz = n + m + 10;
    for (int i = 1; i <= sz; ++i) {
        g[i].clear();
    }
    vector<Rect> a(n + 1);
    for (int i = 1; i <= n; ++i) {
        cin >> a[i].x1 >> a[i].y1 >> a[i].x2 >> a[i].y2;
        for (int x = a[i].x1; x <= a[i].x2; ++x) {
            for (int y = a[i].y1; y <= a[i].y2; ++y) {
                have[x][y].push_back(i);
            }
        }
    }
    for (int i = 1; i <= m; ++i) {
        int ind = i + n;
        Rect cur;
        cin >> cur.x1 >> cur.y1 >> cur.x2 >> cur.y2;
        for (int x = cur.x1; x <= cur.x2; ++x) {
            for (int y = cur.y1; y <= cur.y2; ++y) {
                for (auto &el : have[x][y]) {
                    if (a[el].intersec(cur) * 2 >= a[el].sq()) {
                        g[el + 1].push_back({ind + 1, 1});
                    }
                }
            }
        }
    }
    int ans = 0;
    for (int i = 2; i <= n + 1; ++i) {
        g[1].push_back({i, 1});
    }
    for (int i = n + 2; i <= n + m + 3; ++i) {
        g[i].emplace_back(sz - 1, 1);
    }
    Dinic dinic(sz + 10);
    for (int i = 1; i <= sz; ++i) {
        sort(all(g[i])), g[i].resize(unique(all(g[i])) - g[i].begin());
        for (auto &[to, c] : g[i]) {
            dinic.create(i, to, c);
        }
    }
    cout << dinic.flow(1, sz - 1) << "\n";
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 41ms
memory: 103056kb

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: 36ms
memory: 103052kb

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: 582ms
memory: 313716kb

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: 424ms
memory: 197012kb

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: 307ms
memory: 148940kb

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: 0
Accepted
time: 1471ms
memory: 302584kb

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:

50000
50000
50000
50000
49600

result:

ok 5 lines

Test #7:

score: 0
Accepted
time: 28ms
memory: 103284kb

input:

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

output:

3

result:

ok single line: '3'