[HDU3652]B-number && 数位DP

状态略复杂 d[i, j, k, f]表示i位数 首位为j的数字mod13 == f的数的种类

 k == 1 包含 13 

 k == 0 不包含13

计数的时候要注意如果输入的n某两位是13那么后面的所有数都可以加进来

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int MAXN = 15;
const int p = 13;
int d[MAXN+10][MAXN+10][p+10][2];
int x, Cur[MAXN+10], bit[MAXN+10], Len;
char s[MAXN+10];
void init()
{
	bit[0] = 1; for(int i = 1; i <= MAXN; i++) bit[i] = (bit[i-1] * 10) % p;
	for(int i = 0; i < p; i++) d[1][i][i][0] = 1;
	for(int i = 2; i <= MAXN; i++)
	{
		for(int j = 0; j <= 9; j++)
			for(int k = 0; k <= 9; k++)
				for(int mod = 0; mod < p; mod++)
				{
					int MOD = (j * bit[i-1] + mod) % p;
					if(j == 1 && k == 3)
						d[i][j][MOD][1] += d[i-1][k][mod][0];
					d[i][j][MOD][1] += d[i-1][k][mod][1];
					if(!(j == 1 && k == 3))
						d[i][j][MOD][0] += d[i-1][k][mod][0];
				}
	}
}
int main()
{
	init();
	while(scanf("%d", &x) == 1)
	{
		x++; sprintf(s, "%d", x);
		Len = strlen(s); 
		memset(Cur, 0, sizeof(Cur));
		for(int i = 0; i < Len; i++) Cur[Len-i] = s[i] - '0';
		int ans = 0, mod = 0;
		bool f = false;
		for(int i = Len; i >= 1; i--)
		{
			for(int j = 0; j < Cur[i]; j++)
			{
				ans += d[i][j][mod][1];
				if(f || (Cur[i+1] == 1 && j == 3)) ans += d[i][j][mod][0];
			}
			if(Cur[i+1] == 1 && Cur[i] == 3) f = true;
			mod = (mod - (bit[i-1] * Cur[i]) % p + p) % p;
		}
		printf("%d\n", ans);
	}
}