hdu4352_状压+数位+lis

题目链接

原文：http://www.acmerblog.com/hdu-4352-xhxjs-lis-7363.html

 * 问L到R，各位数字组成的严格上升子序列的长度为K的个数。
 * 0<L<=R<263-1 and 1<=K<=10
 * 注意这里最长上升子序列的定义，和LIS是一样的，不要求是连续的
 * 所以用十位二进制表示0~9出现的情况，和O(nlogn)求LIS一样的方法进行更新

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <ctime>
#include <queue>
#include <list>
#include <set>
#include <map>
using namespace std;
#define INF 0x3f3f3f3f
typedef long long LL;

int bit[35];
LL dp[35][15][1 << 10], k;

int GetNext(int st,int x)  
{  //刚开始用数组存，会超时啊（大哭），还是太嫩了
    //找到第一个大于x的数并且替换他，nlogn求最长递增子序列的思想  
    for(int i=x;i<10;i++)  
        if(st&(1<<i)) return ((st^(1<<i))|(1<<x));  
    return st|(1<<x);  
}  
int GetLen(int st)  
{  
    int cnt=0;  
    while(st)  
    {  
        if(st&1) cnt++;  
        st>>=1;  
    }  
    return cnt;  
}  
LL dfs(int len, int length, int mark, int flag)
{//刚开始一直用dp[len][length][mark], 跟length没有毛线关系啊，wa了n发
    if(len == 0)
        return GetLen(mark) == k;
    if(flag && dp[len][k][mark] >= 0)
        return dp[len][k][mark];
    LL sum = 0;
    int te = flag ? 9 : bit[len];
    for(int i = 0; i <= te; i++)
    {
        int Mark = ((length==0 && i==0) ? 0 : GetNext(mark, i));
        sum += dfs(len - 1, length || i , Mark, flag || (i < te));
    }
    if(flag)
        dp[len][k][mark] = sum;
    return sum;
}
LL solve(LL n)
{
    int len = 0;
    while(n)
    {
        bit[++len] = n % 10;
        n /= 10;
    }
    return dfs(len, 0, 0, 0);
}
int main()
{
    int t;
    LL l, r;
    scanf("%d", &t);
    memset(dp, -1, sizeof(dp));
    for(int i = 1; i <= t; i++)
    {
        cin >> l >> r >> k;
        cout << "Case #"<< i << ": ";
        cout << solve(r) - solve(l - 1) << endl;
    }
    return 0;
}