HDU 4991 Ordered Subsequence （树状数组+离散化）水

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4991

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define read(x,y) scanf("%d%d",&x,&y)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
const int  maxn =2e4+5;
const int mod=123456789;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意：给定一个序列求长度为m的递增子序列的数量。

典型DP+树状数组。
外循环m次，每次把每个位置的长度为m-1的DP值丢到树状数组中，
并查询更新，注意顺序。

*/

int n,m;
ll a[maxn],b[maxn];
int dp[maxn];
///树状数组
int bit[maxn<<1];
int lowbit(int x){return x&(-x);}
void add(int x,int d){for(;x<maxn;(bit[x]+=d)%=mod,x+=lowbit(x));}
int sum(int x){int ret=0;for(;x>0;(ret+=bit[x])%=mod,x-=lowbit(x));return ret;}

int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=0;i<n;i++) dp[i]=1;
        for(int i=0;i<n;i++) scanf("%lld",&a[i]),b[i]=a[i];

        sort(b,b+n);for(int i=0;i<n;i++) a[i]=lower_bound(b,b+n,a[i])-b+1;///离散化

        for(int i=1;i<m;i++)
        {
            memset(bit,0,sizeof(bit));
            for(int j=0;j<n;j++)
            {
                add(a[j]+1,dp[j]);
                dp[j]=sum(a[j]);
            }
        }
        int ans=0;for(int i=0;i<n;i++) (ans+=dp[i])%=mod;
        printf("%d\n",ans);
    }
    return 0;
}