#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
#define ld long double
#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,tmp[1010][1010];
vector<int>a[N];
vector<ii>f[N],sum[N];
void solve(int check)
{
//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[(check==0 ? i : n-i+1)][(check==0 ? j : m-j+1)]+=(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[(check==0 ? i : n-i+1)][(check==0 ? j : m-j+1)]+=(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);
}
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(0);
for(int i=1;i<=q;i++)
{
tv[i].x=n-tv[i].x+1;
tv[i].y=m-tv[i].y+1;
}
solve(1);
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<<(ll)(round((ld)1.0*sum/((ld)1.0*(xe-xs+1))/((ld)1.0*(ye-ys+1))))<<"\n";
}
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1117 ms |
567732 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1121 ms |
567740 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
467 ms |
274952 KB |
Output is correct |
2 |
Correct |
136 ms |
178888 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
668 ms |
284288 KB |
Output is correct |
2 |
Correct |
142 ms |
178832 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1104 ms |
567744 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
740 ms |
335384 KB |
Output is correct |
2 |
Correct |
144 ms |
178800 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
425 ms |
277124 KB |
Output is correct |
2 |
Correct |
146 ms |
179096 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
557 ms |
280592 KB |
Output is correct |
2 |
Correct |
138 ms |
178768 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1122 ms |
574080 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1114 ms |
573956 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
691 ms |
283256 KB |
Output is correct |
2 |
Correct |
708 ms |
283088 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
670 ms |
283108 KB |
Output is correct |
2 |
Correct |
536 ms |
282916 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
802 ms |
285836 KB |
Output is correct |
2 |
Correct |
653 ms |
284224 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
658 ms |
283620 KB |
Output is correct |
2 |
Correct |
673 ms |
478568 KB |
Output is correct |