树状数组专题

基本原理

由于线段树上的任意一个区间(根除外)均可以用父区间减去另一个子区间得到
所以只存左区间也可以维护数的所以信息


如果想求[1,7]的和
由于7=>0111
所以sum[0111]=a[0111]+a[0110]+a[0100]
即：a[7]+a[6]+a[4];

基本工作函数

void add(LL *A,int i,int x){
    while(i<=n)=A[i]+=x,i+=i&-i;
}

LL sum(LL *A,int i){
    LL s=0;
    while(i>0)
    s+=A[i],i-=i&-i;
    return s;
}

1.单点更新，区间查询

【1】把A[i]加x,add(sum_A,i,x);
【2】sum(sum_A,r)-sum(sum_A,l-1);

2.单点查询，区间更新

【1】求A[i],sum(sum_A,i);
【2】把[l,r]加x add(sum_A,l-1,-x),add(sum_A,r,x);

3.区间查询，区间更新

利用A,B bit两个来维护
B[i]=∆A[i]*i;
【1】求[l,r]

l--;
ans=sum(A,r)+sum(B,r)*r//[1,r]
    -sum(A,l)+sum(B,l)*l;[1,l-1]

【2】把[l,r]加x

add(A,l,-x*(l-1));
add(B,l,x);
add(A,r+1,x*r);
add(B,r+1,-x);

例题：

1.MIS

时间限制: 1 Sec 内存限制: 128 MB
提交: 74 解决: 34
[提交][状态][讨论版]
题目描述
长度为m的递增子序列（M length Increase Subsequence），称为MIS。

求长度为n的序列，有多少个MIS。

子序列不需要在原序列中不需要连续，但要保证相对位置。

输入
第一行两个整数n和m

第二行n个整数，表示一个序列。

输出
输出MIS的个数,结果对20140921取余。

样例输入

5 3
1 2 5 3 4

样例输出

5

样例解释：[1 2 5],[1 2 3],[1 2 4],[1 3 4],[2 3 4] 都是MIS。

对于30%的数据，n的范围[1,30],m的范围[1,10];

对于60%的数据，n的范围[1,500],m的的范围[1,50];

对于80%的数据，n的范围[1,1000],m的的范围[1,100]，序列中的元素范围[1,n];

对于100%的数据，n的范围[1,10000],m的的范围[1,100]，序列中的元素范围[1,109];

思路：

对于60%的数据
很容想到dp[i][j]即为sum{dp[1,i-1][j-1];(满足a[x]

for(int i=1;i<=n;i++)dp[i][1]=1;
    for(int i=1;i<=n;i++)
    for(int j=2;j<=m;j++){
        for(int k=1;k<=i-1;k++)if(a[k]<a[i])
        dp[i][j]+=dp[k][j-1];
        if(j==m)ans+=dp[i][j];
    }

所以每次维护sum[j-1,a[i]]表示m=j-1小于等于a[i]的所以dp和
于是我们需要一个数据结构满足单点更新区间查询的功能

void add(LL* A,int i,int x){
    while(i<=n)A[i]+=x,i+=i&-i;
}
int sum(LL* A,int i){
    LL s=0;
    while(i>0)s+=A[i],i-=i&-i;
    return s%P;
}   
int solve(){
    //先离散
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++){
        if(j==1)add(A[1],a[i],1);//dp[1][i]=1;
        else{
            s=sum(A[j-1],a[i]-1);//s=dp[j][i];
            add(A[j],a[i],s);//存下dp[j][i];
        }
    }return sum(A[m],n);
    //因为离散过n即数组中的最大数
}

奶牛来拍照

时间限制: 1 Sec 内存限制: 128 MB

题目描述

农夫约翰很无聊打算去帮奶牛拍照。
由于奶牛很懒它们动也不想动，于是照片就很丑了。
聪明的约翰发现，如果一只奶牛满足下面的条件那么它就会让照片变丑。
即：比左边比它高的奶牛数量li和右边比它高的奶牛数量ri中较大的数比小的数两倍大。
约翰希望我们帮帮它求出有几头奶牛让照片变丑。

输入

第一行 n

之后的n行，每行一个hi 表示奶牛的高度

输出

一个数字，使照片变丑的奶牛数量。

样例输入

7
34
6
23
0
5
99
2

样例输出

3

提示

（1<=n<=100,000,hi<1,000,000,000）

思路：
<1>因为我们要求每个数的比其高的数的数量，首先我们先对h进行排序，之后从大到小来算
<2>我们发现我们现在需要一个数据结构有以下的功能：
1.每次我们可以往里面放一个元素（从大到小放入h的坐标）
2.我们可以快速查寻出里面比它小或大的元素有几个（即下表小的，$L_i$，$R_{i }$）
<3>我们发现我们可以通过树状数组来实现，我们可以通过query（i）表示比 i <script type="math/tex" id="MathJax-Element-3">i</script>小的元素个数

#include<bits/stdc++.h>
using namespace std;
const int M=1e5+5;
int n,num,ans,C[M],sum;
struct node{
    int h,id;
    bool operator<(const node&A)const{return h>A.h;}
    void rd(){scanf("%d",&h);id=++num;}
}A[M];
inline int query(int i){
    int x=0;
    while(i<=n)x+=C[i],i+=i&-i;
    return x;
}
inline void add(int i){
    while(i)C[i]++,i-=i&-i;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)A[i].rd();
    sort(A+1,A+1+n);
    for(int i=1;i<=n;i++){
        int r=query(A[i].id);
        int l=sum-r;
        printf("l=%d r=%d\n",l,r);
        if(l>r*2||r>l*2)ans++;
        add(A[i].id);sum++;
    }printf("%d",ans);
}