QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#180617 | #5668. Cell Nuclei Detection | arseny_y | AC ✓ | 1471ms | 313716kb | C++23 | 4.6kb | 2023-09-16 02:02:37 | 2023-09-16 02:02:37 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
#define all(a) (a).begin(), (a).end()
#define X first
#define Y second
const ll INF = LLONG_MAX;
const int MOD = 1e9 + 7, MAXN = 2e5 + 1337;
const ld EPS = 1e-6;
void solve();
int main() {
#ifdef LOCAL
freopen("../main_cpp/input.txt", "r", stdin);
#else
cin.tie(nullptr)->sync_with_stdio(false);
#endif
int _ = 1;
cin >> _;
while (_--) {
solve();
}
}
vector<pair<int, int>> g[MAXN];
struct Rect {
int x1, y1, x2, y2;
int sq() {
if (min(x2 - x1, y2 - y1) <= 0) return 0;
return (x2 - x1) * (y2 - y1);
}
int intersec(const Rect oth) const {
Rect ans;
ans.x1 = max(x1, oth.x1);
ans.y1 = max(y1, oth.y1);
ans.x2 = min(x2, oth.x2);
ans.y2 = min(y2, oth.y2);
return ans.sq();
}
};
const int MX = 2000 + 10;
vector<int> have[MX][MX];
struct Dinic {
struct edge {
int v, to;
ll f, c;
edge(int a, int b, ll C, ll d) : v(a), to(b), f(C), c(d) {}
edge() = default;
};
vector<vector<int>> g;
vector<edge> e;
vector<int> dist, pnt;
ll low = 1;
int s, t, n;
Dinic() = default;
Dinic(int n, int ss = -1, int tt = -1) : n(n + 2), s(ss == -1 ? n : ss), t(tt == -1 ? n + 1 : tt), dist(n + 2),
pnt(n + 2, 0) {
g.assign(n + 2, {});
}
void init(int nn, int ss, int tt) {
n = nn + 2, s = (ss == -1 ? nn : ss), t = (tt == -1 ? nn + 1 : tt);
dist.resize(n), pnt.assign(n, 0), g.assign(n, {});
}
void create(int v, int to, ll c) {
g[v].push_back(e.size());
e.emplace_back(v, to, 0, c);
g[to].push_back(e.size());
e.emplace_back(to, v, 0, 0);
}
bool bfs() {
fill(dist.begin(), dist.end(), INT_MAX);
dist[s] = 0;
queue<int> q;
q.push(s);
while (!q.empty()) {
int v = q.front();
q.pop();
for (int i: g[v])
if (e[i].c - e[i].f >= low && dist[e[i].to] == INT_MAX) {
dist[e[i].to] = dist[v] + 1;
q.push(e[i].to);
}
}
return dist[t] != INT_MAX;
}
ll dfs(int v, ll minC = INF) {
if (v == t) return minC;
while (pnt[v] < g[v].size()) {
int i = g[v][pnt[v]];
if (dist[e[i].to] == dist[v] + 1 && e[i].c - e[i].f >= low) {
ll delta = dfs(e[i].to, min(minC, e[i].c - e[i].f));
if (delta) {
e[i].f += delta;
e[i ^ 1].f -= delta;
return delta;
}
}
++pnt[v];
}
return 0;
}
ll flow(int S = -1, int Tt = -1) {
low = 1;
s = (S == -1 ? s : S), t = (S == -1 ? t : Tt);
ll ans = 0;
while (bfs()) {
fill(pnt.begin(), pnt.end(), 0);
while (ll y = dfs(s))
ans += y;
}
return ans;
}
};
void solve() {
for (int i = 0; i < MX; ++i) for (int j = 0; j < MX; ++j) have[i][j].clear();
int n, m;
cin >> n >> m;
int sz = n + m + 10;
for (int i = 1; i <= sz; ++i) {
g[i].clear();
}
vector<Rect> a(n + 1);
for (int i = 1; i <= n; ++i) {
cin >> a[i].x1 >> a[i].y1 >> a[i].x2 >> a[i].y2;
for (int x = a[i].x1; x <= a[i].x2; ++x) {
for (int y = a[i].y1; y <= a[i].y2; ++y) {
have[x][y].push_back(i);
}
}
}
for (int i = 1; i <= m; ++i) {
int ind = i + n;
Rect cur;
cin >> cur.x1 >> cur.y1 >> cur.x2 >> cur.y2;
for (int x = cur.x1; x <= cur.x2; ++x) {
for (int y = cur.y1; y <= cur.y2; ++y) {
for (auto &el : have[x][y]) {
if (a[el].intersec(cur) * 2 >= a[el].sq()) {
g[el + 1].push_back({ind + 1, 1});
}
}
}
}
}
int ans = 0;
for (int i = 2; i <= n + 1; ++i) {
g[1].push_back({i, 1});
}
for (int i = n + 2; i <= n + m + 3; ++i) {
g[i].emplace_back(sz - 1, 1);
}
Dinic dinic(sz + 10);
for (int i = 1; i <= sz; ++i) {
sort(all(g[i])), g[i].resize(unique(all(g[i])) - g[i].begin());
for (auto &[to, c] : g[i]) {
dinic.create(i, to, c);
}
}
cout << dinic.flow(1, sz - 1) << "\n";
}
詳細信息
Test #1:
score: 100
Accepted
time: 41ms
memory: 103056kb
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: 36ms
memory: 103052kb
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: 582ms
memory: 313716kb
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: 424ms
memory: 197012kb
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: 307ms
memory: 148940kb
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: 1471ms
memory: 302584kb
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: 28ms
memory: 103284kb
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'