QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#113690#5668. Cell Nuclei DetectionUFRJ#RE 0ms0kbC++202.7kb2023-06-19 04:22:462023-06-19 04:22:47

Judging History

This is the latest submission verdict.

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-06-19 04:22:47]
  • Judged
  • Verdict: RE
  • Time: 0ms
  • Memory: 0kb
  • [2023-06-19 04:22:46]
  • Submitted

answer

#include "bits/stdc++.h"

#define watch(x) cout << (#x) << " is " << (x) << endl

using namespace std;
using lint = int64_t;

const int mx = 2e3 + 1;

vector<vector<set<int>>> grid;

struct bipartite_matching {
    int N, M, T;
    vector<vector<int>> adj;
    vector<int> match, seen;
    bipartite_matching(int a, int b) : N(a), M(a + b), adj(M), 
            match(M, -1), seen(M, -1), T(0) {}
    void add_edge(int a, int b) {
        assert(0 <= a && a < N && b + N < M && N <= b + N);
        adj[a].push_back(b + N);
    }
    void shuffle_edges() { // useful to break some hairy tests
        mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
        for (auto& cur : adj)
            shuffle(cur.begin(), cur.end(), rng);
    } 
    bool dfs(int cur) {
        if (seen[cur] == T) return false;
        seen[cur] = T;
        for (int nxt : adj[cur]) 
            if (match[nxt] == -1) {
                match[nxt] = cur, match[cur] = nxt;
                return true;
            }
        for (int nxt : adj[cur]) 
            if (dfs(match[nxt])) {
                match[nxt] = cur, match[cur] = nxt;
                return true;
            }
        return false;
    }
    int solve() {
        int res = 0;
        while (true) {
            int cur = 0; ++T;
            for (int i = 0; i < N; ++i) 
                if (match[i] == -1) cur += dfs(i);
            if (cur == 0) break;
            else res += cur;
        }
        return res;
    }
};

int main() {
  cin.tie(nullptr)->sync_with_stdio(false); 

  int t, n, m;
  cin >> t;
  while(t--) {
    grid.clear();
    grid.assign(mx, vector<set<int>>(mx));
    map<int,int> area;
    cin >> m >> n;
    bipartite_matching bp(n+1, m+1);
    int x1, x2, y1, y2, seq = 0;
    for(int k = 0; k < n; k++) {
      cin >> x1 >> y1 >> x2 >> y2;
      int id = ++seq;
      int curArea = 0;
      for(int i = x1; i <= x2; i++) {
        for(int j = y1; j <= y2; j++) {
          grid[i][j].insert(id);
          curArea++;
        }
      }
      area[id] = curArea;
    }
    for(int k = 0; k < m; k++) {
      cin >> x1 >> y1 >> x2 >> y2;
      unordered_map<int,int> interArea;
      for(int i = x1; i <= x2; i++) {
        for(int j = y1; j <= y2; j++) {
          for(auto &id : grid[i][j]) {
            interArea[id]++;
          }
        }
      }
      for(auto &id : interArea) {
        int key = id.first, curArea = id.second;
        if(curArea*2 >= area[key]) {
          // cout << "Green: " << k << " Red: " << key << endl;
          bp.add_edge(k, key);
        }
      }
    }
    cout << bp.solve() << endl;
  }
  return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 0
Dangerous Syscalls

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

result: