QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#84133 | #5668. Cell Nuclei Detection | HOLIC# | AC ✓ | 4154ms | 255080kb | C++20 | 3.7kb | 2023-03-05 18:27:04 | 2023-03-05 18:27:07 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
#define endl "\n"
const int N = 2e5 + 1009;
const int M = 8 * N * 40 + N;
int n, m;
namespace Dinic{
const int inf = 0x3f3f3f3f;
int s, t, d[N], cpy[N], head[N], nexts[M], edge[M], ver[M], cnt = 1;
void add(int x, int y, int z) {
ver[++ cnt] = y, nexts[cnt] = head[x], head[x] = cpy[x] = cnt, edge[cnt] = z;
}
void clear() {
for(int i = 1; i <= t; i ++) head[i] = cpy[i] = 0;
cnt = 1;
}
bool bfs() {
memset(d, 0, sizeof(d));
queue<int> q;
q.push(s); d[s] = 1;
while(!q.empty()) {
int x = q.front(); q.pop();
for(int i = head[x]; i; i = nexts[i]) {
int v = ver[i], w = edge[i];
if(w > 0 && !d[v]) d[v] = d[x] + 1, q.push(v);
}
}
if(!d[t]) return false;
return true;
}
int dfs(int x, int flow) {
if(x == t || flow == 0) return flow;
int maxflow = 0, f;
for(int &i = head[x]; i; i = nexts[i]) {
int v = ver[i], w = edge[i];
if(d[v] == d[x] + 1 && w) {
f = dfs(v, min(flow, w));
if(f > 0) {
edge[i] -= f;
edge[i ^ 1] += f;
flow -= f; maxflow += f;
if(!flow) return maxflow;
}
}
}
return maxflow;
}
int dinic() {
int maxflow = 0, flow;
while(bfs()) {
maxflow += dfs(s, inf);
for(int i = 1; i <= t; i ++)
head[i] = cpy[i];
}
return maxflow;
}
}
vector<int> ver[2005][2005], Ver[2005][2005];
int xx[50005], yy[50005], xxx[50005], yyy[50005];
int XX[50005], YY[50005], XXX[50005], YYY[50005];
void work() {
cin >> n >> m;
Dinic::clear();
Dinic::s = n + m + 1;
Dinic::t = n + m + 2;
for(int i = 0; i <= 2000; i ++)
for(int j = 0; j <= 2000; j ++)
ver[i][j].clear(), Ver[i][j].clear();
for(int i = 1; i <= n; i ++) {
cin >> xx[i] >> yy[i] >> xxx[i] >> yyy[i];
for(int I = xx[i]; I <= xxx[i]; I ++) {
for(int J = yy[i]; J <= yyy[i]; J ++) {
ver[I][J].push_back(i);
}
}
Dinic::add(i + m, Dinic::t, 1); Dinic::add(Dinic::t, i + m, 0);
}
for(int i = 1; i <= m; i ++) {
cin >> XX[i] >> YY[i] >> XXX[i] >> YYY[i];
Dinic::add(Dinic::s, i, 1); Dinic::add(i, Dinic::s, 0);
}
for(int i = 1; i <= m; i ++) {
map<int,int> vis;
int x = XX[i], y = YY[i], X = XXX[i], Y = YYY[i];
auto check = [&](int id) {
int lx = xx[id], ly = yy[id];
int rx = xxx[id], ry = yyy[id];
int s = (rx - lx) * (ry - ly);
if(lx < x) lx = x;
if(rx > X) rx = X;
if(ly < y) ly = y;
if(ry > Y) ry = Y;
if(rx - lx <= 0) return false;
if(ry - ly <= 0) return false;
int ss = (rx - lx) * (ry - ly);
return ss * 2 >= s;
};
for(int I = x; I <= X; I ++) {
for(int J = y; J <= Y; J ++) {
for(auto id : ver[I][J]) {
if(check(id)) {
if(vis[id]) continue;
vis[id] = 1;
Dinic::add(i, id + m, 1);
Dinic::add(id + m, i, 0);
}
}
}
}
}
cout << Dinic::dinic() << endl;
}
signed main() {
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
int Case;
cin >> Case;
while(Case --) work();
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 87ms
memory: 192644kb
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: 76ms
memory: 192588kb
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: 1206ms
memory: 243244kb
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: 492ms
memory: 255080kb
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: 434ms
memory: 221916kb
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: 4154ms
memory: 219684kb
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: 47ms
memory: 192716kb
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'