【数位dp && 状态压缩】HDU - 4352 XHXJ's LIS

Problem Description

给你区间范围，给你一个数k。问你区间范围内的数，有多少个满足它的最长上升子序列的个数 等于 k。

思路：

前置技能：O(nlogn)求最长上升子序列的个数。
因为元素不重复，所以最长上升子序列只能由0-9构成，也就是10位，所以可以状态压缩枚举 求最长上升子序列的个数。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll dp[25][1200][15];
int a[25], k;
int Count(int now)//求最长上升子序列的长度
{
    int res = 0;
    for(int i = 0; i < 10; i++)
    {
        if((1<<i) & now) res++;
    }
    return res;
}
int get_now(int now, int x)//插入一个元素后的最长上升子序列的状态（这个方法只适合求个数，不适合求最长上升子序列）
{
    for(int i = x; i < 10; i++)
    {
        if((1<<i) & now) return ((1<<i)^now) | (1<<x);
    }
    return (1<<x)|now;
}
ll dfs(int pos, int now, int sta, int limit)
{
    if(pos < 0) return Count(now) == k;
    if(Count(now) > k) return 0;
    if(!limit && dp[pos][now][k] != -1) return dp[pos][now][k];
    int up = limit ? a[pos] : 9;
    ll res = 0;
    for(int i = 0; i <= up; i++)
    {
        res += dfs(pos-1, (sta || i != 0) ? get_now(now, i):0, sta || i != 0, limit && i == a[pos]);
    }
    if(!limit) return dp[pos][now][k] = res;
    else return res;
}
ll solve(ll x)
{
    int pos = 0;
    while(x)
    {
        a[pos++] = x%10;
        x = x/10;
    }
    return dfs(pos-1, 0, 0, 1);
}
int main()
{
    int T, Case = 1;
    cin >> T;
    ll lt, rt;
    memset(dp, -1, sizeof(dp));
    while(T--)
    {
        cin >> lt >> rt >> k;
        cout << "Case #" << Case++ << ": "<< solve(rt) - solve(lt-1) << endl;
    }
    return 0;
}