【分块】【哈希】bzoj3578 GTY的人类基因组计划2

每个房间用一个集合来维护，具体来说，就是给1-n的数每个数一个long long的hash值，往集合S里insert(i)，就是S^=HASH[i]；erase(i)，也是S^=HASH[i]。

用map/set维护某个集合是否已经做过实验。

分块，对每个块维护一个maxv[i]，代表当前该块内的答案值，要注意，若某个房间的集合的vis[S[i]]==true，则它对其所在块没有贡献。

分块相对于线段树/平衡树/树套树的很大的好处，就是思路简单暴力，代码短，常数小，应对范围广，可以方便地扩张。

像这种区间修改的问题，可以类比线段树打懒标记的思路来进行，使得每次修改是O(sqrt(n))的。

Pushdown函数需要每次对块进行“零散地”操作的之前进行。

 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 #include<cstdlib>
 #include<ctime>
 #include<map>
 using namespace std;
 unsigned long long HASH[];
 struct SET
 {
     unsigned long long S;int size;
     void insert(const int &v){S^=HASH[v];size++;}
     void erase(const int &v){S^=HASH[v];size--;}
 };
 int n,m,q,l[],r[],sumv[],sum,sz,num[],x,y,pos[];
 bool delta[];
 SET sets[];
 map<unsigned long long,bool>vis;
 char op[];
 int Res,Num;char C,CH[];
 inline int G()
 {
     Res=;C='*';
     while(C<''||C>'')C=getchar();
     while(C>=''&&C<=''){Res=Res*+(C-'');C=getchar();}
     return Res;
 }
 inline void P(long long x)
 {
     Num=;if(!x){putchar('');puts("");return;}
     while(x>)CH[++Num]=x%,x/=;
     while(Num)putchar(CH[Num--]+);
     puts("");
 }
 void makeblock()
 {
     srand();
     for(int i=;i<=n;i++) HASH[i]=(unsigned long long)((unsigned long long)rand()<<)|((unsigned long long)rand()<<)|((unsigned long long)rand()<<)|((unsigned long long)rand());
     sz=sqrt((double)m*0.3);
     for(int i=;i<=n;i++) {sets[].insert(i); pos[i]=;}
     sumv[]=n;
     for(sum=;sum*sz<m;sum++)
       {
         l[sum]=(sum-)*sz+;
         r[sum]=sum*sz;
         for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
       }
     l[sum]=sz*(sum-)+; r[sum]=m;
     for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
 }
 inline void Pushdown(const int &p)//对块p下传标记,Pushdown操作需要每次对块进行零散操作之前进行
 {
     if(delta[p])
       {
           for(int i=l[p];i<=r[p];i++)
           if(sets[i].size)
             vis[sets[i].S]=true;
         delta[p]=false;
       }
 }
 inline void Move()
 {
     Pushdown(num[pos[x]]);
     bool flag1=vis[sets[pos[x]].S];
     sets[pos[x]].erase(x);
     bool flag2=vis[sets[pos[x]].S];
     if(flag1&&!flag2) sumv[num[pos[x]]]+=sets[pos[x]].size;
     else if(!flag1&&flag2) sumv[num[pos[x]]]-=(sets[pos[x]].size+);
     else if(!flag1&&!flag2) sumv[num[pos[x]]]--;
     Pushdown(num[y]);
     flag1=vis[sets[y].S];
     sets[y].insert(x);
     flag2=vis[sets[y].S];
     if(flag1&&!flag2) sumv[num[y]]+=sets[y].size;
     else if(!flag1&&flag2) sumv[num[y]]-=(sets[y].size-);
     else if(!flag1&&!flag2) sumv[num[y]]++;
     pos[x]=y;
 }
 inline void Update_Query()
 {
     int ans=;
     Pushdown(num[x]);
     Pushdown(num[y]);
     if(num[x]+>=num[y])
       {
         for(int i=x;i<=y;i++)
           if(sets[i].size)
             if(!vis[sets[i].S])
               {
                   ans+=sets[i].size;
                   sumv[num[i]]-=sets[i].size;
                 vis[sets[i].S]=true;
               }
       }
     else
       {
           for(int i=x;i<=r[num[x]];i++)
             if(sets[i].size)
             if(!vis[sets[i].S])
               {
                   ans+=sets[i].size;
                   sumv[num[x]]-=sets[i].size;
                 vis[sets[i].S]=true;
               }
           for(int i=l[num[y]];i<=y;i++)
             if(sets[i].size)
             if(!vis[sets[i].S])
               {
                   ans+=sets[i].size;
                   sumv[num[y]]-=sets[i].size;
                 vis[sets[i].S]=true;
               }
           for(int i=num[x]+;i<num[y];i++)
           {
             if(!delta[i])
               ans+=sumv[i];
             delta[i]=true;
             sumv[i]=;
           }
       }
     P(ans);
 }
 int main()
 {
     n=G();m=G();q=G();
     makeblock();
     for(;q>;q--)
       {
           scanf("%s",op);x=G();y=G();
           if(op[]=='C') Move();
           else Update_Query();
       }
     return ;
 }