The Battle of Chibi HDU - 5542

https://cn.vjudge.net/problem/HDU-5542

题意：求从长度为n的序列中最多能找出长度为k的严格递增子序列的个数。

思路：dp[i][j],表示第i个数

作为子序列的第j个数的序列有多少个。 

   for (int i = 0; i < n; i++)
        {
            dp[i][1] = 1;
        }
        for (int i = 1; i < n; i++)
        {
            for (int j = 1; j <= min(i+1,m); j++)
            {
                for (int k = i - 1; k >= 0; k--)
                {
                    if (num[k] < num[i])
                        dp[i][j] = (dp[i][j] + dp[k][j - 1]) % MOD;
                }
            }
        }

 会超时，第三重for用树状数组可以优化，二分离散化。

#include <iostream>
#include<algorithm>
#include<queue>
#include<stdio.h>
#include<string.h>
using namespace std;
#define mod 1000000007
int lowbit(int x)
{
    return x&(-x);
}
int n,m,top;
int sum[1003],b[1003];
int dp[1003][1003],a[1003];
void add(int x,int d)
{
    while(x<=n)
    {
        a[x]+=d;
        a[x]%=mod;
        x+=lowbit(x);
    }
}
int queuy(int x)
{
    int ans=0;
    while(x>0)
    {
        ans+=a[x];
        ans%=mod;
        x-=lowbit(x);
    }
    return ans;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&b[i]);
            sum[i]=b[i];
        }
        sort(b+1,b+1+n);
        for(int i=1; i<=n; i++)
            sum[i]=lower_bound(b+1,b+n+1,sum[i])-b;
        for(int i=1; i<=n; i++)
            dp[i][1]=1;
        for(int k=2; k<=m; k++)
        {
            memset(a,0,sizeof(a));
            for(int i=1; i<=n; i++)
            {
                dp[i][k]=queuy(sum[i]-1);
                add(sum[i],dp[i][k-1]);
            }
        }
        int ans=0;
        for(int i=1; i<=n; i++)
        {
            ans+=dp[i][m];
            ans%=mod;
        }
        while(ans<0)
            ans+=mod;
        printf("Case #%d: %d\n",++top,ans);
    }
    return 0;
}