树状数组
树状数组的学习可以看b站董晓算法的讲解(极力推荐)。
大概的思路
无论是往点修往后跳还是求前缀和往前跳都是一次跳2k,k为x二进制最低有效位。
代码模版
template<typename T>
struct Fenwick{
int n;
vector<T> tr;
Fenwick(int n) : n(n), tr(n + 1, 0){}
int lowbit(int x){// 找到x最低有效位
return x & -x;
}
void modify(int x, T c){// 点修
for(int i = x; i <= n; i += lowbit(i)) tr[i] += c;
}
void modify(int l, int r, T c){// 区修
modify(l, c);
if (r + 1 <= n) modify(r + 1, -c);
}
T query(int x){// 前缀和
T res = T();
for(int i = x; i; i -= lowbit(i)) res += tr[i];
return res;
}
T query(int l, int r){// 区间查询
return query(r) - query(l - 1);
}
int find_first(T sum){// 查找第一个满足条件的下标(前缀和
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] < sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans + 1;
}
int find_last(T sum){// 查找最后一个满足条件的下标(前缀和)
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] <= sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans;
}
};
using BIT = Fenwick<int>;
点修+区查
P3374 【模板】树状数组 1
https://www.luogu.com.cn/problem/P3374
点修+区查是操作最简单的,直接用要维持的数据来建树,然后单点修改的话就直接对这个点进行modify,区间查询的话就对这个区间query(l,r);
#include<bits/stdc++.h>
using namespace std;
using u32 = unsigned;
#define i128 __int128;
using ll = long long;
//#define int ll
using u64 = unsigned long long;
const ll inf = 1e9;
const ll INF = 1e18;
template<typename T>
struct Fenwick{
int n;
vector<T> tr;
Fenwick(int n) : n(n), tr(n + 1, 0){}
int lowbit(int x){// 找一个数的最低有效位
return x & -x;
}
void modify(int x, T c){// 点修
for(int i = x; i <= n; i += lowbit(i)) tr[i] += c;
}
void modify(int l, int r, T c){// 区修
modify(l, c);
if (r + 1 <= n) modify(r + 1, -c);
}
T query(int x){// 查询前缀和
T res = T();
for(int i = x; i; i -= lowbit(i)) res += tr[i];
return res;
}
T query(int l, int r){// 查询区间和
return query(r) - query(l - 1);
}
int find_first(T sum){
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] < sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans + 1;
}
int find_last(T sum){
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] <= sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans;
}
};
using BIT = Fenwick<int>;
signed main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
int n,m;
cin>>n>>m;
vector<int>a(n+1);
BIT bit(n);
for(int i=1;i<=n;i++){
int x;
cin>>x;
bit.modify(i,x);
}
for(int i=1;i<=m;i++)
{
int op,x,y;
cin>>op>>x>>y;
if(op==1)
{
bit.modify(x,y);
}
else
{
cout<<bit.query(x,y)<<'\n';
}
}
return 0;
}
区修+点查
树状数组初始值都为0,这是一个差分数组,如果对[l,r]加上某一个数,就modify(l,k),modify(r+1,-k),然后如果要点查的话就直接用a[x]+bit.query(x)就行了
#include<bits/stdc++.h>
using namespace std;
using u32 = unsigned;
#define i128 __int128;
using ll = long long;
#define int ll
using u64 = unsigned long long;
const ll inf = 1e9;
const ll INF = 1e18;
template<typename T>
struct Fenwick{
int n;
vector<T> tr;
Fenwick(int n) : n(n), tr(n + 1, 0){}
int lowbit(int x){
return x & -x;
}
void modify(int x, T c){
for(int i = x; i <= n; i += lowbit(i)) tr[i] += c;
}
void modify(int l, int r, T c){
modify(l, c);
if (r + 1 <= n) modify(r + 1, -c);
}
T query(int x){
T res = T();
for(int i = x; i; i -= lowbit(i)) res += tr[i];
return res;
}
T query(int l, int r){
return query(r) - query(l - 1);
}
int find_first(T sum){
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] < sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans + 1;
}
int find_last(T sum){
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] <= sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans;
}
};
using BIT = Fenwick<int>;
signed main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
int n,m;
cin>>n>>m;
vector<int>a(n+1);
for(int i=1;i<=n;i++)cin>>a[i];
BIT bit(n+1);
for(int i=1;i<=m;i++)
{
int op,x;
cin>>op>>x;
if(op==1)
{
int y,k;
cin>>y>>k;
bit.modify(x,k);
bit.modify(y+1,-k);
}
else
{
cout<<a[x]+bit.query(x)<<'\n';
}
}
return 0;
}
区修+区查
这个是最难的。
P3372 【模板】线段树 1
这一道题的区修就是区间加,区查就是区间和。
#include<bits/stdc++.h>
using namespace std;
using u32 = unsigned;
#define i128 __int128;
using ll = long long;
#define int ll
using u64 = unsigned long long;
const ll inf = 1e9;
const ll INF = 1e18;
template<typename T>
struct Fenwick{
int n;
vector<T> tr;
Fenwick(int n) : n(n), tr(n + 1, 0){}
int lowbit(int x){// 找到x最低有效位
return x & -x;
}
void modify(int x, T c){// 点修
for(int i = x; i <= n; i += lowbit(i)) tr[i] += c;
}
void modify(int l, int r, T c){// 区修
modify(l, c);
if (r + 1 <= n) modify(r + 1, -c);
}
T query(int x){// 前缀和
T res = T();
for(int i = x; i; i -= lowbit(i)) res += tr[i];
return res;
}
T query(int l, int r){// 区间查询
return query(r) - query(l - 1);
}
int find_first(T sum){// 查找第一个满足条件的下标(前缀和
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] < sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans + 1;
}
int find_last(T sum){// 查找最后一个满足条件的下标(前缀和)
int ans = 0; T val = 0;
for(int i = __lg(n); i >= 0; i--){
if ((ans | (1 << i)) <= n && val + tr[ans | (1 << i)] <= sum){
ans |= 1 << i;
val += tr[ans];
}
}
return ans;
}
};
using BIT = Fenwick<int>;
signed main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
int n,m;
cin>>n>>m;
BIT d1(n),d2(n);
vector<int>a(n+1);
for(int i=1;i<=n;i++)cin>>a[i];
for(int i=1;i<=n;i++)
{
d1.modify(i,a[i]);
d1.modify(i+1,-a[i]);
d2.modify(i,i*a[i]);
d2.modify(i+1,-(i+1)*a[i]);
}
for(int i=1;i<=m;i++)
{
int op,x,y;
cin>>op>>x>>y;
if(op==1)
{
int k;
cin>>k;
d1.modify(x,k);
d1.modify(y+1,-k);
d2.modify(x,k*x);
d2.modify(y+1,-k*(y+1));
}
else
{
cout<<(y+1ll)*d1.query(y)-1ll*x*d1.query(x-1)-(d2.query(y)-d2.query(x-1))<<'\n';
}
}
return 0;
}
in>>op>>x>>y;
if(op==1)
{
int k;
cin>>k;
d1.modify(x,k);
d1.modify(y+1,-k);
d2.modify(x,kx);
d2.modify(y+1,-k(y+1));
}
else
{
cout<<(y+1ll)d1.query(y)-1llx*d1.query(x-1)-(d2.query(y)-d2.query(x-1))<<‘\n’;
}
}
return 0;
}
标签:return,树状,--,modify,cin,int,数组,ans,query
From: https://blog.csdn.net/2301_79872524/article/details/143568141