bzoj2738 矩阵乘法 整体二分

链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2738

题意：求出子矩阵$K$小值。

果然除了最小割题面给什么算法不能用什么算法

这个玩意……实际上并没有强制在线的要求……退一步讲这个时空要求也不太适合主席树……然而……$wcx$神犇正面硬上……结果在$MLE$和$RE$间抓狂了一下午……

实际上区间$K$大值离线通解大概就是整体二分吧……二分最大值和最小值，动态维护，如果这个区间有这么多个小于它的数，就设为$mid$，接下来继续递归修改。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=505,maxq=60005;
int C[maxn][maxn],n,q,cnt;
int lowbit(int x)
{
    return x&(-x);
}
int Query(int x,int y)
{
    int res=0;
    for(;x;x-=lowbit(x))
        for(int j=y;j;j-=lowbit(j))res+=C[x][j];
    return res;
}
void Modify(int x,int y,int val)
{
    for(;x<=n;x+=lowbit(x))
        for(int j=y;j<=n;j+=lowbit(j))C[x][j]+=val;
}
struct point
{
    int x,y,val;
    bool operator <(const point &b)const
    {
        return val<b.val;
    }
}P[maxn*maxn];
struct ques
{
    int x1,y1,x2,y2,k;
}Q[maxq];
int id[maxq],tmp[maxq],ans[maxq],T;
bool vis[maxq];
int query(int x)
{
    int x1=Q[x].x1,x2=Q[x].x2,y1=Q[x].y1,y2=Q[x].y2;
    return Query(x2,y2)-Query(x1-1,y2)-Query(x2,y1-1)+Query(x1-1,y1-1);
}
void Solve(int l,int r,int L,int R)
{
    if(l>r)return;if(L==R)return;
    int mid=(L+R)>>1;
    while(P[T+1].val<=mid&&T<cnt)Modify(P[T+1].x,P[T+1].y,1),T++;
    while(P[T].val>mid)Modify(P[T].x,P[T].y,-1),T--;
    int num=0;
    for(int i=l;i<=r;i++)
    {
        if(query(id[i])>=Q[id[i]].k)vis[i]=1,ans[id[i]]=mid,num++;
        else vis[i]=0;
    }
    int l1=l,l2=l+num;
    for(int i=l;i<=r;i++)
        if(vis[i])tmp[l1++]=id[i];else tmp[l2++]=id[i];
    for(int i=l;i<=r;i++)id[i]=tmp[i];
    Solve(l,l1-1,L,mid);Solve(l1,l2-1,mid+1,R);
}
int haha()
{
    scanf("%d%d",&n,&q);int maxx=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            P[++cnt].x=i,P[cnt].y=j;
            scanf("%d",&P[cnt].val);maxx=max(P[cnt].val,maxx);
        }
    sort(P+1,P+cnt+1);
    for(int i=1;i<=q;i++)scanf("%d%d%d%d%d",&Q[i].x1,&Q[i].y1,&Q[i].x2,&Q[i].y2,&Q[i].k);
    for(int i=1;i<=q;i++)id[i]=i;
    Solve(1,q,0,maxx+1);
    for(int i=1;i<=q;i++)printf("%d
",ans[i]);
}
int sb=haha();
int main(){;}