Submission #813446

# Submission time Handle Problem Language Result Execution time Memory
813446 2023-08-07T18:01:57 Z JakobZorz Nuclearia (CEOI15_nuclearia) C++14
55 / 100
1000 ms 22304 KB
#include <iostream>
#include <vector>
#include <queue>
#include <algorithm>
#include <set>
#include <stack>
#include <limits.h>
#include <math.h>
#include <iomanip>
#include <bitset>
#include <unordered_map>
#include <unordered_set>
#include <map>
#include <cstring>
#include <sstream>
 
#pragma GCC target("popcnt")
 
typedef long long ll;
typedef long double ld;
using namespace std;
const int MOD=1e9+7;
typedef pair<ll,ll>point;
//#define x first
//#define y second

int w,h;

class Matrix{
    vector<ll>map_arr;
public:
    void init_map(){
        map_arr.resize(w*h);
    }

    ll get2(int x,int y){
        if(x<0||y<0)
            return 0;
        return map_arr[y*w+x];
    }

    ll&get(int x,int y){
        return map_arr[y*w+x];
    }

    void add_bounded(int x,int y,ll val){
        if(x<0)
            x=0;
        if(y<0)
            y=0;
        if(x>=w)
            return;
        if(y>=h)
            return;
        get(x,y)+=val;
    }
};
 
Matrix mat;

int main(){
    ios::sync_with_stdio(false);
    cout.tie(NULL);
    cin.tie(NULL);
    
    cin>>w>>h;
    mat.init_map();
    int n;
    cin>>n;
    while(n--){
        int x1,y1;
        ll center,step;
        cin>>x1>>y1>>center>>step;
        x1--;y1--;
        
        ll last_step=(center-1)%step+1;
        int num_steps=(int)(center+step-1)/step;
        
        if(x1>y1){
            // hits ceiling
            for(int i=0;i<num_steps-1&&i<y1;i++)
                mat.get(x1-i,y1-i)+=step;
            
            for(int i=y1;i<num_steps-1&&i<x1;i++)
                mat.get(x1-i,0)+=step;
            
            if(x1<num_steps-1)
                mat.get(0,0)+=step*(num_steps-1-x1);
            
        }else{
            // hits left wall
            for(int i=0;i<num_steps-1&&i<x1;i++)
                mat.get(x1-i,y1-i)+=step;
            
            for(int i=x1;i<num_steps-1&&i<y1;i++)
                mat.get(0,y1-i)+=step;
            
            if(y1<num_steps-1)
                mat.get(0,0)+=step*(num_steps-1-y1);
        }
        
        if(x1+y1<w){
            // hits ceiling
            for(int i=0;i<num_steps-1&&i<y1;i++)
                mat.get(x1+i+1,y1-i)-=step;
            
            for(int i=y1;i<num_steps-1&&x1+i+1<w;i++)
                mat.get(x1+i+1,0)-=step;
            
        }else{
            for(int i=0;i<num_steps-1&&x1+i+1<w;i++)
                mat.get(x1+i+1,y1-i)-=step;
        }
        
        if(x1+y1<h){
            // hits left wall
            for(int i=0;i<num_steps-1&&i<x1;i++)
                mat.get(x1-i,y1+i+1)-=step;
            
            for(int i=x1;i<num_steps-1&&y1+i+1<h;i++)
                mat.get(0,y1+i+1)-=step;
            
        }else{
            for(int i=0;i<num_steps-1&&y1+i+1<h;i++)
                mat.get(x1-i,y1+i+1)-=step;
        }
        
        if(w-x1<h-y1){
            for(int i=0;i<num_steps-1&&x1+i+1<w;i++)
                mat.get(x1+i+1,y1+i+1)+=step;
        }else{
            for(int i=0;i<num_steps-1&&y1+i+1<h;i++)
                mat.get(x1+i+1,y1+i+1)+=step;
        }
        
        mat.add_bounded(x1-num_steps+1, y1-num_steps+1, last_step);
        mat.add_bounded(x1-num_steps+1, y1+num_steps,   -last_step);
        mat.add_bounded(x1+num_steps,   y1-num_steps+1, -last_step);
        mat.add_bounded(x1+num_steps,   y1+num_steps,   last_step);
        
        /*for(int x=0;x<w;x++)
            for(int y=0;y<h;y++){
                int dist=max(abs(x1-x),abs(y1-y));
                ll val=max(0LL,center-dist*step);
                get(x,y)+=val;
            }*/
    }
    
    for(int x=0;x<w;x++)
        for(int y=0;y<h;y++)
            mat.get(x,y)+=mat.get2(x-1,y)+mat.get2(x,y-1)-mat.get2(x-1,y-1);
    
    /*for(int y=0;y<h;y++){
        for(int x=0;x<w;x++){
            cout<<mat.get(x,y)<<" ";
        }
        cout<<"\n";
    }*/
    
    for(int x=0;x<w;x++)
        for(int y=0;y<h;y++)
            mat.get(x,y)+=mat.get2(x-1,y)+mat.get2(x,y-1)-mat.get2(x-1,y-1);
    
    int q;
    cin>>q;
    while(q--){
        int x1,y1,x2,y2;
        cin>>x1>>y1>>x2>>y2;
        x1--;y1--;
        
        ll sum=mat.get2(x2-1,y2-1)+mat.get2(x1-1,y1-1)-mat.get2(x1-1,y2-1)-mat.get2(x2-1,y1-1);
        
        ll area=(x2-x1)*(y2-y1);
        ll res=sum/area;
        if((sum%area)*2>=area)
            res++;
        cout<<res<<"\n";
    }

    return 0;
}

/*

4 3
2
1 1 7 3
3 2 4 2
4
1 2 2 3
1 1 4 3
4 2 4 2
1 3 4 3

*/
# Verdict Execution time Memory Grader output
1 Correct 98 ms 19796 KB Output is correct
2 Correct 46 ms 2692 KB Output is correct
3 Correct 44 ms 2284 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 77 ms 19796 KB Output is correct
2 Correct 47 ms 2600 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 23 ms 19796 KB Output is correct
2 Correct 46 ms 2636 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 63 ms 19796 KB Output is correct
2 Correct 47 ms 2620 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 112 ms 22004 KB Output is correct
2 Correct 50 ms 2896 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 126 ms 10500 KB Output is correct
2 Correct 47 ms 2672 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 142 ms 22304 KB Output is correct
2 Correct 47 ms 2856 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 86 ms 10444 KB Output is correct
2 Correct 46 ms 2708 KB Output is correct
# Verdict Execution time Memory Grader output
1 Execution timed out 1072 ms 19880 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1073 ms 19876 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 516 ms 21800 KB Output is correct
2 Correct 852 ms 22060 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 765 ms 21964 KB Output is correct
2 Correct 504 ms 21852 KB Output is correct
# Verdict Execution time Memory Grader output
1 Execution timed out 1069 ms 19936 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Execution timed out 1061 ms 20388 KB Time limit exceeded
2 Halted 0 ms 0 KB -