QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#123413#5668. Cell Nuclei DetectionGenshinImpactsFault#AC ✓593ms135992kbC++143.0kb2023-07-12 15:34:142023-07-12 15:34:17

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-07-12 15:34:17]
  • 评测
  • 测评结果:AC
  • 用时:593ms
  • 内存:135992kb
  • [2023-07-12 15:34:14]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
// #define int long long
#define MOD ((1ll << 31) - 1)
vector<array<int, 3>> pos[2010][2010];
vector<pair<int, int>> cls;

const int N = 100010;
const int INF = N * 10;
int s, t, ans;
int cnt, hed[N], to[N * 100], nxt[N * 100], val[N * 100];
int cur[N], dis[N], g[N];

void init(int n) {
    s = n + 1, t = n + 2, ans = 0;
    for(int i = 0; i <= t; i++) hed[i] = cur[i] = dis[i] = g[i] = 0;
    for(int i = 0; i <= cnt; i++) to[i] = val[i] = nxt[i] = 0;
    cnt = 1;
}
void addedge(int x, int y, int w) {
    // cout << ">>> " << x << " " << y << " " << w << "\n";
    to[++cnt] = y, nxt[cnt] = hed[x], hed[x] = cnt, val[cnt] = w;
    to[++cnt] = x, nxt[cnt] = hed[y], hed[y] = cnt, val[cnt] = 0;
}
int SAP(int u, int flow) {
    if(u == t) return flow;
    int v = flow;
    for(; cur[u]; cur[u] = nxt[cur[u]]) {
        if(dis[to[cur[u]]] == dis[u] - 1 && val[cur[u]]) {
            int f = SAP(to[cur[u]], min(v, val[cur[u]]));
            val[cur[u]] -= f, val[cur[u] ^ 1] += f;
            v -= f;
            if(!v) return flow;
        }
    }
    if(!--g[dis[u]++]) dis[s] = t + 1;
    ++g[dis[u]], cur[u] = hed[u]; return flow - v;
}
void Do() {
    int m, n;
    cin >> m >> n;
    init(n + m);
    for(int i = 1; i <= m; ++ i) {
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        // cerr << x1 << " " << y1 << " " << x2 << " " << y2 << endl;
        pos[x1][y1].push_back({x2, y2, i});
        cls.push_back({x1, y1});
    }
    for(int i = 1; i <= n; ++ i) {
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        for(int x = max(x1 - 3, 0); x <= x2 - 1; ++ x) {
            for(int y = max(y1 - 3, 0); y <= y2 - 1; ++ y) {
                for(auto [x3, y3, id] : pos[x][y]) {
                    // cerr << x3 << " " << y3 << " " << id << endl;
                    // x y x3 y3
                    // x1 y1 x2 y2;
                    int X1 = max(x, x1), Y1 = max(y, y1);
                    int X2 = min(x3, x2), Y2 = min(y3, y2);
                    int S = (X2 - X1) * (Y2 - Y1);
                    if(X1 > X2 || Y1 > Y2) S = 0;
                    int pS = (x3 - x) * (y3 - y);
                    // cerr << pS << " " << S << endl;
                    if(S * 2 >= pS) {
                        // add edge id -> i + m
                        // cerr << id << " " << i + m << endl;
                        addedge(id, i + m, 1);
                    }
                }
            }
        }
    }
    // exit(0);

    for(int i = 1; i <= m; i++) addedge(s, i, 1);
    for(int i = 1; i <= n; i++) addedge(i + m, t, 1);
    for(int i = 1; i <= t; i++) cur[i] = hed[i];
    for(g[0] = t; dis[s] <= t; ans += SAP(s, INF));
    cout << ans << "\n";

    for(auto [x, y] : cls) {
        pos[x][y].clear();
    }
}
signed main() {
	// ios::sync_with_stdio(0);
    // cin.tie(0), cout.tie(0);
	int T;
    cin >> T;
    while(T --) {
        Do();
    }
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 15ms
memory: 102868kb

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: 8ms
memory: 106076kb

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: 285ms
memory: 129424kb

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: 304ms
memory: 135992kb

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: 312ms
memory: 113512kb

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: 593ms
memory: 125784kb

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: 10ms
memory: 102812kb

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'