QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#103485#5668. Cell Nuclei Detectioncardinal_city#AC ✓1982ms15020kbC++173.4kb2023-05-06 05:03:202023-05-06 05:03:41

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-05-06 05:03:41]
  • 评测
  • 测评结果:AC
  • 用时:1982ms
  • 内存:15020kb
  • [2023-05-06 05:03:20]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef long double ld;
typedef vector<int> vi;
typedef pair<int, int> pii;
#define all(x) (x).begin(), (x).end()
#define sz(x) (int)x.size();
#define rep(i,j,k) for (int i = j; i < (k); i++)
#define in(mp, v) (mp.find(v) != mp.end())
#define smx(a, b) a = max(a, b)
#define smn(a, b) a = min(a, b)
#define pb push_back
#define endl '\n'

const ll MOD = 1e9 + 7;
const ld EPS = 1e-9;

bool dfs(int a, int L, vector<vi> &g, vi& btoa, vi& A, vi &B) {
    if(A[a] != L) return 0;
    A[a] = -1;
    for(int b : g[a]) {
        if(B[b] == L+1) {
            B[b] = 0;
            if(btoa[b] == -1 || dfs(btoa[b], L+1, g, btoa, A,B)) {
                return btoa[b] = a, 1;
            }
        }
    }
    return 0;
}

int hopcroftKarp(vector<vi> &g, vi& btoa) {
    int res = 0;
    vi A(g.size()), B(btoa.size()), cur, next;
    for(;;) {
        fill(all(A), 0);
        fill(all(B), 0);
        cur.clear();
        for (int a : btoa) if(a != -1) A[a] = -1;
        for(int a = 0; a < g.size(); a++) {
            if(A[a] == 0) cur.pb(a);
        }
        for(int lay = 1;; lay++) {
            bool islast = 0;
            next.clear();
            for(int a : cur) for(int b : g[a]) {
                if(btoa[b] == -1) {
                    B[b] = lay;
                    islast = 1;
                } else if(btoa[b] != a && !B[b]) {
                    B[b] = lay;
                    next.pb(btoa[b]);
                }
            }
            if(islast) break;
            if(next.empty()) return res;
            for (int a : next) A[a] = lay;
            cur.swap(next);
        }
        for(int a = 0; a < g.size(); a++) {
            res += dfs(a, 0, g, btoa, A, B);
        }
    }
}

struct rect {
    int x1, y1, x2, y2, id;
};

rect inter(rect a, rect b) {
    return {max(a.x1, b.x1), max(a.y1, b.y1), min(a.x2, b.x2), min(a.y2, b.y2)};
}

int area(rect a) {
    if(a.y2 <= a.y1 || a.x2 <= a.x1) return 0;
    return (a.y2 - a.y1) * (a.x2 - a.x1);
}

void print(rect a) {
    cout << a.x1 << " " << a.y1 << " " << a.x2 << " " << a.y2 << endl;
}

void solve() {
    int n,m; cin >> n >> m;

    map<pii, vector<rect>> gt;
    for(int i = 0; i < n; i++) {
        rect a; 
        cin >> a.x1 >> a.y1 >> a.x2 >> a.y2;
        a.id = i;
        gt[{a.x1, a.y1}].pb(a);
    }

    vi btoa(n, -1);
    vector<vi> g(m);

    for(int i = 0; i < m; i++) {
        rect a; 
        cin >> a.x1 >> a.y1 >> a.x2 >> a.y2;
        for(int x = max(0, a.x1 - 3); x < a.x2; x++) {
            for(int y = max(0, a.y1 - 3); y < a.y2; y++) {
                if(gt.count({x,y}) == 0) continue;
                for(auto b : gt[{x,y}]) {
                 //   cout << area(inter(a,b)) << " " << area(b) << endl;
                //    print(inter(a,b));
                    if(area(inter(a,b)) * 2 >= area(b)) {
                        g[i].pb(b.id);
                    }
                }
            }
        }
    }

    // for(int i = 0; i < m; i++) {
    //     cout << i << ": ";
    //     for(auto x : g[i]) cout << x << " ";
    //     cout << endl;
    // }

   cout << hopcroftKarp(g, btoa) << endl;

}

int main() {
    cin.tie(0)->sync_with_stdio(0);

    int t; cin >> t;
    for(int tcase = 0; tcase < t; tcase++) {
        solve();
    }
}

詳細信息

Test #1:

score: 100
Accepted
time: 2ms
memory: 3352kb

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: 2ms
memory: 3356kb

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: 1171ms
memory: 15020kb

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: 753ms
memory: 14628kb

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: 499ms
memory: 11464kb

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: 1982ms
memory: 14988kb

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: 2ms
memory: 3504kb

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'