Submission #565238

# Submission time Handle Problem Language Result Execution time Memory
565238 2022-05-20T14:12:35 Z piOOE Nuclearia (CEOI15_nuclearia) C++17
25 / 100
1000 ms 826944 KB
#include <bits/stdc++.h>

using namespace std;

#define all(x) begin(x), end(x)
#define sz(x) ((int)size(x))
#define trace(x) cout << #x << ": " << (x) << endl;

typedef long long ll;

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int n, m, k;
    cin >> n >> m >> k;
    vector<int> x(k + 1), y(k + 1);
    vector<ll> a(k + 1), b(k + 1);
    for (int i = 1; i <= k; ++i) {
        cin >> x[i] >> y[i] >> a[i] >> b[i];
    }
    vector<vector<ll>>
            v(n + 2, vector<ll>(m + 2)), pref(n + 2, vector<ll>(m + 2)), aa(n + 2, vector<ll>(m + 2)),
            bb(n + 2, vector<ll>(m + 2)), cc(n + 2, vector<ll>(m + 2)), addi(n + 2, vector<ll>(m + 2));

    auto addR = [&](int i, int j, int h, ll x) {
        bool ok1 = 0 <= i && i <= n + 1;
        bool ok2 = 0 <= j && j <= m + 1;
        if (!ok1) {
            while (1) {

            }
        }
        if (!ok2) {
            while (1) {

            }
        }
        bb[i][j] += x;
        int ni = min(n + 1, i + h);
        int nj = min(m + 1, j + h);
        bb[ni][nj] -= x;
    };

    auto addL = [&](int i, int j, int h, ll x) {
        bool ok1 = 0 <= i && i <= n + 1;
        bool ok2 = 0 <= j && j <= m + 1;
        if (!ok1) {
            while (1) {

            }
        }
        if (!ok2) {
            while (1) {

            }
        }
        cc[i][j] += x;
        int ni = min(n + 1, i + h);
        int nj = max(0, j - h);
        cc[ni][nj] -= x;
    };

    for (int i = 1; i <= k; ++i) {
        int h = (a[i] - 1) / b[i];
        ll val = a[i] - h * b[i];
        if (x[i] < h) {
            while (1) {}
        }
        if (y[i] < h) {
            while (1) {}
        }
        if (x[i] + h > n) {
            while (1) {

            }
        }
        if (y[i] + h > m) {
            while (1) {

            }
        }
        aa[x[i] - h][y[i] - h] += val;
        aa[x[i] + h + 1][y[i] + h + 1] += val;
        aa[x[i] - h][y[i] + h + 1] -= val;
        aa[x[i] + h + 1][y[i] - h] -= val;
        if (h) {
            addR(x[i] - h + 1, y[i] - h + 1, h, b[i]);
            addR(x[i] + 1, y[i] + 1, h, b[i]);
            addL(x[i] - h + 1, y[i] + h, h, -b[i]);
            addL(x[i] + 1, y[i], h, -b[i]);
        }
    }
//    trace("bb")
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            bb[i][j] = bb[i][j] + bb[i - 1][j - 1];
//            cout << bb[i][j] << " ";
        }
//        cout << endl;
    }
//    trace("cc")
    for (int j = m; j > 0; --j) {
        for (int i = 1; i <= n; ++i) {
            cc[i][j] += cc[i - 1][j + 1];
        }
    }
//    for (int i = 1; i <= n; ++i) {
//        for (int j = 1; j <= m; ++j) {
//            cout << cc[i][j] << " ";
//        }
//        cout << endl;
//    }
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            v[i][j] += bb[i][j] + cc[i][j] + aa[i][j];
            v[i][j] += v[i - 1][j];
            v[i][j] += v[i][j - 1];
            v[i][j] -= v[i - 1][j - 1];
//            cout << v[i][j] << " ";
        }
//        cout << endl;
    }
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            pref[i][j] = v[i][j];
            pref[i][j] += pref[i - 1][j];
            pref[i][j] += pref[i][j - 1];
            pref[i][j] -= pref[i - 1][j - 1];
        }
    }
    auto get = [&pref](int x1, int y1, int x2, int y2) {
        ll ans = pref[x2][y2];
        if (x1) ans -= pref[x1 - 1][y2];
        if (y1) ans -= pref[x2][y1 - 1];
        if (x1 && y1) ans += pref[x1 - 1][y1 - 1];
        return ans;
    };
    int q;
    cin >> q;
    for (int i = 1; i <= q; ++i) {
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        ll sum = get(x1, y1, x2, y2);
//        trace(sum)
        ll S = (x2 - x1 + 1) * (ll) (y2 - y1 + 1);
        assert(S);
        cout << (sum * 2 / S + 1) / 2 << '\n';
    }
    return 0;
}
# Verdict Execution time Memory Grader output
1 Execution timed out 1151 ms 822212 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1139 ms 822184 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1100 ms 119348 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1074 ms 130536 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1128 ms 822148 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1110 ms 329012 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1101 ms 118276 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1098 ms 211568 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1124 ms 826944 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1134 ms 826904 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 269 ms 126108 KB Output is correct
2 Correct 306 ms 132616 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 290 ms 125492 KB Output is correct
2 Correct 270 ms 132516 KB Output is correct
# Verdict Execution time Memory Grader output
1 Execution timed out 1101 ms 124928 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1089 ms 124108 KB Time limit exceeded
2 Halted 0 ms 0 KB -