HDU3652：B-number——题解

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

题目大意：给一个数n，求1～n所有满足下列条件的数的个数：

1.包含一个子串为“13”

2.能被13整除。

————————————————————

数位dp……然而不像各位大佬记忆化搜索的写法那样清真。

欢迎访问旁边友链的学姐来获取这个我抄的代码

设dp[i][j][k][0/1][0/1]表示第i位填j，现在该数被13除余数为k，无/有“13”，前i位小于等于/大于n的前i位。

然后推的方程……太长了，看代码吧。

虽然很长但是好理解不是吗……

#include<cstdio>
#include<cstring>
using namespace std;
int a[],b[],f[][][][][];
int dp(int x){
    int len=;
    while(x)a[++len]=x%,x/=;
    if(len==)a[++len]=;
    memset(f,,sizeof(f));
    for(int i=;i<=;i++){
    if(i<=a[])f[][i][i][][]=;
    else f[][i][i][][]=;
    }
    for(int i=;i<=len;i++){
    for(int j=;j<=;j++){
        for(int k=;k<=;k++){
        for(int l=;l<;l++){
            bool t=;
            if(j==&&k==)t=;
            if(j<a[i])
            if(t)
                f[i][j][(l+j*b[i-])%][t][]+=f[i-][k][l][][]+f[i-][k][l][][]
                +f[i-][k][l][][]+f[i-][k][l][][];
            else{
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][]+f[i-][k][l][][];
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][]+f[i-][k][l][][];
            }
            else if(j==a[i])
            if(t){
                f[i][j][(l+j*b[i-])%][t][]+=f[i-][k][l][][]+f[i-][k][l][][];
                f[i][j][(l+j*b[i-])%][t][]+=f[i-][k][l][][]+f[i-][k][l][][];
            }
            else{
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][];
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][];
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][];
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][];
            }
            else
            if(t)
                f[i][j][(l+j*b[i-])%][t][]+=f[i-][k][l][][]+f[i-][k][l][][]
                +f[i-][k][l][][]+f[i-][k][l][][];
            else{
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][]+f[i-][k][l][][];
                f[i][j][(l+j*b[i-])%][][]+=f[i-][k][l][][]+f[i-][k][l][][];
            }
        }
        }
    }
    }
    int ans=;
    for(int i=;i<=a[len];i++)ans+=f[len][i][][][];
    for(int i=len-;i;i--){
    for(int j=;j<=;j++){
        ans+=f[i][j][][][]+f[i][j][][][];
    }
    }
    return ans;
}
int main(){
    b[]=;
    for(int i=;i<=;i++)b[i]=b[i-]*%;
    int n;
    while(scanf("%d",&n)!=EOF&&n){
    printf("%d\n",dp(n));
    }
    return ;
}