[Sdoi2010]粟粟的书架
首先看到这种另外\(50\%\)的数据,就知道肯定是强行2合1的题目了。
如果\(R,C \le 200\),显然可以二分,二分最小的书籍的厚度,然后显然如果最小的越大,用的书就越少。然后直接\(check\)即可。
考虑当\(R=1\)时,在主席树上二分即可。
/*
mail: mleautomaton@foxmail.com
author: MLEAutoMaton
This Code is made by MLEAutoMaton
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define ll long long
#define re register
#define REP(a,b,c) for(int a=b;a<=c;a++)
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
inline int gi()
{
int f=1,sum=0;char ch=getchar();
while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
return f*sum;
}
#define ERROR puts("Poor QLW");
int r,c,m,h,x1,x2,y1,y2;
const int N=210,NN=2000010;
int p[N][N],sum[N][N][1010],num[N][N][1010];
int Num(int x1,int y1,int x2,int y2,int mid){return num[x2][y2][mid]-num[x2][y1-1][mid]-num[x1-1][y2][mid]+num[x1-1][y1-1][mid];}
int Sum(int x1,int y1,int x2,int y2,int mid){return sum[x2][y2][mid]-sum[x2][y1-1][mid]-sum[x1-1][y2][mid]+sum[x1-1][y1-1][mid];}
int number(int x1,int y1,int x2,int y2,int h,int mid){return Num(x1,y1,x2,y2,mid)-(Sum(x1,y1,x2,y2,mid)-h)/mid;}
bool check(int mid){return Sum(x1,y1,x2,y2,mid)>=h;}
void work1()
{
int mx=0;
REP(i,1,r)REP(j,1,c)p[i][j]=gi(),mx=max(mx,p[i][j]);
REP(k,0,mx)REP(i,1,r)REP(j,1,c)
{
num[i][j][k]=num[i-1][j][k]+num[i][j-1][k]-num[i-1][j-1][k]+(p[i][j]>=k);
sum[i][j][k]=sum[i-1][j][k]+sum[i][j-1][k]-sum[i-1][j-1][k]+((p[i][j]>=k)?p[i][j]:0);
}
while(m--)
{
x1=gi(),y1=gi(),x2=gi(),y2=gi(),h=gi();
int l=0,r=mx,ans=-1;
while(l<=r)
{
int mid=(l+r)>>1;
if(check(mid)){ans=mid;l=mid+1;}
else r=mid-1;
}
if(ans!=-1)printf("%d\n",number(x1,y1,x2,y2,h,ans));
else ERROR;
}
}
struct node
{
int ls,rs,sum,siz;
}t[NN*5];
int a[NN],n,cnt,rt[NN];
void modify(int x,int &o,int l,int r,int pos)
{
o=++cnt;t[o]=t[x];
t[o].sum+=pos;t[o].siz++;if(l==r)return;
int mid=(l+r)>>1;
if(pos<=mid)modify(t[x].ls,t[o].ls,l,mid,pos);
else modify(t[x].rs,t[o].rs,mid+1,r,pos);
}
int query(int x,int y,int l,int r,int sum)
{
if(l==r)return sum/l+(sum%l>0);
int mid=(l+r)>>1,S=(t[t[y].rs].sum-t[t[x].rs].sum);
if(S>=sum)return query(t[x].rs,t[y].rs,mid+1,r,sum);
else return query(t[x].ls,t[y].ls,l,mid,sum-S)+t[t[y].rs].siz-t[t[x].rs].siz;
}
void work2()
{
n=c;
for(int i=1;i<=n;i++)a[i]=gi();
for(int i=1;i<=n;i++)
modify(rt[i-1],rt[i],1,1000,a[i]),a[i]+=a[i-1];
while(m--)
{
x1=gi(),y1=gi(),x2=gi(),y2=gi();h=gi();
if(a[y2]-a[y1-1]<h){ERROR;}
else printf("%d\n",query(rt[y1-1],rt[y2],1,1000,h));
}
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.in","r",stdin);
#endif
r=gi();c=gi();m=gi();
if(r!=1 || c<=200)work1();
else work2();
return 0;
}