結果

問題 No.1490 スライムと爆弾
ユーザー 👑 emthrmemthrm
提出日時 2021-04-23 23:10:57
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 608 ms / 2,000 ms
コード長 3,715 bytes
コンパイル時間 2,227 ms
コンパイル使用メモリ 211,496 KB
実行使用メモリ 130,560 KB
最終ジャッジ日時 2024-07-04 08:42:06
合計ジャッジ時間 7,359 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 1 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 1 ms
5,376 KB
testcase_13 AC 194 ms
58,880 KB
testcase_14 AC 223 ms
76,288 KB
testcase_15 AC 129 ms
64,384 KB
testcase_16 AC 140 ms
39,680 KB
testcase_17 AC 111 ms
74,880 KB
testcase_18 AC 232 ms
76,416 KB
testcase_19 AC 197 ms
70,272 KB
testcase_20 AC 146 ms
58,752 KB
testcase_21 AC 93 ms
56,576 KB
testcase_22 AC 171 ms
45,952 KB
testcase_23 AC 2 ms
6,940 KB
testcase_24 AC 2 ms
6,944 KB
testcase_25 AC 327 ms
130,560 KB
testcase_26 AC 323 ms
130,432 KB
testcase_27 AC 325 ms
130,560 KB
testcase_28 AC 288 ms
128,640 KB
testcase_29 AC 346 ms
128,640 KB
testcase_30 AC 608 ms
130,560 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 1000000007;
// constexpr int MOD = 998244353;
constexpr int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1};
constexpr int dy8[] = {1, 1, 0, -1, -1, -1, 0, 1}, dx8[] = {0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U> inline bool chmax(T &a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U> inline bool chmin(T &a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
  IOSetup() {
    std::cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << fixed << setprecision(20);
  }
} iosetup;

template <typename Abelian>
struct BIT2DRangeAdd {
  BIT2DRangeAdd(int height_, int width_, const Abelian ID = 0) : height(height_), width(width_), ID(ID) {
    ++height; ++width;
    dat_const.assign(height, std::vector<Abelian>(width, ID));
    dat_linear[0].assign(height, std::vector<Abelian>(width, ID));
    dat_linear[1].assign(height, std::vector<Abelian>(width, ID));
    dat_quadratic.assign(height, std::vector<Abelian>(width, ID));
  }

  void add(int y1, int x1, int y2, int x2, const Abelian val) {
    ++y1; ++x1; ++y2; ++x2;
    for (int i = y1; i < height; i += i & -i) for (int j = x1; j < width; j += j & -j) {
      dat_const[i][j] += val * (y1 - 1) * (x1 - 1);
      dat_linear[0][i][j] -= val * (x1 - 1);
      dat_linear[1][i][j] -= val * (y1 - 1);
      dat_quadratic[i][j] += val;
    }
    for (int i = y1; i < height; i += i & -i) for (int j = x2 + 1; j < width; j += j & -j) {
      dat_const[i][j] -= val * (y1 - 1) * x2;
      dat_linear[0][i][j] += val * x2;
      dat_linear[1][i][j] += val * (y1 - 1);
      dat_quadratic[i][j] -= val;
    }
    for (int i = y2 + 1; i < height; i += i & -i) for (int j = x1; j < width; j += j & -j) {
      dat_const[i][j] -= val * y2 * (x1 - 1);
      dat_linear[0][i][j] += val * (x1 - 1);
      dat_linear[1][i][j] += val * y2;
      dat_quadratic[i][j] -= val;
    }
    for (int i = y2 + 1; i < height; i += i & -i) for (int j = x2 + 1; j < width; j += j & -j) {
      dat_const[i][j] += val * y2 * x2;
      dat_linear[0][i][j] -= val * x2;
      dat_linear[1][i][j] -= val * y2;
      dat_quadratic[i][j] += val;
    }
  }

  Abelian sum(int y, int x) const {
    ++y; ++x;
    Abelian quad = ID, cons = ID, line[2] = {ID, ID};
    for (int i = y; i > 0; i -= i & -i) for (int j = x; j > 0; j -= j & -j) {
      quad += dat_quadratic[i][j];
      line[0] += dat_linear[0][i][j];
      line[1] += dat_linear[1][i][j];
      cons += dat_const[i][j];
    }
    return quad * y * x + line[0] * y + line[1] * x + cons;
  }

  Abelian sum(int y1, int x1, int y2, int x2) const {
    return y1 <= y2 && x1 <= x2 ? sum(y2, x2) - sum(y2, x1 - 1) - sum(y1 - 1, x2) + sum(y1 - 1, x1 - 1) : ID;
  }

private:
  int height, width;
  const Abelian ID;
  std::vector<std::vector<Abelian>> dat_const, dat_quadratic;
  std::vector<std::vector<Abelian>> dat_linear[2];
};

int main() {
  int h, w, n, m; cin >> h >> w >> n >> m;
  vector<int> t(n), u(n), l(n), r(n), a(n);
  REP(i, n) cin >> t[i] >> u[i] >> l[i] >> r[i] >> a[i], --t[i], --u[i], --l[i], --r[i];
  BIT2DRangeAdd<ll> bit(h, w);
  while (m--) {
    int x, y, b, c; cin >> x >> y >> b >> c; --x; --y;
    bit.add(max(x - b, 0), max(y - b, 0), min(x + b, h - 1), min(y + b, w - 1), c);
  }
  int ans = 0;
  REP(i, n) ans += bit.sum(t[i], l[i], u[i], r[i]) < a[i];
  cout << ans << '\n';
  return 0;
}
0