2019南昌网络赛-I. Yukino With Subinterval 线段树套树状数组

TMD。。。这题卡内存卡的真优秀。。。

所以以后还是别用主席树的写法。。。不然怎么死的都不知道。。。

树套树中，主席树方法开权值线段树。。。会造成空间的浪费。。。这道题内存卡的很紧。。。

由于树套树已经不需要持久化了，直接动态开点就完事了。。。用主席树方法开过不去，要么超内存，要么越界。。。

大概思路。。。这题要求的[L,R]区间内，满足x<=a[i]<=y的连续的段数，

这题大概是个套路题，我们很容易想到，我们把连续的区间变成单点，把一段区间，类似1 1 1 3 5 变成 1 0 0 3 5 这样我们

只需要维护区间内部，在某个范围内数字的个数，这样的求法有个很显然的弊端如果1 0 0 3 5 序列， 我们求的区间是 2到5

这样求出来的答案是2，但是答案是3，究其原因，是我们b[l]端点出了问题，如果b[l]是0，我们求出来的就比答案少一个，那

么如何求出答案呢？？？很简单，可以先求出a[l+1]-a[r]的个数，这个算出来的是肯定不准确的，在b[l+1]=0的时候，在其他的

情况下正确的，如果b[l+1]=0 或者b[l+1]!=0，我们只要判断a[l]是否满足，如果满足条件就+1，这样就保证b[l+1]的情况，并且

在b[l+1]!=0的情况下，也是正确的。

考虑修改，单点更新，如果当前点是a[l]==a[l-1]那么a[l]将变成一个新的左端点。

如果修改后,a[l]==a[l+1]那么a[l+1]将不再是一个左端点，需要舍去。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxx = 2e5+6;
struct node{
   int l,r;
   int cnt;
}tree[maxx*220];
int cnt,n,cnt1,cnt2;
int root[maxx],a[maxx],trl[maxx],trr[maxx];
void inserts(int &now,int l,int r,int pos,int w){
   if(!now)now=++cnt;
   tree[now].cnt+=w;
   if (l==r){
    return ;
   }
   int mid=(l+r)>>1;
   if (pos<=mid){
    inserts(tree[now].l,l,mid,pos,w);
   }else{
    inserts(tree[now].r,mid+1,r,pos,w);
   }
}
int lowbit(int x){
   return x&(-x);
}
void add(int x,int w){
   for(int i=x;i<=n;i+=lowbit(i)){
       inserts(root[i],1,n,a[x],w);
   }
}
int query(int l,int r,int ql,int qr){
    if (r<ql || l>qr){
        return 0;
    }
    int tmpl[1000],tmpr[1000];
    int s=0,mid=(l+r)>>1,sum=0;
    for(int i=1;i<=cnt1;i++)s-=tree[trl[i]].cnt,tmpl[i]=trl[i];
    for(int i=1;i<=cnt2;i++)s+=tree[trr[i]].cnt,tmpr[i]=trr[i];
    if (ql<=l && r<=qr){
         return s;
    }
    if (mid>=ql){
        for (int i=1;i<=cnt1;i++){
            trl[i]=tree[tmpl[i]].l;
        }
        for (int i=1;i<=cnt2;i++){
            trr[i]=tree[tmpr[i]].l;
        }
        sum+=query(l,mid,ql,qr);
    }
    if (mid<qr){
        for (int i=1;i<=cnt1;i++){
            trl[i]=tree[tmpl[i]].r;
        }
        for (int i=1;i<=cnt2;i++){
            trr[i]=tree[tmpr[i]].r;
        }
        sum+=query(mid+1,r,ql,qr);
    }
    return sum;
}
int main(){
  int m;
  scanf("%d%d",&n,&m);
    cnt=cnt1=cnt2=0;
  for (int i=1;i<=n;i++){
    scanf("%d",&a[i]);
    if (a[i]!=a[i-1])add(i,1);
  }
  int op;
  int pos,v,l,r,x,y;
  while(m--){
    scanf("%d",&op);
    if(op==1){
        scanf("%d%d",&pos,&v);
        if (a[pos]==v)continue;
        if(a[pos]!=a[pos-1]){
            add(pos,-1);
        }
        if(a[pos]==a[pos+1]){
            add(pos+1,1);
        }else if (a[pos+1]==v){
            add(pos+1,-1);
        }
        a[pos]=v;
        if (a[pos]!=a[pos-1]){
            add(pos,1);
        }
    }else {
        scanf("%d%d%d%d",&l,&r,&x,&y);
        cnt1=cnt2=0;
        int f=0;
        for (int i=l;i;i-=lowbit(i)){
            trl[++cnt1]=root[i];
        }
        for (int i=r;i;i-=lowbit(i)){
            trr[++cnt2]=root[i];
        }
        if (a[l]>=x && a[l]<=y)f++;
        printf("%d\n",query(1,n,x,y)+f);
    }
  }
  return 0;
}

原文地址：https://www.cnblogs.com/bluefly-hrbust/p/11516694.html