hdu4352 数位dp+状态压缩+一个tip

按照nlogn求lis的方法，把lis的状态压缩了，每次新加一个数就把它右边第一个数的位置置为0，然后把这个数加进去

一个需要注意的地方，如果前面都是0，那么状态s中代表0的位置不可以是1，因为这种情况下0不可以被算作是lis里的一位

/*
dp[i][j][k]表示后面还有i位，前面的lis状态是j，要求lis长度为k的个数
*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll n,k,a[],dp[][<<][];
ll update(ll x,ll s){
    for(int i=x;i<;i++)
        if(s&(<<i))
            return (s^(<<i))|(<<x);
    return s|(<<x);
}
ll calc(ll s){
    ll res=;
    for(int i=;i<=;i++)
        if(s&(<<i))res++;
    return res;
}
ll dfs(ll pos,ll s,ll lim,int zero){//zero=1表示前面都是0
    if(pos<=)return calc(s)==k;
    if(!lim && dp[pos][s][k]!=-)return dp[pos][s][k];

    ll res=,num=lim?a[pos]:;
    for(int i=;i<=num;i++){
        int z=zero&&(i==);
        res+=dfs(pos-,z?:update(i,s),lim&&i==num,z);
    }
    if(!lim)dp[pos][s][k]=res;
    return res;
}
ll solve(ll x){
    if(x<)return ;
    int len=;
    memset(a,,sizeof a);
    while(x){
        a[++len]=x%;
        x/=;
    }
    return dfs(len,,,);
}
int main(){
    int t;
    cin>>t;
    memset(dp,-,sizeof dp);
    ll a,b;
    for(int tt=;tt<=t;tt++){
        cin>>a>>b>>k;
        printf("Case #%d: ",tt);
        cout<<solve(b)-solve(a-)<<endl;
    }
}