QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#84807 | #5668. Cell Nuclei Detection | runewrz | AC ✓ | 2257ms | 132624kb | C++17 | 3.3kb | 2023-03-06 19:16:09 | 2023-03-06 19:16:44 |
Judging History
answer
#include <bits/stdc++.h>
#define fo(i, x, y) for (int i = int(x); i <= int(y); ++i)
#define fd(i, x, y) for (int i = int(x); i >= int(y); --i)
#define fi first
#define se second
using namespace std;
using pii = pair<int, int>;
using i64 = long long;
constexpr int X = 2e3, N = 1e5 + 5, inf = 1e9;
struct Dinic
{
int n, m, s, t;
struct Edge
{
int from, to, cap;
};
vector<Edge> edges;
vector<int> G[N];
int d[N];
int cur[N];
void init(int n)
{
this->n = n;
fo(i, 0, n - 1)
G[i].clear();
edges.clear();
}
void AddEdge(int from, int to, int cap)
{
edges.push_back({from, to, cap});
edges.push_back({to, from, 0});
m = edges.size();
G[from].push_back(m - 2);
G[to].push_back(m - 1);
}
bool BFS()
{
memset(d, 0, sizeof(d));
queue<int> q;
q.push(s);
d[s] = 1;
while (q.size())
{
int x = q.front();
q.pop();
for (int i = 0; i < G[x].size(); i++)
{
Edge &e = edges[G[x][i]];
if (!d[e.to] && e.cap > 0)
{
d[e.to] = d[x] + 1;
q.emplace(e.to);
}
}
}
return d[t] > 0;
}
int DFS(int x, int a)
{
if (x == t || !a) return a;
int flow = 0, f;
for (int &i = cur[x]; i < G[x].size(); i++)
{
Edge &e = edges[G[x][i]], &re = edges[G[x][i] ^ 1];
if (d[x] + 1 == d[e.to] && (f = DFS(e.to, min(a, e.cap))) > 0)
{
e.cap -= f;
re.cap += f;
flow += f;
a -= f;
if (!a) break;
}
}
return flow;
}
int MaxFlow(int s, int t)
{
this->s = s; this->t = t;
int flow = 0;
while (BFS())
{
memset(cur, 0, sizeof(cur));
flow += DFS(s, inf);
}
return flow;
}
}fl;
void work()
{
int m, n;
cin >> m >> n;
vector<array<int, 4>> a(m), b(n);
for (auto& [x, y, xx, yy]: a)
{
cin >> x >> y >> xx >> yy;
xx--; yy--;
}
for (auto& [x, y, xx, yy]: b)
{
cin >> x >> y >> xx >> yy;
xx--; yy--;
}
fl.init(n + m + 2);
fl.s = n + m; fl.t = n + m + 1;
map<int, vector<int>> id;
fo(i, 0, n - 1)
{
auto &[x, y, xx, yy] = b[i];
fo(p, x, xx)
fo(q, y, yy)
id[p * (X + 1) + q].emplace_back(i);
}
fo(i, 0, m - 1)
fl.AddEdge(fl.s, i, 1);
fo(i, 0, n - 1)
fl.AddEdge(i + m, fl.t, 1);
fo(ia, 0, m - 1)
{
auto &[x, y, xx, yy] = a[ia];
vector<int> tid;
fo(p, x, xx)
fo(q, y, yy)
{
if (!id.count(p * (X + 1) + q)) continue;
for (auto& i: id[p * (X + 1) + q])
{
constexpr auto area = [](const array<int, 4>& a)
{
auto& [x, y, xx, yy] = a;
return (xx - x + 1) * (yy - y + 1);
};
constexpr auto ad = [](const array<int, 4>& a, const array<int, 4>& b) -> array<int, 4>
{
auto& [x, y, xx, yy] = a;
auto& [bx, by, bxx, byy] = b;
return { max(x, bx), max(y, by), min(xx, bxx), min(yy, byy) };
};
auto check = [&](const array<int, 4>& a, const array<int, 4>& b) -> bool
{
int s1 = area(a), s2 = area(ad(a, b));
return 2 * s2 >= s1;
};
if (check({ x, y, xx, yy }, b[i]))
tid.emplace_back(i);
}
}
sort(tid.begin(), tid.end());
tid.erase(unique(tid.begin(), tid.end()), tid.end());
for (auto &ib : tid)
fl.AddEdge(ia, ib + m, 1);
}
int ans = fl.MaxFlow(fl.s, fl.t);
cout << ans << '\n';
}
int main()
{
ios::sync_with_stdio(0), cin.tie(0);
#ifdef LC
freopen("t.in", "r", stdin);
freopen("t.out", "w", stdout);
#endif
int T;
cin >> T;
while (T--)
work();
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 6492kb
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: 3ms
memory: 6500kb
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: 2257ms
memory: 128960kb
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: 1687ms
memory: 132624kb
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: 1297ms
memory: 48544kb
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: 1468ms
memory: 47592kb
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: 4ms
memory: 6476kb
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'