$$OI中的数据结构$$
$$By;TYQ$$
线性结构
大部略
单调栈
栈 , 但是push的时候要弹出所有比他小/大的(多用于优化Dp)
单调队列
队列 , 同单调栈
树状结构
树状数组
核心:lowbit(x) = (x) & (-x)
...其实lowbit(x) = 2^x的最低非0位
PION8012初赛中考了...但只涉及正数...
- 为什么lowbit(x) = (x) & (-x)
考虑x>0时-x等于多少:-x在二进制中的意义为x所有位取反后+1 , 那么他的第一个非0位以前的都是0 , &后结果为0 , 在第一个非0位时-x发生进位使哪位为1 , 那么这位为1 , 再以后呢?为0 , 所以&也为0
- 树状数组每个节点的值
(C_{i} = sum_{j=i}^{j!=0 , j=lowbit(j)} C_{j})
- 树状数组怎么求和 :
我们为什么要设计前面的(C_{i})为这个奇怪的数?
答案揭晓!为了查询!当查询时我们只需要遍历查询数的每个lowbit , 将值加上(C_{x})
查询时间复杂度为(O(logN))
那么这为什么对呢?因为:
.
.
.
画出这颗树发现的确是对的QAQ其实是我不会证
关于代码:
void modify(int x , int y)/*modify a[x] to a[x]+y*/{for(int i=x;i<=n;C[i]+=y,i+=(i)&(-i)) ; }
int query(int x)/*query a[1]+...+a[x]*/{int ret = 0 ;for(int i=x;i;i-=(i)&(-i))ret+=C[i] ; return ret ;}
线段树
其实一开始你觉得比较难 , 但其实很基础
线段树每次将区间分成两个小区间 , 到底层递归 ;
用指针写伪代码就是
Query()
int sum = 0 ;
if(!this->ls) sum+=this->ls.Query()
if(!this->rs)sum+=this->rs.Query()
return sum ;
//什么?你问我为什么要写指针?代码短!
但是你代码肯定不能这么写...线段树维护的不只和 , 需要zici区间查询 :
关于区间查询
Query(i,j,ni,nj,p)
if(i<=ni && j<=nj) return Value[p] ;
int m = (s + t) >> 2, sum = 0;
if (l <= m) sum += Query(l, r, s, m, p * 2);
if (r > m) sum += Query(l, r, m + 1, t, p * 2 + 1);
return sum;
pushdown
pushdown(p,m,s,t)
d[p * 2] += b[p] * (m - s + 1), d[p * 2 + 1] += b[p] * (t - m),b[p * 2] += b[p], b[p * 2 + 1] += b[p];
区间修改:
每次执行单点修改并给节点打上LazyTag , 需要时下放
板子
#include<iostream>
#include<cstdio>
#define ll long long
#define MAXN 100000
int p , arr[MAXN+5];
using namespace std ;
struct node{
ll v,mul,add;
}tr[4*MAXN+5];
void build(int root,int l,int r){
tr[root].mul=1,tr[root].add=0;
if(l==r) tr[root].v=arr[l];
else {
int m=(l+r)/2;
build(root*2,l,m) ;
build(root*2+1,m+1,r) ;
tr[root].v=tr[root*2].v+tr[root*2+1].v ;
}
tr[root].v%=p ;
}
void pushdown(int root,int l,int r){
int mid=(l+r)/2 ;
tr[root*2].v=(tr[root*2].v*tr[root].mul+tr[root].add*(mid-l+1))%p,
tr[root*2+1].v=(tr[root*2+1].v*tr[root].mul+tr[root].add*(r-mid))%p,
tr[root*2].mul=(tr[root*2].mul*tr[root].mul)%p,
tr[root*2+1].mul=(tr[root*2+1].mul*tr[root].mul)%p,
tr[root*2].add=(tr[root*2].add*tr[root].mul+tr[root].add)%p,
tr[root*2+1].add=(tr[root*2+1].add*tr[root].mul+tr[root].add)%p,
tr[root].mul=1,tr[root].add=0;
}
void U_M(int root,int il,int ir,int l,int r,ll k){
if(r<il||ir<l) return ;
if(l<=il&&ir<=r) {
tr[root].v=(tr[root].v*k)%p,
tr[root].mul=(tr[root].mul*k)%p,
tr[root].add=(tr[root].add*k)%p;
return ;
}
pushdown(root,il,ir) ;
int mid=(il+ir)/2 ;
U_M(root*2,il,mid,l,r,k) ;
U_M(root*2+1,mid+1,ir,l,r,k) ;
tr[root].v=(tr[root*2].v+tr[root*2+1].v)%p ;
return ;
}
void U_A(int root,int il,int ir,int l,int r,ll k){
if(r<il||ir<l) return ;
if(l<=il&&ir<=r) {
tr[root].add=(tr[root].add+k)%p,
tr[root].v=(tr[root].v+k*(ir-il+1))%p;
return ;
}
pushdown(root,il,ir) ;
int mid=(il+ir)/2 ;
U_A(root*2,il,mid,l,r,k) ;
U_A(root*2+1,mid+1,ir,l,r,k) ;
tr[root].v=(tr[root*2].v+tr[root*2+1].v)%p ;
return ;
}
ll query(int root,int il,int ir,int l,int r){
if(r<il||ir<l) return 0;
if(l<=il&&ir<=r) {return tr[root].v;}
pushdown(root,il,ir) ;
int mid=(il+ir)/2 ;
return (query(root*2,il,mid,l,r)+query(root*2+1,mid+1,ir,l,r))%p;
}
int main(){
int n,m ;
scanf("%d%d%d",&n,&m,&p);
for(int i=1;i<=n;++i) scanf("%d",&arr[i]) ; build(1,1,n) ;
//for(int i=1;i<=n;++i) cout<<tr[i].v<<" " ; cout<<endl ;
while(m--){
int x,y,z; ll k;
scanf("%d",&x) ;
if(x==1)scanf("%d%d%lld",&y,&z,&k) , U_M(1,1,n,y,z,k) ;
else if(x==2)scanf("%d%d%lld",&y,&z,&k) , U_A(1,1,n,y,z,k) ;
else scanf("%d%d",&y,&z),printf("%lld
",query(1,1,n,y,z)) ;
//for(int i=1;i<=n;++i) cout<<tr[i].v<<" " ; cout<<endl ;
}
return 0 ;
}