uva 11361 Investigating Div-Sum Property

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2346

　　不知道算什么题，好像是数位dp，不过在切递推的题，所以这题应该算是递推。

　　首先要预处理出dp[i][j][k][p]，表示一共有i位数字，其中数字和mod p余j而且这个数mod p余k的个数。然后就是枚举每一位，分段处理出结果。

代码如下：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;

typedef long long LL;
const int N = 12;
const int M = 55;
LL dp[N][M][M][M], ep[N];

void PRE() {
    for (int i = 0; i < M; i++) dp[0][0][0][i] = 1;
    for (int i = 0; i < N - 1; i++) {
        for (int j = 0; j < M; j++) {
            for (int k = 0; k < M; k++) {
                for (int p = 1; p < M; p++) {
                    for (int q = 0; q < 10; q++) {
                        dp[i + 1][(j + q) % p][(k * 10 + q) % p][p] += dp[i][j][k][p];
                    }
                }
            }
        }
    }
    ep[0] = 1;
    for (int i = 0; i < N - 1; i++) ep[i + 1] = ep[i] * 10;
    //cout << dp[1][0][0][1] << endl;
    //cout << dp[1][0][0][2] << endl;
    //cout << dp[1][0][0][3] << endl;
    //cout << dp[3][1][1][4] << endl;
}

LL cal(int x, int k) {
    int p = 0, dgs = 0;
    if (k >= M) return 0;
    LL ret = 0;
    while (x >= ep[p]) p++;
    while (p > 0) {
        p--;
        int cdg = x / ep[p] % 10;
        for (int i = 0; i < cdg; i++) {
            ret += dp[p][(k - (dgs + i) % k) % k][(k - (x / ep[p + 1] * 10 + i) * ep[p] % k) % k][k];
        }
        dgs += cdg;
    }
    return ret + dp[0][(k - dgs % k) % k][(k - x % k) % k][k];
}

int main() {
    PRE();
    int T, n, m, k;
    cin >> T;
    while (T-- && cin >> n >> m >> k) cout << cal(m, k) - cal(n - 1, k) << endl;
    return 0;
}

——written by Lyon