引入:用数组模拟树形结构,可以解决一些区间更新,求和的问题。比线段树更简单一些的数据结构,树状数组能够解决的问题,线段树(还没学)都可以解决。

原理啥的,我也讲不明白,看大佬博客

构建树状数组

单点更新,区间查询

  1. int n,a[1000010],c[1000010];
  2. int lowbit(int x){
  3. return x&-x;//返回二进制最后一位1
  4. }
  5. void updata(int i,int k){//在i的位置上加上k
  6. for(;i<=n;i+=lowbit(i)) c[i]+=k;
  7. }
  8. int getsum(int i){//求a[1]-a[i]的和
  9. int ans=0;
  10. for(;i>0;i-=lowbit(i)) ans+=c[i];
  11. return ans;
  12. }

来一个模板题

HDU1166

代码:

  1. #include<bits/stdc++.h>
  2. using namespace std;
  3. int n,a[1000010],c[1000010];
  4. int lowbit(int x){
  5. return x&-x;
  6. }
  7. void updata(int i,int k){
  8. for(;i<=n;i+=lowbit(i)) c[i]+=k;
  9. }
  10. int getsum(int i){
  11. int ans=0;
  12. for(;i>0;i-=lowbit(i)) ans+=c[i];
  13. return ans;
  14. }
  15. int main(){
  16. int t;
  17. cin>>t;
  18. for(int tot=1;tot<=t;tot++){
  19. cout<<"Case "<<tot<<":"<<endl;
  20. memset(a,0,sizeof a);
  21. memset(c,0,sizeof c);
  22. cin>>n;
  23. for(int i=1;i<=n;i++) {
  24. cin>>a[i];
  25. updata(i,a[i]);
  26. }
  27. string s;
  28. int x,y;
  29. while(cin>>s&&s[0]!='E'){
  30. cin>>x>>y;
  31. if(s[0]=='Q'){
  32. int sum=getsum(y)-getsum(x-1);
  33. cout<<sum<<endl;
  34. }
  35. else if(s[0]=='S'){
  36. updata(x,-y);
  37. }
  38. else if(s[0]=='A'){
  39. updata(x,y);
  40. }
  41. }
  42. }
  43. }

我的板子好像有问题,以后都不敢用了,TLE一下午不知道咋回事,世界未解之谜,代码还是写简单一些不加乱七八糟的东西,不然debug浪费时间

上面这种情况是单点更新、区间查询

单点更新,单点查询

普通数组就可以了。

区间更新、单点查询

给某一段区间同时加上一个数,最后加之后某个点的值,我们知道在某一个区间加上一个数,可以用差分,这里我们用差分数组建立树状数组。

int n,a[1000010],c[1000010];//a原数组,c差分数组

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

void updata(int i,int k){
    for(;i<=n;i+=lowbit(i)) c[i]+=k;
}

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

void add(int l,int r,int k){//在区间[l-r]加上一个数k
    updata(l,k);
    updata(r+1,-k);
}

洛谷P4939

模板题,给你两个操作,0,在给你区间+1,1输出某个点的值。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=1000010;
int n,a[10000010],c[10000010];//a原数组,c差分数组

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

void updata(int i,int k){
    for(;i<=n;i+=lowbit(i)) c[i]+=k;
}

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

void add(int l,int r,int k){//在区间[l-r]加上一个数k
    updata(l,k);
    updata(r+1,-k);
}

signed main(){
    int m;
    scanf("%lld%lld",&n,&m);
    while(m--){
        int op,l,r;
        scanf("%lld",&op);
        if(op&1){
            scanf("%lld",&l);
            int ans=getsum(l);
            printf("%lld\n",ans);
        }
        else {
            scanf("%lld%lld",&l,&r);
            add(l,r,1);
        }
    }
}

洛谷P5057

跟上一题类似,判断最后的奇偶就行了。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=1000010;
int n,a[10000010],c[10000010];//a原数组,c差分数组

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

void updata(int i,int k){
    for(;i<=n;i+=lowbit(i)) c[i]+=k;
}

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

void add(int l,int r,int k){//在区间[l-r]加上一个数k
    updata(l,k);
    updata(r+1,-k);
}

signed main(){
    int m;
    scanf("%lld%lld",&n,&m);
    while(m--){
        int op,l,r;
        scanf("%lld",&op);
        if(op==2){
            scanf("%lld",&l);
            int ans=getsum(l)&1;
            printf("%lld\n",ans);
        }
        else {
            scanf("%lld%lld",&l,&r);
            add(l,r,1);
        }
    }
}

一个简单的整数问题

给了原始数组的模板题,树状数组,区间更新,单点查询。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=1000010;
int n,a[10000010],c[10000010];//a原数组,c差分数组

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

void updata(int i,int k){
    for(;i<=n;i+=lowbit(i)) c[i]+=k;
}

int getsum(int i){
    int ans=a[i];
    for(;i>0;i-=lowbit(i)) ans+=c[i];
    return ans;
}

void add(int l,int r,int k){//在区间[l-r]加上一个数k
    updata(l,k);
    updata(r+1,-k);
}

signed main(){
    int m;
    scanf("%lld%lld",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%lld",&a[i]);
    }
    while(m--){
        int l,r,d;
        char op;
        cin>>op;//用scanf,多了一个空格,要用getchar吃掉,靠
        if(op=='Q'){
            scanf("%lld",&l);
            int ans=getsum(l);
            printf("%lld\n",ans);
        }
        else if(op=='C'){
            scanf("%lld%lld%lld",&l,&r,&d);
            add(l,r,d);
        }
    }
}

区间更新,区间查询

const int N=1000010;
int n,a[1000010],c[1000010];//a原数组,c差分数组
int b1[N],b2[N];
inline int lowbit(int x){return x&-x;}

void updata(int i,int k){
    int tmp=i*k;
    for(;i<=n;i+=lowbit(i)) b1[i]+=k,b2[i]+=tmp;
}

int query(int b[],int i){
    int ans=0;
    for(;i>0;i-=lowbit(i)) ans+=b[i];
    return ans;
}

void add(int l,int r,int k){//在区间[l-r]加上一个数k
    updata(l,k);
    updata(r+1,-k);
}

洛谷P2357

做了这一题发现,区间更新,区间查询也可以用于做单点更新单点查询,区间更新,单点查询。
这一题5个操作,1给区间都加一个数,2给 树状数组 - 图1加上一个数,3给树状数组 - 图2减一个数,4求一个区间的和,5求树状数组 - 图3的值。

代码:

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=1000010;
int n,a[1000010],c[1000010];//a原数组,c差分数组
int b1[N],b2[N];
inline int lowbit(int x){return x&-x;}

void updata(int i,int k){
    int tmp=i*k;
    for(;i<=n;i+=lowbit(i)) b1[i]+=k,b2[i]+=tmp;
}

int query(int b[],int i){
    int ans=0;
    for(;i>0;i-=lowbit(i)) ans+=b[i];
    return ans;
}

void add(int l,int r,int k){//在区间[l-r]加上一个数k
    updata(l,k);
    updata(r+1,-k);
}

int query(int l,int r){
    return (r+1)*query(b1,r)-query(b2,r)-(l*query(b1,l-1)-query(b2,l-1));
}

signed main(){
    int f;
    cin>>n>>f;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        b1[i]+=a[i]-a[i-1];
        b2[i]+=i*(a[i]-a[i-1]);
        int j=i+lowbit(i);
        if(j<=n) b1[j]+=b1[i],b2[j]+=b2[i];
    }
    while(f--){
        int op,l,r,k;
        cin>>op;
        if(op==1){
            cin>>l>>r>>k;
            add(l,r,k);
        }
        else if(op==2){
            cin>>k;
            add(1,1,k);
        }
        else if(op==3){
            cin>>k;
            add(1,1,-k);
        }
        else if(op==4){
            cin>>l>>r;
            int ans=query(l,r);
            cout<<ans<<endl;
        }
        else{
            int ans=query(1,1);
            cout<<ans<<endl;
        }
    }
    return 0;
}

洛谷P1908

求逆序对。

树状数组+离散化

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define sc(n) scanf("%lld",&n)
#define SC(a,b) scanf("%lld%lld",&a,&b)
#define pr(a) printf("%lld",a)
const int N=1000010;
int n,a[1000010],c[1000010],b[N];

int lowbit(int x){
    return x&-x;//返回二进制最后一位1
}

void updata(int i,int k){//在i的位置上加上k
    for(;i<=N;i+=lowbit(i)) c[i]+=k;
}

int getsum(int i){//求a[1]-a[i]的和
    int ans=0;
    for(;i>0;i-=lowbit(i)) ans+=c[i];
    return ans;
}


signed main(){
    sc(n);
    for(int i=1;i<=n;i++) sc(a[i]),b[i]=a[i];//离散化操作
    sort(b+1,b+1+n);
    int cnt=unique(b+1,b+1+n)-(b+1);
    for(int i=1;i<=n;i++){
        a[i]=lower_bound(b+1,b+1+n,a[i])-b;
    }
    int ans=0;
    for(int i=1;i<=n;i++){
        ans+=i-1-getsum(a[i]);//getsum(a[i]) 求得是a[i]之前有多少个比a[i]更小的数,
        updata(a[i],1);
    }
    pr(ans);
    return 0;
}

洛谷P1637 三元上升子序列

题意:Erwin 最近对一种叫 thair 的东西巨感兴趣。。。在含有 n 个整数的序列$ a_1,a_2,…,a_n 树状数组 - 图4i<j<k$ 且 树状数组 - 图5;求一个序列中 thair 的个数。

输入

4
2 1 3 4
5
1 2 2 3 4

输出

2
7

DP+动态规划+离散化

  • 假设树状数组 - 图6表示以a[j]结尾的长度为 i 的上升子序列的个数。

  • 动态转移方程树状数组 - 图7 , 树状数组 - 图8;

  • 要满足两个性质:

    • 树状数组 - 图9;
    • 树状数组 - 图10

求和我们可以用树状数组+离散化来维护,

代码:

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define sc(n) scanf("%lld",&n)
#define SC(a,b) scanf("%lld%lld",&a,&b)
#define pr(a) printf("%lld",a)
const int N=1000010;
int n,a[1000010],c[1000010],b[N],m,ans,f[4][300010];
//f[i][j]表示以a[j]结尾的长度为i,上升子序列的个数
int lowbit(int x){
    return x&-x;//返回二进制最后一位1
}

void updata(int i,int k){//在i的位置上加上k
    for(;i<=n;i+=lowbit(i)) c[i]+=k;
}

int getsum(int i){//求a[1]-a[i]的和
    int ans=0;
    for(;i>0;i-=lowbit(i)) ans+=c[i];
    return ans;
}


signed main(){
    sc(n);
    for(int i=1;i<=n;i++) sc(a[i]),b[i]=a[i];//离散化操作
    sort(b+1,b+1+n);
    int cnt=unique(b+1,b+1+n)-(b+1);
    for(int i=1;i<=n;i++){
        f[1][i]=1;//长度为1的个数初始化为1
        a[i]=lower_bound(b+1,b+1+n,a[i])-b;//离散化
    }
    for(int i=2;i<=3;i++){
        memset(c,0,sizeof(c));
        for(int j=1;j<=n;j++){
            f[i][j]=getsum(a[j]-1);//动态转移方程
            updata(a[j],f[i-1][j]);
        }
    }
    for(int i=1;i<=n;i++){
        ans+=f[3][i];
    }
    pr(ans);
    return 0;
}

力扣1904

题意:给一个序列a,求出以每一个树状数组 - 图11结尾的最长上升子序列的长度。

树状数组+离散化

class Solution {
public:
    vector<int> tree;
    int n;
    int lowbit(int i) {
        return i & -i;
    }
    void update(int id, int x) { // 将id更新为x
        while (id <= n) {
            tree[id] = max(tree[id], x);
            id += lowbit(id);
        }
    }
    int query(int id) {
        int ans = 0;
        while (id) {
            ans = max(ans, tree[id]);
            id -= lowbit(id);
        }
        return ans;
    }
    vector<int> longestObstacleCourseAtEachPosition(vector<int>& obstacles) {
        vector<int> diff = obstacles; // 离散化数组
        sort(diff.begin(), diff.end());
        n = unique(diff.begin(), diff.end()) - diff.begin();
        tree.resize(n + 1, 0);
        vector<int> ans;
        for (auto& i : obstacles) {
            int idx = lower_bound(diff.begin(), diff.begin() + n, i) - diff.begin() + 1;
            int res = query(idx);
            ans.push_back(res + 1);
            update(idx, res + 1);
        }
        return ans;
    }
};