上一题的加强版~
就是问有几个长度为m的严格上升的子序列~
一样~
排序之后往里插
然后用dp[i][j-1]的和算dp[i][j]的
dp就是一共插入i个数严格上升j个数的序列的个数~
#include <iostream>
#include <cstdio>
#include <string.h>
#include <algorithm>
using namespace std;
/*
还是树状数组应该
逆序对数
选n个数
*/
int dp[1010][1010],c[1010][1010];
const int mod = 1000000007;
int n,m;
struct node
{
long long int zhi;
int pos;
}a[1010];
int lowbit(int x)
{
return x&-x;
}
long long int sum(int x,int j)
{
long long int s=0;
while(x>0)
{
s+=c[x][j];
s%=mod;
x-=lowbit(x);
}
return s;
}
bool cmp(node a,node b)
{
if(a.zhi == b.zhi) return a.pos > b.pos;
return a.zhi < b.zhi;
}
void add(int x,int j,long long int k)
{
while(x<=n)
{
c[x][j]+=k;
c[x][j]%=mod;
x+=lowbit(x);
}
}
int main()
{
int cases;
scanf("%d",&cases);
for(int y=1;y<=cases;y++)
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%lld",&a[i].zhi);
a[i].pos = i;
}
sort(a+1,a+n+1,cmp);
memset(c,0,sizeof(c));
add(1,0,1);
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++)
{
for(int j=1;j<=i;j++)
{
dp[i][j]=sum(a[i].pos,j-1);
dp[i][j]%=mod;
add(a[i].pos+1,j,dp[i][j]);
if(dp[i][j]==0) break; //不加这句优化就会超时= =
}
}
long long int aaa=0;
for(int i=m;i<=n;i++)
{
aaa+=dp[i][m];
aaa%=mod;
}
printf("Case #%d: ",y);
printf("%lld\n",aaa);
}
return 0;
}