HDU 4734 F(x)（数位DP）

最新推荐文章于 2020-09-19 00:23:43 发布

glq007

最新推荐文章于 2020-09-19 00:23:43 发布

阅读量662

点赞数

分类专栏： DP 文章标签： dp

本文链接：https://blog.csdn.net/glqac/article/details/41130065

版权

DP 专栏收录该内容

66 篇文章 0 订阅

订阅专栏

题意很简单， F(x) = A_n * 2^n-1 + A_n-1 * 2^n-2 + ... + A₂ * 2 + A₁ * 1.，给出A，B，求出0到B内f（x）<= f(A)的个数。

时间给的太少了，0.5s，所以必须预处理，dp[i][j]代表位数为i的数字（首位可以为0），f(x)等于j的数字个数。

转移方程就是dp[i][j] = dp[i-1][j-k*jie[i-1]]，jie[]是2的阶乘。然后再对dp[i][j]求前缀和，那么dp[i][j]就变为位数为i，f(x)小于等于j的数字个数。

预处理完成后，从B的高位开始枚举，比如最高位是第5位枚举为4时，那么答案就加上dp[4][f( A ) - 4*2^4]，枚举到B[5]-1，然后转移到第4位的时候把f（A）- B[5]*2^4。以此类推。最后还要判断一下是否剩下的大于等于0，这意味着f(B) >= f(A)，说明这也要加到答案里。

AC代码：

//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<ctype.h>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
#include<cstdlib>
#include<stack>
#include<cmath>
#include<queue>
#include<set>
#include<map>
#include<ctime>
#include<string.h>
#include<string>
#include<sstream>
#include<bitset>
using namespace std;
#define ll long long
#define eps 1e-13
#define NMAX 200005
#define MOD 1000000009
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define PI acos(-1)
template<class T>
inline void scan_d(T &ret)
{
    char c;
    int flag = 0;
    ret=0;
    while(((c=getchar())<'0'||c>'9')&&c!='-');
    if(c == '-')
    {
        flag = 1;
        c = getchar();
    }
    while(c>='0'&&c<='9') ret=ret*10+(c-'0'),c=getchar();
    if(flag) ret = -ret;
}
int dp[11][4599+10];//位数等于第一维f(x)值为第二维的个数
int jie[10],shu[10],nct,pp;
void init()
{
    for(int k = 1; k <= 9; k++)
    {
        for(int i = 0; i <= 9; i++)
            for(int j = i*jie[k-1]; j <= 4600; j++)
                dp[k][j] += dp[k-1][j - i*jie[k-1]];
//        for(int i = 1; i <= 4599; i++) dp[k][i] += dp[k][i-1];
    }
}

int main()
{
#ifdef GLQ
    freopen("input.txt","r",stdin);
//    freopen("o4.txt","w",stdout);
#endif // GLQ
    jie[0] = 1;
    for(int i = 1; i <= 8; i++) jie[i] = jie[i-1]*2;
    dp[0][0] = 1;
    init();
    for(int i = 0; i <= 9; i++)
        for(int j = 1; j <= 4599; j++) dp[i][j] += dp[i][j-1];
    int t,a,b,cas=1;
    scanf("%d",&t);
//    cout<<dp[5][123]<<endl;
    while(t--)
    {
        scanf("%d%d",&a,&b);
        int tmp=b,ans = 0;
        nct= pp = 0;
        while(a)
        {
            pp += (a%10)*jie[nct++];
            a /= 10;
        }
        nct = 0;
        while(tmp)
        {
            shu[nct+1] = tmp%10;
            tmp /= 10;
            nct++;
        }
//        ans += dp[nct-1][pp];
        for(int i = nct; i >= 1; i--)
        {
            for(int j = shu[i]-1; j >= 0; j--)
                if(pp >= j*jie[i-1]) ans += dp[i-1][pp-j*jie[i-1]];
            pp -= shu[i]*jie[i-1];
        }
        if(pp>=0) ans++;
        printf("Case #%d: %d\n",cas++,ans);
    }
    return 0;
}