【HDU4991】Ordered Subsequence(离散化+dp+树状数组)

记录一个菜逼的成长。。

题目链接
题目大意:
给你n个数的序列，问长度为m的严格上升子序列有多少个。答案对123456789取模
$(n\leq10000 , m\leq100)$

有一个暴力dp
$dp[i][j]$:=表示以第i个数结尾，上升子序列长度为j的个数
如果a[k] < a[i]
$dp[i][j] += dp[k][j-1]$
复杂度$O(n^2*m)$显然会T

这里要用树状数组优化。
首先，先对序列n离散化，相同的数离散化后的值要一样，这里可以用map记录。
其实就是对上升子序列长度为1~m的都建一颗树
然后$dp[i][j]$就是$\sum dp[k][j-1]$，用树状数组来维护这个和

#include <cstdio>
#include <map>
#include <cstring>
#include <algorithm>
using namespace std;
#define clr clear()
#define cl(a,b) memset(a,b,sizeof(a))
typedef long long LL;
const int maxn = 10000 + 1;
const LL mod = 123456789;
int a[maxn],b[maxn],n,m;
LL dp[maxn][101],c[101][maxn];
int lowbit(int x)
{
    return x&(-x);
}
void modify(int x,int y,int v)
{
    while(x < maxn){
        c[y][x] = (c[y][x] + v) % mod;
        x += lowbit(x);
    }
}
LL getSum(int x,int y)
{
    LL ans = 0;
    while(x > 0){
        ans = (ans + c[y][x]) % mod;
        x -= lowbit(x);
    }
    return ans;
}
map<int,int>ma;
int main()
{
    while(~scanf("%d%d",&n,&m)){
        cl(dp,0);
        cl(c,0);
        ma.clr;
        for( int i = 1; i <= n; i++ ){
            scanf("%d",&a[i]);
            b[i] = a[i];
        }
        sort(b+1,b+1+n);
        for( int i = 1; i <= n; i++ ){
            if(!ma[b[i]])ma[b[i]] = i;
        }
        LL ans = 0;
        for( int i = 1; i <= n; i++ ){
            //cout<<ma[a[i]]<<' ';
            dp[i][1] = 1;
            modify(ma[a[i]],1,1);
            for( int j = 1; j < m; j++ ){
                dp[i][j+1] = getSum(ma[a[i]]-1,j);
                modify(ma[a[i]],j+1,dp[i][j+1]);
            }
        }

        for( int i = 1; i <= n; i++ ){
            ans = (ans + dp[i][m]) % mod;
        }
        printf("%lld\n",ans);
    }
    return 0;
}