题目传送：https://www.luogu.org/problemnew/show/P2468
这是一个二合一的题目，前50% \(n!=1\)的分数中，我们考虑用动态规划来做。
设\(sum[i][j][k]\)表示大于等于k的\([1,1]\)到\([i,j]\)范围内的数的总和，然后用\(num[i][j][k]\)表示数的个数。
之后就可以用二维前缀和来维护了，具体见代码subtask1
后50%的题目因为数据范围很大，所以我们只能考虑使用log复杂度的数据结构，因为还涉及到区间排名的问题，所以自然查询的时候要换成主席树。
不需要离散化。在查询的时候因为是查找大于等于当前查找值的数的个数，所以我们要从右子树开始查找（显然是要从大到小），之后如果右子树的总和大于等于当前查找值，那么就进入右子树，如果不是，那么加上右子树的数的个数然后查找左子树。之后要注意的是如果到了叶子节点，我们要求的个数等于还需要的总和除以该节点值的向上取整。
具体见代码subtask2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 1000*10+500000*15+10
int n,m,kkk;
namespace subtask1
{
int maxx=0;
int sum[210][210][1010],num[210][210][1010],a[1010][1010];
inline int max(int x,int y){return x>y?x:y;}
inline int getsum(int x1,int x2,int y1,int y2,int k){return sum[x2][y2][k]-sum[x1-1][y2][k]-sum[x2][y1-1][k]+sum[x1-1][y1-1][k];}
inline int getnum(int x1,int x2,int y1,int y2,int k){return num[x2][y2][k]-num[x1-1][y2][k]-num[x2][y1-1][k]+num[x1-1][y1-1][k];}
inline void init()
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]),maxx=max(maxx,a[i][j]);
for(int k=0;k<=maxx;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
sum[i][j][k]=sum[i][j-1][k]+sum[i-1][j][k]-sum[i-1][j-1][k]+(a[i][j]>=k?a[i][j]:0);
num[i][j][k]=num[i][j-1][k]+num[i-1][j][k]-num[i-1][j-1][k]+(a[i][j]>=k?1:0);
}
}
inline void work(int x1,int x2,int y1,int y2,int h)
{
if(getsum(x1,x2,y1,y2,0)<h)
{
printf("Poor QLW\n");
return;
}
int l=0,r=maxx+1,ans=-0x3f;
while(l<r)
{
int mid=(l+r)>>1;
if(getsum(x1,x2,y1,y2,mid)>=h) ans=mid,l=mid+1;
else r=mid;
}
if(ans==-0x3f)
{
printf("Poor QLW\n");
return;
}
printf("%d\n",getnum(x1,x2,y1,y2,ans)-(getsum(x1,x2,y1,y2,ans)-h)/ans);
}
}
namespace subtask2
{
int cnt;
int sum[MAXN],siz[MAXN],ls[MAXN],rs[MAXN],rt[MAXN],a[500010];
inline int max(int x,int y){return x>y?x:y;}
inline void build(int &x,int l,int r)
{
x=++cnt;
if(l==r) return;
int mid=(l+r)>>1;
build(ls[x],l,mid);
build(rs[x],mid+1,r);
}
inline int modify(int x,int l,int r,int p)
{
int to=++cnt;
ls[to]=ls[x],rs[to]=rs[x],siz[to]=siz[x]+1,sum[to]=sum[x]+p;
if(l==r) return to;
int mid=(l+r)>>1;
if(mid<p) rs[to]=modify(rs[to],mid+1,r,p);
else ls[to]=modify(ls[to],l,mid,p);
return to;
}
inline int query(int u,int v,int l,int r,int k)
{
if(l==r) return (k+l-1)/l;
int mid=(l+r)>>1;
int x=sum[rs[v]]-sum[rs[u]];
if(k<=x) return query(rs[u],rs[v],mid+1,r,k);
else return siz[rs[v]]-siz[rs[u]]+query(ls[u],ls[v],l,mid,k-x);
}
inline void work()
{
int maxx=0;
for(int i=1;i<=m;i++) scanf("%d",&a[i]),maxx=max(maxx,a[i]);
build(rt[0],1,m);
for(int i=1;i<=m;i++) rt[i]=modify(rt[i-1],1,maxx,a[i]);
for(int i=1;i<=kkk;i++)
{
int cur1,cur2,u,v,h;
scanf("%d",&cur1),scanf("%d",&u);
scanf("%d",&cur2),scanf("%d",&v);
scanf("%d",&h);
if(sum[rt[v]]-sum[rt[u-1]]<h)
{
printf("Poor QLW\n");
continue;
}
printf("%d\n",query(rt[u-1],rt[v],1,maxx,h));
}
}
}
using namespace std;
using namespace subtask1;
using namespace subtask2;
int main()
{
scanf("%d%d%d",&n,&m,&kkk);
if(n!=1)
{
subtask1::init();
for(int i=1;i<=kkk;i++)
{
int xx1,xx2,yy1,yy2;
int val;
scanf("%d%d%d%d%d",&xx1,&yy1,&xx2,&yy2,&val);
subtask1::work(xx1,xx2,yy1,yy2,val);
}
}
else
subtask2::work();
return 0;
}