Submission #274870

#	Time	Username	Problem	Language	Result	Execution time	Memory
274870		MarcoMeijer	Nuclearia (CEOI15_nuclearia)	C++14	30 / 100	901 ms	554360 KiB

nuclearia

#include <bits/stdc++.h>

using namespace std;

#define REP(a,b,c) for(int a=int(b); a<int(c); a++)     
#define REV(a,b,c) for(int a=int(c-1); a>=int(b); a--)     
#define RE(a,b) REP(a,0,b)
#define RE1(a,b) REP(a,1,b+1)
#define REI(a,b,c) REP(a,b,c+1)
#define FOR(a,b) for(auto& a : b)
#define all(a) a.begin(), a.end()
#define pb push_back
#define sz size()
#define fi first
#define se second

typedef long long ll;
typedef pair<int,int> ii;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector<ii> vii;

template<class T> void IN(T& a) {cin>>a;}
template<class T, class... L> void IN(T& a, L&... b) {IN(a); IN(b...);}
template<class T> void OUT(const T& a) {cout<<a;}
template<class T, class... L> void OUT(const T& a, const L&... b) {OUT(a); OUT(b...);}
template<class... T> void OUTL(const T&... a) {OUT(a..., '\n');}

const int MX = 1e5;
ll w, h, n;
ll x[MX], y[MX], a[MX], b[MX];
vector<vll> c, d, V, SUM;

bool inside(ll i, ll j) {return i>=0 && i<=w && j>=0 && j<=h;}
ll getV(ll i, ll X, ll Y) {
    return a[i] - b[i]*max(0ll, max(abs(x[i]-X), abs(y[i]-Y)));
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    // input
    IN(w,h,n);
    RE(i,n) IN(x[i],y[i],a[i],b[i]), x[i]--, y[i]--;

    c.assign(w+1,vll(h+1,0));
    d.assign(w+1,vll(h+1,0));
    V.assign(w+1,vll(h+1,0));
    SUM.assign(w+1,vll(h+1,0));

    RE(i,n) {
        ll r = (a[i]-1)/b[i];
        ll mnj=max(y[i]-r,0ll);
        ll mxj=min(h,y[i]+r);
        REI(j,mnj,mxj) {
            ll cr = abs(j-y[i]);

            // increasing
            ll bg = max(0ll, x[i] - r);
            ll ed = x[i] - cr - 1;
            if(ed >= 0) {
                ll v = getV(i, bg, j);
                d[bg][j] += b[i];
                c[bg][j] += v - b[i]*bg;
                d[ed+1][j] -= b[i];
                c[ed+1][j] -= v - b[i]*bg;
            }

            // same
            bg = max(0ll, x[i] - cr);
            ed = min(w-1, x[i] + cr);
            ll v = getV(i, bg, j);
            c[bg][j] += v;
            c[ed+1][j] -= v;

            // decreasing
            bg = x[i] + cr + 1;
            ed = min(w-1, x[i] + r);
            if(bg < w) {
                ll v = getV(i, bg, j);
                d[bg][j] -= b[i];
                c[bg][j] += v + b[i]*bg;
                d[ed+1][j] += b[i];
                c[ed+1][j] -= v + b[i]*bg;
            }
        }
    }
    REP(i,1,w+1) RE(j,h+1) {
        c[i][j] += c[i-1][j];
        d[i][j] += d[i-1][j];
    }

    RE(i,w) RE(j,h) V[i][j] = c[i][j] + (ll)i*d[i][j];
    RE(i,w) SUM[i][0] = 0;
    RE(i,h) SUM[0][i] = 0;
    RE(i,w) RE(j,h) SUM[i+1][j+1] = SUM[i+1][j] + SUM[i][j+1] - SUM[i][j] + V[i][j];

    // anser queries
    int q;
    IN(q);
    RE(i,q) {
        ll x1, x2, y1, y2;
        IN(x1,y1,x2,y2); x1--; y1--;
        ll sm = SUM[x2][y2] - SUM[x1][y2] - SUM[x2][y1] + SUM[x1][y1];
        ll area = abs(x2-x1)*abs(y2-y1);
        OUTL((sm+area/2)/area);
    }
}

• Subtask 1: 3 / 3
• Subtask 2: 2 / 2
• Subtask 3: 3 / 3
• Subtask 4: 2 / 2
• Subtask 5: 6 / 6
• Subtask 6: 4 / 4
• Subtask 7: 6 / 6
• Subtask 8: 4 / 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