The 2019 Asia Nanchang First Round Online Programming Contest I. Yukino With Subinterval(cdq分治)

最新推荐文章于 2020-10-30 22:55:27 发布

hongism

最新推荐文章于 2020-10-30 22:55:27 发布

阅读量96

点赞数

分类专栏： cdq分治文章标签： cdq

本文链接：https://blog.csdn.net/qq_40851016/article/details/101346122

版权

cdq分治专栏收录该内容

2 篇文章 0 订阅

订阅专栏

I. Yukino With Subinterval

题意：

输入 $n, m (2 e 5)$ ;
接下来输入 $a_1,a_2,\dots,a_n(1\leq a_i\leq n)$ ；
接下来 $m$ 行表示有 $m$ 个操作；两种操作：
1、 $1, p o s, v$ 将 $a_{pos}$ 修改为 $v$ ;
2、 $2, l, r, x, y$ 输出在 $l 到 r$ 范围内数值大于等于 $x$ ，小于等于 $y$ ，相等数的段数。

题解：

可以用BIT套线段树来写，~~然而我太菜了，不会啊！~~
用cdq分治来写，来分析一下：
因为要求段数，就可以直接保留最左边断点的值，这样看 $l - r$ 区间内有多少个 $x\leq a_i\leq y$ ,这样有可能会忽视掉 $l$ 即左端的那一段，因此要判断一下 $a_l$ 和 $a_{l-1}$ 的关系，修改操作，照这种思路修改就好，查询的话，就和找矩形内有多少点是一样的。
注意：加入点不能太多，否则会超时，修改时尽量加入比较少的点来维护正确答案，查询时可以不加入四个点，可以只加入两个点，加的时候只加 $x - y$ 范围内的就好，这样复杂度降了二，不过这里我没有写这个优化。

代码：

//#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+9;
int n,m,a[N];
struct T{
    int t,x,y,v,z;
}b[N*4];
int nb;
int ans[N];int nans;
int f[N];
void add(int x,int y){for(;x<N;x+=(x&-x))f[x]+=y;}
inline int query(int x){int ans=0;while(x)ans+=f[x],x-=(x&-x);return ans;}
bool cmp(T t1,T t2){return t1.x<t2.x;}
T tmp[N*4];
void cdq(int l,int r){
    if(l==r)return;
    int mid=(l+r)>>1;
    cdq(l,mid),cdq(mid+1,r);
    sort(b+l,b+mid+1,cmp),sort(b+mid+1,b+r+1,cmp);
    int p=l;
    for(int i=mid+1;i<=r;i++){
        while(p<=mid&&b[p].x<=b[i].x){if(b[p].z==0)add(b[p].y,b[p].v);p++;}
        if(b[i].z)ans[b[i].v]+=b[i].z*query(b[i].y);
    }
    for(int i=l;i<p;i++)if(b[i].z==0)add(b[i].y,-b[i].v);
}
int main(){
    //freopen("tt.in","r",stdin),freopen("tt.out","w",stdout);
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    cin>>n>>m;
   //scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        cin>>a[i];
        //scanf("%d",&a[i]);
        if(a[i]!=a[i-1])b[++nb]=(T){0,i,a[i],1,0};
    }
    int opt,l,r,x,y,p,v;
    for(int i=1;i<=m;i++){
        cin>>opt;
     //   scanf("%d",&opt);
        if(opt==1){
            cin>>p>>v;
            //scanf("%d%d",&p,&v);
            if(a[p]==v)continue;
            if(a[p-1]!=a[p])b[++nb]=(T){i,p,a[p],-1,0};
            if(a[p+1]!=a[p]&&a[p+1]==v)b[++nb]=(T){i,p+1,a[p+1],-1,0};
            a[p]=v;
            if(a[p-1]!=a[p])b[++nb]=(T){i,p,a[p],1,0};
            //if(a[p+1]!=a[p])b[++nb]=(T){i,p+1,a[p+1],1,0};
        }else{
            cin>>l>>r>>x>>y;
           //scanf("%d%d%d%d",&l,&r,&x,&y);
            ++nans;
            if(a[l]==a[l-1]&&a[l]>=x&&a[l]<=y)ans[nans]++;
            b[++nb]=(T){i,r,y,nans,1};
            b[++nb]=(T){i,l-1,y,nans,-1};
            b[++nb]=(T){i,r,x-1,nans,-1};
            b[++nb]=(T){i,l-1,x-1,nans,1};
        }
    }
  //  for(int i=1;i<=nans;i++)cout<<ans[i]<<" ";cout<<endl;
  //  for(int i=1;i<=nb;i++)cout<<b[i].t<<" "<<b[i].x<<" "<<b[i].y<<" "<<b[i].v<<" "<<b[i].z<<endl;
    cdq(1,nb);
    for(int i=1;i<=nans;i++)cout<<ans[i]<<endl;
    return 0;
}