結果

問題 No.1194 Replace
ユーザー Mayimg
提出日時 2020-08-25 02:50:46
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,128 ms / 2,000 ms
コード長 3,051 bytes
コンパイル時間 3,053 ms
コンパイル使用メモリ 217,340 KB
最終ジャッジ日時 2025-01-13 13:45:39
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;

map<int, int> zip;
vector<int> unzip;
template <typename T> void compress(vector<T>& v) {
  for (int i = 0; i < (int) v.size(); i++) unzip.push_back(v[i]);
  sort(unzip.begin(), unzip.end());
  unzip.erase(unique(unzip.begin(), unzip.end()), unzip.end());
  for (int i = 0; i < (int) unzip.size(); i++) zip[unzip[i]] = i;
}

class StronglyConnectedComponents {
private:
  int siz;
  int itr = 0;
  vector<bool> done;
  vector<int> order, comp;
  void dfs (int cur, const vector<vector<int>>& g) {
    done[cur] = true;
    for (int nxt : g[cur]) {
      if (!done[nxt]) dfs(nxt, g);
    }
    order.push_back(cur);
  }
  void rdfs (int cur, const vector<vector<int>>& rg) {
    done[cur] = true;
    comp[cur] = itr;
    for (int nxt : rg[cur]) {
      if (!done[nxt]) rdfs(nxt, rg);
    }
  }
  void build (const vector<vector<int>>& g, const vector<vector<int>>& rg) {
    for (int i = 0; i < siz; ++i) if (!done[i]) dfs (i, g);
    reverse(order.begin(), order.end());
    fill(done.begin(), done.end(), false);
    for (int i = 0; i < siz; ++i) if (!done[order[i]]) rdfs (order[i], rg), ++itr;
  }

public:
  StronglyConnectedComponents(const vector<vector<int>>& g, const vector<vector<int>>& rg, const int& n) : siz(n), done(siz, false), comp(siz, -1) {
    build(g, rg);
  }
  bool  has_cycle() { return *max_element(comp.begin(), comp.end()) < siz - 1; }
  vector<int> get_comp() { return comp;}
};
using SCC = StronglyConnectedComponents;

int dfs (int cur, vector<set<int>>& g, vector<int>& val, vector<int>& vis) {
  if (vis[cur]) return val[cur];
  vis[cur] = 1;
  for (int nxt : g[cur]) {
    val[cur] = max(val[cur], dfs(nxt, g, val, vis));
  }
  return val[cur];
}

signed main() { 
  ios::sync_with_stdio(false); cin.tie(0);
  int n, m;
  cin >> n >> m;
  vector<int> all;
  vector<int> b(m), c(m);
  for (int i = 0; i < m; i++) {
    cin >> b[i] >> c[i];
    all.push_back(b[i]);
    all.push_back(c[i]);
  }
  compress(all);
  for (int i = 0; i < m; i++) {
    b[i] = zip[b[i]];
    c[i] = zip[c[i]];
  }
  int nn = zip.size();
  vector<vector<int>> g(nn), rg(nn);
  for (int i = 0; i < m; i++) {
    g[b[i]].push_back(c[i]);
    rg[c[i]].push_back(b[i]);
  }
  SCC scc(g, rg, nn);
  vector<int> comp = scc.get_comp();
  int nnn = *max_element(comp.begin(), comp.end()) + 1;
  vector<set<int>> gg(nnn);
  vector<int> val(nnn);
  for (int i = 0; i < m; i++) {
    assert(comp[b[i]] < nnn);
    assert(comp[c[i]] < nnn);
    val[comp[b[i]]] = max(val[comp[b[i]]], unzip[b[i]]);
    val[comp[c[i]]] = max(val[comp[c[i]]], unzip[c[i]]);
    if (comp[b[i]] == comp[c[i]]) continue;
    gg[comp[b[i]]].insert(comp[c[i]]);
  }
  vector<int> vis(nnn);
  for (int i = 0; i < nnn; i++) dfs (i, gg, val, vis);
  long long ans = 1LL * n * (n + 1) / 2;
  set<int> app;
  for (auto x : b) {
    if (app.count(x)) continue;
    app.insert(x);
    int y = unzip[x];
    ans -= y;
    assert(y <= val[comp[x]]);
    ans += val[comp[x]];
  }
  cout << ans << endl;
  return 0;
}
0