HDU 3652 B-number（数位DP）

#include<stdio.h>
#include<string.h>
#define LL long long
const int MN=100;
LL dp[MN][20][10];//位置，余数，状态，刚开始余数那里数组开小了
int digit[MN];

LL DFS(int pos,int mod,int have,int doing)
{
if(pos==-1) return ((have==2) && (mod==0));
if(!doing && dp[pos][mod][have]!=-1) return dp[pos][mod][have];
int end=doing?digit[pos]:9;
LL ans=0;
for(int i=0;i<=end;i++)
{
int MOD=(mod*10+i)%13;
int nhave=have;
if(have==0 && i==1)
nhave=1;
if(have==1 && i==3)
nhave=2;
if(have==1 && i!=3)
nhave=0;
if(have==1 && i==1)
nhave=1;
ans+=DFS(pos-1,MOD,nhave,(doing && i==end));
}
if(!doing)
{
dp[pos][mod][have]=ans;
}
return ans;
}

LL cal(LL n)
{
int len=0;
while(n)
{
digit[len++]=n%10;
n/=10;
}
memset(dp,-1,sizeof(dp));
return DFS(len-1,0,0,1);
}

int main()
{
int i,j;
LL n;
while(scanf("%I64d",&n)!=EOF)
{
printf("%I64d\n",cal(n));
}
return 0;
}