HDU 4352

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cmath>
 4 #include <string.h>
 5 #include <algorithm>
 6 #define LL __int64
 7 using namespace std;
 8 
 9 LL dp[25][1<<10][11];
10 LL aa,bb;
11 int a[30],k;
12 /*这道题的难点就在于DP的状态表示了，确实是这样的。DP【pos】【state】【K】,pos表示当前的数位，state表示当前的LIS的状态，K
13 表示要求的K长度。那个DFS就不解释了，都成模板了。而对于state的转移，是用二进制表示当前状态，因为最长的
14 序列就是0~9十个数表示，所以有十位的二进制。例如12456的上升序列，下一次出现3，则可表示为12356，至于转移，要联系nlogn解法的
15 LIS来看了 */ 
16 int getnews(int x,int s){
17     for(int i=x;i<10;i++)
18         if(s&(1<<i))return (s^(1<<i))|(1<<x);
19     return s|(1<<x);
20 }
21 int getnum(int s){
22     int ret=0;
23     while(s)
24     {
25         if(s&1)ret++;
26         s>>=1;
27     }
28     return ret;
29 }
30 
31 LL dfs(int len,int st,bool zero,bool flag){
32     if(len==0) return getnum(st)==k;
33     if(!flag&&dp[len][st][k]!=-1)
34     return dp[len][st][k];
35     int en=flag?a[len]:9;
36     LL ans=0;
37     for(int i=0;i<=en;i++){    
38         ans+=dfs(len-1,(zero&&i==0)?0:getnews(i,st),zero&&i==0,flag&&i==en);
39     }
40     if(!flag) dp[len][st][k]=ans;
41     return ans;
42 }
43 
44 LL cal(LL n){
45     LL tmp=n;
46     int len=0;
47     while(tmp){
48         a[++len]=(int)(tmp%10);
49         tmp/=10;
50     }
51     return dfs(len,0,true,true);
52 }
53 
54 int main(){
55     memset(dp,-1,sizeof(dp));
56     int T,t=0;
57     scanf("%d",&T);
58     while(++t<=T){
59         scanf("%I64d%I64d%d",&aa,&bb,&k);
60         printf("Case #%d: %I64d\n",t,cal(bb)-cal(aa-1LL));
61     }
62     return 0;
63 }