#pragma GCC optimize("Ofast")
#pragma GCC target("avx,avx2,fma")
#include <iostream>
#include <algorithm>
#include <climits>
#include <cmath>
#include <map>
#include <set>
#include <iomanip>
#include <vector>
#include <fstream>
using namespace std;
#define vll vector<ll>
#define sll set<ll>
#define vstr vector<string>
#define ll long long
#define ld long double
#define supra main
#define pb push_back
#define add insert
#define rall(x) rbegin(x),rend(x)
#define all(x) (x).begin(),(x).end()
#define I ios_base::sync_with_stdio(false);
#define Hear cin.tie(NULL);
#define Shots cout.tie(NULL);
#define Ratatatata
ll powmod(ll a,ll b,ll modulo)
{
if(b==0){
return 1;
}
ll temp=powmod(a,b/2,modulo);
if(b%2==0){
return (temp*temp)%modulo;
}
else{
return (a*((temp*temp)%modulo))%modulo;
}
}
bool Prime(ll n){
if (n <= 1)
return false;
for (ll i = 2; i <= sqrt(n); i++)
if (n % i == 0)
return false;
return true;
}
#define int long long
struct SegmentTree
{
int sum=0,lazy=0;
SegmentTree* next[2];
int s,e;
SegmentTree(int l,int r)
{
s=l;
e=r;
if(s!=e)
{
int mid=(s+e)/2;
next[0]=new SegmentTree(s,mid);
next[1]=new SegmentTree(mid+1,e);
}
}
void push()
{
if(!lazy)
return;
sum+=(e-s+1)*lazy;
if(s!=e)
{
next[0]->lazy+=lazy;
next[1]->lazy+=lazy;
}
lazy=0;
}
void Update(int& l,int& r,int& val)
{
push();
if(e<l or r<s)
return;
if(l<=s and e<=r)
{
lazy+=val;
push();
return;
}
next[0]->Update(l,r,val);
next[1]->Update(l,r,val);
sum=next[0]->sum+next[1]->sum;
}
int get(int& l,int& r)
{
push();
if(e<l or r<s)
{
return 0;
}
if(l<=s and e<=r)
{
return sum;
}
int ans=next[0]->get(l,r)+next[1]->get(l,r);
sum=next[0]->sum+next[1]->sum;
return ans;
}
};
SegmentTree* dif;
ll n,m;
void Call_Me_When_You_Want(int& l,int& r,int& cur,int& add)
{
if(r<l)
{
return;
}
// dif[l]+=cur;
// dif[l+1]+=add;
// dif[r]+=add;
// dif[r+1]-=(cur+(add*(r-l)))
// a[l]+=cur;
// a[l+1]+=cur+add;
// | += |
// a[r]+=cur+(add*(r-l))
// cout<<"Adding "<<cur<<" To range "<<l<<' '<<l<<endl;
dif->Update(l,l,cur);
// for(int i=1;i<=m;i++)
// {
// cout<<dif->get(i,i)<<' ';;
// }
// cout<<endl;
// cout<<"Adding "<<add<<" To range "<<l+1<<' '<<r<<endl;
l++;
dif->Update(l,r,add);
l--;
// for(int i=1;i<=m;i++)
// {
// cout<<dif->get(i,i)<<' ';;
// }
// cout<<endl;
// cout<<"Adding "<<-1*(cur+(r-l)*add)<<" To range "<<r+1<<' '<<r+1<<endl;
r++;
cur+=(r-l-1)*add;
cur*=-1;
dif->Update(r,r,cur);
cur*=-1;
cur-=(r-l-1)*add;
r--;
// for(int i=1;i<=m;i++)
// {
// cout<<dif->get(i,i)<<' ';;
// }
// cout<<endl;
}
void solve()
{
cin>>m>>n;
ll p;
cin>>p;
dif=new SegmentTree(1,m+1);
for(int asp=0;asp<p;asp++)
{
ll i,j,final,add;
cin>>j>>i>>final>>add;
ll l=max(1ll,j-(final/add));
ll cur=final-(add*(j-l));
j--;
Call_Me_When_You_Want(l,j,cur,add);
j++;
ll r=min(m,j+(final/add));
add*=-1;
Call_Me_When_You_Want(j,r,final,add);
}
vll pre={0};
int one=1;
for(int i=1;i<=m;i++)
{
pre.pb(pre.back()+dif->get(one,i));
}
ll q;
cin>>q;
while(q--)
{
ll y1,x1,y2,x2;
cin>>y1>>x1>>y2>>x2;
ll sum=0,cell=(y2-y1+1)*(x2-x1+1);
sum=pre[y2]-pre[y1-1];
ld cur=((ld)sum/(ld)cell);
cur+=0.5;
cout<<(ll)cur<<'\n';
}
}
signed supra(){
I Hear Shots Ratatatata
// cout<<setprecision(1000);
// freopen("input.txt", "r", stdin);
// freopen("output.txt", "w", stdout);
ll tqwertyuiop=1;
// cin>>tqwertyuiop;
for(int tp=1;tp<=tqwertyuiop;tp++)
{
// cout<<"Case #"<<tp<<": ";
solve();
}
return 0;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
522 ms |
346784 KB |
Output is correct |
2 |
Correct |
45 ms |
2640 KB |
Output is correct |
3 |
Correct |
43 ms |
2388 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
503 ms |
347816 KB |
Output is correct |
2 |
Correct |
46 ms |
2644 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
1 ms |
1116 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
0 ms |
344 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
596 ms |
347104 KB |
Output is correct |
2 |
Correct |
56 ms |
3152 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
276 ms |
135864 KB |
Output is correct |
2 |
Correct |
45 ms |
2640 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
55 ms |
2128 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
159 ms |
69296 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
1051 ms |
317752 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
1049 ms |
318396 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
286 ms |
2844 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
274 ms |
2536 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
161 ms |
2128 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
279 ms |
2972 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |