QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#548994 | #5668. Cell Nuclei Detection | skrghariapa# | WA | 569ms | 7592kb | C++14 | 4.3kb | 2024-09-05 23:01:58 | 2024-09-05 23:01:58 |
Judging History
answer
#include <bits/stdc++.h>
#define NIL 0
#define INF INT_MAX
using namespace std;
class BipartiteMatching {
int U, V;
vector<vector<int>> adj;
vector<int> pairU, pairV, dist;
public:
BipartiteMatching(int U, int V) {
this->U = U;
this->V = V;
adj.resize(U + 1);
pairU.resize(U + 1, NIL);
pairV.resize(V + 1, NIL);
dist.resize(U + 1, INF);
}
void addEdge(int u, int v) {
adj[u].push_back(v);
}
bool bfs() {
queue<int> q;
for (int u = 1; u <= U; u++) {
if (pairU[u] == NIL) {
dist[u] = 0;
q.push(u);
} else {
dist[u] = INF;
}
}
dist[NIL] = INF;
while (!q.empty()) {
int u = q.front();
q.pop();
if (dist[u] < dist[NIL]) {
for (int v : adj[u]) {
if (dist[pairV[v]] == INF) {
dist[pairV[v]] = dist[u] + 1;
q.push(pairV[v]);
}
}
}
}
return dist[NIL] != INF;
}
bool dfs(int u) {
if (u != NIL) {
for (int v : adj[u]) {
if (dist[pairV[v]] == dist[u] + 1 && dfs(pairV[v])) {
pairV[v] = u;
pairU[u] = v;
return true;
}
}
dist[u] = INF;
return false;
}
return true;
}
int hopcroftKarp() {
int matching = 0;
while (bfs()) {
for (int u = 1; u <= U; u++) {
if (pairU[u] == NIL && dfs(u)) {
matching++;
}
}
}
return matching;
}
};
int areaIntersect(array<int, 4> p1, array<int, 4> p2) {
int x_overlap = max(0, min(p1[2], p2[2]) - max(p1[0], p2[0]));
int y_overlap = max(0, min(p1[3], p2[3]) - max(p1[1], p2[1]));
return x_overlap * y_overlap;
}
int area(array<int, 4> p) {
int x = p[2] - p[0];
int y = p[3] - p[1];
return x * y;
}
void solve() {
int m, n;
cin >> m >> n;
BipartiteMatching bm(m, n);
vector<array<int, 4>> base, target;
map<pair<float, float>, int> mp;
for (int i = 0; i < m; i++) {
int a, b, c, d;
cin >> a >> b >> c >> d;
base.push_back({a, b, c, d});
}
for (int i = 0; i < n; i++) {
int a, b, c, d;
cin >> a >> b >> c >> d;
target.push_back({a, b, c, d});
float midx = (float) (a + c)/ 2.0;
float midy = (float) (b + d)/ 2.0;
mp[{midx, midx}] = i;
}
for (int i = 0; i < m; i++) {
float midx = (float) (base[i][0] + base[i][2]) / 2.0;
float midy = (float) (base[i][1] + base[i][3]) / 2.0;
for (float j = 0; j <= 2.0; j += 0.5) {
for (float k = 0; k <= 2.0; k += 0.5) {
if (mp.find({midx + j, midy + k}) != mp.end() && area(base[i]) <= 2 * areaIntersect(base[i], target[mp[{midx + j, midy + k}]])) {
// 2 * area(base[i]) <= areaIntersect(base[i], target[mp[{midx + j, midy + k}]])
// cout << i << ": " << area(base[i]) << ", " << mp[{midx + j, midy + k}] << ": " << areaIntersect(base[i], target[mp[{midx + j, midy + k}]]) << '\n';
bm.addEdge(i + 1, mp[{midx + j, midy + k}] + 1);
}
if (mp.find({midx + j, midy - k}) != mp.end() && area(base[i]) <= 2 * areaIntersect(base[i], target[mp[{midx + j, midy - k}]])) {
// cout << i << ": " << area(base[i]) << ", " << mp[{midx + j, midy - k}] << ": " << areaIntersect(base[i], target[mp[{midx + j, midy - k}]]) << '\n';
bm.addEdge(i + 1, mp[{midx + j, midy - k}] + 1);
}
if (mp.find({midx - j, midy + k}) != mp.end() && area(base[i]) <= 2 * areaIntersect(base[i], target[mp[{midx - j, midy + k}]])) {
// cout << i << ": " << area(base[i]) << ", " << mp[{midx - j, midy + k}] << ": " << areaIntersect(base[i], target[mp[{midx - j, midy + k}]]) << '\n';
bm.addEdge(i + 1, mp[{midx - j, midy + k}] + 1);
}
if (mp.find({midx - j, midy - k}) != mp.end() && area(base[i]) <= 2 * areaIntersect(base[i], target[mp[{midx - j, midy - k}]])) {
// cout << i << ": " << area(base[i]) << ", " << mp[{midx - j, midy - k}] << ": " << areaIntersect(base[i], target[mp[{midx - j, midy - k}]]) << '\n';
bm.addEdge(i + 1, mp[{midx - j, midy - k}] + 1);
}
}
}
}
int res = bm.hopcroftKarp();
cout << res << '\n';
}
int main() {
int t = 1;
cin >> t;
while (t--) {
solve();
}
}
详细
Test #1:
score: 100
Accepted
time: 0ms
memory: 3620kb
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: 3616kb
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: -100
Wrong Answer
time: 569ms
memory: 7592kb
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:
0 0 0 0 0
result:
wrong answer 1st lines differ - expected: '50000', found: '0'