【分块】bzoj1593 [Usaco2008 Feb]Hotel 旅馆

分块，记录每个块内包括左端点的最大连续白段的长度，

整个块内的最大连续白段的长度，

和包括右端点的最大连续白段的长度。

Because 是区间染色，所以要打标记。

至于怎样在O(sqrt(n))的时间内找到最左的白色段呢？非常恶心…… 要考虑跨块的白段和块内的白段，并且顺序不能反。(代码中YuDing()函数)

这题完全体现不出分块编程复杂度低的优势，完全逊色于线段树。

So 综上，对于要进行复杂的成段修改的题目，分块的编程复杂度较高，几乎没有优势，不推荐写。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int len[300],lenl[300],lenr[300],n,m,l[300],r[300],num[50001],sum,sz,x,y;
bool a[50001];
int delta[300],op;
int Res,Num;char C,CH[12];
inline int G()
{
    Res=0;C='*'; 
    while(C<'0'||C>'9')C=getchar();
    while(C>='0'&&C<='9'){Res=Res*10+(C-'0');C=getchar();}
    return Res;
}
inline void P(long long x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
void makeblock()
{
    memset(delta,-1,sizeof(delta));
    sz=sqrt(n);
    for(sum=1;sum*sz<n;sum++)
      {
        l[sum]=(sum-1)*sz+1;
        r[sum]=sum*sz;
        for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
        len[sum]=lenl[sum]=lenr[sum]=sz;
      }
    l[sum]=sz*(sum-1)+1; r[sum]=n;
    for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
    len[sum]=lenl[sum]=lenr[sum]=r[sum]-l[sum]+1;
}
inline void Pushdown(const int &p)
{
    if(delta[p]!=-1)
      {for(int i=l[p];i<=r[p];i++) a[i]=delta[p];
      delta[p]=-1;}
}
inline void Work(const int &Lb,const int &Rb,const int &sym)
{
    Pushdown(num[Lb]);
    for(int i=Lb;i<=Rb;i++) a[i]=sym;
    int cnt=0;
    for(int i=l[num[Lb]];i<=r[num[Lb]];i++) {if(a[i]) break; cnt++;}
    lenl[num[Lb]]=cnt; cnt=0;
    for(int i=r[num[Lb]];i>=l[num[Lb]];i--) {if(a[i]) break; cnt++;}
    lenr[num[Lb]]=cnt; cnt=0;
    int Longest=0;
    for(int i=l[num[Lb]];i<=r[num[Lb]];i++)
      {
          if(a[i]) cnt=0; else cnt++;
          if(cnt>Longest) Longest=cnt;
      }
    len[num[Lb]]=Longest;
}
inline void Update(const int &L,const int &R,const bool &sym)
{
    if(num[L]==num[R]) Work(L,R,sym);
    else
      {
          Work(L,r[num[L]],sym);
          Work(l[num[R]],R,sym);
          for(int i=num[L]+1;i<num[R];i++)
            {
                delta[i]=sym;
                len[i]=lenl[i]=lenr[i]=sym ? 0 : r[i]-l[i]+1;
            }
      }
}
inline void YuDing()
{
    int kua=0,kuasta;
    for(int i=1;i<=sum;i++)
      {
          if(kua) kua+=lenl[i];
        if(kua>=x)//预定sta~sta+x-1 
          {
              Update(kuasta,kuasta+x-1,true);
              P(kuasta);
              return;
          }
          if(len[i]!=r[i]-l[i]+1) kua=0;
        if(!kua&&lenr[i])
            {
                kua=lenr[i];
                kuasta=r[i]-lenr[i]+1;
            }
          if(len[i]>=x)
          {
              Pushdown(i);
              int cnt=0,Nowleft;
            for(int j=l[i];j<=r[i];j++)
              {
                  if(a[j]) cnt=0; else cnt++;
                  if(cnt==1) Nowleft=j;
                  if(cnt==x)
                    {
                        Update(Nowleft,Nowleft+x-1,true);
                        P(Nowleft);
                        return;
                    }
              }
          }
      }
    puts("0");
}
int main()
{
    n=G();m=G();
    makeblock();
    for(int i=1;i<=m;i++)
      {
          op=G();x=G();
          if(op==1) YuDing();
          else {y=G(); Update(x,x+y-1,false);}
      }
    return 0;
}