#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
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 |
1121 ms |
567712 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1119 ms |
567748 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
464 ms |
274956 KB |
Output is correct |
2 |
Correct |
145 ms |
178920 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
469 ms |
284288 KB |
Output is correct |
2 |
Correct |
148 ms |
178948 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1124 ms |
567796 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
725 ms |
335412 KB |
Output is correct |
2 |
Correct |
160 ms |
178764 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
459 ms |
277056 KB |
Output is correct |
2 |
Correct |
142 ms |
179036 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
567 ms |
280620 KB |
Output is correct |
2 |
Correct |
136 ms |
178856 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1117 ms |
573972 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Execution timed out |
1123 ms |
574040 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
708 ms |
283116 KB |
Output is correct |
2 |
Correct |
696 ms |
283028 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
696 ms |
283112 KB |
Output is correct |
2 |
Correct |
573 ms |
282924 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
622 ms |
285804 KB |
Output is correct |
2 |
Correct |
714 ms |
284188 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
720 ms |
283640 KB |
Output is correct |
2 |
Correct |
687 ms |
478640 KB |
Output is correct |