zoukankan      html  css  js  c++  java
  • 线段树练习

    试题传送:http://codevs.cn/problem/1080/

    1080 线段树练习

     

     时间限制: 1 s
     空间限制: 128000 KB
     题目等级 : 钻石 Diamond
     
     
     
    题目描述 Description

    一行N个方格,开始每个格子里都有一个整数。现在动态地提出一些问题和修改:提问的形式是求某一个特定的子区间[a,b]中所有元素的和;修改的规则是指定某一个格子x,加上或者减去一个特定的值A。现在要求你能对每个提问作出正确的回答。1≤N<100000,,提问和修改的总数m<10000条。

    输入描述 Input Description

    输入文件第一行为一个整数N,接下来是n行n个整数,表示格子中原来的整数。接下一个正整数m,再接下来有m行,表示m个询问,第一个整数表示询问代号,询问代号1表示增加,后面的两个数x和A表示给位置X上的数值增加A,询问代号2表示区间求和,后面两个整数表示a和b,表示要求[a,b]之间的区间和。

    输出描述 Output Description

    共m行,每个整数

    样例输入 Sample Input

    6

    3

    4

    1 3 5

    2 1 4

    1 1 9

    2 2 6

    样例输出 Sample Output

    22

    22

    数据范围及提示 Data Size & Hint

    1≤N≤100000, m≤10000 。

    分类标签 Tags 

    代码

    #include<cstdio>
    #include<cstdlib>
    #include<iostream>
    #include<cstring>
    #include<algorithm>
    #include<vector>
    #define MAXN 101000
    #define MAXT MAXN*4
    #define lch now<<1
    #define rch ((now<<1)+1)
    #define smid ((l+r)>>1)
    int a[MAXN];
    using namespace std;
    struct node
    {
            int sum;
            int l,r;
    }sgt[MAXT];
    int read(){
        register int f=1,x=0;
        register char ch=getchar();
        while(ch>'9'||ch<'0'){if(ch=='-') f=-1;ch=getchar();}
        while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
        return x*f;
    }
    void build(int now,int l,int r){
        sgt[now].l=l;
        sgt[now].r=r;
        if(l==r) {sgt[now].sum=a[l];return;}
        build(lch,l,smid);
        build(rch,smid+1,r);
        sgt[now].sum=sgt[lch].sum+sgt[rch].sum;
    }
    void modify(int now,int pos,int l,int r,int v){
        if(l==r){
            sgt[now].sum+=v;return;
        }
        if(pos<=smid)
           modify(lch,pos,l,smid,v);
        else 
           modify(rch,pos,smid+1,r,v);
        sgt[now].sum=sgt[lch].sum+sgt[rch].sum;
        
    }
    int query(int now,int l,int r,int x,int y){
        if(l==x&&r==y) return sgt[now].sum;
        if(y<=smid) 
           return query(lch,l,smid,x,y);
        else if(smid<x)
           return query(rch,smid+1,r,x,y);
        else
           return query(lch,l,smid,x,smid)+query(rch,smid+1,r,smid+1,y);
    }
    int main(){
        int n=read();
        for(int i=1;i<=n;i++)
           a[i]=read();
        build(1,1,n);
        int m=read();
        for(int i=1;i<=m;i++){
            int opt=read();
            if(opt==1){
               int pos=read(),v=read();
               modify(1,pos,1,n,v);
            }else{
                int x=read(),y=read();
                printf("%d
    ",query(1,1,n,x,y));
            }    
        }
        return 0;
    }

     更新一下代码风格

    #include<bits/stdc++.h>
    using namespace std;
    #define N 100010
    #define lc k<<1
    #define rc k<<1|1
    #define mid (l+r>>1)
    int tmp[N],a[N<<2];
    void build(int k,int l,int r){
        if(l==r){
            a[k]=tmp[l];return ;
        }
        build(lc,l,mid);
        build(rc,mid+1,r);
        a[k]=a[lc]+a[rc];
    }
    void add(int k,int l,int r,int pos,int y){
        if(l==r){
            a[k]+=y;return ;
        }
        if(pos<=mid) 
            add(lc,l,mid,pos,y);
        else 
            add(rc,mid+1,r,pos,y);
        a[k]=a[lc]+a[rc];
    }
    int query(int k,int l,int r,int x,int y){
        if(l==x&&r==y) return a[k];
        if(y<=mid) return query(lc,l,mid,x,y);
        else if(x>mid) return query(rc,mid+1,r,x,y);
        else return query(lc,l,mid,x,mid)+query(rc,mid+1,r,mid+1,y);
    }
    int main(){
        int n,m,opt,x,y;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d",tmp+i);
        build(1,1,n);
        scanf("%d",&m);
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&opt,&x,&y);
            if(opt==1) add(1,1,n,x,y);
            else printf("%d
    ",query(1,1,n,x,y));
        }
        return 0;
    }

    树状数组版

    #include<bits/stdc++.h>
    using namespace std;
    #define N 100010
    int n,m,opt,x,y,c[N];
    int lowbit(int x){
        return x&-x;
    }
    void updata(int p,int v){
        for(int i=p;i<=n;i+=lowbit(i))
            c[i]+=v;
    }
    int query(int p){
        int ans=0;
        for(int i=p;i>=1;i-=lowbit(i))
            ans+=c[i];
        return ans;
    }
    int main(){
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d",&x),updata(i,x);//初始化树状数组 
        scanf("%d",&m);
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&opt,&x,&y);
            if(opt==1) updata(x,y);
            else printf("%d
    ",query(y)-query(x-1));//类似一个前缀和 
        }
        return 0;
    }

    -----------------------------------------------------------------------------------------------------------------------------------------------

    华丽的分割线

    ------------------------------------------------------------------------------------------------------------------------------------------------

    ps:补一下树状数组

    如果给定一个数组,要你求里面所有数的和,一般都会想到累加。但是当那个数组很大的时候,累加就显得太耗时了,时间复杂度为O(n),并且采用累加的方法还有一个局限,那就是,当修改掉数组中的元素后,仍然要你求数组中某段元素的和,就显得麻烦了。所以我们就要用到树状数组,他的时间复杂度为O(lgn),相比之下就快得多。下面就讲一下什么是树状数组:

             一般讲到树状数组都会少不了下面这个图:

             

             

    下面来分析一下上面那个图看能得出什么规律:

    据图可知:c1=a1,c2=a1+a2,c3=a3,c4=a1+a2+a3+a4,c5=a5,c6=a5+a6,c7=a7,c8=a1+a2+a3+a4+a5+a6+a7+a8,c9=a9,c10=a9+a10,c11=a11........c16=a1+a2+a3+a4+a5+.......+a16。

    分析上面的几组式子可知,当 i 为奇数时,ci=ai ;当 i 为偶数时,就要看 i 的因子中最多有二的多少次幂,例如,6 的因子中有 2 的一次幂,等于 2 ,所以 c6=a5+a6(由六向前数两个数的和),4 的因子中有 2 的两次幂,等于 4 ,所以 c4=a1+a2+a3+a4(由四向前数四个数的和)。

    (一)有公式:cn=a(n-a^k+1)+.........+an(其中 k 为 n 的二进制表示中从右往左数的 0 的个数)。

    那么,如何求 a^k 呢?求法如下:

    int lowbit(int x){
        return x&-x;
    }

    lowbit()的返回值就是 2^k 次方的值。

    求出来 2^k 之后,数组 c 的值就都出来了,接下来我们要求数组中所有元素的和。

     

    (二)求数组的和的算法如下:

    (1)首先,令sum=0,转向第二步;

    (2)接下来判断,如果 n>0 的话,就令sum=sum+cn转向第三步,否则的话,终止算法,返回 sum 的值;

    (3)n=n - lowbit(n)(将n的二进制表示的最后一个零删掉),回第二步。

    代码实现:

    int query(int p){
        int ans=0;
        for(int i=p;i>=1;i-=lowbit(i))
            ans+=c[i];
        return ans;
    }

     

    (三)当数组中的元素有变更时,树状数组就发挥它的优势了,算法如下(修改为给某个节点 i 加上 x ):

    (1)当 i<=n 时,执行下一步;否则的话,算法结束;

    (2)ci=ci+x ,i=i+lowbit(i)(在 i 的二进制表示的最后加零),返回第一步。

    代码实现:

    void updata(int p,int v){
        for(int i=p;i<=n;i+=lowbit(i))
            c[i]+=v;
    }
  • 相关阅读:
    Django进阶(三)
    Django进阶(二)
    全局变量初始化顺序探究
    调试实战 —— dll 加载失败之全局变量初始化篇
    10 个实验搞懂命令提示符
    帮 vs2019 找回丢失的 SDK
    从堆里找回“丢失”的代码
    排错实战 —— 解决 c++ 工程编译错: error C2059 'string' illegal token on right side of xxx
    善用 vs 中的错误列表和输出窗口,高效查找 C++ 多工程编译错误
    使用 VMware + win10 + vs2019 从零搭建双机内核调试环境
  • 原文地址:https://www.cnblogs.com/shenben/p/5459703.html
Copyright © 2011-2022 走看看