Submission #564391

#	Time	Username	Problem	Language	Result	Execution time	Memory
564391		piOOE	Nuclearia (CEOI15_nuclearia)	C++17	0 / 100	1098 ms	165472 KiB

nuclearia

#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;

const int N = 200000;

int x[N], y[N], a[N], b[N];

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int nn, mm, n;
    cin >> mm >> nn >> n;
    for (int i = 0; i < n; ++i) {
        cin >> x[i] >> y[i] >> a[i] >> b[i];
        swap(x[i], y[i]);
        --x[i], --y[i];
    }
    vector<vector<ll>> v(nn, vector<ll>(mm)), pref(nn, vector<ll>(mm));
    if (nn == 1) {
        vector<ll> added(mm + 1), bb(mm + 1);
        auto add = [&](int l, int r, ll x, ll y) {
            if (l > r) return;
            added[l] += x;
            bb[l + 1] += y;
            bb[r + 1] -= y;
            added[r + 1] -= x + y * (r - l);
        };
        for (int i = 0; i < n; ++i) {
            int u = a[i] / b[i];
            int L = max(0, y[i] - u);
            add(L, y[i], a[i] - b[i] * 1LL * (y[i] - L), b[i]);
            int R = min(mm - 1, y[i] + u);
            add(y[i] + 1, R, a[i] - b[i], -b[i]);
        }
        ll ad = 0;
        ll sum = 0;
        for (int i = 0; i < mm; ++i) {
            sum += added[i];
            ad += bb[i];
            sum += ad;
            v[0][i] = sum;
        }
    } else if (n * nn * 1LL * mm <= 100000000) {
        for (int k = 0; k < n; ++k) {
            for (int i = 0; i < nn; ++i) {
                for (int j = 0; j < mm; ++j) {
                    v[i][j] += max(0ll, a[k] - b[k] * (ll) max(abs(i - x[k]), abs(j - y[k])));
                }
            }
        }
    }
    for (int i = 0; i < nn; ++i) {
        for (int j = 0; j < mm; ++j) {
            pref[i][j] = v[i][j];
            if (i) pref[i][j] += pref[i - 1][j];
            if (j) pref[i][j] += pref[i][j - 1];
            if (i && j) 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 = 0; i < q; ++i) {
        int x1, y1, x2, y2;
        swap(x1, y1);
        swap(x2, y2);
        cin >> x1 >> y1 >> x2 >> y2;
        --x1, --y1, --x2, --y2;
        ll sum = get(x1, y1, x2, y2);
        ll S = (x2 - x1 + 1) * (ll) (y2 - y1 + 1);
        cout << (sum * 2 / S + 1) / 2 << '\n';
    }
    return 0;
}

• Subtask 1: 0 / 3
• Subtask 2: 0 / 2
• Subtask 3: 0 / 3
• Subtask 4: 0 / 2
• Subtask 5: 0 / 6
• Subtask 6: 0 / 4
• Subtask 7: 0 / 6
• Subtask 8: 0 / 4
• Subtask 9: 0 / 15
• Subtask 10: 0 / 10
• Subtask 11: 0 / 15
• Subtask 12: 0 / 10
• Subtask 13: 0 / 12
• Subtask 14: 0 / 8