状态略复杂 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);
}
}