QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#95712#5668. Cell Nuclei DetectionSGColin#AC ✓948ms138852kbC++202.6kb2023-04-11 15:29:082023-04-11 15:29:12

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-04-11 15:29:12]
  • 评测
  • 测评结果:AC
  • 用时:948ms
  • 内存:138852kb
  • [2023-04-11 15:29:08]
  • 提交

answer

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

inline int rd() {
	int x = 0;
	bool f = 0;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) f |= (c == '-');
	for (; isdigit(c); c = getchar()) x = x * 10 + (c ^ 48);
	return f ? -x : x;
}

inline bool getmin(int &a, int b) {return a > b ? a = b, true : false;}

template<const int V, const int inf = 1000000000>
struct Hopcraft { 
	bool vis[V];
	vector<int> e[V];
	int nl, nr, ml[V], mr[V], dl[V], dr[V];
	inline void reset(int _nl, int _nr) {
		nl = _nl; nr = _nr;
		for (int i = 1; i <= nl; ++i) e[i].clear();
	}
	inline bool bfs() {
		static int q[V], hd, tl, dT;
		hd = 1; tl = 0; dT = inf;
		memset(dl, -1, sizeof(int) * (nl + 1));
		memset(dr, -1, sizeof(int) * (nr + 1));
		for (int i = 1; i <= nl; ++i) 
			if (!ml[i]) {dl[i] = 0; q[++tl] = i;}
		for (int u; hd <= tl;) {
			if (dl[u = q[hd++]] >= dT) break;
			for (auto v : e[u])
				if (dr[v] == -1) {
					dr[v] = dl[u] + 1;
					if (!mr[v]) getmin(dT, dr[v] + 1);
					else {dl[mr[v]] = dr[v] + 1; q[++tl] = mr[v];}
				}
		}
		return dT != inf;
	}
	bool dfs(int u) {
		for (auto v : e[u]) {
			if (vis[v] || dl[u] + 1 != dr[v]) continue;
			vis[v] = true;
			if (!mr[v] || dfs(mr[v])) {mr[v] = u; ml[u] = v; return true;}
		}
		return false;
	}
	inline void add(int u, int v) {e[u].push_back(v);}
	inline int max_matching() {
		int ans = 0;
		memset(ml, 0, sizeof(ml));
		memset(mr, 0, sizeof(mr));
		while (bfs()) {
			memset(vis, 0, sizeof(bool) * (nr + 1));
			for (int i = 1; i <= nl; ++i)
				if (!ml[i]) ans += dfs(i);
		}
		return ans;
	}
};

#define N 50007

Hopcraft<N> E;

vector<int> id[2007][2007];

set<int> tmp;

int lx[N], ly[N], rx[N], ry[N];

inline void work() {
	for (int i = 0; i <= 2000; ++i)
		for (int j = 0; j <= 2000; ++j) id[i][j].clear();
	int m = rd(), n = rd();
	E.reset(m, n);
	for (int i = 1; i <= m; ++i) {
		int ax = rd(), ay = rd(), bx = rd(), by = rd();
		lx[i] = ax; ly[i] = ay; rx[i] = bx; ry[i] = by;
		for (int x = ax; x < bx; ++x)
			for (int y = ay; y < by; ++y) id[x][y].push_back(i);
	}
	for (int i = 1; i <= n; ++i) {
		int ax = rd(), ay = rd(), bx = rd(), by = rd();
		tmp.clear();
		for (int x = ax; x < bx; ++x)
			for (int y = ay; y < by; ++y) 
				for (auto t : id[x][y])
					if (!tmp.count(t)) {
						tmp.insert(t);
						if (max(0, (min(bx, rx[t]) - max(ax, lx[t]))) * max(0, (min(by, ry[t]) - max(ay, ly[t]))) * 2 >= (rx[t] - lx[t]) * (ry[t] - ly[t])) {
							E.add(t, i);
						}
					}
	}
	printf("%d\n", E.max_matching());
}

int main() {
	for (int t = rd(); t; --t) work();
	return 0;
}

详细

Test #1:

score: 100
Accepted
time: 34ms
memory: 100064kb

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: 37ms
memory: 99584kb

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: 520ms
memory: 132424kb

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: 289ms
memory: 138852kb

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: 247ms
memory: 121452kb

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: 948ms
memory: 109000kb

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: 24ms
memory: 99632kb

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'