[BZOJ2738]矩阵乘法（整体二分）

题目：

我是超链接

题解：

又是一个整体二分了，因为本人的手残导致现在才调出来，具体来说就是
首先把权值离散
将所有的点看成是单个插入操作，所有的询问看成是查询操作
二分第k小的数mid，将插入操作按照权值排序了之后就是在一段区间里
将<=mid的数加入到bit里，然后对于每一个询问在bit中查询
如果不够的话就往大的走，如果小的话就往小的走
注意如果向mid,r走的话要将k减去当前的答案
要加快读卡常才能过Orz

代码：

#include <map>
#include <cstdio>
#include <algorithm>
using namespace std;
const int M=60005;
const int N=505;
const int pf=255025;
struct po{int x,y,val;}zb[pf];
struct hh{int x1,y1,x2,y2,id,k,num;}q[M];
int c[N][N],a[N][N],n,ans[M],s[pf];
map<int,int>hash;
int read()
{
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x;
}
void add(int x,int y,int v)
{
    for (int i=x;i<=n;i+=i&(-i))
      for (int j=y;j<=n;j+=j&(-j)) c[i][j]+=v;
}
int qurry(int x,int y)
{
    int ans=0;
    for (int i=x;i>=1;i-=i&(-i))
      for (int j=y;j>=1;j-=j&(-j))
        ans+=c[i][j];
    return ans;
}
int cmp(hh a,hh b){return a.num<b.num;}
void ef(int l,int r,int a,int b,int x,int y)
{
    if (a>b) return;
    int mid=(l+r)>>1;
    if (l==r)
    {
        for (int i=a;i<=b;i++) ans[q[i].id]=l;
        return;
    }
    int pm=0;
    for (int i=x;i<=y;i++)
      if (zb[i].val<=mid) pm=i,add(zb[i].x,zb[i].y,1);
      else break;
    int pa=0,pb=b-a+1;
    for (int i=a;i<=b;i++)
    {
        int t=qurry(q[i].x2,q[i].y2)-qurry(q[i].x1-1,q[i].y2)-qurry(q[i].x2,q[i].y1-1)+qurry(q[i].x1-1,q[i].y1-1);
        if (t>=q[i].k) q[i].num=++pa;
        else q[i].num=++pb,q[i].k-=t;
    }
    for (int i=x;i<=pm;i++) add(zb[i].x,zb[i].y,-1);
    sort(q+a,q+b+1,cmp);
    ef(l,mid,a,a+pa-1,x,pm);
    ef(mid+1,r,a+pa,b,pm+1,y);
}
int cmp1(po a,po b){return a.val<b.val;}
int main()
{
    int Q;n=read();Q=read();int num=0;
    for (int i=1;i<=n;i++) 
      for (int j=1;j<=n;j++) a[i][j]=read(),s[++num]=a[i][j];
    sort(s+1,s+num+1);
    int ha=unique(s+1,s+num+1)-s-1;
    for (int i=1;i<=ha;i++) hash[s[i]]=i;
    num=0;
    for (int i=1;i<=n;i++)
      for (int j=1;j<=n;j++) zb[++num].val=hash[a[i][j]],zb[num].x=i,zb[num].y=j;
    for (int i=1;i<=Q;i++) q[i].x1=read(),q[i].y1=read(),q[i].x2=read(),q[i].y2=read(),q[i].k=read(),q[i].id=i;
    sort(zb+1,zb+num+1,cmp1);
    ef(1,ha,1,Q,1,num);
    for (int i=1;i<=Q;i++) printf("%d\n",s[ans[i]]);
}