HDU5542 The Battle of Chibi （DP）

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

题目大意：给出n个数让你求长度恰好为m的递增子序列的方案数。

解题思路：dp[i][j]表示第i个位置取长度为j，且子序列最后一个是a[i]的方案数。

因为n最多1e3个，但大小范围是1e9，所以先离散化。

然后利用树状数组来找到长度为j-1的并且最后一个数a[k]<a[i]。

然后还有些细节方面需要优化，不然会TLE。我的也就3000+ms卡过。

/* ***********************************************
┆  ┏┓　　　┏┓ ┆
┆┏┛┻━━━┛┻┓ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　━　　　┃ ┆
┆┃　┳┛　┗┳　┃ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　┻　　　┃ ┆
┆┗━┓　马　┏━┛ ┆
┆　　┃　勒　┃　　┆　　　　　　
┆　　┃　戈　┗━━━┓ ┆
┆　　┃　壁　　　　　┣┓┆
┆　　┃　的草泥马　　┏┛┆
┆　　┗┓┓┏━┳┓┏┛ ┆
┆　　　┃┫┫　┃┫┫ ┆
┆　　　┗┻┛　┗┻┛ ┆
************************************************ */
//#pragma comment(linker, "/STACK:102400000,102400000")
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <bitset>
using namespace std;

#define rep(i,a,b) for (int i=(a),_ed=(b);i<=_ed;i++)
#define per(i,a,b) for (int i=(b),_ed=(a);i>=_ed;i--)
#define pb push_back
#define mp make_pair
const int inf_int = 2e9;
const long long inf_ll = 2e18;
#define inf_add 0x3f3f3f3f
#define mod 1000000007
#define LL long long
#define ULL unsigned long long
#define MS0(X) memset((X), 0, sizeof((X)))
#define SelfType int
SelfType Gcd(SelfType p,SelfType q){return q==0?p:Gcd(q,p%q);}
SelfType Pow(SelfType p,SelfType q){SelfType ans=1;while(q){if(q&1)ans=ans*p;p=p*p;q>>=1;}return ans;}
#define Sd(X) int (X); scanf("%d", &X)
#define Sdd(X, Y) int X, Y; scanf("%d%d", &X, &Y)
#define Sddd(X, Y, Z) int X, Y, Z; scanf("%d%d%d", &X, &Y, &Z)
#define reunique(v) v.resize(std::unique(v.begin(), v.end()) - v.begin())
#define all(a) a.begin(), a.end()
#define   mem(x,v)      memset(x,v,sizeof(x))
typedef pair<int, int> pii;
typedef pair<long long, long long> pll;
typedef vector<int> vi;
typedef vector<long long> vll;
inline int read(){int ra,fh;char rx;rx=getchar(),ra=0,fh=1;while((rx<'0'||rx>'9')&&rx!='-')rx=getchar();if(rx=='-')fh=-1,rx=getchar();while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra*fh;}

const int N = 1005;
int dp[N][N],a[N],t[N][N];
int top;

inline int lowbit(int x) {return x&-x;}

void add(int i,int j,int c)
{
    while(j<top)
    {
        t[i][j] += c;
        if(t[i][j]>=mod) t[i][j] %= mod;
        j += lowbit(j);
    }
}

int sum(int i,int j)
{
    int res = 0;
    while(j)
    {
        res += t[i][j];
        if(res>=mod) res %= mod;
        j -= lowbit(j);
    }
    return res;
}





int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	//ios::sync_with_stdio(0);
	//cin.tie(0);
	int tt = read();
	int cas = 1;
	while(tt--)
    {
        int n = read(),m = read();
        MS0(dp);
        MS0(t);
        set<int>s;
        map<int,int>mp;
        for(int i=1;i<=n;i++)
        {
            a[i] = read();
            s.insert(a[i]);
        }
        top = 1;
        for(set<int>::iterator it=s.begin();it!=s.end();it++)
        {
            mp[*it] = top++;
        }
        for(int i=1;i<=n;i++)
        {
            int p = min(i,m);
            for(int j=1;j<=p;j++)
            {
                if(j==1) dp[i][j] = 1;
                else dp[i][j] += sum(j-1,mp[a[i]]-1);
                if(dp[i][j]>=mod) dp[i][j] %= mod;
                add(j,mp[a[i]],dp[i][j]);
            }
        }
        int ans = 0;
        for(int i=m;i<=n;i++)
        {
            ans += dp[i][m];
            if(ans>=mod) ans %= mod;
        }
        printf("Case #%d: %d\n",cas++,ans);
    }


	return 0;
}