QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#578387 | #5668. Cell Nuclei Detection | CelestialCoder# | TL | 4266ms | 30124kb | C++20 | 3.4kb | 2024-09-20 18:47:08 | 2024-09-20 18:47:08 |
Judging History
answer
#include <bits/stdc++.h>
#ifdef SHARAELONG
#include "../../cpp-header/debug.hpp"
#endif
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
struct HopcroftKarp {
int n, m;
vector<vector<int>> g;
vector<int> dst, le, ri;
vector<char> visit, track;
HopcroftKarp(int n, int m): n(n), m(m), g(n), dst(n), le(n, -1), ri(m, -1), visit(n), track(n+m) {}
void add_edge(int s, int e) {
g[s].push_back(e);
}
bool bfs() {
bool res = false;
queue<int> que;
fill(dst.begin(), dst.end(), 0);
for (int i=0; i<n; ++i) {
if (le[i] == -1) que.push(i), dst[i] = 1;
}
while (!que.empty()) {
int v = que.front();
que.pop();
for (auto i: g[v]) {
if (ri[i] == -1) res = true;
else if (!dst[ri[i]]) dst[ri[i]] = dst[v]+1, que.push(ri[i]);
}
}
return res;
}
bool dfs(int v) {
if (visit[v]) return false;
visit[v] = 1;
for (auto i: g[v]) {
if (ri[i] == -1 || !visit[ri[i]] && dst[ri[i]] == dst[v]+1 && dfs(ri[i])) {
le[v] = i;
ri[i] = v;
return true;
}
}
return false;
}
int maximum_matching() {
int res = 0;
fill(le.begin(), le.end(), -1);
fill(ri.begin(), ri.end(), -1);
while (bfs()) {
fill(visit.begin(), visit.end(), 0);
for (int i=0; i<n; ++i) {
if (le[i] == -1) res += dfs(i);
}
}
return res;
}
};
void solve() {
int m, n;
cin >> m >> n;
vector<array<int, 4>> gr(m), de(n);
for (int i=0; i<m; ++i) {
for (int j=0; j<4; ++j) cin >> gr[i][j];
}
for (int i=0; i<n; ++i) {
for (int j=0; j<4; ++j) cin >> de[i][j];
}
map<pii, vector<int>> mp;
for (int i=0; i<n; ++i) {
auto&[x1, y1, x2, y2] = de[i];
mp[{ x1,y1 }].push_back(i);
mp[{ x1,y2 }].push_back(i);
mp[{ x2,y1 }].push_back(i);
mp[{ x2,y2 }].push_back(i);
}
auto inter = [](const array<int, 4>& g, const array<int, 4>& d) {
if (min(g[2], d[2]) < max(g[0], d[0])) return false;
if (min(g[3], d[3]) < max(g[1], d[1])) return false;
int sz = (min(g[2], d[2]) - max(g[0], d[0])) * (min(g[3], d[3]) - max(g[1], d[1]));
return (2*sz >= (g[2]-g[0])*(g[3]-g[1]));
};
auto search = [&](int i, pii c, set<int>& s) {
for (int x=c.first-4; x<=c.first+4; ++x) {
for (int y=c.second-4; y<=c.second+4; ++y) {
if (mp.find({ x,y }) != mp.end()) {
auto& v = mp[{ x,y }];
for (int j: v) {
if (inter(gr[i], de[j])) s.insert(j);
}
}
}
}
};
HopcroftKarp g(m, n);
for (int i=0; i<m; ++i) {
auto[x1, y1, x2, y2] = gr[i];
set<int> s;
search(i, { x1,y1 }, s);
search(i, { x1,y2 }, s);
search(i, { x2,y1 }, s);
search(i, { x2,y2 }, s);
for (int j: s) g.add_edge(i, j);
}
cout << g.maximum_matching() << '\n';
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
int cases;
cin >> cases;
while (cases--) solve();
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 3564kb
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: 0ms
memory: 3488kb
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: 3703ms
memory: 16676kb
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: 3672ms
memory: 30124kb
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: 4266ms
memory: 29516kb
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: -100
Time Limit Exceeded
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...