【dp】B-number

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3652

题解：先预处理([0,0][1,1],[2,2]....[0,9],[10, 19],[20,29]......[900000000, 1000000000]  区域中 有子串13 （用1表示）余数为0,1,2.。。12的个数。以及 无子串13（用0表示）余数为0,1,2,3.。。12 的个数。

四维数组dp【一共多少位数】【最高位的数】【是否含有子串13】【余数】 例如  dp[4][1][1][0] 表示 （【1000,2000)）（四位， 最高位数为1）， 子串中含有13， 且余数为0的个数。

dp[a][b][c][d]

当求a = 2时必须 a=1的元素都已知，例如 求[20, 29）中含有13， 余数为3的个数。 即 dp【2】【2】【0】【3】 ，由于 [20, 29）可看做 20 + x（0—9）。而20 % 13 = 7， 只需找余数为9. 则 sum（ dp【1】【j】【0】【9】）（j= 0,1,2.。。。9）。

需注意的是 当b = 1 时， 由于在低一位上有3，  所以这种情况要分开求。

/***Good Luck***/
#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <algorithm>
#include <stack>
#include <map>
#include <queue>
#include <vector>
#include <set>
#include <functional>
#include <cmath>
#include <numeric>

#define Zero(a) memset(a, 0, sizeof(a))
#define Neg(a)  memset(a, -1, sizeof(a))
#define All(a) a.begin(), a.end()
#define PB push_back
#define inf 0x3f3f3f3f
#define inf2 0x7fffffffffffffff
#define ll long long
using namespace std;
//#pragma comment(linker, "/STACK:102400000,102400000")
void get_val(int &a) {
    int value = , s = ;
    char c;
    while ((c = getchar()) == ' ' || c == '\n');
    if (c == '-') s = -s; else value = c - ;
    while ((c = getchar()) >= '' && c <= '')
        value = value *  + c - ;
    a = s * value;
}
const int maxn = ;
int dp[maxn][maxn][][maxn];
int n;

void change1(int i, int j, int m) {
    for (int k = ; k <= ; ++k)
    for (int kk = ; kk <= ; ++kk) {
        dp[i][j][][(kk - m + ) % ] += dp[i - ][][k][kk];
    }
}

void change2(int i, int m,int j, int jj) {
    for (int k = ; k <= ; ++k)
    for (int kk = ; kk <= ; ++kk) {
        dp[i][j][k][(kk + m) % ] += dp[i - ][jj][k][kk];
    }
}
void init() {
    Zero(dp);
    int w = , m;
    dp[][][][] = ;
    for (int i = ; i <= ; ++i) dp[][i][][i] = ;
    for (int i = ; i <= ; ++i) {
        for (int j = ; j <= ; ++j) {

            for (int jj = ; jj <= ; ++jj) {
                m = (j * w) % ;
                if (j ==  && jj == ) {
                    m = (jj * w /  )% ;
                    change1(i, j, m);
                } else
                    change2(i, m, j, jj);
            }
        }
        w *= ;
    }
}

int cal(int n1) {
    if (n1 == ) return ;
    char ch[];
    Zero(ch);
    int i = ;
    int flag = -;
    int ret = ;
    while (n1) {
        ch[i++] = n1 %  + '';
        n1 /= ;
    }
    for (int j = i - ; j > ; --j) {
        if (ch[j] == '' && ch[j - ] == '') {
            flag = j - ;
            break;
        }
    }
    n1 = n / ;
    int w = ;
    for (int j = ; j <= i; ++j) {
        for (int jj = ; jj < ch[j] - ''; ++jj) {
            int t = (n1 *  ) * w % ;
            if (jj ==  && ch[j + ] == '') {
                t = (n1 /  * ) * w % ;
                ret += dp[j][][][( - t) % ] + dp[j ][][][( - t) % ];
                continue;
            }
            if (flag > j) {
                ret += dp[j][jj][][( - t) % ] + dp[j][jj][][( - t) % ];
            } else {
                ret += dp[j][jj][][( - t) % ];
            }
        }
        w *= ;
        n1 /= ;
    }
    if (n %  ==  &&  flag != -) ret++;
    return ret;
}
int main() {
    //freopen("data.out", "w", stdout);
    //freopen("data.in", "r", stdin);
    //cin.sync_with_stdio(false);
    init();
    while (cin >> n) {
        cout << cal(n) << endl;
    }
    return ;
}