HDU 4991 Ordered Subsequence -- DP+树状数组

题意：求数组中长度为k的上升子序列的个数。

思路：dp。 dp[i][j] 表示以第i个为结尾，长度为j的上升子序列的个数。则 dp[i][j] = sigma(dp[k][j-1]) (a[k]<ai]) 。由于n和k较大，需要用树状数组优化。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<algorithm>
#include<iostream>
#include<string>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
const int N = 1005;
typedef long long ll;
int n;
int mod = 123456789,cnt;
int c[105][10005];
inline int lowbit(int t)
{
    return t&(-t);
}
void update(int x,int y,int v)
{
    for (int j=y; j<=cnt; j+=lowbit(j))
        c[x][j] =(c[x][j] +v)%mod;
}
int query(int x,int y)
{
    ll s=0;
    for (int j=y; j>0; j-=lowbit(j))
        s = (s+c[x][j])%mod;
    return s;
}

int a[10005],b[10005];
int main ()
{
    int m;
    while(~scanf("%d%d",&n,&m))
    {
        memset(c,0,sizeof(c));
        cnt = 0;
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&a[i]);
            b[cnt++] = a[i];
        }
        sort(b,b+cnt);
        cnt = unique(b,b+cnt) - b;
        for(int i=1; i<=n; i++)  a[i] = lower_bound(b,b+cnt,a[i]) - b + 1;
        for(int i=1; i<=n; i++)
        {
            for(int j = 2; j<=m; j++)
            {
                    int tot = query(j-1,a[i]-1);
                    update(j,a[i],tot);
            }
            update(1,a[i],1);
        }
        int ans = query(m,cnt) ;
        printf("%d\n",ans);
    }

    return 0;
}