题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4352
题目大意:题意:介绍了电子科大的一个传奇学姐,最后几句话才是题意,这题意思就是给你一个LL范围内的区间[L, R],问你在这个区间内最长递增子序列长度恰为K的数有多少个。
在dp过程中,我们要维护LIS。而且数字只有0-9,如果你知道怎么在O(n * logn)的时间复杂度维护LIS,那么就好办了。但是这个方法用维护一个栈,每次替换栈中第一个大于当前值的数。但是我们怎么把一个大小最大为10的栈表示为状态。就是状态压缩:用一个int来表示栈0000000000。这表示一个空栈。
举例来说:如果sta按照数字13425来更新。
首先遇到1,变成 0100000000 (或者0000000010,其实这是完全一样的,只要保证不同状态的sta不一样就行了)
然后遇到3,很明显,之前没有比3更大的数字,然后变成0101000000
遇到4,sta变成0101100000
在这里打断一下,不能看出,sta中1的个数,就是LIS的长度。
然后遇到2,这时大于等于2的有一个3.于是把3的二进制1交给2,sta变成0110100000 实现替换。遇到4,sta变成0101100000
所以:sta中1的个数,就是LIS的长度。
当然要判断前导0,例如:1 10 2,如果没有判断前导0,那么pos=1时,当前数为0, 栈变1000000000, 再枚举pos=0时,那么1 2 3 4 5 6 7 8 9都满足条件了。所以必须考虑前导0的影响。
#include <bits/stdc++.h>
#define LL long long
using namespace std;
LL a[20];
LL dp[20][1024][15];
LL Size(LL x)//二进制0的个数
{
LL ans=0;
while(x)
{
ans+=(x&1);
x>>=1;
}
return ans;
}
LL Lis(LL st, LL k)//维护LIS
{
for(LL i=k;i<10;i++)//第一个>k的数
{
if(st&(1<<i))
{
return (st^(1<<i))|(1<<k);
}
}
return st|(1<<k);
}
LL dfs(LL pos, LL st, LL Lead, LL mt, LL k)
{
if(pos==-1)
{
if(Size(st)==k)//LIS==k
{
return 1;
}
else
{
return 0;
}
}
if(!Lead&&!mt&&dp[pos][st][k]!=-1)
{
return dp[pos][st][k];
}
LL Len=mt?a[pos]:9;
LL ans=0;
for(LL i=0;i<=Len;i++)
{
ans+=dfs(pos-1, (i==0&&Lead)?st:Lis(st, i), Lead&&(i==0), mt&&(i==a[pos]), k);//考虑前导0
}
if(!mt&&!Lead)
{
dp[pos][st][k]=ans;
}
return ans;
}
LL slove(LL x, LL k)
{
LL pos=0;
while(x)
{
a[pos++]=x%10;
x/=10;
}
return dfs(pos-1, 0, 1, 1, k);
}
int main()
{
LL T, CUT=1;
memset(dp, -1, sizeof(dp));
scanf("%lld",&T);
while(T--)
{
LL L, R, k;
scanf("%lld%lld%lld",&L,&R,&k);
printf("Case #%lld: %lld\n",CUT++, slove(R, k)-slove(L-1, k));
}
return 0;
}