QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#650603#5668. Cell Nuclei DetectionXfJbUhpyzgaW#AC ✓967ms77004kbC++142.8kb2024-10-18 15:48:052024-10-18 15:48:05

Judging History

你现在查看的是最新测评结果

  • [2024-10-18 15:48:05]
  • 评测
  • 测评结果:AC
  • 用时:967ms
  • 内存:77004kb
  • [2024-10-18 15:48:05]
  • 提交

answer

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+5,maxn=1e7;
struct mat{int x1,y1,x2,y2;}a[N],b[N];
int head[maxn],go[maxn],nxt[maxn],flow[maxn],dis[maxn],cur[maxn];
vector<int> mp[505][505];
int n=0,m,T,tot,St,Ed;
bool check(mat a,mat b){
    int x1=max(a.x1,b.x1),x2=min(a.x2,b.x2),y1=max(a.y1,b.y1),y2=min(a.y2,b.y2);
    if (x1>x2 || y1>y2) return 0;
    return 2*(x2-x1)*(y2-y1)>=(a.x2-a.x1)*(a.y2-a.y1);
}
void add(int x,int y,int z){
    // printf("%d %d %d\n",x,y,z);
    ++tot;go[tot]=y;nxt[tot]=head[x];head[x]=tot;flow[tot]=z;
    ++tot;go[tot]=x;nxt[tot]=head[y];head[y]=tot;flow[tot]=0;
}
vector<pair<int,int>> edge;
void link(int x,int y,int id){
    for (int i:mp[x][y])
        if (check(a[i],b[id]))
            edge.emplace_back(i,id+n);
}
void clear(){
    for (int i=0;i<=500;++i)
        for (int j=0;j<=500;++j)
            mp[i][j].clear();
    tot=1;St=n+m+1;Ed=n+m+2;
    for (int i=1;i<=n+m+2;++i) head[i]=0;
    edge.clear();
}
queue<int> q;
bool bfs(){
    for (int i=1;i<=n+m+2;++i) dis[i]=-1;
    dis[St]=1;q.emplace(St);
    while (!q.empty()){
        int x=q.front();q.pop();
        for (int i=head[x];i;i=nxt[i])
            if (flow[i] && dis[go[i]]==-1){
                dis[go[i]]=dis[x]+1;
                q.emplace(go[i]);
            }
    }
    return dis[Ed]!=-1;
}
int dfs(int x,int flw){
    if (x==Ed) return flw;
    int sum=0;
    for (int &i=cur[x];i;i=nxt[i])
        if (dis[go[i]]==dis[x]+1 && flow[i]){
            int t=dfs(go[i],min(flw-sum,flow[i]));
            flow[i]-=t,flow[i^1]+=t,sum+=t;
            if (sum==flw) return flw;
        }
    return sum;
}
int dinic(){
    int res=0;
    while (bfs()){
        for (int i=1;i<=n+m+2;++i) cur[i]=head[i];
        res+=dfs(St,1e8);
        // printf("%d\n",res);
    }
    return res;
}
void work(){
    scanf("%d%d",&n,&m);
    clear();
    for (int i=1;i<=n;++i){
        scanf("%d%d%d%d",&a[i].x1,&a[i].y1,&a[i].x2,&a[i].y2);
        mp[a[i].x1/4][a[i].y1/4].emplace_back(i);
        mp[a[i].x1/4][a[i].y2/4].emplace_back(i);
        mp[a[i].x2/4][a[i].y1/4].emplace_back(i);
        mp[a[i].x2/4][a[i].y2/4].emplace_back(i);
    }
    for (int i=1;i<=m;++i){
        scanf("%d%d%d%d",&b[i].x1,&b[i].y1,&b[i].x2,&b[i].y2);
        link(b[i].x1/4,b[i].y1/4,i);
        link(b[i].x1/4,b[i].y2/4,i);
        link(b[i].x2/4,b[i].y1/4,i);
        link(b[i].x2/4,b[i].y2/4,i);
    }
    sort(edge.begin(),edge.end());
    edge.erase(unique(edge.begin(),edge.end()),edge.end());
    for (auto e:edge) add(e.first,e.second,1);
    for (int i=1;i<=n;++i) add(St,i,1);
    for (int i=1;i<=m;++i) add(i+n,Ed,1);
    printf("%d\n",dinic());    
}
int main(){
    // freopen("a.in","r",stdin);
    for (scanf("%d",&T);T;--T) work();
    return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 23968kb

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: 0ms
memory: 23612kb

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: 288ms
memory: 63300kb

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: 317ms
memory: 74868kb

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: 122ms
memory: 35624kb

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: 967ms
memory: 77004kb

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: 20756kb

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'