oj.uz

Submission #276534

# Submission time Handle Problem Language Result Execution time Memory
276534 2020-08-20T13:34:00 Z MarcoMeijer Nuclearia (CEOI15_nuclearia) C++14
61 / 100
 1000 ms 243960 KB 
#pragma GCC optimize ("O3")
#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 = 5e6+1e5;
int w, h, n, w1;
ll x[MX], y[MX], a[MX], b[MX];
ll c[MX], d[MX], e[MX], g[MX], v[MX], temp[MX];
bool firstTime = 1;

int f(int i, int j) {return i+j*w1;}

void rotate() {
    // [x,y] => [h-1-y,x]
    RE(i,n) {
        y[i] = h-1-y[i];
        swap(x[i], y[i]);
    }
    int h1=h+1;
    RE(i,w) RE(j,h) temp[h-1-j + i*h1] = v[f(i,j)];
    swap(w,h);
    w1 = h1;
    RE(i,w) RE(j,h) v[f(i,j)] = temp[f(i,j)];
}
inline void addArea(int bx, int by, int ex, int ey, ll val) {
    bx = max(0,bx);
    by = max(0,by);
    ex = min((int)w,ex+1);
    ey = min((int)h,ey+1);
    if(bx > ex) return;
    if(by > ey) return;
    c[f(bx,by)] += val;
    c[f(ex,by)] -= val;
    c[f(bx,ey)] -= val;
    c[f(ex,ey)] += val;
}
inline void add(int x, int y, ll del) {
    if(x >= w || y >= h) return;
    if(x < 0 && y < 0) {
        int add = min(-x, -y);
        x += add;
        y += add;
        c[0] += del*add;
    }
    if(x < 0) {
        g[f(0,y)] += del;
        g[f(0,y-x)] -= del;
        y -= x;
        x = 0;
    }
    if(y < 0) {
        d[f(x,0)] += del;
        d[f(x-y,0)] -= del;
        x -= y;
        y = 0;
    }
    e[f(x,y)] += del;
}
void solve() {
    RE(i,w1*(h+1)) c[i]=d[i]=e[i]=g[i]=0;
    RE(i,n) {
        ll r = (a[i]-1)/b[i];
        add(x[i]-r, y[i]-r, b[i]);
        add(x[i]+r+1, y[i]+r+1, -b[i]);
        addArea(x[i]+r+1,y[i]+r+1,w,h,-ll(r*2+1)*b[i]);
        addArea(x[i]-r,y[i]+r+1,x[i]+r,h,-ll(r*2+2)*b[i]);
        if(firstTime)
            addArea(x[i]-r, y[i]-r, x[i]+r, y[i]+r, (-b[i]*r+a[i])*2ll -ll(r*2+4)*b[i]);
    }
    RE(i,w) RE(j,h) {
        int idx = f(i,j);
        int idxl = f(i-1,j);
        int idxu = f(i,j-1);
        int idxul = f(i-1,j-1);
        d[idx] += (i ? d[idxl] : 0ll);
        g[idx] += (j ? g[idxu] : 0ll);
        e[idx] += (i&&j ? e[idxul] : 0ll);
        c[idx] += d[idx]+e[idx]+g[idx];
        if(i) c[idx] += c[idxl];
        if(j) c[idx] += c[idxu];
        if(i&&j) c[idx] -= c[idxul];
        v[idx] += c[idx];
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
 
    // input
    IN(w,h,n); w1 = w+1;
    RE(i,n) IN(x[i],y[i],a[i],b[i]), x[i]--, y[i]--;

    RE(i,4) {
        solve();
        rotate();
        firstTime = 0;
    }

    RE(i,w) RE(j,h) v[f(i,j)] >>= 1;

    // create sums
    REV(i,0,w) REV(j,0,h) v[f(i+1,j+1)] = v[f(i,j)];
    RE(i,w) v[f(i,0)]=0;
    RE(i,h) v[f(0,i)]=0;
    RE(i,w) RE(j,h) v[f(i+1,j+1)] += v[f(i,j+1)]+v[f(i+1,j)]-v[f(i,j)];

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

Subtask #1
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 290 ms 235256 KB Output is correct
2 Correct 74 ms 2808 KB Output is correct
3 Correct 69 ms 2424 KB Output is correct

Subtask #2
2.0 / 2.0

# Verdict Execution time Memory Grader output
1 Correct 282 ms 235256 KB Output is correct
2 Correct 76 ms 2864 KB Output is correct

Subtask #3
3.0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 903 ms 118172 KB Output is correct
2 Correct 79 ms 2808 KB Output is correct

Subtask #4
2.0 / 2.0

# Verdict Execution time Memory Grader output
1 Correct 694 ms 123796 KB Output is correct
2 Correct 75 ms 2808 KB Output is correct

Subtask #5
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 435 ms 237492 KB Output is correct
2 Correct 80 ms 3192 KB Output is correct

Subtask #6
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 307 ms 96848 KB Output is correct
2 Correct 86 ms 2936 KB Output is correct

Subtask #7
0 / 6.0

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

Subtask #8
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 224 ms 73464 KB Output is correct
2 Correct 77 ms 2932 KB Output is correct

Subtask #9
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 565 ms 243960 KB Output is correct

Subtask #10
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 576 ms 243960 KB Output is correct

Subtask #11
0 / 15.0

# Verdict Execution time Memory Grader output
1 Execution timed out 1056 ms 124420 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 1099 ms 124408 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #13
12.0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 683 ms 127384 KB Output is correct
2 Correct 863 ms 126984 KB Output is correct

Subtask #14
0 / 8.0

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