二合一题....
前面 $50$ 分：
考虑取书显然优先取厚的，所以答案满足单调性
发现 $P_{i,j}$ 不大，所以考虑二分最小厚度 $mid$，把大于等于 $mid$ 的书取走
维护 $cnt[i][j][k]$ 表示位置 $i,j$ 为右下角一直到 $1,1$ 的矩形内厚度大于等于 $k$ 的书的数量
维护 $sum[i][j][k]$ 表示位置 $i,j$ 为右下角一直到 $1,1$ 的矩形内厚度大于等于 $k$ 的书的总厚度
然后就可以愉快地二分答案了，注意最后要特判一下厚度为 $ans$ 的书不完全取完的情况！
后面 $50$ 分：
因为只有一行，考虑用主席树维护每个区间的所有值，然后处理询问在主席树上走就好了
同样注意特判最后一种厚度的书不完全取完的情况
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
typedef long long ll;
inline int read()
{
int x=,f=; char ch=getchar();
while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
return x*f;
}
const int N=,M=5e5+,INF=1e9+;
int n,m,Q; int cnt[N][N][],sum[N][N][];//sum是总和，cnt是数量
int xa,ya,xb,yb,H,mx;
inline int clac(int p) { return sum[xb][yb][p]-sum[xa-][yb][p]-sum[xb][ya-][p]+sum[xa-][ya-][p]; }//求一个矩形的sum
inline int CNT(int p) { return cnt[xb][yb][p]-cnt[xa-][yb][p]-cnt[xb][ya-][p]+cnt[xa-][ya-][p]; }//求一个矩形的cnt
inline int query()//二分答案
{
int res=-,l=,r=mx+,mid=l+r>>;
while(l<=r)
{
mid=l+r>>;
if(clac(mid)>=H) l=mid+,res=mid;
else r=mid-;
}
return res;
}
int bk[N][N];
void solve1()//前面50分
{
for(int i=;i<=n;i++)
for(int j=;j<=m;j++) bk[i][j]=read(),mx=max(mx,bk[i][j]);
for(int k=;k<=mx;k++)
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
sum[i][j][k]=sum[i-][j][k]+sum[i][j-][k]-sum[i-][j-][k]+(bk[i][j]>=k ? bk[i][j] : );
cnt[i][j][k]=cnt[i-][j][k]+cnt[i][j-][k]-cnt[i-][j-][k]+(bk[i][j]>=k ? : );
}
while(Q--)
{
xa=read(),ya=read(),xb=read(),yb=read(),H=read();
int t=query();
if(t==-) printf("Poor QLW\n");
else printf("%d\n", CNT(t) - (clac(t)-H)/t );//注意特判厚度为t的书不全取
}
} int t[M<<],sz[M<<],L[M<<],R[M<<],rt[M],tot;
int res,val;
inline void ins(int &o,int l,int r,int pre)
{
o=++tot; t[o]=t[pre]+val; sz[o]=sz[pre]+;
int mid=l+r>>;
if(l==r) return;
if(val<=mid) ins(L[o],l,mid,L[pre]),R[o]=R[pre];
else ins(R[o],mid+,r,R[pre]),L[o]=L[pre];
}
inline void query(int hea,int las,int l,int r)//在主席树上走
{
int mid=l+r>>;
if(l==r)
{
res+= (val/l+(val%l!=)) ;//注意特判
return;
}
if(t[R[las]]-t[R[hea]]<=val)//优先取大的
{
res+=sz[R[las]]-sz[R[hea]]; val-=t[R[las]]-t[R[hea]];
query(L[hea],L[las],l,mid);
}
else query(R[hea],R[las],mid+,r);
}
inline void solve2()//后50分
{
int l,r,a,b;
for(int i=;i<=m;i++) val=read(),ins(rt[i],,,rt[i-]);
while(Q--)
{
a=read(),l=read(),b=read(),r=read(),val=read();
if(t[rt[r]]-t[rt[l-]]<val) { printf("Poor QLW\n"); continue; }//特判
res=; query(rt[l-],rt[r],,);
printf("%d\n",res);
}
} int main()
{
n=read(),m=read(),Q=read();
if(n>) solve1();
else solve2();
return ;
}