[BZOJ 1926] 粟粟的书架

#include <bits/stdc++.h>

using namespace std;

inline int read()
{
    char ch;int num,f=0;
    while(!isdigit(ch=getchar())) f|=(ch==‘-‘);
    num=ch-‘0‘;
    while(isdigit(ch=getchar())) num=num*10+ch-‘0‘;
    return f?-num:num;
}

const int MAXN=205,MAXM=1005;
const int N=5e5+10,M=1e7+10;

struct FunSeg
{
    int ls,rs,sum,cnt;
}seg[M];
int root[N],S[N],cnt;

int n,m,q,dat[MAXN][MAXN],f[MAXN][MAXN][MAXM],g[MAXN][MAXN][MAXM];

int cal_f(int x1,int y1,int x2,int y2,int k)
{
    return f[x2][y2][k]-f[x1-1][y2][k]-f[x2][y1-1][k]+f[x1-1][y1-1][k];
}

int cal_g(int x1,int y1,int x2,int y2,int k)
{
    return g[x2][y2][k]-g[x1-1][y2][k]-g[x2][y1-1][k]+g[x1-1][y1-1][k];
}

void solve1()
{
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) dat[i][j]=read();
        
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            for(int k=1;k<=1000;k++)
                f[i][j][k]=f[i-1][j][k]+f[i][j-1][k]-f[i-1][j-1][k],
                g[i][j][k]=g[i-1][j][k]+g[i][j-1][k]-g[i-1][j-1][k];
            f[i][j][dat[i][j]]+=dat[i][j];g[i][j][dat[i][j]]++;
        }
                
    for(int k=999;k>=1;k--)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                f[i][j][k]+=f[i][j][k+1],g[i][j][k]+=g[i][j][k+1];
    
    while(q--)
    {
        int x1=read(),y1=read(),x2=read(),y2=read(),h=read(),l=1,r=1000;
        
        while(l<=r)
        {
            int mid=(l+r)>>1;
            if(cal_f(x1,y1,x2,y2,mid)>=h) l=mid+1;
            else r=mid-1;
        }
        
        if(!r){puts("Poor QLW");continue;};
        cout << cal_g(x1,y1,x2,y2,r)-(cal_f(x1,y1,x2,y2,r)-h)/r << endl;  //将这个厚度中多余的书减掉
    }
}

void Update(int x,int &y,int val,int l,int r)
{
    y=++cnt;seg[y]=seg[x];  //不用改变的节点沿用上一次已经完全构造完的
    seg[y].sum+=val;seg[y].cnt++;
    if(l==r) return;
    
    int mid=(l+r)>>1;
    if(val<=mid) Update(seg[x].ls,seg[y].ls,val,l,mid);
    else Update(seg[x].rs,seg[y].rs,val,mid+1,r);
}

int Query(int x,int y,int val,int l,int r)
{
    if(l==r) return val/l+(val%l>0);  //对整除情况特殊处理
    
    int mid=(l+r)>>1;
    int t=seg[seg[y].rs].sum-seg[seg[x].rs].sum;
    if(t>=val) return Query(seg[x].rs,seg[y].rs,val,mid+1,r);  //厚度范围是(mid+1,r)时足够
    else return Query(seg[x].ls,seg[y].ls,val-t,l,mid)+seg[seg[y].rs].cnt-seg[seg[x].rs].cnt;  //不够时
}

void solve2()
{
    for(int i=1;i<=m;i++)
        S[i]=read(),Update(root[i-1],root[i],S[i],1,1000),S[i]+=S[i-1];   //以1~i-1的前缀树为模板建树
        
    while(q--)
    {
        int x1=read(),y1=read(),x2=read(),y2=read(),h=read();
        if(S[y2]-S[y1-1]<h) puts("Poor QLW");
        else cout << Query(root[y1-1],root[y2],h,1,1000) << endl;
    }
}

int main()
{
    n=read();m=read();q=read();
    if(n>1) solve1();
    else solve2();
    return 0;
}