CQOI2016 Number - 数位dp

题目描述

分析:

dp[哪一位][前一个数的倒数第二位][前一个数的倒数最后一位][0/1表示是否满足了至少有三个相邻数字的要求][0/1/2 (0:没有8和4;1:有8;2:有4)][0/1表示最高位到当前位小于或者等于R]
具体转移都在代码里了。

#include<cstdio>
#include<cstring>
#define MAXL 20
typedef long long LL;

int lmt[MAXL+10];
LL dp[MAXL+10][20][20][5][5][5];

int Getorg(LL num){
    int ret=0;
    for(LL x=num;x;x/=10)
        ret++;
    memset(lmt,0,sizeof lmt);
    int k=0;
    for(LL x=num;x;x/=10,k++)
        lmt[ret-k]=x%10;
    return ret;
}
inline int f(int x){
    if(x==8) return 1;
    else if(x==4) return 2;
    else return 0;
}
inline int h(int x,int y,int z){
    if(x==y&&y==z) return 1;
    else return 0;
}
LL DP(LL R)
{
    int n=Getorg(R);
    if(R==10000000000LL-1) return 0;
    memset(dp,0,sizeof dp);
    for(int i=1;i<lmt[1];i++)
        dp[1][0][i][0][f(i)][0]=1;
    dp[1][0][lmt[1]][0][f(lmt[1])][1]=1;
    for(int i=2;i<=n;i++){
        for(int j=0;j<=9;j++)
            for(int k=0;k<=9;k++){
                // < R
                for(int p=0;p<=9;p++)
                    dp[i][k][p][h(j,k,p)][f(p)][0]+=dp[i-1][j][k][0][0][0];
                for(int p=0;p<=9;p++){
                    if(p==4) continue;
                    dp[i][k][p][h(j,k,p)][1][0]+=dp[i-1][j][k][0][1][0];
                }
                for(int p=0;p<=9;p++){
                    if(p==8) continue;
                    dp[i][k][p][h(j,k,p)][2][0]+=dp[i-1][j][k][0][2][0];
                }
                for(int p=0;p<=9;p++)
                    dp[i][k][p][1][f(p)][0]+=dp[i-1][j][k][1][0][0];
                for(int p=0;p<=9;p++){
                    if(p==4) continue;
                    dp[i][k][p][1][1][0]+=dp[i-1][j][k][1][1][0];
                }
                for(int p=0;p<=9;p++){
                    if(p==8) continue;
                    dp[i][k][p][1][2][0]+=dp[i-1][j][k][1][2][0];
                }
                // < R
                for(int p=0;p<lmt[i];p++)
                    dp[i][k][p][h(j,k,p)][f(p)][0]+=dp[i-1][j][k][0][0][1];
                for(int p=0;p<lmt[i];p++){
                    if(p==4) continue;
                    dp[i][k][p][h(j,k,p)][1][0]+=dp[i-1][j][k][0][1][1];
                }
                for(int p=0;p<lmt[i];p++){
                    if(p==8) continue;
                    dp[i][k][p][h(j,k,p)][2][0]+=dp[i-1][j][k][0][2][1];
                }
                for(int p=0;p<lmt[i];p++)
                    dp[i][k][p][1][f(p)][0]+=dp[i-1][j][k][1][0][1];
                for(int p=0;p<lmt[i];p++){
                    if(p==4) continue;
                    dp[i][k][p][1][1][0]+=dp[i-1][j][k][1][1][1];
                }
                for(int p=0;p<lmt[i];p++){
                    if(p==8) continue;
                    dp[i][k][p][1][2][0]+=dp[i-1][j][k][1][2][1];
                }

                // == R
                for(int p=lmt[i];p<=lmt[i];p++)
                    dp[i][k][p][h(j,k,p)][f(p)][1]+=dp[i-1][j][k][0][0][1];
                for(int p=lmt[i];p<=lmt[i];p++){
                    if(p==4) continue;
                    dp[i][k][p][h(j,k,p)][1][1]+=dp[i-1][j][k][0][1][1];
                }
                for(int p=lmt[i];p<=lmt[i];p++){
                    if(p==8) continue;
                    dp[i][k][p][h(j,k,p)][2][1]+=dp[i-1][j][k][0][2][1];
                }
                for(int p=lmt[i];p<=lmt[i];p++)
                    dp[i][k][p][1][f(p)][1]+=dp[i-1][j][k][1][0][1];
                for(int p=lmt[i];p<=lmt[i];p++){
                    if(p==4) continue;
                    dp[i][k][p][1][1][1]+=dp[i-1][j][k][1][1][1];
                }
                for(int p=lmt[i];p<=lmt[i];p++){
                    if(p==8) continue;
                    dp[i][k][p][1][2][1]+=dp[i-1][j][k][1][2][1];
                }
            }
    }
    LL ret=0;
    for(int j=0;j<=9;j++)
        for(int k=0;k<=9;k++)
            for(int t=0;t<3;t++)
                for(int d=0;d<2;d++)
                    ret+=dp[n][j][k][1][t][d];
    return ret;
}
int main()
{
    //freopen("number.in","r",stdin);
    //freopen("number.out","w",stdout);

    LL L,R;
    scanf("%lld%lld",&L,&R);
    printf("%lld\n",DP(R)-DP(L-1));
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值