递增子序列
题目链接:ybt金牌导航1-2-1
题目大意
给定一个长度为 n 的序列,求有多少个子序列满足长度为 m,并且递增。
思路
首先,看到数这么大,然后又只有这几个数,还是要递增,那我们直接先离散化一波。
然后我们因为递增,我们先考虑 dp。
就设
f
i
,
j
f_{i,j}
fi,j 为最后一个数是
i
i
i,长度为
j
j
j 的递增子序列有多少个。
可以想到要先枚举到第几个数
k
k
k,然后再枚举
i
i
i 和
j
j
j 进行转移,会超时。
那我们看看转移方程: f i , j = ∑ k = 1 i − 1 f k , j − 1 × ( a k < a i ) f_{i,j}=\sum\limits_{k=1}^{i-1}f_{k,j-1}\times(a_k<a_i) fi,j=k=1∑i−1fk,j−1×(ak<ai)。
我们会发现就是找顺序对,那可以用树状数组来做。
那就可以了。
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#define mo 123456789
using namespace std;
struct node {
int x, num;
}a[10001], b[10001];
int n, m, tot;
long long tree[101][100001], re, ans;
bool cmp1(node x, node y) {
return x.x < y.x;
}
bool cmp2(node x, node y) {
return x.num < y.num;
}
long long ask(int x, int y) {//树状数组
re = 0;
for (int i = x; i; i -= i & (-i))
re = (tree[y][i] + re) % mo;
return re;
}
void build(int x, int y, long long z) {
for (int i = x; i <= 10000; i += i & (-i))
tree[y][i] = (tree[y][i] + z) % mo;
}
int main() {
while (scanf("%d %d", &n, &m) != EOF) {
tot = 0;
ans = 0;
memset(tree, 0, sizeof(tree));
for (int i = 1; i <= n; i++) {//离散化
scanf("%d", &a[i].x);
a[i].x++;
a[i].num = i;
}
sort(a + 1, a + n + 1, cmp1);
for (int i = 1; i <= n; i++) {
if (a[i].x != a[i - 1].x) {
b[i].x = ++tot;
}
else {
b[i].x = tot;
}
b[i].num = a[i].num;
}
sort(b + 1, b + n + 1, cmp2);
for (int i = 1; i <= n; i++) {
build(b[i].x, 1, 1ll);//单个
for (int j = 2; j <= m; j++) {
re = ask(b[i].x - 1, j - 1);//多个,类似dp的思想
if (!re) break;
build(b[i].x, j, re);
}
}
printf("%lld\n", ask(10000, m));
}
return 0;
}