QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#269968 | #5668. Cell Nuclei Detection | STnofarjo# | AC ✓ | 803ms | 106740kb | C++20 | 3.2kb | 2023-11-30 12:50:50 | 2023-11-30 12:50:50 |
Judging History
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;
}
詳細信息
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'