HDU4352 XHXJ's LIS(LIS 状压)

时间:2021-06-02 16:24:34

题意

题目链接

Sol

刚开始的思路是\(f[i][j]\)表示到第\(i\)位,LIS长度为\(j\)的方案。

然而发现根本不能转移,除非知道了之前的状态然后重新dp一遍。。

题解,,,挺暴力的把,直接把求LIS过程中的单调栈当成一个状态压进去了。。

自己真是不长记性,明明已经被这个单调栈坑过一次了。。

考虑到\(k\)非常小,于是直接对\(k\)进行状压

设\(f[i][sta][j]\)表示长度为\(i\),单调栈内状态为\(sta\), LIS长度为\(k\)的方案数

最后一维如果是单组数据的话是不必要的。但是考虑到时多组数据,就一起加进来吧。

转移的时候枚举一下这一位放了什么,然后暴力的改一下LIS的状态。

#include<bits/stdc++.h>
#define LL long long
#define int long long
using namespace std;
const int MAXN = 1e5 + 10;
LL T, l, r, K;
int f[64][1 << 10][11];//长度为i,lis为k的方案数
int a[MAXN], num;
int change(int S, int x) {//向状态s中加入一个数x
for(int i = x; i <= 9; i++)
if(S & (1 << i)) {S ^= (1 << i); break;}
return S | (1 << x);
}
int dfs(int now, int lim, int lead, int s) {
// printf("%d %d %d %d\n", now, lim, lead, s);
if(!now) return (__builtin_popcount(s) == K);
if(!lim && f[now][s][K] != -1) return f[now][s][K];
int ans = 0;
for(int i = 0; i <= (lim ? a[now] : 9); i++)
ans += dfs(now - 1, lim && i == a[now], lead && (!i), (lead && (!i)) ? 0 : change(s, i));
if(!lim) f[now][s][K] = ans;
return ans;
}
LL solve(LL x) {
num = 0;
while(x) a[++num] = x % 10, x /= 10;
// cout << num << endl;
dfs(num, 1, 1, 0);
}
main() {
memset(f, -1, sizeof(f));
cin >> T;
for(int i = 1; i <= T; i++) {
cin >> l >> r >> K;
printf("Case #%d: ", i);
cout << solve(r) - solve(l - 1) << endl;
}
return 0;
}