oj.uz

Submission #570281

# Submission time Handle Problem Language Result Execution time Memory
570281 2022-05-29T06:57:38 Z 600Mihnea Nuclearia (CEOI15_nuclearia) C++17
30 / 100
 1000 ms 315628 KB 
#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

signed main() {
  ios::sync_with_stdio(0); cin.tie(0);

 /// freopen ("input.txt", "r", stdin);

  bool Inverse = 0;

  int n, m;
  cin >> n >> m;
  if (n > m) {
    Inverse = 1;
    swap(n, m);
  }
  assert(n <= m);
  vector<vector<ll>> radiation(n + 2, vector<ll> (m + 2, 0)),
                          deri(n + 2, vector<ll> (m + 2, 0)),
                          jmen(n + 2, vector<ll> (m + 2, 0)),
                          difs(n + 2, vector<ll> (m + 2, 0)),
                             v(n + 2, vector<ll> (m + 2, 0));

  function<void(int, int, int, int, int)> proc = [&] (int dist, int r, int c, int a, int b) {
    int stop = min(max(m - c, c - 1), a / b);
    assert(a - dist * b >= 1);
    int L = max(1, c - dist), R = min(m, c + dist);
    jmen[r][L] += (a - dist * b);
    jmen[r][R + 1] -= (a - dist * b);
    if (dist + 1 <= stop) {
      if (dist + 1 <= min(stop, m - c)) {
        L = c + dist + 1;
        R = min(c + stop, m);
        assert(L <= R);
        jmen[r][L] += (a + 1LL * b * c);
        jmen[r][R + 1] -= (a + 1LL * b * c);
        difs[r][L] -= 1LL * b * L;
        deri[r][L + 1] -= b;
        deri[r][R + 1] += b;
        difs[r][R + 1] += 1LL * b * R;
      }
      if (dist + 1 <= min(stop, c - 1)) {
        L = max(c - stop, 1);
        R = c - dist - 1;
        assert(L <= R);
        jmen[r][L] += (a - 1LL * b * c);
        jmen[r][R + 1] -= (a - 1LL * b * c);
        difs[r][L] += 1LL * b * L;
        deri[r][L + 1] += b;
        deri[r][R + 1] -= b;
        difs[r][R + 1] -= 1LL * b * R;
      }
    }
  };

  int ops;
  cin >> ops;
  for (int i = 1; i <= ops; i++) {
    int r, c, a, b;
    cin >> r >> c >> a >> b;
    if (Inverse) {
      swap(r, c);
    }
    for (int r2 = 1; r2 <= n; r2++) {
      for (int c2 = 1; c2 <= m; c2++) {
        int dist = max(abs(r - r2), abs(c - c2));
        radiation[r2][c2] += max(0LL, a - 1LL * b * dist);
      }
    }

    for (int dist = 0; (r + dist <= n || r - dist >= 1) && a - dist * b >= 1; dist++) {
      if (r + dist <= n) {
        proc(dist, r + dist, c, a, b);
      }
      if (r - dist >= 1 && dist >= 1) {
        proc(dist, r - dist, c, a, b);
      }
    }
  }
  for (int r = 1; r <= n; r++) {
    for (int c = 1; c <= m; c++) {
      deri[r][c] += deri[r][c - 1];
      difs[r][c] += deri[r][c];
      difs[r][c] += difs[r][c - 1];
      jmen[r][c] += jmen[r][c - 1];
      v[r][c] += jmen[r][c];
      v[r][c] += difs[r][c];
    }
  }
  assert(v == radiation);

  {
    /// transform the radiation in the form of prefix rectangle sum

    for (int i = 1; i <= n; i++) {
      ll cur = 0;
      for (int j = 1; j <= m; j++) {
        cur += radiation[i][j];
        radiation[i][j] = radiation[i - 1][j] + cur;
      }
    }
  }

  int q;
  cin >> q;
  while (q--) {
    int r1, c1, r2, c2;
    cin >> r1 >> c1 >> r2 >> c2;
    if (Inverse) {
      swap(r1, c1);
      swap(r2, c2);
    }
    ll total = radiation[r2][c2] - radiation[r1 - 1][c2] - radiation[r2][c1 - 1] + radiation[r1 - 1][c1 - 1];
    ll area = (r2 - r1 + 1) * (c2 - c1 + 1);

    cout << total / area + (total % area >= area - total % area) << "\n";
  }
}

Subtask #1
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 365 ms 313304 KB Output is correct
2 Correct 57 ms 2636 KB Output is correct
3 Correct 64 ms 2252 KB Output is correct

Subtask #2
2.0 / 2.0

# Verdict Execution time Memory Grader output
1 Correct 343 ms 313424 KB Output is correct
2 Correct 63 ms 2664 KB Output is correct

Subtask #3
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 239 ms 98644 KB Output is correct
2 Correct 54 ms 2636 KB Output is correct

Subtask #4
2.0 / 2.0

# Verdict Execution time Memory Grader output
1 Correct 229 ms 109048 KB Output is correct
2 Correct 71 ms 2636 KB Output is correct

Subtask #5
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 306 ms 315628 KB Output is correct
2 Correct 63 ms 3048 KB Output is correct

Subtask #6
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 351 ms 127964 KB Output is correct
2 Correct 58 ms 2704 KB Output is correct

Subtask #7
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 369 ms 101136 KB Output is correct
2 Correct 61 ms 2900 KB Output is correct

Subtask #8
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 309 ms 84872 KB Output is correct
2 Correct 57 ms 2636 KB Output is correct

Subtask #9
0 / 15.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1110 ms 313300 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #10
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1104 ms 313440 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #11
0 / 15.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1092 ms 98644 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #12
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1097 ms 98544 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #13
0 / 12.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1094 ms 100320 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #14
0 / 8.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1101 ms 98684 KB Time limit exceeded
2 Halted 0 ms 0 KB -