【HDU 3652】 B-numbers

【题目链接】

点击打开链接

【算法】

数位DP

f[i][j][k][l]表示i位数，第一位为j，除以13的余数为k，是/否包括子串“13”的方案数

当然，我们也可以先打表，然后，对于每次询问，二分即可

【代码】

#include<bits/stdc++.h>
using namespace std;
#define MAXL 13

int i,n;
int a[MAXL],r[MAXL];
int f[MAXL][][][];

inline void dp()
{
    int i,j,k,x,m;
    for (i = ; i <= ; i++) f[][i][i][] = ;
    for (i = ; i < MAXL; i++)
    {
        for (j = ; j <= ; j++)
        {
            for (k = ; k <= ; k++)
            {
                m = (r[i-] * j) % ;
                for (x = ; x <= ; x++)
                {
                    if (j ==  && x == ) continue;
                    f[i][j][k][] += f[i-][x][(k-m+)%][];
                }
                for (x = ; x <= ; x++)
                {
                    f[i][j][k][] += f[i-][x][(k-m+)%][];
                }
                if (j == ) f[i][j][k][] += f[i-][][(k-m+)%][];
            }
        }
    }
}
inline int calc(int n)
{
    int i,j,len = ,ans = ,m = ,flag = false;
    memset(a,,sizeof(a));
    while (n)
    {
        a[++len] = n % ;
        n /= ;
    }
    for (i = len; i >= ; i--)
    {
        if (flag)
        {
            for (j = ; j < a[i]; j++)
            {
                ans = ans + f[i][j][(-m)%][] + f[i][j][(-m)%][];
            }
        } else
        {
            for (j = ; j < a[i]; j++)
            {
                ans += f[i][j][(-m)%][];
            }
            if (a[i+] ==  && a[i] > ) ans += f[i][][(-m)%][];
        }
        m = (m + a[i] * r[i-]) % ;
        if (a[i] ==  && a[i+] == ) flag = true;
    }
    return ans;
}

int main()
{

    r[] = ;
    for (i = ; i < MAXL; i++) r[i] = (r[i-] * ) % ;
    dp();
    while (scanf("%d",&n) != EOF) printf("%d\n",calc(n+));    

    return ;
}