hdu 4352 XHXJ's LIS 状压数位dp

求a - b 的按每位值来最长lis长度为k的个数
当求lis时 每加入一个数看此时的lis中有没有数大于等于它且前一个数小于它,有的话将这个位置替换成这个数,没有的话把它填入到末尾 lis+1
所以dp[i][sym][len] 代表 推到第i位,此时组成状态为sym, 最后lis长度为len的方案数
注意处理前导0

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#define LL long long
using namespace std;
LL dp[21][1026][23];
LL l, x[100];
LL update(LL n, LL va)
{
    for(int i = va; i<= 9; i++)
        if(n & (1 << i)) return n ^ (1 << i) | (1 << va);
    return n | (1 << va);
}
LL getNum(LL n)
{
   int sum = 0;
   while(n)
   {
       if(n & 1) sum++;
       n >>= 1;
   }
   return sum;
}
LL dfs(LL i, bool e, LL pre, LL len, bool sym) //状态
{
   if(i == -1) return getNum(pre) == len;
   if(e && dp[i][pre][len] != -1) return dp[i][pre][len];
   int Max = e ? 9 : x[i];
   LL ans = 0;
   for(int j = 0; j <= Max; j++)
       ans += dfs(i - 1, !(!e && j == x[i]),  sym||j? update(pre, j) : 0, len, sym || j);
   if(e) dp[i][pre][len] = ans;
   return ans;
}
LL cal(LL n, LL len)
{
    l = 0;
    while(n)
    {
       x[l++] = n % 10;
       n /= 10;
    }
    return dfs(l - 1, 0, 0, len, 0);
}
int main()
{
    memset(dp, -1, sizeof(dp));
    int t, i1 = 1;
    scanf("%d", &t);
    while(t--)
    {
         LL a, b, k;
         scanf("%I64d %I64d %I64d", &a, &b, &k);
         printf("Case #%d: %I64d\n", i1++, cal(b, k) - cal(a - 1, k));
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值