BZOJ3591: 最长上升子序列 状压DP

Description
给出1~n的一个排列的一个最长上升子序列,求原排列可能的种类数。


Sample Input
5
3
1 3 4


Sample Output
11


挺好的DP题。。。
我们先考虑普通LIS,我们维护一个序列,每一次尽量去替换序列中的值。
对于这个,可以用一个三进制维护,0表示未进队列,1表示进了队列且在序列中,2表示进了队列但已被踢出队列。
考虑题目的限制条件,那么我们只要每次转移时使他前面的数都填满了且最长上升子序列不超过给定长度即可。


#include 
#include 

using namespace std;

bool v[16];
int tp, sta[16], pos[16];
int a[16], d[16], f[14348907];

int main() {
	int n, k; scanf("%d%d", &n, &k);
	for(int i = 1; i <= k; i++) {
		scanf("%d", &a[i]); pos[a[i]] = i;
		if(a[i - 1] > a[i]) {printf("0\n"); return 0;}
	}
	d[0] = 1; for(int i = 1; i <= n; i++) d[i] = d[i - 1] * 3;
	f[0] = 1; int ans =0;
	for(int i = 0; i < d[n]; i++) {
		if(!f[i]) continue; tp = 0;
		int u = i, y = 0; memset(v, 0, sizeof(v));
		for(int j = 1; j <= n; j++) {
			int x = u % 3; u /= 3;
			if(!x) v[j] = 1; else y++;
			if(x == 1) sta[++tp] = j;
		} if(y == n) {ans += f[i]; continue;}
		int hh = 1;
		for(int j = 1; j <= n; j++) if(v[j]){
			if(pos[j] > 1 && v[a[pos[j] - 1]]) continue;
			int yy = i;
			while(j > sta[hh] && hh <= tp) hh++;
			if(hh == k + 1) continue;
			yy += d[j - 1];
			if(tp >= hh) yy += d[sta[hh] - 1];
			f[yy] += f[i];
		}
	} printf("%d\n", ans);
	return 0;
}

你可能感兴趣的:(状压DP)