HDU 4352 XHXJ's LIS

数位DP，以dp[index][sta][k]来记录搜索到当前层时，最长序列的状态为sta，要求的是长度为k的序列的个数。状态sta中1的个数即是当前这个数的最长上子升序列的长度，

这个状态的求法和nlogn求最长上升子序列的求法是一样的。

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <map>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <cstdlib>
#include <algorithm>
using namespace std;
typedef __int64 int64;
typedef long long ll;
#define M 100005
#define N 1000005
#define max_inf 0x7f7f7f7f
#define min_inf 0x80808080
#define mod 1000000007

int k;
int dig[30];
int64 dp[30][1<<11][11];

int Getlen (int sta)//求出最长上升子序列的长度
{
	int ret = 0;
	while (sta)
	{
		ret += (sta&1);
		sta >>= 1;
	}
	return ret;
}

int Getnew(int x , int sta)//更新当前数的最长上升子序列的状态
{
	int i;
	for (i = x ; i <= 9 ; i++)
		if ((sta & (1<<i)))
			return (sta-(1<<i))|(1<<x);
	return sta|(1<<x);
}

int64 Dfs(int index , int fg , int sta , int lim)//fg标记前面是否出现过非零位
{
	if (!index)return Getlen(sta) == k;
	if (!lim && dp[index][sta][k] != -1)return dp[index][sta][k];
	int i , up = lim ? dig[index] : 9;
	int64 ret = 0;
	for (i = 0 ; i <= up ; i++)
	{
		ret += Dfs(index-1 , fg||i!=0 , (!fg&&i==0) ? 0 : Getnew(i,sta) , lim&&i==up);
	}
	if (!lim)dp[index][sta][k] = ret;
	return ret;

}

int64 Solve(int64 num)
{
	int len = 0;
	while (num)
	{
		dig[++len] = num%10;
		num /= 10;
	}
	return Dfs(len , 0 , 0 , 1);
}

int main()
{
	int t , tcase = 1;
	int64 l , r;
	scanf("%d",&t);
	memset(dp , -1 , sizeof dp);
	while (t--)
	{
		scanf("%I64d%I64d%d",&l,&r,&k);
		printf("Case #%d: %I64d\n",tcase++,Solve(r)-Solve(l-1));
	}
	return  0;
}