#include<bits/stdc++.h>
using namespace std;
#define N 2500010
#define ll long long
#define ii pair<int,int>
#define fs first
#define sc second
#define pb push_back
#define iii pair<int,ii>
#define int ll
#pragma GCC optimize("Ofast")
#pragma GCC optimize("inline")
#pragma GCC optimize("omit-frame-pointer")
#pragma GCC optimize("unroll-loops")
struct haha
{
int x,y,a,b;
}tv[200010];
int n,m,q;
vector<int>a[N];
vector<ii>f[N],sum[N],cs[N];
void solve()
{
//duoi
//cout<<"cc\n";
for(int i=0;i<=n+1;i++)
for(int j=0;j<=m+1;j++)
f[i][j]=sum[i][j]={0,0};
for(int i=1;i<=q;i++)
{
f[tv[i].x][tv[i].y].fs+=tv[i].a-tv[i].b*tv[i].x;
f[tv[i].x][tv[i].y].sc+=tv[i].b;
int L=(tv[i].a+tv[i].b-1)/tv[i].b;
if(tv[i].x-L>=1)
{
f[tv[i].x-L][max(1ll,tv[i].y-L)].fs-=(tv[i].a-tv[i].b*tv[i].x);
f[tv[i].x-L][max(1ll,tv[i].y-L)].sc-=tv[i].b;
}
}
for(int i=n;i>=1;i--)
for(int j=m;j>=1;j--)
{
sum[i][j].fs+=f[i][j].fs;
sum[i][j].sc+=f[i][j].sc;
if(i>1)
{
f[i-1][max(1ll,j-1)].fs+=f[i][j].fs;
f[i-1][max(1ll,j-1)].sc+=f[i][j].sc;
}
}
for(int i=0;i<=n+1;i++)
for(int j=0;j<=m+1;j++)
f[i][j]={0,0};
for(int i=1;i<=q;i++)
{
f[tv[i].x][tv[i].y].fs-=(tv[i].a-tv[i].b*tv[i].x);
f[tv[i].x][tv[i].y].sc-=(tv[i].b);
int L=(tv[i].a+tv[i].b-1)/tv[i].b;
if(tv[i].x-L>=1)
{
// cout<<i<<" "<<L<<endl;
f[tv[i].x-L][min(m+1,tv[i].y+L)].fs+=(tv[i].a-tv[i].b*tv[i].x);
f[tv[i].x-L][min(m+1,tv[i].y+L)].sc+=tv[i].b;
}
}
for(int i=n;i>=1;i--)
for(int j=m;j>=1;j--)
{
sum[i][j].fs+=f[i][j].fs;
sum[i][j].sc+=f[i][j].sc;
if(i>1)
{
f[i-1][min(m+1,j+1)].fs+=f[i][j].fs;
f[i-1][min(m+1,j+1)].sc+=f[i][j].sc;
}
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
sum[i][j].fs+=sum[i][j-1].fs;
sum[i][j].sc+=sum[i][j-1].sc;
a[cs[i][j].fs][cs[i][j].sc]+=(sum[i][j].fs+sum[i][j].sc*i);
}
//trai
for(int i=0;i<=n+1;i++)
for(int j=0;j<=m+1;j++)
f[i][j]=sum[i][j]={0,0};
for(int i=1;i<=q;i++)
{
f[tv[i].x][tv[i].y].fs-=(tv[i].a-tv[i].b*tv[i].y);
f[tv[i].x][tv[i].y].sc-=(tv[i].b);
int L=(tv[i].a+tv[i].b-1)/tv[i].b;
// cout<<tv[i].y<<" "<<L<<endl;
if(tv[i].y-L>=1)
{
// cout<<i<<endl;
f[max(0ll,tv[i].x-L)][tv[i].y-L].fs+=(tv[i].a-tv[i].b*tv[i].y);
f[max(0ll,tv[i].x-L)][tv[i].y-L].sc+=tv[i].b;
}
}
for(int j=m;j>=1;j--)
for(int i=n;i>=1;i--)
{
sum[i][j].fs+=f[i][j].fs;
sum[i][j].sc+=f[i][j].sc;
if(j>1)
{
f[max(0ll,i-1)][j-1].fs+=f[i][j].fs;
f[max(0ll,i-1)][j-1].sc+=f[i][j].sc;
}
}
for(int i=0;i<=n+1;i++)
for(int j=0;j<=m+1;j++)
f[i][j]={0,0};
for(int i=1;i<=q;i++)
{
f[tv[i].x][tv[i].y].fs+=(tv[i].a-tv[i].b*tv[i].y);
f[tv[i].x][tv[i].y].sc+=tv[i].b;
int L=(tv[i].a+tv[i].b-1)/tv[i].b;
if(tv[i].y-L>=1)
{
f[min(n,tv[i].x+L)][tv[i].y-L].fs-=(tv[i].a-tv[i].b*tv[i].y);
f[min(n,tv[i].x+L)][tv[i].y-L].sc-=tv[i].b;
}
}
for(int j=m;j>=1;j--)
for(int i=n;i>=1;i--)
{
sum[i][j].fs+=f[i][j].fs;
sum[i][j].sc+=f[i][j].sc;
if(j>1)
{
f[min(n,i+1)][j-1].fs+=f[i][j].fs;
f[min(n,i+1)][j-1].sc+=f[i][j].sc;
}
}
for(int i=n;i>=1;i--)
for(int j=1;j<=m;j++)
{
sum[i][j].fs+=sum[i+1][j].fs;
sum[i][j].sc+=sum[i+1][j].sc;
a[cs[i][j].fs][cs[i][j].sc]+=(sum[i][j].fs+sum[i][j].sc*j);
}
}
int32_t main()
{
// freopen("A.inp","r",stdin);
// freopen("A.out","w",stdout);
ios::sync_with_stdio(0);
cin.tie(NULL);
cout.tie(NULL);
cin>>n>>m>>q;
for(int i=0;i<=n+1;i++)
{
a[i].resize(m+2,0);
f[i].resize(m+2);
sum[i].resize(m+2);
cs[i].resize(m+2);
for(int j=1;j<=m;j++)
cs[i][j]={i,j};
}
for(int i=1;i<=q;i++)
{
cin>>tv[i].x>>tv[i].y>>tv[i].a>>tv[i].b;
a[tv[i].x][tv[i].y]+=tv[i].a;
}
solve();
for(int i=1;i<=q;i++)
{
tv[i].x=n-tv[i].x+1;
tv[i].y=m-tv[i].y+1;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
cs[i][j]={n-i+1,m-j+1};
solve();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
a[i][j]+=a[i-1][j];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
a[i][j]+=a[i][j-1];
cin>>q;
while(q--)
{
int xs,ys,xe,ye;
cin>>xs>>ys>>xe>>ye;
int sum=a[xe][ye]+a[xs-1][ys-1]-a[xe][ys-1]-a[xs-1][ye];
cout<<round(1.0*sum/(1.0*(xe-xs+1)*(ye-ys+1)))<<"\n";
}
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1124 ms |
783040 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1112 ms |
783068 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
561 ms |
373044 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
680 ms |
386112 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1073 ms |
783092 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1012 ms |
460136 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
599 ms |
378452 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
766 ms |
381780 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1059 ms |
793164 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1114 ms |
793100 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
878 ms |
389120 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
874 ms |
388652 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
966 ms |
391736 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
831 ms |
388692 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |