HDU 4352 XHXJ's LIS (数位DP+状压DP+LIS)好题*

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4352

题目大意

给定l和r，要求在这个区间内，
，满足要求的数字的个数，
要求是：数位上最长上升子序列长度为k。

题目分析 

数位DP+状压DP。
思维很脑洞，主要是用状压的思想去搞
最长上升子序列，首先观察到其最大长度不超过10 ，
且每位数字0到9，那么我们可以用1<<10来压缩状态，
我们用位标记哪些数构成了这个子序列，那么子序列的长度
可以用二进制数值上的1的个数来得出了。
假设我们已经维护出x1,x2,...xi位，其构成了最长序列，
那么新遇到x，如果这个x是最大的直接添加标记状态，
如果不是，我们用弥补的策略，把x到10上最近一次出现的位删去并把
x放进去，首先长度上答案不造成影响的，然后如果以后的最大长度是由x
产生的，那么在x到当前最大长度中肯定新添了更多的数值，这样上面的策略
也能应对，所以可行性得到了保证。。。T_T我是第一次接触这样的技巧，，
还是太弱了。。。
下面状态找好了，按照数位DP的套路，要唯一标记状态了，
dp[pos][num][k],这里新增一位k，因为我们初始化和记忆化针对的是
全部的查询，这也算是一个常用的技巧吧。
然后按照正常的套路走就可以了，注意前导零。
 

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=5e2+5;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意：
给定l和r，要求在这个区间内，
，满足要求的数字的个数，
要求是：数位上最长上升子序列长度为k。

题目分析：
数位DP+状压DP。
思维很脑洞，主要是用状压的思想去搞
最长上升子序列，首先观察到其最大长度不超过10 ，
且每位数字0到9，那么我们可以用1<<10来压缩状态，
我们用位标记哪些数构成了这个子序列，那么子序列的长度
可以用二进制数值上的1的个数来得出了。
假设我们已经维护出x1,x2,...xi位，其构成了最长序列，
那么新遇到x，如果这个x是最大的直接添加标记状态，
如果不是，我们用弥补的策略，把x到10上最近一次出现的位删去并把
x放进去，首先长度上答案不造成影响的，然后如果以后的最大长度是由x
产生的，那么在x到当前最大长度中肯定新添了更多的数值，这样上面的策略
也能应对，所以可行性得到了保证。。。T_T我是第一次接触这样的技巧，，
还是太弱了。。。
下面状态找好了，按照数位DP的套路，要唯一标记状态了，
dp[pos][num][k],这里新增一位k，因为我们初始化和记忆化针对的是
全部的查询，这也算是一个常用的技巧吧。
然后按照正常的套路走就可以了，注意前导零。
*/
ll l,r;
int nxt[1<<10][10],cnt[1<<10];
int idx=0,dig[30],k;
ll dp[20][1<<10][11];
int getnxt(int x,int y){
    rep(i,y,10) if(x&(1<<i)){
        x^=(1<<i);
        break;
    }
    return x|(1<<y);
}
void init(){
    mst(dp,-1);
    rep(i,0,(1<<10)){
        cnt[i]=0;
        rep(j,0,10){
            if(i&(1<<j)) cnt[i]++;
            nxt[i][j]=getnxt(i,j);
        }
    }
}

ll dfs(int pos,int num,int fg,int zero){
    if(pos<0) return cnt[num]==k;
    if(!fg&&dp[pos][num][k]!=-1) return dp[pos][num][k];
    if(cnt[num]>k) return dp[pos][num][k]=0;
    int s=fg?dig[pos]:9;
    ll ret=0;
    for(int i=0;i<=s;i++){
        ret+=dfs(pos-1,(zero&&i==0)?num:nxt[num][i],fg&&(i==s),zero&&(i==0));
    }
    if(!fg) dp[pos][num][k]=ret;
    return ret;
}

ll solve(ll x){
    idx=0;while(x) dig[idx++]=x%10,x/=10;
    return dfs(idx-1,0,1,1);
}

int main(){
    int t;scanf("%d",&t);
    init();
    for(int ca=1;ca<=t;ca++){
       scanf("%lld%lld%d",&l,&r,&k);
       printf("Case #%d: %lld\n",ca,solve(r)-solve(l-1));
    }
    return 0;
}