HDU - 4352（数位dp+状态压缩+LIS）

题意

求区间范围内最长上升子序列数位k的数字有多少个

解释

首先由数位dp引入思路，dp转移方程为：
dfs(pos,…status.) = 累加dfs(pos-1,.stautselse…)
可备忘录记录的dp数组为 dp[pos][status][k] ,之所以要备忘录k,是为了满足不同k的访问
status记录的是开头数字到pos位置所选择的数字组合，通过二进制压缩状态为一个数字，二进制位数 i 对应1代表存在，0代表不存在。
得到stautus的方式是一种等效的LIS的方式,注意是等效，访问具体位数上的1存在可能并不是正确的，但对最长上升子序列的长度是没变的。
由此数位dp中最后判断一个数是结果方式就变成了-1位情况下取出了一个确定数，该数的最长上升子序列的长度是否等于k。其它按照数位dp正常方式来做，注意前导0，因为该题使用的LIC算法是没有考虑x=0的情况，会造成二进制0位上也出现1，导致最长上升子序列结果+1。

#include <cstdio> 
#include <cstring>
#define ll long long
using namespace std;
int k;
int dig[25];
ll dp[25][1<<10][11];
int changes(int s, int x){
	for(int i = x; i < 10; i++){
		if(s&(1<<i))
			return (s^(1<<i))|(1<<x);
	}
	return s|(1<<x);
	
} 
int getnum(int s){
	int ans = 0;
	while(s){
		if(s&1)
			ans++;
		s>>=1;
	}
	return ans;
}
ll dfs(int pos,int sta, int limit, int pre){
	if(pos == -1) return getnum(sta) == k;
	if(!limit && dp[pos][sta][k]!= -1) return dp[pos][sta][k];
	int top = limit? dig[pos]:9;
	ll ans = 0;
	for(int i = 0; i <= top; i++){
		ans += dfs(pos-1, (pre&&(i == 0))?0:changes(sta, i),limit&&(i == top), pre&&(i == 0));
	} 
	if(!limit) dp[pos][sta][k] = ans;
	return ans;
}

ll slove(ll x){
	int pos = 0;
	while(x){
		dig[pos++] = x%10;
		x /= 10; 
	}
	ll ans = dfs(pos-1,0,1,1);
	return ans;
}
int main(){
	int t;
	ll n, m; 
	memset(dp, -1, sizeof(dp));	
	scanf("%d", &t);
	int cont = 1;

	while(t--){
		scanf("%lld %lld %d", &n, &m, &k);
		printf("Case #%d: %lld\n", cont++, slove(m) - slove(n-1));
	}
	
	return 0;
}