QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#657653#5668. Cell Nuclei Detectionucup-team5234#WA 1044ms144992kbC++234.2kb2024-10-19 15:09:552024-10-19 15:09:58

Judging History

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

  • [2024-10-19 15:09:58]
  • 评测
  • 测评结果:WA
  • 用时:1044ms
  • 内存:144992kb
  • [2024-10-19 15:09:55]
  • 提交

answer

#include "bits/stdc++.h"
using namespace std;
using ll = long long;
using vvi = vector<vector<int>>;
using vb = vector<bool>;
using vi = vector<int>;

namespace atcoder {
    template <class Cap> struct mf_graph {
        public:
        mf_graph(int n) : _n(n), g(n) {}

        int add_edge(int from, int to, Cap cap){
            int m = int(pos.size());
            pos.push_back({from, int(g[from].size())});
            int from_id = int(g[from].size());
            int to_id = int(g[to].size());
            if(from == to)to_id++;
            g[from].push_back(_edge{to, to_id, cap});
            g[to].push_back(_edge{from, from_id, 0});
            return m;
        }

        struct edge {
            int from, to;
            Cap cap, flow;
        };

        Cap flow(int s, int t) {
            vi level(_n), iter(_n);
            queue<int> que;
            auto bfs = [&](){
                fill(level.begin(),level.end(), -1);
                level[s] = 0;
                while(!que.empty())que.pop();
                que.push(s);
                while(!que.empty()){
                    int v = que.front();
                    que.pop();
                    for(auto e : g[v]){
                        if(e.cap == 0 || level[e.to] >= 0)continue;
                        level[e.to] = level[v] + 1;
                        if(e.to== t)return;
                        que.push(e.to);
                    }
                }
            };
            auto dfs = [&](auto self, int v, Cap up){
                if(v == s)return up;
                Cap res = 0;
                int level_v = level[v];
                for(int& i = iter[v]; i < int(g[v].size()); i++){
                    _edge& e = g[v][i];
                    if(level_v <= level[e.to] || g[e.to][e.rev].cap == 0)continue;
                    Cap d = self(self, e.to, min(up - res, g[e.to][e.rev].cap));
                    if(d <= 0)continue;
                    g[v][i].cap += d;
                    g[e.to][e.rev].cap -= d;
                    res += d;
                    if(res == up)break;
                }
                return res;
            };
            Cap flow = 0;
            while(true){
                bfs();
                if(level[t] == -1)break;
                fill(iter.begin(), iter.end(), 0);
                while(true) {
                    Cap f = dfs(dfs, t, numeric_limits<Cap>::max());
                    if(!f)break;
                    flow += f;
                }
            } 
            return flow;
        }
        private:
    int _n;
    struct _edge {
        int to, rev;
        Cap cap;
    };
    vector<pair<int,int>> pos;
    vector<vector<_edge>> g;
    };
};

using namespace atcoder;

int main(){
    int t;
    cin >> t;
    while(t -- ){
        int m,n;
        cin >> m >> n;
        vector<vector<vi>> v(2010,vvi(2010));
        vector<int> sz(m);
        for(int i = 0; i < m; i++){
            int x1,y1,x2,y2;
            cin >> x1 >> y1 >> x2 >> y2;
            sz[i] = (x2-x1 + 1) * (y2-y1 + 1);
            for(int x = x1; x <= x2; x++){
                for(int y = y1; y <= y2; y++){
                    v[x][y].push_back(i);
                }
            }
        }
        vector<int> cnt(m);
        mf_graph<ll> g(m + n + 2);
        int s = m + n;
        int t = s + 1;
        for(int i = 0; i < m; i++)g.add_edge(s, i, 1);
        for(int i = 0; i < n; i++)g.add_edge(i + m, t, 1);
        for(int i = 0; i < n; i++){
            int x1,y1,x2,y2;
            cin >> x1 >> y1 >> x2 >> y2;
            set<int> S;
            for(int x = x1; x <= x2; x++){
                for(int y = y1; y <= y2; y++){
                    for(auto e: v[x][y]){
                        S.insert(e);
                        cnt[e]++;
                    }
                }
            }
            for(auto e: S){
                // cout << e << " " << i << " " << cnt[e] << " " << sz[e] << endl;
                if(sz[e] <= 2 * cnt[e]){
                    // cout << "!" << e << " " << i << endl;
                    g.add_edge(e, m + i, 1);
                }
                cnt[e] = 0;
            }
        }
        cout << g.flow(s, t) << endl;
    }
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 28ms
memory: 97880kb

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: 52ms
memory: 97948kb

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: -100
Wrong Answer
time: 1044ms
memory: 144992kb

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
3213

result:

wrong answer 5th lines differ - expected: '3150', found: '3213'