HDU 4734 - F(x)

http://acm.hdu.edu.cn/showproblem.php?pid=4734

For a decimal number x with n digits (A_nA_n-1A_n-2 ... A₂A₁), we define its weight as F(x) = A_n * 2^n-1 + A_n-1 * 2^n-2 + ... + A₂ * 2 + A₁ * 1. Now you are given two numbers A and B, please calculate how many numbers are there between 0 and B, inclusive, whose weight is no more than F(A).

Input

The first line has a number T (T <= 10000) , indicating the number of test cases.
For each test case, there are two numbers A and B (0 <= A,B < 10⁹)

Output

For every case,you should output "Case #t: " at first, without quotes. The t is the case number starting from 1. Then output the answer.

Sample Input

3
0 100
1 10
5 100

Sample Output

Case #1: 1
Case #2: 2
Case #3: 13

求∑Bi=0[F(i)<=F(A)]

https://blog.csdn.net/wust_zzwh/article/details/52100392
题目给了个f(x)的定义：F(x) = An * 2n-1 + An-1 * 2n-2 + ... + A2 * 2 + A1 * 1，Ai是十进制数位，然后给出a,b求区间[0,b]内满足f(i)<=f(a)的i的个数。
常规想：这个f(x)计算就和数位计算是一样的，就是加了权值，所以dp[pos][sum]，这状态是基本的。a是题目给定的，f(a)是变化的不过f(a)最大好像是4600的样子。如果要memset优化就要加一维存f(a)的不同取值，那就是dp[10][4600][4600]，这显然不合法。
这个时候就要用减法了，dp[pos][sum]，sum不是存当前枚举的数的前缀和(加权的)，而是枚举到当前pos位，后面还需要凑sum的权值和的个数，
也就是说初始的是时候sum是f(a),枚举一位就减去这一位在计算f(i)的权值，那么最后枚举完所有位 sum>=0时就是满足的，后面的位数凑足sum位就可以了。
仔细想想这个状态是与f(a)无关的(新手似乎很难理解)，一个状态只有在sum>=0时才满足，如果我们按常规的思想求f(i)的话，那么最后sum>=f(a)才是满足的条件。
int dp[12][N];
int f(int x)
{
    if(x==0) return 0;
    int ans=f(x/10);
    return ans*2+(x%10);
}
int all;
int a[12];
int dfs(int pos,int sum,bool limit)
{
    if(pos==-1) {return sum<=all;}
    if(sum>all) return 0;
    if(!limit && dp[pos][all-sum]!=-1) return dp[pos][all-sum];
    int up=limit ? a[pos] : 9;
    int ans=0;
    for(int i=0;i<=up;i++)
    {
        ans+=dfs(pos-1,sum+i*(1<<pos),limit && i==a[pos]);
    }
    if(!limit) dp[pos][all-sum]=ans;
    return ans;
}
int solve(int x)
{
    int pos=0;
    while(x)
    {
        a[pos++]=x%10;
        x/=10;
    }
    return dfs(pos-1,0,true);
}
int main()
{
    int a,ri;
    int T_T;
    int kase=1;
    scanf("%d",&T_T);
    memset(dp,-1,sizeof dp);
    while(T_T--)
    {
        scanf("%d%d",&a,&ri);
        all=f(a);
        printf("Case #%d: %d\n",kase++,solve(ri));
    }
    return 0;
}
https://blog.csdn.net/feng_zhiyu/article/details/78796133
dp[i][j]表示i位数比j小的数的个数， 另外注意状态转移方程： 
dp[i][j]+=dp[i-1][j-k*(1<<(i-1))]；
https://blog.csdn.net/deepseazbw/article/details/78247376
http://www.voidcn.com/article/p-kjzblhod-bmt.html
https://github.com/AplusB/ACEveryDay/blob/master/QAsQ/Blog/HDU%204734%20F(x).md

常规想：这个f(x)计算就和数位计算是一样的，就是加了权值，所以dp[pos][sum]，这状态是基本的。a是题目给定的，f(a)是变化的不过f(a)最大好像是4600的样子。如果要memset优化就要加一维存f(a)的不同取值，那就是dp[10][4600][4600]，这显然不合法。

这个时候就要用减法了，dp[pos][sum]，sum不是存当前枚举的数的前缀和(加权的)，而是枚举到当前pos位，后面还需要凑sum的权值和的个数，

也就是说初始的是时候sum是f(a),枚举一位就减去这一位在计算f(i)的权值，那么最后枚举完所有位 sum>=0时就是满足的，后面的位数凑足sum位就可以了。

仔细想想这个状态是与f(a)无关的(新手似乎很难理解)，一个状态只有在sum>=0时才满足，如果我们按常规的思想求f(i)的话，那么最后sum>=f(a)才是满足的条件。

1. const int N =1e4+5;
2. int dp[12][N];
3. int all;
4. int a[12];
5. int f(int x)
6. {
7. if(x==0)
8. return 0;
9. int ans=f(x/10);
10. return ans*2+(x%10);
11. }
12. int dfs(int pos,int sum,bool limit)
13. {
14. if(pos==-1) return sum<=all;
15. if(sum>all) return 0;
16. if(!limit&&dp[pos][all-sum]!=-1) return dp[pos][all-sum];
17. int up=limit?a[pos]:9;
18. int ans=0;
19. for(int i=0;i<=up;i++)
20. ans+=dfs(pos-1,sum+i*(1<<pos),limit&&i==a[pos]);
21. if(!limit) dp[pos][all-sum]=ans;
22. return ans;
24. }
25. int solve(int x)
26. {
27. int pos=0;
28. while(x)
29. {
30. a[pos++]=x%10;
31. x/=10;
32. }
33. return dfs(pos-1,0,true);
34. }
35. int main()
36. {
37. int a,ri,icase,kase=1;
38. scanf("%d",&icase);
39. while(icase--)
40. {
41. scanf("%d%d",&a,&ri);
42. all=f(a);
43. printf("Case #%d: %d\n",kase++,solve(ri)-0);
44. }
45. return 0;
46. }

http://semon0x29a.cn/2018/09/12/HDU-4734/

一开始我的做法是基于dp[u][prefix]$dp[u][prefix]$的数位dp， 第一维表示当前数位， 第二维表示前部的级数和， 但是数据组数太多，会TLE。
于是考虑将第二维倒过来，变为dp[u][F(A)−prefix]$dp[u][F(A)-prefix]$， 这样就可以复用前面数据计算过的dp数组， 于是O(Accepted)$O(Accepted)$

int n, a, b;
int v[12];
int dp[12][9300];

int dfs(int u, bool limits, int pre) {
    if (pre < 0) return 0;
    if (u == -1) return 1;
    if (!limits && dp[u][pre] != -1) {
        return dp[u][pre];
    }
    int res = 0, top = limits ? v[u] : 9;
    for (int i = 0; i <= top; i++) {
        res += dfs(u - 1, limits && (i == top), pre - i * (1 << u));
    }
    if (!limits) dp[u][pre] = res;
    return res;
}

int calc() {
    if (b == 0) return 1;
    for (n = 0; b > 0; b /= 10) {
        v[n++] = b % 10;
    }
    int sum = 0, mul = 1;
    while (a > 0) {
        sum += (a % 10) * mul;
        mul <<= 1;
        a /= 10;
    }
    return dfs(n - 1, true, sum);
}

http://vawait.com/2013/09/hdu-4734/
f[i][j]表示有i位数随意填F()值不超过j的方案数。

void init()

{

    b[1] = 1;

    rep(i,2,12) b[i] = b[i-1] * 2;

    clr(f,0);

    f[0][0] = 1;

    rep(i,1,10)

        rep(j,0,5000)

            rep(k,0,9) if ( k * b[i] <= j ) f[i][j] += f[i-1][j-k*b[i]];

    rep(i,0,10)

        rep(j,1,5000) f[i][j] += f[i][j-1];

}

int calc(int k)

{

    s = 0;

    while ( k ) {

        a[++s] = k % 10;

        k /= 10;

    }

    int sum=0;

    rep(i,1,s) sum += a[i] * b[i];

    return sum;

}

void dfs(int t,int s)

{

    int k;

    if ( s < 1 || t > m ) return;

    rep(i,0,a[s] - 1) {

        k = t + i * b[s];

        if ( k > m ) continue;

        ans += f[s-1][m-k];

    }

    dfs( t + a[s] * b[s], s - 1);

}