QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#84807#5668. Cell Nuclei DetectionrunewrzAC ✓2257ms132624kbC++173.3kb2023-03-06 19:16:092023-03-06 19:16:44

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-06 19:16:44]
  • 评测
  • 测评结果:AC
  • 用时:2257ms
  • 内存:132624kb
  • [2023-03-06 19:16:09]
  • 提交

answer

#include <bits/stdc++.h>

#define fo(i, x, y) for (int i = int(x); i <= int(y); ++i)
#define fd(i, x, y) for (int i = int(x); i >= int(y); --i)
#define fi first
#define se second

using namespace std;
using pii = pair<int, int>;
using i64 = long long;

constexpr int X = 2e3, N = 1e5 + 5, inf = 1e9;

struct Dinic
{
	int n, m, s, t;
	struct Edge
	{
		int from, to, cap;
	};
	vector<Edge> edges;
	vector<int> G[N];
	int d[N];
	int cur[N];
	void init(int n)
	{
		this->n = n;
		fo(i, 0, n - 1)
			G[i].clear();
		edges.clear();
	}
	void AddEdge(int from, int to, int cap)
	{
		edges.push_back({from, to, cap});
		edges.push_back({to, from, 0});
		m = edges.size();
		G[from].push_back(m - 2);
		G[to].push_back(m - 1);
	}
	bool BFS()
	{
		memset(d, 0, sizeof(d));
		queue<int> q;
		q.push(s);
		d[s] = 1;
		while (q.size())
		{
			int x = q.front();
			q.pop();
			for (int i = 0; i < G[x].size(); i++)
			{
				Edge &e = edges[G[x][i]];
				if (!d[e.to] && e.cap > 0)
				{
					d[e.to] = d[x] + 1;
					q.emplace(e.to);
				}
			}
		}
		return d[t] > 0;
	}
	int DFS(int x, int a)
	{
		if (x == t || !a) return a;
		int flow = 0, f;
		for (int &i = cur[x]; i < G[x].size(); i++)
		{
			Edge &e = edges[G[x][i]], &re = edges[G[x][i] ^ 1];
			if (d[x] + 1 == d[e.to] && (f = DFS(e.to, min(a, e.cap))) > 0)
			{
				e.cap -= f;
				re.cap += f;
				flow += f;
				a -= f;
				if (!a) break;
			}
		}
		return flow;
	}
	int MaxFlow(int s, int t)
	{
		this->s = s; this->t = t;
		int flow = 0;
		while (BFS())
		{
			memset(cur, 0, sizeof(cur));
			flow += DFS(s, inf);
		}
		return flow;
	}
}fl;

void work()
{
	int m, n;
	cin >> m >> n;
	vector<array<int, 4>> a(m), b(n);
	for (auto& [x, y, xx, yy]: a)
	{
		cin >> x >> y >> xx >> yy;
		xx--; yy--;
	}
	for (auto& [x, y, xx, yy]: b)
	{
		cin >> x >> y >> xx >> yy;
		xx--; yy--;
	}
	fl.init(n + m + 2);
	fl.s = n + m; fl.t = n + m + 1;
	map<int, vector<int>> id;
	fo(i, 0, n - 1)
	{
		auto &[x, y, xx, yy] = b[i];
		fo(p, x, xx)
			fo(q, y, yy)
				id[p * (X + 1) + q].emplace_back(i);
	}
	fo(i, 0, m - 1)
		fl.AddEdge(fl.s, i, 1);
	fo(i, 0, n - 1)
		fl.AddEdge(i + m, fl.t, 1);
	fo(ia, 0, m - 1)
	{
		auto &[x, y, xx, yy] = a[ia];
		vector<int> tid;
		fo(p, x, xx)
			fo(q, y, yy)
			{
				if (!id.count(p * (X + 1) + q)) continue;
				for (auto& i: id[p * (X + 1) + q])
				{
					constexpr auto area = [](const array<int, 4>& a)
					{
						auto& [x, y, xx, yy] = a;
						return (xx - x + 1) * (yy - y + 1);
					};
					constexpr auto ad = [](const array<int, 4>& a, const array<int, 4>& b) -> array<int, 4>
					{
						auto& [x, y, xx, yy] = a;
						auto& [bx, by, bxx, byy] = b;
						return { max(x, bx), max(y, by), min(xx, bxx), min(yy, byy) };
					};
					auto check = [&](const array<int, 4>& a, const array<int, 4>& b) -> bool
					{
						int s1 = area(a), s2 = area(ad(a, b));
						return 2 * s2 >= s1;
					};
					if (check({ x, y, xx, yy }, b[i]))
						tid.emplace_back(i);
				}
			}
		sort(tid.begin(), tid.end());
		tid.erase(unique(tid.begin(), tid.end()), tid.end());
		for (auto &ib : tid)
			fl.AddEdge(ia, ib + m, 1);
	}
	int ans = fl.MaxFlow(fl.s, fl.t);
	cout << ans << '\n';
}

int main()
{
	ios::sync_with_stdio(0), cin.tie(0);
#ifdef LC
	freopen("t.in", "r", stdin);
	freopen("t.out", "w", stdout);
#endif
	int T;
	cin >> T;
	while (T--)
		work();
	return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 6492kb

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: 3ms
memory: 6500kb

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: 2257ms
memory: 128960kb

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: 1687ms
memory: 132624kb

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: 1297ms
memory: 48544kb

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: 1468ms
memory: 47592kb

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: 4ms
memory: 6476kb

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'