結果
問題 | No.1479 Matrix Eraser |
ユーザー | SHIJOU |
提出日時 | 2021-04-16 22:38:10 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 205 ms / 3,000 ms |
コード長 | 3,237 bytes |
コンパイル時間 | 2,423 ms |
コンパイル使用メモリ | 220,104 KB |
実行使用メモリ | 23,904 KB |
最終ジャッジ日時 | 2023-09-16 01:40:53 |
合計ジャッジ時間 | 8,472 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 8 ms
15,044 KB |
testcase_01 | AC | 8 ms
15,156 KB |
testcase_02 | AC | 8 ms
15,024 KB |
testcase_03 | AC | 8 ms
15,072 KB |
testcase_04 | AC | 8 ms
15,116 KB |
testcase_05 | AC | 8 ms
15,188 KB |
testcase_06 | AC | 8 ms
15,068 KB |
testcase_07 | AC | 29 ms
16,040 KB |
testcase_08 | AC | 47 ms
16,664 KB |
testcase_09 | AC | 90 ms
18,252 KB |
testcase_10 | AC | 164 ms
20,424 KB |
testcase_11 | AC | 105 ms
18,852 KB |
testcase_12 | AC | 37 ms
16,280 KB |
testcase_13 | AC | 47 ms
16,560 KB |
testcase_14 | AC | 38 ms
16,228 KB |
testcase_15 | AC | 14 ms
15,312 KB |
testcase_16 | AC | 41 ms
16,436 KB |
testcase_17 | AC | 199 ms
21,404 KB |
testcase_18 | AC | 202 ms
21,292 KB |
testcase_19 | AC | 201 ms
21,296 KB |
testcase_20 | AC | 201 ms
21,472 KB |
testcase_21 | AC | 202 ms
21,388 KB |
testcase_22 | AC | 199 ms
21,360 KB |
testcase_23 | AC | 197 ms
21,412 KB |
testcase_24 | AC | 200 ms
21,472 KB |
testcase_25 | AC | 205 ms
21,332 KB |
testcase_26 | AC | 201 ms
21,372 KB |
testcase_27 | AC | 166 ms
17,936 KB |
testcase_28 | AC | 163 ms
17,964 KB |
testcase_29 | AC | 165 ms
18,088 KB |
testcase_30 | AC | 164 ms
17,956 KB |
testcase_31 | AC | 167 ms
18,020 KB |
testcase_32 | AC | 78 ms
23,864 KB |
testcase_33 | AC | 78 ms
23,820 KB |
testcase_34 | AC | 78 ms
23,904 KB |
testcase_35 | AC | 78 ms
23,828 KB |
testcase_36 | AC | 77 ms
23,708 KB |
testcase_37 | AC | 25 ms
16,888 KB |
testcase_38 | AC | 145 ms
17,344 KB |
testcase_39 | AC | 186 ms
22,892 KB |
testcase_40 | AC | 8 ms
15,132 KB |
ソースコード
//#define _GLIBCXX_DEBUG #include <bits/stdc++.h> using namespace std; #define rep(i, n) for(int i=0; i<n; ++i) #define all(v) v.begin(), v.end() #define rall(v) v.rbegin(), v.rend() using ll = int64_t; using ull = uint64_t; using ld = long double; using P = pair<int, int>; using vs = vector<string>; using vi = vector<int>; using vvi = vector<vi>; template<class T> using PQ = priority_queue<T>; template<class T> using PQG = priority_queue<T, vector<T>, greater<T>>; const int INF = 0xccccccc; const ll LINF = 0xcccccccccccccccLL; template<typename T1, typename T2> inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);} template<typename T1, typename T2> inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);} template<typename T1, typename T2> istream &operator>>(istream &is, pair<T1, T2> &p) { return is >> p.first >> p.second;} template<typename T1, typename T2> ostream &operator<<(ostream &os, const pair<T1, T2> &p) { return os << p.first << ' ' << p.second;} struct edge{ int to, rev; int64_t cap; edge(int t=0, int64_t c=0, int r=0):to(t), cap(c), rev(r) {} }; struct max_flow { vector<vector<edge> > G; vector<int> level, iter; int n; max_flow(int n=0):n(n), G(n), level(n), iter(n) {} inline int add(int from, int to, int64_t cap) { G[from].emplace_back(to, cap, G[to].size()); G[to].emplace_back(from, 0, G[from].size()-1); return G[from].size()-1; } void bfs(int s) { for(int i = 0; i < n; i++) level[i] = -1; queue<int> que; level[s] = 0; que.emplace(s); while(!que.empty()) { int now = que.front(); que.pop(); for(edge &e:G[now]) { if(e.cap > 0 and level[e.to] < 0) { level[e.to] = level[now] + 1; que.emplace(e.to); } } } } int64_t dfs(int v, int t, int64_t f) { if(v == t) return f; for(int &i = iter[v]; i < G[v].size(); i++) { edge &e = G[v][i]; if(e.cap > 0 and level[v] < level[e.to]) { int64_t d = dfs(e.to, t, min(f, e.cap)); if(d > 0) { e.cap -= d; G[e.to][e.rev].cap += d; return d; } } } return 0; } int64_t Dinic(int s, int t) { int64_t flow = 0; while(true) { bfs(s); if(level[t] < 0) return flow; for(int i = 0; i < n; i++) iter[i] = 0; int64_t f; while((f = dfs(s, t, INT64_MAX)) > 0) flow += f; } } }; //head #define N 500010 vector<P> _[N]; int main() { ios::sync_with_stdio(false); cin.tie(0); int h, w; cin >> h >> w; rep(i, h) rep(j, w) { int a; cin >> a; _[a].emplace_back(i, j); } int ans = 0; for(int i = N-1; i > 0; i--) if(not _[i].empty()) { vi lst; for(auto [a, b]:_[i]) { lst.emplace_back(a); lst.emplace_back(b+N); } sort(all(lst)); lst.erase(unique(all(lst)), lst.end()); const int s = lst.size(), t = s+1; max_flow mfl(t+1); for(auto [a, b]:_[i]) { int ia = lower_bound(all(lst), a)-lst.begin(); int ib = lower_bound(all(lst), b+N)-lst.begin(); mfl.add(ia, ib, 1); } rep(i, s) { if(lst[i] >= N) mfl.add(i, t, 1); else mfl.add(s, i, 1); } ans += mfl.Dinic(s, t); } cout << ans << endl; }