HDU4352 XHXJ's LIS(LIS 状压)

题意

题目链接

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;
}