很容易想到n^3的dp,,但是要超时,用二维树状数组优化一下枚举上一个点的第三重循环就好了。
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn = 1005;
const int mod = 1000000007;
typedef long long ll;
int n;
int num[maxn];
int tmp[maxn];
ll tree[maxn][maxn];
inline int lowbit(int i)
{
return i&(-i);
}
inline void add(int i,int j,int x)
{
while(i<=n)
{
tree[i][j] += x%mod;
tree[i][j]%=mod;
i += lowbit(i);
}
}
inline ll getsum(int i,int j)
{
ll sum = 0;
while(i)
{
sum += tree[i][j]%mod;
sum %= mod;
i -= lowbit(i);
}
return sum;
}
int main()
{
int cases,m,t=1;
scanf("%d",&cases);
while(cases--)
{
memset(tree,0,sizeof(tree));
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d",&num[i]),tmp[i] = num[i];
sort(tmp+1,tmp+1+n);
for(int i=1;i<=n;i++) num[i] = lower_bound(tmp+1,tmp+1+n,num[i])-tmp;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=min(i,m);j++)
{
if(j==1) add(num[i],1,1);
else
{
ll temp = getsum(num[i]-1,j-1);
add(num[i],j,temp);
}
// for(int k=1;k<i;k++)
// if(num[k]<num[i])
// dp[i][j] += dp[k][j-1];
}
}
printf("Case #%d: ",t++);
printf("%lld\n",getsum(n,m));
}
return 0;
}