QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#123413 | #5668. Cell Nuclei Detection | GenshinImpactsFault# | AC ✓ | 593ms | 135992kb | C++14 | 3.0kb | 2023-07-12 15:34:14 | 2023-07-12 15:34:17 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
// #define int long long
#define MOD ((1ll << 31) - 1)
vector<array<int, 3>> pos[2010][2010];
vector<pair<int, int>> cls;
const int N = 100010;
const int INF = N * 10;
int s, t, ans;
int cnt, hed[N], to[N * 100], nxt[N * 100], val[N * 100];
int cur[N], dis[N], g[N];
void init(int n) {
s = n + 1, t = n + 2, ans = 0;
for(int i = 0; i <= t; i++) hed[i] = cur[i] = dis[i] = g[i] = 0;
for(int i = 0; i <= cnt; i++) to[i] = val[i] = nxt[i] = 0;
cnt = 1;
}
void addedge(int x, int y, int w) {
// cout << ">>> " << x << " " << y << " " << w << "\n";
to[++cnt] = y, nxt[cnt] = hed[x], hed[x] = cnt, val[cnt] = w;
to[++cnt] = x, nxt[cnt] = hed[y], hed[y] = cnt, val[cnt] = 0;
}
int SAP(int u, int flow) {
if(u == t) return flow;
int v = flow;
for(; cur[u]; cur[u] = nxt[cur[u]]) {
if(dis[to[cur[u]]] == dis[u] - 1 && val[cur[u]]) {
int f = SAP(to[cur[u]], min(v, val[cur[u]]));
val[cur[u]] -= f, val[cur[u] ^ 1] += f;
v -= f;
if(!v) return flow;
}
}
if(!--g[dis[u]++]) dis[s] = t + 1;
++g[dis[u]], cur[u] = hed[u]; return flow - v;
}
void Do() {
int m, n;
cin >> m >> n;
init(n + m);
for(int i = 1; i <= m; ++ i) {
int x1, y1, x2, y2;
cin >> x1 >> y1 >> x2 >> y2;
// cerr << x1 << " " << y1 << " " << x2 << " " << y2 << endl;
pos[x1][y1].push_back({x2, y2, i});
cls.push_back({x1, y1});
}
for(int i = 1; i <= n; ++ i) {
int x1, y1, x2, y2;
cin >> x1 >> y1 >> x2 >> y2;
for(int x = max(x1 - 3, 0); x <= x2 - 1; ++ x) {
for(int y = max(y1 - 3, 0); y <= y2 - 1; ++ y) {
for(auto [x3, y3, id] : pos[x][y]) {
// cerr << x3 << " " << y3 << " " << id << endl;
// x y x3 y3
// x1 y1 x2 y2;
int X1 = max(x, x1), Y1 = max(y, y1);
int X2 = min(x3, x2), Y2 = min(y3, y2);
int S = (X2 - X1) * (Y2 - Y1);
if(X1 > X2 || Y1 > Y2) S = 0;
int pS = (x3 - x) * (y3 - y);
// cerr << pS << " " << S << endl;
if(S * 2 >= pS) {
// add edge id -> i + m
// cerr << id << " " << i + m << endl;
addedge(id, i + m, 1);
}
}
}
}
}
// exit(0);
for(int i = 1; i <= m; i++) addedge(s, i, 1);
for(int i = 1; i <= n; i++) addedge(i + m, t, 1);
for(int i = 1; i <= t; i++) cur[i] = hed[i];
for(g[0] = t; dis[s] <= t; ans += SAP(s, INF));
cout << ans << "\n";
for(auto [x, y] : cls) {
pos[x][y].clear();
}
}
signed main() {
// ios::sync_with_stdio(0);
// cin.tie(0), cout.tie(0);
int T;
cin >> T;
while(T --) {
Do();
}
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 15ms
memory: 102868kb
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: 8ms
memory: 106076kb
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: 285ms
memory: 129424kb
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: 304ms
memory: 135992kb
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: 312ms
memory: 113512kb
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: 593ms
memory: 125784kb
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: 10ms
memory: 102812kb
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'