题意:[l,r]区间之内,有多少数字构成了长度为k的严格上升子序列。
首先我们需要确定一个状态来表示各个数位之间的共性,由于数字数目很小,只有10个,所以可以考虑状压。如果用dp[i][state]表示考虑到第i位,0-9哪些数字选或不选,并不能直接的表示出当前的状态,即使再加一维[j]表示最长上升子序列长度,每次添加的时候也不知道会加多少。所以我们可以稍微改变state的定义,用第i个为1的数来表示长度为i的最长上升子序列的最后一位最小是多少,就如nlogn最长上升子序列求法定义的那样。(最后多加了一维k,表示输入的那个k下的dp[i][state],因为询问很多,k很少,所以每次都用预处理出来的可以减小常数)。
接着考虑转移,有了上文nlogn最长上升子序列的借鉴,我们就可以类比转移了,枚举下一个填1-9哪个数字,看能放在当前state下的第几位,再把原来比它大的表示长度为i的最长上升子序列的那一位删掉即可。
最后边界条件即为递归到最后一层,且最长上升子序列长度为k则答案是1,否则答案是0。
下附AC代码。
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define maxn 65
using namespace std;
typedef long long ll;
ll l,r,k,tot;
ll bit[maxn];
ll dp[maxn][(1<<10)][11];
ll lowbit(ll now)
{
return (now&(-now));
}
ll calbit(ll now)
{
ll ans=0;
while(now)
ans++,now-=lowbit(now);
return ans;
}
ll getstate(ll state,int pos)
{
for(int i=pos;i<=9;i++)
if((state&(1<<i)))
return ((state^(1<<i))|(1<<pos));
return state|(1<<pos);
}
ll dfs(ll pos,ll state,int iszero,int isin)
{
if(pos==0) return calbit(state)==k;
if(!isin && dp[pos][state][k]!=-1) return dp[pos][state][k];
int maxx=isin?bit[pos]:9;
ll ans=0;
for(int i=0;i<=maxx;i++)
ans+=dfs(pos-1,(iszero&&i==0)?0:getstate(state,i),(iszero&&i==0),(isin&&i==bit[pos]));
if(!isin) dp[pos][state][k]=ans;
return ans;
}
ll solve(ll now)
{
tot=0;
while(now)
{
bit[++tot]=now%10;
now/=10;
}
return dfs(tot,0,1,1);
}
int main()
{
memset(dp,-1,sizeof(dp));
int _,kase=0;;
scanf("%d",&_);
while(_--)
{
scanf("%lld%lld%lld",&l,&r,&k);
printf("Case #%d: %lld\n",++kase,solve(r)-solve(l-1));
}
}