QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#84133#5668. Cell Nuclei DetectionHOLIC#AC ✓4154ms255080kbC++203.7kb2023-03-05 18:27:042023-03-05 18:27:07

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-03-05 18:27:07]
  • 评测
  • 测评结果:AC
  • 用时:4154ms
  • 内存:255080kb
  • [2023-03-05 18:27:04]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
#define endl "\n"
const int N = 2e5 + 1009;
const int M = 8 * N * 40 + N;
int n, m;
namespace Dinic{
    const int inf = 0x3f3f3f3f;
    int s, t,  d[N], cpy[N], head[N], nexts[M], edge[M], ver[M], cnt = 1;
    void add(int x, int y, int z) {
        ver[++ cnt] = y, nexts[cnt] = head[x], head[x] = cpy[x] = cnt, edge[cnt] = z;
    }
    void clear() {
        for(int i = 1; i <= t; i ++) head[i] = cpy[i] = 0;
        cnt = 1;
    }
    bool bfs() {
        memset(d, 0, sizeof(d));
        queue<int> q;
        q.push(s);  d[s] = 1;
        while(!q.empty()) {
            int x = q.front(); q.pop();
            for(int i = head[x]; i; i = nexts[i]) {
                int v = ver[i], w = edge[i];
                if(w > 0 && !d[v]) d[v] = d[x] + 1, q.push(v);
            }
        }
        if(!d[t]) return false;
        return true;
    }
    int dfs(int x, int flow) {
        if(x == t || flow == 0) return flow;
        int maxflow = 0, f;
        for(int &i = head[x]; i; i = nexts[i]) {
            int v = ver[i], w = edge[i];
            if(d[v] == d[x] + 1 && w) {
                f = dfs(v, min(flow, w));
                if(f > 0) {
                    edge[i] -= f;
                    edge[i ^ 1] += f;
                    flow -= f; maxflow += f;
                    if(!flow) return maxflow;
                }
            }
        }
        return maxflow;
    }
    int dinic() {
        int maxflow = 0, flow;
        while(bfs()) {
            maxflow += dfs(s, inf);
            for(int i = 1; i <= t; i ++)
                head[i] = cpy[i];
        }
        return maxflow;
    }
}
vector<int> ver[2005][2005], Ver[2005][2005];
int xx[50005], yy[50005], xxx[50005], yyy[50005];
int XX[50005], YY[50005], XXX[50005], YYY[50005];
void work() {
    cin >> n >> m;
    Dinic::clear();
    Dinic::s = n + m + 1;
    Dinic::t = n + m + 2;
    for(int i = 0; i <= 2000; i ++)
        for(int j = 0; j <= 2000; j ++)
            ver[i][j].clear(), Ver[i][j].clear();

    for(int i = 1; i <= n; i ++) {
        cin >> xx[i] >> yy[i] >> xxx[i] >> yyy[i];
        for(int I = xx[i]; I <= xxx[i]; I ++) {
            for(int J = yy[i]; J <= yyy[i]; J ++) {
                ver[I][J].push_back(i);
            }
        }
        Dinic::add(i + m, Dinic::t, 1); Dinic::add(Dinic::t, i + m, 0);
    }

    for(int i = 1; i <= m; i ++) {
        cin >> XX[i] >> YY[i] >> XXX[i] >> YYY[i];
        Dinic::add(Dinic::s, i, 1); Dinic::add(i, Dinic::s, 0);
    }
    for(int i = 1; i <= m; i ++) {
        map<int,int> vis;
        int x = XX[i], y = YY[i], X = XXX[i], Y = YYY[i];
        auto check = [&](int id) {

            int lx = xx[id], ly = yy[id];
            int rx = xxx[id], ry = yyy[id];
            int s = (rx - lx) * (ry - ly);

            if(lx < x) lx = x;
            if(rx > X) rx = X;
            if(ly < y) ly = y;
            if(ry > Y) ry = Y;
            if(rx - lx <= 0) return false;
            if(ry - ly <= 0) return false;
            int ss = (rx - lx) * (ry - ly);

            return ss * 2 >= s;
        };
        for(int I = x; I <= X; I ++) {
            for(int J = y; J <= Y; J ++) {
                for(auto id : ver[I][J]) {
                    if(check(id)) {
                        if(vis[id]) continue;
                        vis[id] = 1;
                        Dinic::add(i, id + m, 1);
                        Dinic::add(id + m, i, 0);
                    }
                }
            }
        }
    }
    cout << Dinic::dinic() << endl;
}
signed main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int Case;
    cin >> Case;
    while(Case --) work();
    return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 87ms
memory: 192644kb

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: 76ms
memory: 192588kb

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: 1206ms
memory: 243244kb

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: 492ms
memory: 255080kb

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: 434ms
memory: 221916kb

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: 4154ms
memory: 219684kb

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: 47ms
memory: 192716kb

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'