【分块】【链表】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[];
 bool operator < (const Val3 &a,const Val3 &b){return a.v<b.v;}
 int n,m,sumv[][],q,anss[];
 bool vis[][];
 void makeblock()
 {
     int sz=sqrt(q),sum=; if(!sz) sz=;
     for(;sum*sz<q;++sum)
       {
           int l=(sum-)*sz+; int r=sum*sz;
           for(int i=l;i<=r;++i) vis[Ins[i].x][Ins[i].y]=;
           for(int i=;i<=n;++i)
             for(int j=;j<=n;++j)
               sumv[i][j]=sumv[i-][j]+sumv[i][j-]-sumv[i-][j-]+vis[i][j];
           ITER it=Q.begin();
           while(it!=Q.end())
             {
                 int hav=sumv[(*it).x2][(*it).y2]-sumv[(*it).x1-][(*it).y2]-
                 sumv[(*it).x2][(*it).y1-]+sumv[(*it).x1-][(*it).y1-];
                 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-)*sz+; int r=q;
     for(int i=l;i<=r;++i) vis[Ins[i].x][Ins[i].y]=;
     for(int i=;i<=n;++i)
       for(int j=;j<=n;++j)
           sumv[i][j]=sumv[i-][j]+sumv[i][j-]-sumv[i-][j-]+vis[i][j];
     ITER it=Q.begin();
     while(it!=Q.end())
       {
           int hav=sumv[(*it).x2][(*it).y2]-sumv[(*it).x1-][(*it).y2]-
             sumv[(*it).x2][(*it).y1-]+sumv[(*it).x1-][(*it).y1-];
           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=;i<=n;++i)
       for(int j=;j<=n;++j)
         {
             scanf("%d",&Ins[++q].v);
             Ins[q].x=i; Ins[q].y=j;
         }
     sort(Ins+,Ins+q+);
     for(int i=;i<=m;++i) Q.push_back(Ask(i));
     makeblock();
     for(int i=;i<=m;++i) printf("%d\n",anss[i]);
     return ;
 }