結果

問題 No.1920 Territory
ユーザー 👑 hos.lyrichos.lyric
提出日時 2022-04-29 22:48:13
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 718 ms / 5,000 ms
コード長 5,036 bytes
コンパイル時間 1,426 ms
コンパイル使用メモリ 121,276 KB
実行使用メモリ 60,620 KB
最終ジャッジ日時 2023-09-11 14:12:33
合計ジャッジ時間 13,327 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 34 ms
28,884 KB
testcase_04 AC 30 ms
28,780 KB
testcase_05 AC 31 ms
28,888 KB
testcase_06 AC 32 ms
28,784 KB
testcase_07 AC 31 ms
28,716 KB
testcase_08 AC 7 ms
4,376 KB
testcase_09 AC 7 ms
4,376 KB
testcase_10 AC 7 ms
4,376 KB
testcase_11 AC 6 ms
4,384 KB
testcase_12 AC 7 ms
4,380 KB
testcase_13 AC 665 ms
60,484 KB
testcase_14 AC 681 ms
60,620 KB
testcase_15 AC 659 ms
60,448 KB
testcase_16 AC 677 ms
60,584 KB
testcase_17 AC 718 ms
60,556 KB
testcase_18 AC 607 ms
41,772 KB
testcase_19 AC 578 ms
41,856 KB
testcase_20 AC 580 ms
41,868 KB
testcase_21 AC 571 ms
41,764 KB
testcase_22 AC 563 ms
41,756 KB
testcase_23 AC 401 ms
43,468 KB
testcase_24 AC 222 ms
40,508 KB
testcase_25 AC 172 ms
47,488 KB
testcase_26 AC 171 ms
47,516 KB
testcase_27 AC 181 ms
49,176 KB
testcase_28 AC 177 ms
49,176 KB
testcase_29 AC 321 ms
55,128 KB
testcase_30 AC 319 ms
55,284 KB
testcase_31 AC 323 ms
42,984 KB
testcase_32 AC 323 ms
43,500 KB
testcase_33 AC 303 ms
48,536 KB
testcase_34 AC 339 ms
48,876 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }


int uf[200'010];
int root(int u) {
  return (uf[u] < 0) ? u : (uf[u] = root(uf[u]));
}
bool connect(int u, int v) {
  u = root(u);
  v = root(v);
  if (u == v) return false;
  if (uf[u] > uf[v]) swap(u, v);
  uf[u] += uf[v];
  uf[v] = u;
  return true;
}


template <class T> void bAdd(vector<T> &bit, int pos, const T &val) {
  const int bitN = bit.size();
  for (int x = pos; x < bitN; x |= x + 1) bit[x] += val;
}
template <class T> T bSum(const vector<T> &bit, int pos) {
  T ret = 0;
  for (int x = pos - 1; x >= 0; x = (x & (x + 1)) - 1) ret += bit[x];
  return ret;
}
template <class T> T bSum(const vector<T> &bit, int pos0, int pos1) {
  return bSum(bit, pos1) - bSum(bit, pos0);
}


int M, N;
vector<int> A, B, C;
vector<int> P, Q, R;

int main() {
  for (; ~scanf("%d%d", &M, &N); ) {
    A.resize(M);
    B.resize(M);
    C.resize(M);
    for (int i = 0; i < M; ++i) scanf("%d%d%d", &A[i], &B[i], &C[i]);
    P.resize(N);
    Q.resize(N);
    R.resize(N);
    for (int j = 0; j < N; ++j) scanf("%d%d%d", &P[j], &Q[j], &R[j]);
    
    int lim = 0;
    for (int i = 0; i < M; ++i) {
      chmax(lim, A[i]);
      chmax(lim, B[i]);
      chmax(lim, C[i]);
    }
    for (int j = 0; j < N; ++j) {
      chmax(lim, P[j]);
      chmax(lim, Q[j]);
      chmax(lim, R[j]);
    }
    lim += 5;
    
    Int numComps = 0;
    {
      int segN;
      for (segN = 1; segN < lim; segN <<= 1) {}
      vector<vector<int>> iss(segN << 1), jss(segN << 1);
      for (int i = 0; i < M; ++i) {
        // x \in [B[i], C[i]]
        int a = B[i], b = C[i] + 1;
        for (a += segN, b += segN; a < b; a >>= 1, b >>= 1) {
          if (a & 1) iss[a++].push_back(i);
          if (b & 1) iss[--b].push_back(i);
        }
      }
      for (int j = 0; j < N; ++j) {
        int a = P[j];
        for (a += segN; a; a >>= 1) {
          jss[a].push_back(j);
        }
      }
      
      fill(uf, uf + (M + N), -1);
      for (int a = 1; a < segN << 1; ++a) if (!iss[a].empty() && !jss[a].empty()) {
// cerr<<"is = ";pv(iss[a].begin(),iss[a].end());
// cerr<<"js = ";pv(jss[a].begin(),jss[a].end());
        vector<pair<pair<int, int>, int>> events;
        for (const int i : iss[a]) {
          events.emplace_back(make_pair(A[i], 1), i);
        }
        for (const int j : jss[a]) {
          events.emplace_back(make_pair(Q[j], 0), j);
        }
        sort(events.begin(), events.end());
        using Entry = pair<int, int>;
        priority_queue<Entry, vector<Entry>, greater<Entry>> que;
        for (const auto &e : events) {
          if (e.first.second) {
            const int i = e.second;
            for (; !que.empty() && que.top().first < A[i]; que.pop()) {}
            if (!que.empty()) {
              for (; ; ) {
                const auto p = que.top();
                que.pop();
                connect(i, M + p.second);
                if (que.empty()) {
                  que.push(p);
                  break;
                }
              }
            }
          } else {
            const int j = e.second;
            que.emplace(R[j], j);
          }
        }
      }
      for (int u = 0; u < M + N; ++u) if (uf[u] < 0) {
        ++numComps;
      }
    }
// cerr<<"uf = ";pv(uf,uf+(M+N));
cerr<<"numComps = "<<numComps<<endl;
    
    Int numPoints = 0;
    {
      vector<vector<int>> adds(lim), rems(lim), jss(lim);
      for (int i = 0; i < M; ++i) {
        adds[B[i]].push_back(i);
        rems[C[i]].push_back(i);
      }
      for (int j = 0; j < N; ++j) {
        jss[P[j]].push_back(j);
      }
      vector<int> bit(lim, 0);
      for (int x = 0; x < lim; ++x) {
        for (const int i : adds[x]) {
          bAdd(bit, A[i], +1);
        }
        for (const int j : jss[x]) {
          Int num = 0;
          num += bSum(bit, R[j] + 1);
          num -= bSum(bit, Q[j]);
          numPoints += num;
        }
        for (const int i : rems[x]) {
          bAdd(bit, A[i], -1);
        }
      }
    }
cerr<<"numPoints = "<<numPoints<<endl;
    
    Int ans = 0;
    ans += numPoints;
    ans -= (M + N);
    ans += numComps;
    printf("%lld\n", ans);
  }
  return 0;
}
0