数位DP HDU-3652

题目大意：计算0~n内（含n）含有子串“13”且可被13整除的数的个数

数位DP   写的记忆化搜索（参照了某大牛的代码）

dp[i][j][k]: i:长度 j:对13的余数 k：末位1/末位非1/含有13；

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include <algorithm>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define ll long long
#define eps 1e-8
#define ms(x,y) (memset(x,y,sizeof(x)))
#define fr(i,x,y) for(int i=x;i<=y;i++)
using namespace std;
int dp[13][13][3];
// 位数 %13 尾1/尾非1/含13
int bit[15];
int dfs(int len,int mod,int is,int lim)
{
    //cout<<len<<' '<<mod<<' '<<is<<' '<<lim<<endl;
    int s=0;
    if(len==0)return (mod==0)&&(is==2);//余数是零且含13
    if((!lim)&&dp[len][mod][is]!=-1)return dp[len][mod][is];
    int num=lim?bit[len]:9;
    fr(i,0,num)
    {
        int t=0;
        if(is==2)t=2;
        else if(is==1&&i==3)t=2;
        else if(i==1)t=1;
        s+=dfs(len-1,(mod*10+i)%13,t,(i==num)&&lim);
    }
    if(!lim)dp[len][mod][is]=s;
    return s;
}
int main()
{
    ms(dp,-1);
    ms(bit,0);
    int n;
    while(cin>>n)
    {
        int l=0;
        while(n)
        {
            bit[++l]=n%10;
            n/=10;
        }
        cout<<dfs(l,0,0,1)<<endl;
    }
}