hdu 4352 XHXJ's LIS (数位dp+状态压缩)

Description

#define xhxj (Xin Hang senior sister(学姐))
If you do not know xhxj, then carefully reading the entire description is very important.
As the strongest fighting force in UESTC, xhxj grew up in Jintang, a border town of Chengdu.
Like many god cattles, xhxj has a legendary life:
2010.04, had not yet begun to learn the algorithm, xhxj won the second prize in the university contest. And in this fall, xhxj got one gold medal and one silver medal of regional contest. In the next year's summer, xhxj was invited to Beijing to attend the astar onsite. A few months later, xhxj got two gold medals and was also qualified for world's final. However, xhxj was defeated by zhymaoiing in the competition that determined who would go to the world's final(there is only one team for every university to send to the world's final) .Now, xhxj is much more stronger than ever,and she will go to the dreaming country to compete in TCO final.
As you see, xhxj always keeps a short hair(reasons unknown), so she looks like a boy( I will not tell you she is actually a lovely girl), wearing yellow T-shirt. When she is not talking, her round face feels very lovely, attracting others to touch her face gently。Unlike God Luo's, another UESTC god cattle who has cool and noble charm, xhxj is quite approachable, lively, clever. On the other hand,xhxj is very sensitive to the beautiful properties, "this problem has a very good properties",she always said that after ACing a very hard problem. She often helps in finding solutions, even though she is not good at the problems of that type.
Xhxj loves many games such as,Dota, ocg, mahjong, Starcraft 2, Diablo 3.etc,if you can beat her in any game above, you will get her admire and become a god cattle. She is very concerned with her younger schoolfellows, if she saw someone on a DOTA platform, she would say: "Why do not you go to improve your programming skill". When she receives sincere compliments from others, she would say modestly: "Please don’t flatter at me.(Please don't black)."As she will graduate after no more than one year, xhxj also wants to fall in love. However, the man in her dreams has not yet appeared, so she now prefers girls.
Another hobby of xhxj is yy(speculation) some magical problems to discover the special properties. For example, when she see a number, she would think whether the digits of a number are strictly increasing. If you consider the number as a string and can get a longest strictly increasing subsequence the length of which is equal to k, the power of this number is k.. It is very simple to determine a single number’s power, but is it also easy to solve this problem with the numbers within an interval? xhxj has a little tired,she want a god cattle to help her solve this problem,the problem is: Determine how many numbers have the power value k in [L,R] in O(1)time.
For the first one to solve this problem,xhxj will upgrade 20 favorability rate。

Input

First a integer T(T<=10000),then T lines follow, every line has three positive integer L,R,K.(
0<L<=R<2 63-1 and 1<=K<=10).

Output

For each query, print "Case #t: ans" in a line, in which t is the number of the test case starting from 1 and ans is the answer.

Sample Input

1
123 321 2

Sample Output

Case #1: 139 

题目的大概意思就是让你统计在给定区间内,符合要求的数的个数。 一个数如果它的各个数位的最长上升子序列(LIS)长度为k,那么它就是符合要求的。
这题分为三个点。
1.首先这个题符合区间减法,我们只需要求出0~l-1和0~r的合法数的个数,再做减法即可。
2.对LIS的处理我们采用状态压缩来处理
LIS状压:
这个数的二进制的第i个1的位置表示当前序列长度为i的LIS最后一位最小是多少。这样1的个数就是LIS的长度
比如010101101(从左到右,每一位代表0~9)表示这个序列以9结尾(最后一个1在9的位置)
上升子序列长度分别为1,2,3,4,5的子序列的结尾(最后一位)最小分别是2,4,6,7,9。
状态的更新:
假设现在状态为0100100110,最长序列为1 4 7 8,如果我们下一个dp位的值为6,那么长度为3的上升子序列就由原来的1 4 7,变为1 4 6。
相应的我们更新后状态为0100101010,相当于6把7在二进制数上替换了(红色)。
3.剩下的交给深搜...
 1 #include <bits/stdc++.h>
 2 
 3 using namespace std;
 4 long long dp[25][1<<10][11];
 5 //dp[i][j][k] ,i为当前进行到的数位,j状态压缩,k要求的上升子序列的度
 6 int k,bit[25];
 7 int getnews(int x,int s)//更新状态
 8 {
 9     for (int i=x;i<10;++i)
10     if (s&1<<i) return (s^(1<<i)|(1<<x));//让x替换x位后的第一个数
11     return s|(1<<x);//如果x最大,将x对应位变为1
12 }
13 int getnum(int s)//求出当前状态LIS的长度也就是1的个数
14 {
15     int ret=0;
16     while (s)
17     {
18         if (s&1)
19         ret++;
20         s>>=1;
21     }
22     return ret;
23 }
24 long long dfs (int pos,int s,bool e,bool z)
25 //pos表示处理过的位置,s表示LIS的状态,e表示是否到达边界,z表示前面的所有位是否全为0
26 {
27     if (pos==-1) return getnum(s)==k;
28     if (!e&&dp[pos][s][k]!=-1) return dp[pos][s][k];
29     long long ans=0;
30     int endd=e?bit[pos]:9;
31     for (int i=0;i<=endd;++i)
32     ans+=dfs(pos-1,(z&&i==0)?0:getnews(i,s),e&&i==endd,z&&(i==0));
33     if (!e) dp[pos][s][k]=ans;
34     return ans;
35 }
36 long long calc (long long n)
37 {
38     int len=0;
39     while (n)
40     {
41         bit[len++]=n%10;
42         n/=10;
43     }
44     return dfs(len-1,0,1,1);
45 }
46 int main()
47 {
48     //freopen("de.txt","r",stdin);
49     int t;
50     long long l,r;
51     memset(dp,-1,sizeof dp);
52     scanf("%d",&t);
53     int casee=0;
54     while (t--)
55     {
56         scanf("%I64d%I64d%d",&l,&r,&k);
57         printf("Case #%d: ",++casee);
58         printf("%I64d
",calc(r)-calc(l-1));
59     }
60     return 0;
61 }

结尾自满一下,这个题网上应该没有人比我写的题解详细,虽然我也是基本上照搬kuangbin大神的模板,再看看网上的题解综合自己的想法写出来的。

 
原文地址:https://www.cnblogs.com/agenthtb/p/5910483.html