QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#88412#5668. Cell Nuclei DetectionCCSU_WineAC ✓3808ms127712kbC++203.1kb2023-03-16 10:27:272023-03-16 10:28: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-16 10:28:07]
  • 评测
  • 测评结果:AC
  • 用时:3808ms
  • 内存:127712kb
  • [2023-03-16 10:27:27]
  • 提交

answer

#include <queue>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 2010, M = 50010, inf = 1e5, SIZ = 2000;

int n, m;
struct Box{
    int x1, y1, x2, y2, area;
    void input(){ cin >> x1 >> y1 >> x2 >> y2; }
    int get_id(){ return x1 * SIZ + y1; } //返回左下角坐标对应编号
    void get_area(){
        area = (x2 - x1) * (y2 - y1);
    }
}tr[M], de[M];

struct edge{
    int to, nex, w;
}e[M * 200];
int tot, s, t, head[M * 2];
void add(int from, int to, int w){
    e[++ tot].w = w;
    e[tot].to = to;
    e[tot].nex = head[from];
    head[from] = tot;
}

int dep[M * 2];
bool bfs(){
    queue<int>q;
    q.push(s);
    for(int i = s; i <= t; i ++) dep[i] = 0;
    dep[s] = 1;
    while(!q.empty()){
        int u = q.front(); q.pop();
        for(int i = head[u]; i; i = e[i].nex){
            int v = e[i].to;
            if(e[i].w && !dep[v]){
                dep[v] = dep[u] + 1;
                if(v == t) return true;
                q.push(v);
            }
        }
    }
    return dep[t];
}

int dfs(int u, int inflow){
    if(u == t) return inflow;
    int outflow = 0;
    for(int i = head[u]; i && inflow; i = e[i].nex){
        int v = e[i].to;
        if(e[i].w && dep[v] == dep[u] + 1){
            int flow = dfs(v, min(e[i].w, inflow));
            e[i].w -= flow;
            e[i ^ 1].w += flow;
            inflow -= flow;
            outflow += flow;
        }
    }
    if(outflow == 0) dep[u] = 0;
    return outflow;
}

vector<int>g[N * N];//存坐标所有的矩形左上角编号
void init(){
    tot = 1;
    for(int i = 1; i <= m; i ++){
        g[tr[i].get_id()].clear();
    }
    for(int i = s; i <= t; i ++) head[i] = 0;
}


int overlap_area(Box &a, Box &b){//计算线段交
    int x = max(0, min(a.x2, b.x2) - max(a.x1, b.x1));
    int y = max(0, min(a.y2, b.y2) - max(a.y1, b.y1));
    return x * y;
}

void build_edge(int u){
    Box &A = de[u];
    for(int i = max(A.x1 - 2, 0); i <= A.x2 + 2; i ++){
        for(int j = max(A.y1 - 2, 0); j <= A.y2 + 2; j ++){
            int idx = i * SIZ + j;
            for(int id : g[idx]){
                if(overlap_area(A, tr[id]) * 2 >= tr[id].area){
                    add(u, id + n, 1);
                    add(id + n, u, 0);
                }
            }
        }
    }
}

void solve(){
    init();
    cin >> m >> n;
    for(int i = 1; i <= m; i ++){//真实框
        tr[i].input();
        tr[i].get_area();
        g[tr[i].get_id()].push_back(i);
    }
    for(int i = 1; i <= n; i ++){// 搜索框
        de[i].input();
        de[i].get_area();
    }
    
    s = 0, t = n + m + 1;
    for(int i = 1; i <= n; i ++){
        add(s, i, 1);
        add(i, s, 0);
        build_edge(i);
    }
    for(int i = 1; i <= m; i ++){
        add(i + n, t, 1);
        add(t, i + n, 0);
    }

    int ans = 0;
    while(bfs()){
        ans += dfs(s, inf);
    }
    cout << ans << "\n";
}
int main(){
    ios::sync_with_stdio(false);
    cout.tie(NULL);
    int t;
    cin >> t;
    while(t --){
        solve();
    }
    return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 24ms
memory: 98064kb

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: 15ms
memory: 98220kb

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: 2459ms
memory: 121524kb

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: 273ms
memory: 127712kb

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: 254ms
memory: 105164kb

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: 3808ms
memory: 120192kb

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: 16ms
memory: 98060kb

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'