bzoj5157 [Tjoi2014]上升子序列（dp，树状数组优化）

设dp[i]为以a[i]结尾的上升子序列个数，则显然有

$dp[i]=\sum\limits_{j=1}^{i-1}dp[j]*[a[j]<a[i]]+1$。

由于要求本质不同，所以如果a[i]还在前面出现过，要把那个位置的dp值设为0，避免重复。我们显然可以以权值建树状数组来优化一下。复杂度$O(nlogn)$

答案就是所有的dp值的和，因为要求长度>=2，所以我们最后还要减掉长度为1的。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
#define N 100010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,a[N],aa[N],m,c[N],last[N];
inline void add(int x,int val){
    for(;x<=m;x+=x&-x) c[x]=(c[x]+val)%mod;
}
inline int ask(int x){
    int res=0;for(;x;x-=x&-x) res=(res+c[x])%mod;return res;
}
int main(){
    n=read();
    for(int i=1;i<=n;++i) a[i]=aa[i]=read();
    sort(aa+1,aa+n+1);m=unique(aa+1,aa+n+1)-aa-1;
    for(int i=1;i<=n;++i){
        int x=lower_bound(aa+1,aa+m+1,a[i])-aa;
        int res=ask(x-1);
        add(x,(res+1-last[x]+mod)%mod);last[x]=(res+1)%mod;
    }printf("%d\n",ask(m)-m);
    return 0;
}