QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#269968#5668. Cell Nuclei DetectionSTnofarjo#AC ✓803ms106740kbC++203.2kb2023-11-30 12:50:502023-11-30 12:50:50

Judging History

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

  • [2023-11-30 12:50:50]
  • 评测
  • 测评结果:AC
  • 用时:803ms
  • 内存:106740kb
  • [2023-11-30 12:50:50]
  • 提交

answer

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

#define x1 lqceopqiwecm
#define x2 poqwiecqweco
#define y1 qpwoeicmpoiq
#define y2 weopcmrqvoiu

#define rep(i,a,b) for (int i=(a); i<(b); ++i)
#define sz(x) (int)(x).size()
#define all(x) (x).begin(), (x).end()
typedef vector<int> vi;

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;
		rep(a,0,sz(g)) if (A[a] == 0) cur.push_back(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.push_back(btoa[b]);
				}
			}
			if (islast) break;
			if (next.empty()) return res;
			for (int a : next) A[a] = lay;
			cur.swap(next);
		}
		rep(a,0,sz(g))
			res += dfs(a, 0, g, btoa, A, B);
	}
}

const int maxx = 2003, maxd = 4;

void solve() {
	int m, n;
	cin >> m >> n;
	vector<pair<pair<int, int>, pair<int, int>>> ra(m), rb(n);
	for (int i=0; i<m; i++) {
		int x1, y1, x2, y2;
		cin >> x1 >> y1 >> x2 >> y2;
		ra[i] = {{x1, y1}, {x2, y2}};
	}
	for (int i=0; i<n; i++) {
		int x1, y1, x2, y2;
		cin >> x1 >> y1 >> x2 >> y2;
		rb[i] = {{x1, y1}, {x2, y2}};
	}
	
	vector<vector<vector<int>>> who(maxx, vector<vector<int>>(maxx));
	for (int i=0; i<m; i++) {
		auto [p1, p2] = ra[i];
		auto [x1, y1] = p1;
		auto [x2, y2] = p2;
		// for (int x=x1; x<=x2; x++) {
		// 	for (int y=y1; y<=y2; y++) {
		// 		who[x][y].push_back(i);
		// 	}
		// }
		who[x1][y1].push_back(i);
	}
	vector<vi> g(n);
	for (int i=0; i<n; i++) {
		auto [p1, p2] = rb[i];
		auto [x1, y1] = p1;
		auto [x2, y2] = p2;
		int l = max(0, x1 - maxd), r = min(maxx - 1, x2);
		int d = max(0, y1 - maxd), u = min(maxx - 1, y2);
		set<int> cands;
		for (int x=l; x<=r; x++) {
			for (int y=d; y<=u; y++) {
				for (int j : who[x][y]) cands.insert(j);
			}
		}

		for (int j : cands) {
			auto [p3, p4] = ra[j];
			auto [x3, y3] = p3;
			auto [x4, y4] = p4;
			int mxl = max(x1, x3), mnr = min(x2, x4);
			int mxd = max(y1, y3), mnu = min(y2, y4);
			int iarea = (mnr - mxl) * (mnu - mxd);
			int barea = (x4 - x3) * (y4 - y3);
			if (mxl <= mnr && mxd <= mnu &&  iarea * 2 >= barea) {
				// cout << "tebel " << i << " tipis " << j << '\n';
				// cout << "isect: " << mxl << ' ' << mxd << ' ' << mnr << ' ' << mnu << ' ' << iarea << '\n';
				// cout << "tebel: " << x1 << ' ' << y1 << ' ' << x2 << ' ' << y2 << '\n';
				// cout << "tipis: " << x3 << ' ' << y3 << ' ' << x4 << ' ' << y4 << ' ' << barea << '\n';
				g[i].push_back(j);
			}
		}
	}

	vi btoa(m, -1); 
	int ans = hopcroftKarp(g, btoa);
	cout << ans << '\n';
}

int main() {
	ios_base::sync_with_stdio(0); cin.tie(0);
	int tcs;
	cin >> tcs;
	while (tcs--) {
		solve();
	}
	return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 48ms
memory: 97348kb

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: 60ms
memory: 97440kb

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: 416ms
memory: 106500kb

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: 381ms
memory: 106304kb

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: 245ms
memory: 103304kb

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: 803ms
memory: 106740kb

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: 21ms
memory: 97200kb

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'