結果

問題 No.1479 Matrix Eraser
ユーザー SHIJOUSHIJOU
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

//#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;
}
0