HDU-4352 XHXJ's LIS

数位DP+状态压缩

d[len][num][K]保存len位数，状态为num，且LIS长度为K的数字个数

num为最长上升子序列的二进制压缩，大小为1024

如果当前LIS为 2 5 8 9，i为3，LIS可更新为 2 3 8 9（从i位置向后寻找第一个比它大的替换，没有就直接加到后边）

注意，这里的num就不是真正的LIS压缩了，但是长度还是一样的，并且这样可以向后寻找最优解

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=20;
ll d[N][1024+10][10];
int dis[N],K;
int getnews(int x,int s)
{
    for(int i=x;i<10;i++)
        if(s&(1<<i)) return (s^(1<<i))|(1<<x);
    return s|(1<<x);
}
int getnum(int num)
{
	int t=0;
	while(num)
	{
		t+=num&1;
		num/=2;
	}
	return t;
}
ll dp(int len,int num,bool z,bool flag)
{
	if(len==-1) return getnum(num)==K;
	if(!flag&&d[len][num][K]!=-1) return d[len][num][K];
	ll ans=0;
	int end=flag?dis[len]:9;
	for(int i=0;i<=end;i++)
		ans+=dp(len-1,(z&&i==0)?0:getnews(i,num),z&&(i==0),flag&&(i==end));
	if(!flag) d[len][num][K]=ans;
	return ans;
}
ll solve(ll n)
{
	int t=0;
	while(n)
	{
		dis[t++]=n%10;
		n/=10;
	}
	int ans=0;
	return dp(t-1,0,true,true);
}
int main()
{
	memset(d,-1,sizeof(d));
	int T,kase=0;
	ll l,r;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lld%lld%d",&l,&r,&K);
		printf("Case #%d: %lld\n",++kase,solve(r)-solve(l-1));
	}	
	return 0;
}