bzoj 1926: [Sdoi2010]粟粟的书架

 #include<cstdio>
 #include<iostream>
 #define N 201
 #define M 500008
 using namespace std;
 int cnt,r,c,m,ls[*M],rs[*M],root[M],ans,sum[*M],sum1[*M],shu,num[][N][N],num1[][N][N];
 void jia(int l,int r,int x,int &y,int v)
 {
     y=++cnt;
     sum[y]=sum[x]+v;
     sum1[y]=sum1[x]+;
     if(l==r)
       return;
     ls[y]=ls[x];
     rs[y]=rs[x];
     int mid=(l+r)>>;
     if(v<=mid)
       jia(l,mid,ls[x],ls[y],v);
     else
       jia(mid+,r,rs[x],rs[y],v);
     return;
 }
 int xun(int a1,int a2,int l,int r,int v)
 {
     if(!a1)
       return ;
     if(l>=v)
       {
         shu-=sum[a1]-sum[a2];
         return sum1[a1]-sum1[a2];
       }
     int s=,mid=(l+r)>>;
     if(mid>=v)
       s+=xun(ls[a1],ls[a2],l,mid,v);
     s+=xun(rs[a1],rs[a2],mid+,r,v);
     return s;
 }
 void solve1()
 {
     for(int i=;i<=c;i++)
       {
         int a1;
         scanf("%d",&a1);
         jia(,,root[i-],root[i],a1);
       }
     for(int i=;i<=m;i++)
       {
         int x0,y0,x1,y1,v;
         scanf("%d%d%d%d%d",&x0,&y0,&x1,&y1,&v);
         int l=,r=,mid=(l+r)>>,b1,md,md1;
         ans=;
         for(;l<=r;)
           {
             shu=v;
             b1=xun(root[y1],root[y0-],,,mid);
             if(shu<=)
             {
                 md1=v-shu;
                 md=mid;
                 ans=b1;
                 l=mid+;
             }
             else
               r=mid-;
             mid=(l+r)>>;
           }
         if(!ans)
           printf("Poor QLW\n");
         else
           {
             int a1;
             r=xun(root[y1],root[y0-],,,md+);
             r=ans-r;
             l=;
             mid=(l+r)>>;
             for(;l<=r;)
               {
                 if(md1-mid*md>=v)
                   {
                     a1=ans-mid;
                     l=mid+;
                   }
                 else
                   r=mid-;
                 mid=(l+r)>>;
               }
             printf("%d\n",a1);
           }
       }
     return;
 }
 void solve2()
 {
     for(int i=;i<=r;i++)
       for(int j=;j<=c;j++)
         {
             int a1;
             scanf("%d",&a1);
             for(int k=;k<=;k++)
               {
                 num[k][i][j]=(num[k][i][j-]+num[k][i-][j]-num[k][i-][j-]);
                 num1[k][i][j]=(num1[k][i][j-]+num1[k][i-][j]-num1[k][i-][j-]);
                 if(a1>=k)
                   {
                     num[k][i][j]+=a1;
                     num1[k][i][j]++;
                   }
               }
         }
     for(int i=;i<=m;i++)
       {
         int x0,y0,x1,y1,v;
         scanf("%d%d%d%d%d",&x0,&y0,&x1,&y1,&v);
         int l=,r=,mid=(l+r)>>,tm=,tm1=,tm2=;
         for(;l<=r;)
           {
             if(num[mid][x1][y1]-num[mid][x1][y0-]-num[mid][x0-][y1]+num[mid][x0-][y0-]>=v)
               {
                 tm=num1[mid][x1][y1]-num1[mid][x1][y0-]-num1[mid][x0-][y1]+num1[mid][x0-][y0-];
                 tm1=mid;
                 l=mid+;
               }
             else
               r=mid-;
             mid=(l+r)>>;
           }
         if(!tm)
           printf("Poor QLW\n");
         else
           {
             tm2=num[tm1][x1][y1]-num[tm1][x1][y0-]-num[tm1][x0-][y1]+num[tm1][x0-][y0-];
             l=,r=num1[tm1+][x1][y1]-num1[tm1+][x1][y0-]-num1[tm1+][x0-][y1]+num1[tm1+][x0-][y0-];
             r=tm-r;
             mid=(l+r)>>;
             for(;l<=r;)
               {
                  if(tm2-mid*tm1>=v)
                    {
                      ans=tm-mid;
                      l=mid+;
                    }
                 else
                   r=mid-;
                 mid=(l+r)>>;
               }
             printf("%d\n",ans);
           }
       }
     return;
 }
 int main()
 {
     scanf("%d%d%d",&r,&c,&m);
     if(r==)
       solve1();
     else
       solve2();
     return ;
 }

两种情况 对于R，C<=200,用前缀和暴力二分，对于另外的数据，在主席树上二分。