树状数组

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

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

构建树状数组

单点更新,区间查询

int n,a[1000010],c[1000010];

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;
}

来一个模板题

HDU1166

代码:

#include<bits/stdc++.h>
using namespace std;

int n,a[1000010],c[1000010];

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;
}


int main(){
	int t;
	cin>>t;
	for(int tot=1;tot<=t;tot++){
		cout<<"Case "<<tot<<":"<<endl;
		memset(a,0,sizeof a);
		memset(c,0,sizeof c);
		cin>>n;
		for(int i=1;i<=n;i++) {
			cin>>a[i];
			updata(i,a[i]);
		}
		string s;
		int x,y;
		while(cin>>s&&s[0]!='E'){
			cin>>x>>y;
			if(s[0]=='Q'){
				int sum=getsum(y)-getsum(x-1);
				cout<<sum<<endl;
			}
			else if(s[0]=='S'){
				updata(x,-y);
			}
			else if(s[0]=='A'){
				updata(x,y);
			}
		}
	}
	
}

我的板子好像有问题,以后都不敢用了,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);
		}
	}
}

区间更新,区间查询

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给 a 1 a_1 a1加上一个数,3给 a 1 a_{1} a1减一个数,4求一个区间的和,5求 a 1 a_{1} a1的值。

代码:

#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]);
		updata(a[i],1);
	}
	pr(ans);
	return 0;
}

洛谷P1637 三元上升子序列

题意:Erwin 最近对一种叫 thair 的东西巨感兴趣。。。在含有 n 个整数的序列$ a_1,a_2,…,a_n ​ ​ 中 , 三 个 数 被 称 作 ‘ t h a i r ‘ 当 且 仅 当 ​​中,三个数被称作`thair`当且仅当 thairi<j<k$ 且 a i < a j < a k a_i<a_j<a_k ai<aj<ak;求一个序列中 thair 的个数。

输入

4
2 1 3 4
5
1 2 2 3 4    

输出

2
7

DP+动态规划+离散化

  • 假设 f [ i ] [ j ] f[i][j] f[i][j]表示以 a [ j ] a[j] a[j]​​结尾的长度为 i 的上升子序列的个数。

  • 动态转移方程 f [ i ] [ j ] = ∑ k j f [ i − 1 ] [ k ] f[i][j]=\sum_k^j f[i-1][k] f[i][j]=kjf[i1][k] , a [ k ] < a [ j ] a[k]<a[j] a[k]<a[j];

  • 要满足两个性质:

    • k < j k<j k<j;
    • a [ k ] < a [ j ] ; a[k]<a[j]; a[k]<a[j];

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

代码:

#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,求出以每一个 a i a_i ai结尾的最长上升子序列的长度。

树状数组+离散化

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;
    }
};

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

万伏小太阳

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值