hdu4352 (数位dp+状态压缩)

HDU 4352 XHXJ’s LIS

传送门
题目描述：
给定一个区间中，将区间的每一个数看成一个字符串，求这个区间内每个字符串的最大上升
子序列等于k的个数。
思路：
dp[pos][sta][k]
pos:搜索到第几位
sta：之前的最长上升子序列的状态。
k：就是输入的那个k
牛逼的地方：
sta储存为0-9每个数字的状态 类似最长上升子序列更新方式。sta二进制下1的个数就是当前最长上升子序列的长度
比如：1245->12453 LIS都为3 状态sta变化为 00000110110->00000101110
即寻找3在原数组lower_bound的位置为4 把4位置0，3位置1

#pragma GCC optimize(3,"Ofast","inline")  	//G++
#include<bits/stdc++.h>
#define mem(a,x) memset(a,x,sizeof(a))
#define debug(x) cout << #x << ": " << x << endl;
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define fcout cout<<setprecision(4)<<fixed
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pii;

const ll inf=0x3f3f3f3f;
const ll mod=1e9+7;
const ll maxn = 1e6+5;

template<typename T> void read(T &x){
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
template<typename T, typename... Args> void read(T &first, Args& ... args) {
    read(first);
    read(args...);
}
ll dp[20][1<<10][11];
ll a[20];

ll get_sta(ll sta,ll k){
    for(ll i=k;i<=9;i++){
        if(sta&(1<<i)) return (sta^(1<<i))|(1<<k);
    }
    return sta|(1<<k);
}
ll dfs(ll pos,ll sta,ll k,bool limit,bool lead){
    if(!limit&&dp[pos][sta][k]!=-1)  return dp[pos][sta][k];
    if(pos==0) return __builtin_popcount(sta)==k;
    ll up=limit?a[pos]:9;
    ll res=0;
    for(ll i=0;i<=up;i++){
        ll t=sta;
        if(!(lead&&i==0)){
            t=get_sta(sta,i);
        }
        res+=dfs(pos-1,t,k,limit&&a[pos]==i,lead&&i==0);
    }
    if(!limit) dp[pos][sta][k]=res;
    return res;
}

ll solve(ll n,ll k){
    ll pos=0;
    while(n){
        a[++pos]=n%10;
        n/=10;
    }
    return dfs(pos,0,k,1,1);
}
 main()
{
    int Tim=0;
    mem(dp,-1);
    ll T;
    read(T);
    while(T--){
        ll l,r,k;
        read(l,r,k);
        printf("Case #%d: ",++Tim);
        cout<<solve(r,k)-solve(l-1,k)<<"\n";
    }
}

__builtin_popcount 返回值为unsigned int