[bzoj 5064]B-number

王女士是一个数学爱好者，这些天她对B数很感兴趣。 B数的定义：能被13整除且本身包含字符串”13”的数。例如：130和2613是B数，但是143和2639不是B数。 你的任务是计算1到n之间有多少个数是B数。

这道题可以算是一道较水的数位dp，可以说是比较套路了，把握好f数组就好了。定义f[i][j][k1][k2],i为枚举到的位数（从后往前），j为枚举到的前一个的数，k1表示这个数是否包含“13”，k2为目前枚举到的数部分%13的数。之后就不详解了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
int a[20];
long long f[18][10][2][15],s[18];
long long dfs(int pos,int sta,long long ss,bool limt,bool lead)
{
    if(pos==0)
    {
        if(ss==0 && lead==true)return 1;
        else return 0;
    }
    if(limt==false && f[pos][sta][lead][ss]!=-1)return f[pos][sta][lead][ss];
    int up=9;long long ans=0;
    if(limt==true)up=a[pos];
    for(int i=0;i<=up;i++)
    {
        bool bk1=false,bk2=lead;
        if(limt==true && i==a[pos])bk1=true;
        if(sta==1 && i==3)bk2=true;
        ans+=dfs(pos-1,i,(ss+i*s[pos])%13,bk1,bk2);
    }
    if(limt==false)f[pos][sta][lead][ss]=ans;
    return ans;
}
long long solve(long long x)
{
    int pos=0;
    while(x)
    {
        a[++pos]=x%10;
        x/=10;
    }
    memset(f,-1,sizeof(f));
    return dfs(pos,0,0,true,false);
}
int main()
{
    s[1]=1;
    for(int i=2;i<=15;i++)s[i]=s[i-1]*10;
    long long n;
    scanf("%lld",&n);
    printf("%lld\n",solve(n));
    return 0;
}