이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
#pragma GCC optimize("O3")
#define int long long
using namespace std;
typedef long long ll;
typedef long double ld;
vector<vector<ll>> v;
int n,m,k,q;
int nradds=0;
void add(int x1,int y1,int x2,int y2,ll val)
{
x1=max(x1,1LL);
y1=max(y1,1LL);
x2=min(x2,n);
y2=min(y2,m);
v[x1][y1]+=val;
v[x2+1][y1]-=val;
v[x1][y2+1]-=val;
v[x2+1][y2+1]+=val;
}
ll getsum(int x1,int y1,int x2,int y2)
{
return v[x2][y2]-v[x1-1][y2]-v[x2][y1-1]+v[x1-1][y1-1];
}
ll put[2][2500005];
ll s[2500005];
signed main()
{
ios_base::sync_with_stdio(false);
cin.tie(0);
cin>>n>>m>>k;
if(m==1)
{
/*for(int z=1;z<=k;z++)
{
ll x,y,a,b;
cin>>x>>y>>a>>b;
assert(y==1);
ll nr=a/b;
put[0][x]+=a;
put[0][x+1]-=a;
ll st=x-nr;
ll dr=x-1;
if(dr>0)
{
st=max(st,1LL);
ll lg=dr-st+2;
put[0][st]+=a-lg*b;
put[0][dr+1]-=a-lg*b;
put[0][dr+1]-=(lg-1)*b;
put[1][st]+=b;
put[1][dr+1]-=b;
}
st=x+1;
dr=x+nr;
if(st<=n)
{
dr=min(dr,n*1LL);
ll lg=dr-st+2;
put[0][st]+=a;
put[0][dr+1]-=a;
put[0][dr+1]-=(lg-1)*(-b);
put[1][st]-=b;
put[1][dr+1]+=b;
}
}
ll suma=0;
for(int i=1;i<=n;i++)
{
suma+=put[0][i];
s[i]=suma;
}
ll si=0;
suma=0;
for(int i=1;i<=n;i++)
{
si+=put[1][i];
suma+=si;
s[i]+=suma;
}*/
for(int z=1;z<=k;z++)
{
ll x,y,a,b;
cin>>x>>y>>a>>b;
s[x]+=a;
ll nr=a-b;
for(int i=x-1;i>=1&&nr>0;i--,nr-=b)
s[i]+=nr;
nr=a-b;
for(int i=x+1;i<=n&&nr>0;i++,nr-=b)
s[i]+=nr;
}
for(int i=1;i<=n;i++)
s[i]+=s[i-1];
cin>>q;
while(q--)
{
int x1,y1,x2,y2;
cin>>x1>>y1>>x2>>y2;
assert(y1==1&&y2==1);
cout<<s[x2]-s[x1-1]<<'\n';
}
return 0;
}
v.resize(n+2);
for(int i=0;i<=n+1;i++)
v[i].resize(m+2);
for(int z=1;z<=k;z++)
{
int x,y,a,b;
cin>>x>>y>>a>>b;
int nr=a;
int x1=x,y1=y,x2=x,y2=y;
int cnt=0;
while(nr>0&&cnt<=max(n,m))
{
add(x1,y1,x2,y2,b);
nr-=b;
cnt++;
x1--;
x2++;
y1--;
y2++;
if(x1<=0&&y1<=0&&x2>n&&y2>m)
break;
}
//cout<<v[1][1]<<' ';
x1++;
x2--;
y1++;
y2--;
add(x1,y1,x2,y2,a);
add(x1,y1,x2,y2,-cnt*b);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
v[i][j]+=v[i-1][j]+v[i][j-1]-v[i-1][j-1];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
v[i][j]=v[i][j]+v[i-1][j]+v[i][j-1]-v[i-1][j-1];
//cout<<s[i][j]<<' ';
}
//assert(nradds<=k*max(n,m));
cin>>q;
while(q--)
{
int x1,y1,x2,y2;
cin>>x1>>y1>>x2>>y2;
ll suma=getsum(x1,y1,x2,y2);
ll cells=(y2-y1+1)*(x2-x1+1);
ll rez=round(ld(suma)/ld(cells));
cout<<rez<<'\n';
}
return 0;
}
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |