HDU 4352 XHXJ's LIS(数位dp)

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

 

题目大意：求l到r中有多少数字满足LIS长度=k

 

题目思路：感觉这个题特别有意思，结合了lis nlogn做法和数位dp，首先得知道一个事情，就是nlogn的做法是找到大于等于当前数的第一个位置，把那个位置上的数改成当前数(不会的去找博客学一下)，在这题里，直接就是要插入的数字往后找第一个是1的位，把它变成0，然后插入自己，如果没有1的位直接插入自己，status用一个数字表示，每一位表示维护的序列的情况，最后有多少1就是长度多少

数位dp常见问题，千万注意前导0

 

以下是代码：

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define rep(i,a,b) for(ll i=a;i<=b;i++)
#define per(i,a,b) for(ll i=a;i>=b;i--)
#define ll long long
const ll MAXN = 100+5;
const ll MOD =1e9+7;
ll t,l,r,k,a[MAXN],dp[MAXN][1200][15];
ll check(ll x){
    ll num=0;
    while(x){
        if(x&1)num++;
        x>>=1;
    }
    return num==k;
}
ll getsta(ll sta,ll pos){
    rep(i,pos,9){
        if(sta&(1<<i))return (sta^(1<<i))|(1<<pos);
    }
    return sta|(1<<pos);
}
ll dfs(ll pos,ll sta,ll flag,ll limit){
    if(!pos) return check(sta);
    if(!limit&&dp[pos][sta][k]!=-1)return dp[pos][sta][k];
    ll up=limit?a[pos]:9,ans=0;
    rep(i,0,up){
        ans+=dfs(pos-1,(!flag&&!i)?0:getsta(sta,i),flag||i,limit&&i==a[pos]);
    }
    if(!limit)dp[pos][sta][k]=ans;
    return ans;
}
ll solve(ll x){
    ll pos=0;
    while(x){
        a[++pos]=x%10;
        x/=10;
    }
    return dfs(pos,0,0,1);
}
int main(){
    scanf("%lld",&t);
    memset(dp,-1,sizeof(dp));
    rep(_,1,t){
        scanf("%lld%lld%lld",&l,&r,&k);
        printf("Case #%lld: %lld\n",_,solve(r)-solve(l-1));
    }
    return 0;
}