QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#186230#5668. Cell Nuclei DetectionBUET_POTATOES#TL 2335ms554276kbC++203.4kb2023-09-23 14:27:192023-09-23 14:27:20

Judging History

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

  • [2023-09-23 14:27:20]
  • 评测
  • 测评结果:TL
  • 用时:2335ms
  • 内存:554276kb
  • [2023-09-23 14:27:19]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int nmax = 2e3+2;

vector<int>bam[nmax][nmax];
vector<int>dan[nmax][nmax];

const int N = 5e4+5;
int X1[2][N], X2[2][N], Y1[2][N], Y2[2][N];

int interSecc(int L1, int R1, int L2, int R2){
    int L = max(L1, L2), R = min(R1, R2);
    if(L<=R) return R-L;
    return 0;
}

bool hasEdge(int i, int j){
    int aria = (X2[0][i] - X1[0][i]) * (Y2[0][i] - Y1[0][i]);
    int interAria = interSecc(X1[0][i], X2[0][i], X1[1][j], X2[1][j]) * interSecc(Y1[0][i], Y2[0][i], Y1[1][j], Y2[1][j]);
    if(interAria*2 >= aria){
        return true;
    }
    return false;
}

namespace HopcroftKarp { // 1-indexed. L,R indep
const int maxN=1e6+7,maxM=1e6+7; int n, m, match;
int vis[maxN], level[maxN], ml[maxN], mr[maxM];
vector<int> edge[maxN];
void init(int N, int M) { //N=left nodes, M = right
  n = N, m = M;
  for (int i=1;i<=n;i++) edge[i].clear(),ml[i]=-1;
  for (int i=1;i<=m;i++)  mr[i] = -1;
  match = 0;
}
void add(int u, int v) {
  edge[u].push_back(v);
  if(ml[u]==-1&&mr[v]==-1)ml[u]=v,mr[v]=u,match++;
}
bool dfs(int u) {
  vis[u] = true;
  for (int x: edge[u]) {
    int v = mr[x];
    if (v == -1 || (!vis[v]
              && level[u] < level[v] && dfs(v))) {
      ml[u] = x; mr[x] = u; return true;
    }
  }
  return false;
}
int matching() {
  while (true) {
    queue<int> q;
    for (int i = 1; i <= n; ++i) {
      if (ml[i] == -1)  level[i] = 0, q.push(i);
      else        level[i] = -1;
    }
    while (!q.empty()) {
      int u = q.front(); q.pop();
      for (int x: edge[u]) {
        int v = mr[x];
        if (v != -1 && level[v] < 0) {
          level[v] = level[u] + 1; q.push(v);
        }
      }
    }
    for (int i = 1; i <= n; ++i) vis[i] = false;
    int d = 0;
    for(int i=1;i<=n;++i)if(ml[i]==-1&&dfs(i))++d;
    if (d == 0) return match; match += d;
  }
} }

void solv(int ces){
    int m,n;
    cin>>m>>n;
    for(int i=1; i<=m; i++){
        int x1,y1,x2,y2;
        cin>>x1>>y1>>x2>>y2;
        X1[0][i] = x1;
        X2[0][i] = x2;
        Y1[0][i] = y1;
        Y2[0][i] = y2;

        for(int x=x1; x<=x2; x++){
            for(int y=y1; y<=y2; y++){
                bam[x][y].push_back(i);
            }
        }
    }

    for(int i=1; i<=n; i++){
        int x1,y1,x2,y2;
        cin>>x1>>y1>>x2>>y2;
        X1[1][i] = x1;
        X2[1][i] = x2;
        Y1[1][i] = y1;
        Y2[1][i] = y2;

        for(int x=x1; x<=x2; x++){
            for(int y=y1; y<=y2; y++){
                dan[x][y].push_back(i);
            }
        }
    }
    vector<pii>allP;
    for(int i=0; i<nmax; i++){
        for(int j=0; j<nmax; j++){
            for(int a : bam[i][j]){
                for(int b : dan[i][j]){
                    allP.emplace_back(a, b);
                }
            }
            bam[i][j].clear();
            dan[i][j].clear();
        }
    }
    sort(allP.begin(), allP.end());
    allP.erase(unique(allP.begin(), allP.end()), allP.end());

    HopcroftKarp::init(m, n);

    for(auto [i,j] : allP){
        if(hasEdge(i, j)){
            HopcroftKarp::add(i, j);
        }
    }

    cout<< HopcroftKarp::matching() <<"\n";
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(NULL);

    int tc = 1;
    cin>>tc;
    for(int ces=1; ces<=tc; ces++){
        solv(ces);
    }

    return 0;
}




Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 68ms
memory: 223828kb

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: 67ms
memory: 223828kb

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: 2335ms
memory: 554276kb

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: 856ms
memory: 327112kb

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: 493ms
memory: 290160kb

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: -100
Time Limit Exceeded

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:


result: