bzoj1926: [Sdoi2010]粟粟的书架

题目链接

传送门

题解

大佬好强

看了看数据范围，不是很明白为什么明明两道题却硬要挤在一道题里

0.5倍经验orz……

前一半的数据$R,C<=200$

前缀和+二分

$num[i][j][k]表示(1,1)到(i,j)中比大于等于k的数有几个$
$sum[i][j][k]表示(1,1)到(i,j)中比大于等于k的数的总和$

二分找k+1，对k直接暴力填

后一半的数据就是个序列

直接开个主席树

二分一下取多少本书

然后直接树上跑一跑就行

ps：主席树查询时应优先挑选高度更高的书，原因应该不用解释

// luogu-judger-enable-o2
//minamoto
#include<bits/stdc++.h>
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
inline int read(){
    #define num ch-'0'
    char ch;bool flag=0;int res;
    while(!isdigit(ch=getc()))
    (ch=='-')&&(flag=true);
    for(res=num;isdigit(ch=getc());res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
int n,m,q;
namespace solve1{
    #define N 205
    int a[N][N],num[N][N][1005],sum[N][N][1005];
    int mx=0;
    int getsum(int a,int b,int c,int d,int h){
        return sum[c][d][h]-sum[c][b][h]-sum[a][d][h]+sum[a][b][h];
    }
    int getnum(int a,int b,int c,int d,int h){
        return num[c][d][h]-num[c][b][h]-num[a][d][h]+num[a][b][h];
    }
    int getans(int a,int b,int c,int d,int k,int need){
        int res=getnum(a,b,c,d,k+1);
        need-=getsum(a,b,c,d,k+1);
        if(need>=0){
            int num=need/k;
            res+=num;
            need-=num*k;
            if(need>0) ++res,need-=k;
        }
        return res;
    }
    bool check(int a,int b,int c,int d,int mid,int h){
        return getsum(a,b,c,d,mid)>h?1:0;
    }
    void main(){
        for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j)
        cmax(mx,a[i][j]=read());
        for(int k=0;k<=mx;++k)
        for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j){
            num[i][j][k]=num[i-1][j][k]+num[i][j-1][k]-num[i-1][j-1][k]+(a[i][j]>=k?1:0);
            sum[i][j][k]=sum[i-1][j][k]+sum[i][j-1][k]-sum[i-1][j-1][k]+(a[i][j]>=k?a[i][j]:0);
        }
        while(q--){
            int a,b,c,d,h;
            a=read(),b=read(),c=read(),d=read(),h=read();
            --a,--b;
            if(getsum(a,b,c,d,0)<h){
                puts("Poor QLW");continue;
            }
            int l=1,r=mx,ans=0;
            while(l<=r){
                int mid=(l+r)>>1;
                if(check(a,b,c,d,mid,h)) ans=mid,l=mid+1;
                else r=mid-1;
            }
            if(!ans) ans=1;
            printf("%d
",getans(a,b,c,d,ans,h));
        }
    }
    #undef N
}
namespace solve2{
    #define N 500005
    int sum[N<<4],L[N<<4],R[N<<4],sz[N<<4];
    int rt[N];
    int mx=0,cnt=0;
    void update(int last,int &now,int l,int r,int x){
        sz[now=++cnt]=sz[last]+1,sum[now]=sum[last];
        if(l==r){sum[now]+=x;return;}
        int mid=(l+r)>>1;
        if(x<=mid) R[now]=R[last],update(L[last],L[now],l,mid,x);
        else L[now]=L[last],update(R[last],R[now],mid+1,r,x);
        sum[now]=sum[L[now]]+sum[R[now]];
    }
    int query(int u,int v,int l,int r,int k){
        if(sz[v]-sz[u]<=k) return sum[v]-sum[u];
        if(l==r) return (sum[v]-sum[u])/(sz[v]-sz[u])*k;
        int mid=(l+r)>>1;
        if(sz[R[v]]-sz[R[u]]>=k) return query(R[u],R[v],mid+1,r,k);
        else return query(L[u],L[v],l,mid,k-(sz[R[v]]-sz[R[u]]))+sum[R[v]]-sum[R[u]];
    }
    void main(){
        mx=1000;
        for(int i=1;i<=m;++i){
            int x=read();
            update(rt[i-1],rt[i],1,mx,x);
        }
        while(q--){
            int a=read(),b=read(),c=read(),d=read(),h=read();
            --b;
            int l=1,r=d-b+1;
            while(l<r){
                int mid=(l+r)>>1;
                if(query(rt[b],rt[d],1,mx,mid)>=h) r=mid;
                else l=mid+1;
            }
            if(r==d-b+1) puts("Poor QLW");
            else printf("%d
",l);
        }
    }
}
int main(){
    //freopen("testdata.in","r",stdin);
    n=read(),m=read(),q=read();
    if(n!=1) solve1::main();
    else solve2::main();
    return 0;
}