QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#88412 | #5668. Cell Nuclei Detection | CCSU_Wine | AC ✓ | 3808ms | 127712kb | C++20 | 3.1kb | 2023-03-16 10:27:27 | 2023-03-16 10:28:07 |
Judging History
answer
#include <queue>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 2010, M = 50010, inf = 1e5, SIZ = 2000;
int n, m;
struct Box{
int x1, y1, x2, y2, area;
void input(){ cin >> x1 >> y1 >> x2 >> y2; }
int get_id(){ return x1 * SIZ + y1; } //返回左下角坐标对应编号
void get_area(){
area = (x2 - x1) * (y2 - y1);
}
}tr[M], de[M];
struct edge{
int to, nex, w;
}e[M * 200];
int tot, s, t, head[M * 2];
void add(int from, int to, int w){
e[++ tot].w = w;
e[tot].to = to;
e[tot].nex = head[from];
head[from] = tot;
}
int dep[M * 2];
bool bfs(){
queue<int>q;
q.push(s);
for(int i = s; i <= t; i ++) dep[i] = 0;
dep[s] = 1;
while(!q.empty()){
int u = q.front(); q.pop();
for(int i = head[u]; i; i = e[i].nex){
int v = e[i].to;
if(e[i].w && !dep[v]){
dep[v] = dep[u] + 1;
if(v == t) return true;
q.push(v);
}
}
}
return dep[t];
}
int dfs(int u, int inflow){
if(u == t) return inflow;
int outflow = 0;
for(int i = head[u]; i && inflow; i = e[i].nex){
int v = e[i].to;
if(e[i].w && dep[v] == dep[u] + 1){
int flow = dfs(v, min(e[i].w, inflow));
e[i].w -= flow;
e[i ^ 1].w += flow;
inflow -= flow;
outflow += flow;
}
}
if(outflow == 0) dep[u] = 0;
return outflow;
}
vector<int>g[N * N];//存坐标所有的矩形左上角编号
void init(){
tot = 1;
for(int i = 1; i <= m; i ++){
g[tr[i].get_id()].clear();
}
for(int i = s; i <= t; i ++) head[i] = 0;
}
int overlap_area(Box &a, Box &b){//计算线段交
int x = max(0, min(a.x2, b.x2) - max(a.x1, b.x1));
int y = max(0, min(a.y2, b.y2) - max(a.y1, b.y1));
return x * y;
}
void build_edge(int u){
Box &A = de[u];
for(int i = max(A.x1 - 2, 0); i <= A.x2 + 2; i ++){
for(int j = max(A.y1 - 2, 0); j <= A.y2 + 2; j ++){
int idx = i * SIZ + j;
for(int id : g[idx]){
if(overlap_area(A, tr[id]) * 2 >= tr[id].area){
add(u, id + n, 1);
add(id + n, u, 0);
}
}
}
}
}
void solve(){
init();
cin >> m >> n;
for(int i = 1; i <= m; i ++){//真实框
tr[i].input();
tr[i].get_area();
g[tr[i].get_id()].push_back(i);
}
for(int i = 1; i <= n; i ++){// 搜索框
de[i].input();
de[i].get_area();
}
s = 0, t = n + m + 1;
for(int i = 1; i <= n; i ++){
add(s, i, 1);
add(i, s, 0);
build_edge(i);
}
for(int i = 1; i <= m; i ++){
add(i + n, t, 1);
add(t, i + n, 0);
}
int ans = 0;
while(bfs()){
ans += dfs(s, inf);
}
cout << ans << "\n";
}
int main(){
ios::sync_with_stdio(false);
cout.tie(NULL);
int t;
cin >> t;
while(t --){
solve();
}
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 24ms
memory: 98064kb
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: 15ms
memory: 98220kb
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: 2459ms
memory: 121524kb
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: 273ms
memory: 127712kb
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: 254ms
memory: 105164kb
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: 3808ms
memory: 120192kb
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: 16ms
memory: 98060kb
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'