欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  IT编程

HDU4352 XHXJ's LIS(LIS 状压)

程序员文章站 2022-04-09 18:40:23
题意 "题目链接" Sol 刚开始的思路是$f[i][j]$表示到第$i$位,LIS长度为$j$的方案。 然而发现根本不能转移,除非知道了之前的状态然后重新dp一遍。。 题解,,,挺暴力的把,直接把求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;
}