J - Junior Mathematician（数位DP)

链接：https://codeforces.com/gym/102452/problem/J
题意:
给你L,R，m，问

满足这个公式的有多少个数。
L< R<10^5000；
其中k为x的位数。
思路：
数位DP,
$dp[pos][wsum][f[x]-x$%$m]$
pos表示当前位数。
wsum表示之前位数和
f[x]-x%m表示差值，为0就计算一次。
这样开数组不会炸掉。如果把f[x]和x单独开个数组也可以计算，但是数组会炸。
dfs(pos，limit，wsum，fsum)

注意两点：
1，左端点要减一。
2，计算答案取mod的时候不能放在循环内，不然会T35，要放到循环外。

#include <algorithm>
#include <iostream>
#include <map>
#include <math.h>
#include <queue>
#include <set>
#include <stack>
#include <stdio.h>
#include <string.h>
#include <string>
#include <vector>
#define emplace_back push_back
#define pb push_back
using namespace std;
typedef long long LL;
const int mod = 1e9 + 7;
const int N = 5e3 + 10;
char s1[N], s2[N];
int dp[N][61][61];
int m;
int p[N];
int a[N];
int dfs(int pos, int limit, int wsum, int fsum)
{
    if (!pos)
        return fsum == 0;
    if (!limit && (dp[pos][wsum][fsum] != -1))
    {
        return dp[pos][wsum][fsum];
    }
    int up = limit ? a[pos] : 9;
    LL val = 0;
    for (int i = 0; i <=up;i++)
    {
        int _fsum=((fsum+wsum*i-i*p[pos-1])%m+m)%m;
        val += dfs(pos - 1, limit && i == a[pos],(wsum +i)%m,_fsum);
    }
    val %= mod;
    if(!limit)
    {
        dp[pos][wsum][fsum] = val;
    }
    return val;
}
int solve(char *s, int l)
{
    for (int i = 0; i < l; i++)
    {
        a[l - i] = s[i] - '0';
    }
    return dfs(l, 1, 0, 0);
}
int main()
{
    int t;
    scanf("%d", &t);
    while (t--)
    {
        scanf("%s%s", s1, s2);
        scanf("%d", &m);
        int l = strlen(s1);
        for (int i = l - 1; i >= 0; i--)
        {
            if (s1[i] >= '0')
            {
                s1[i]--;
                break;
            }
            else
            {
                s1[i] = '9';
            }
        }
        int l2 = strlen(s2);
        p[0] = 1;
        for (int i = 1; i <= l2; i++)
        {
            p[i] = p[i - 1] * 10 % m;
        }
        for (int i = 1; i <= l2; i++)
        {
            for (int j = 0; j <= m; j++)
            {
                for (int k = 0; k <= m; k++)
                {
                    dp[i][j][k] = -1;
                }
            }
        }
        printf("%d\n", (solve(s2, l2) - solve(s1, l) + mod) % mod);
    }
    return 0;
}