【分块】【链表】bzoj2738 矩阵乘法

http://www.cnblogs.com/jianglangcaijin/p/3460012.html

首先将矩阵的数字排序。设置size，每次将size个数字插入。插入时，我们用h[i][j]记录该位置的数字是否已经插入；用sum[i][j]表示子矩阵(1,1)到(i,j)已经插入的数字个数总和。每次插入后，暴力扫一次询问，若查询子矩阵的数字个数大于等于K则答案就在此次插入的数字中①；然后将该询问从链表中删除。

①怎样知道是插入的哪个数字呢？暴力O(sqrt(n))地扫描当前这段区间，依次判断这个数插入时是否恰好达到k个。

复杂度真心不知道怎么分析了……跪了……

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<list>
using namespace std;
struct Ask{int x1,y1,x2,y2,k,p;
Ask(const int &x){scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&k); p=x;}Ask(){}};
list<Ask>Q;
typedef list<Ask>::iterator ITER;
struct Val3{int v,x,y;}Ins[250001];
bool operator < (const Val3 &a,const Val3 &b){return a.v<b.v;}
int n,m,sumv[501][501],q,anss[60001];
bool vis[501][501];
void makeblock()
{
    int sz=sqrt(q),sum=1; if(!sz) sz=1; 
    for(;sum*sz<q;++sum)
      {
          int l=(sum-1)*sz+1; int r=sum*sz;
          for(int i=l;i<=r;++i) vis[Ins[i].x][Ins[i].y]=1;
          for(int i=1;i<=n;++i)
            for(int j=1;j<=n;++j)
              sumv[i][j]=sumv[i-1][j]+sumv[i][j-1]-sumv[i-1][j-1]+vis[i][j];
          ITER it=Q.begin();
          while(it!=Q.end())
            {
                int hav=sumv[(*it).x2][(*it).y2]-sumv[(*it).x1-1][(*it).y2]-
                sumv[(*it).x2][(*it).y1-1]+sumv[(*it).x1-1][(*it).y1-1];
                if(hav>=(*it).k)
              {
                for(int i=r;i>=l;--i)
                  if(Ins[i].x>=(*it).x1 && Ins[i].x<=(*it).x2 &&
                     Ins[i].y>=(*it).y1 && Ins[i].y<=(*it).y2)
                       {
                           if(hav==(*it).k)
                             {
                                 anss[(*it).p]=Ins[i].v;
                                 break;
                             } --hav;
                       }
                ITER t=it; ++it; Q.erase(t);
              }
            else ++it;
            }
      }
    int l=(sum-1)*sz+1; int r=q;
    for(int i=l;i<=r;++i) vis[Ins[i].x][Ins[i].y]=1;
    for(int i=1;i<=n;++i)
      for(int j=1;j<=n;++j)
          sumv[i][j]=sumv[i-1][j]+sumv[i][j-1]-sumv[i-1][j-1]+vis[i][j];
    ITER it=Q.begin();
    while(it!=Q.end())
      {
          int hav=sumv[(*it).x2][(*it).y2]-sumv[(*it).x1-1][(*it).y2]-
            sumv[(*it).x2][(*it).y1-1]+sumv[(*it).x1-1][(*it).y1-1];
          if(hav>=(*it).k)
          {
            for(int i=r;i>=l;--i)
              if(Ins[i].x>=(*it).x1 && Ins[i].x<=(*it).x2 &&
                 Ins[i].y>=(*it).y1 && Ins[i].y<=(*it).y2)
                   {
                     if(hav==(*it).k)
                       {
                           anss[(*it).p]=Ins[i].v;
                           break;
                       } --hav;
                  }
            ITER t=it; ++it; Q.erase(t);
          }
        else ++it;
      }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i)
      for(int j=1;j<=n;++j)
        {
            scanf("%d",&Ins[++q].v);
            Ins[q].x=i; Ins[q].y=j;
        }
    sort(Ins+1,Ins+q+1);
    for(int i=1;i<=m;++i) Q.push_back(Ask(i));
    makeblock();
    for(int i=1;i<=m;++i) printf("%d\n",anss[i]);
    return 0;
}