小朋友 center temp tmp namespace out 當前 要求 n)

1926: [Sdoi2010]粟粟的書架

Time Limit: 30 Sec Memory Limit: 552 MB
Submit: 1013 Solved: 396
[Submit][Status][Discuss]

Description

幸福幼兒園 B29 班的粟粟是一個聰明機靈、乖巧可愛的小朋友，她的愛好是畫畫和讀書，尤其喜歡 Thomas H. Co rmen 的文章。粟粟家中有一個 R行C 列的巨型書架，書架的每一個位置都擺有一本書，上數第i 行、左數第j 列 擺放的書有Pi,j頁厚。粟粟每天除了讀書之外，還有一件必不可少的工作就是摘蘋果，她每天必須摘取一個指定的 蘋果。粟粟家果樹上的蘋果有的高、有的低，但無論如何憑粟粟自己的個頭都難以摘到。不過她發現， 如果在腳 下放上幾本書，就可以夠著蘋果；她同時註意到，對於第 i 天指定的那個蘋果，只要她腳下放置書的總頁數之和 不低於Hi，就一定能夠摘到。由於書架內的書過多，父母擔心粟粟一天內就把所有書看完而耽誤了上幼兒園，於是 每天只允許粟粟在一個特定區域內拿書。這個區域是一個矩形，第 i 天給定區域的左上角是上數第 x1i行的左數 第 y1i本書，右下角是上數第 x2i行的左數第y2i本書。換句話說，粟粟在這一天，只能在這﹙x2i－x1i＋1﹚×﹙ y2i－y1i＋1﹚本書中挑選若幹本墊在腳下，摘取蘋果。粟粟每次取書時都能及時放回原位，並且她的書架不會再 撤下書目或換上新書，摘蘋果的任務會一直持續 M天。給出每本書籍的頁數和每天的區域限制及采摘要求，請你告 訴粟粟，她每天至少拿取多少本書，就可以摘到當天指定的蘋果。

Input

第一行是三個正整數R，C，M。 接下來是一個R行C列的矩陣，從上到下、從左向右依次給出了每本書的頁數Pi，j。 接下來M行，第i行給出正整數x1i，y1i，x2i，y2i，Hi，表示第i天的指定區域是﹙x1i，y1i﹚與﹙x2i，y2i﹚間 的矩形，總頁數之和要求不低於Hi。 保證1≤x1i≤x2i≤R，1≤y1i≤y2i≤C。

Output

有M行，第i 行回答粟粟在第 i 天時為摘到蘋果至少需要 拿取多少本書。如果即使取走所有書都無法摘到蘋果， 則在該行輸出“Poor QLW” （不含引號）。

Sample Input

5 5 7
14 15 9 26 53
58 9 7 9 32
38 46 26 43 38
32 7 9 50 28
8 41 9 7 17
1 2 5 3 139
3 1 5 5 399
3 3 4 5 91
4 1 4 1 33
1 3 5 4 185
3 3 4 3 23
3 1 3 3 108

Sample Output

6
15
2
Poor QLW
9
1
3

HINT

對於 10%的數據，滿足 R, C≤10；

對於 20%的數據，滿足 R, C≤40；

對於 50%的數據，滿足 R, C≤200，M≤200,000；

另有 50%的數據，滿足 R＝1，C≤500,000，M≤20,000；

對於 100%的數據，滿足 1≤Pi,j≤1,000，1≤Hi≤2,000,000,000 思路{ 　　第二問主席樹裸題。第一問差不多是一樣的。 　　主席樹上的節點維護$(1,1)$到$(a,b)$的信息，每次就直接查矩形中的數的和的形式查矩形。 　　構建主席樹的時候考場上是基於$(a-1,b)$，暴力把$(a,1)$到$(a,b)$插入主席樹。 　　然後$MLE$了T_T.看著$200$隨便搞啊，算都沒算空間。。。。 　　在$HN01---cjk$指點下發現可以每次建一行的主席樹，然後把這個版本的主席樹和當前$(a-1,b)$版本的主席樹合並一下就是答案了。 　　這樣空間優秀得多。 }

#include<bits/stdc++.h>
#define ll long long
#define RG register
#define il inline
#define N 1000010
#define maxx 20000010
#define mid ((l+r)>>1)
#define db double
using namespace std;
int r,c,q;
ll sum[maxx];int cnt[maxx],idn,ls[maxx],rs[maxx];
int rt[500010],a[500010];
int P(int x,int y){return (x-1)*c+y;}
void Insert(int &x,int l,int r,ll num){
  sum[++idn]=sum[x]+num,cnt[idn]=cnt[x]+1;
  rs[idn]=rs[x],ls[idn]=ls[x];
  x=idn;if(l==r)return;
  if(mid<num)Insert(rs[x],mid+1,r,num);
  else Insert(ls[x],l,mid,num);
}
int merge(int x,int y){
  if(!x||!y)return x+y;
  int p=++idn;
  sum[p]=sum[x]+sum[y];cnt[p]=cnt[x]+cnt[y];
  rs[p]=merge(rs[x],rs[y]);
  ls[p]=merge(ls[x],ls[y]);
  return p;
}
void query(int x,int y,int a,int b,int l,int r,ll num,int tmp){
  if(l==r){
    int cnt=(ceil)((db)num/l);
    printf("%d\n",cnt+tmp);
    return;
  }
  ll temp=sum[rs[x]]+sum[rs[y]]-sum[rs[a]]-sum[rs[b]];
  if(temp<num)query(ls[x],ls[y],ls[a],ls[b],l,mid,num-temp,tmp+cnt[rs[x]]+cnt[rs[y]]-cnt[rs[a]]-cnt[rs[b]]);
  else query(rs[x],rs[y],rs[a],rs[b],mid+1,r,num,tmp);
}
void work(){
  for(int i=1;i<=r;++i)
    for(int j=1;j<=c;++j)
      scanf("%d",&a[P(i,j)]);
  for(int i=1;i<=r;++i){
    int RT=0;
    for(int j=1;j<=c;++j){
      if(i!=1)rt[P(i,j)]=merge(rt[P(i,j)],rt[P(i-1,j)]);
      Insert(RT,1,1000,a[P(i,j)]);
      rt[P(i,j)]=merge(rt[P(i,j)],RT);
    }
  }
  for(int i=1;i<=q;++i){
    int x,y,xx,yy;ll H;
    scanf("%d%d%d%d%lld",&x,&y,&xx,&yy,&H);
    int X=rt[P(xx,yy)],Y=rt[P(x-1,y-1)],XX=rt[P(x-1,yy)],YY=rt[P(xx,y-1)];
    if(x==1)Y=0,XX=0;
    if(y==1)YY=0,Y=0;
    if(sum[X]+sum[Y]-sum[YY]-sum[XX]<H)cout<<"Poor QLW\n";
    else query(X,Y,XX,YY,1,1000,H,0);
  }
}
int main(){
  freopen("susu.in","r",stdin);
  freopen("susu.out","w",stdout);
  scanf("%d%d%d",&r,&c,&q);work();
  return 0;
}

BZOJ1926: [Sdoi2010]粟粟的書架