【题解】AHOI2009同类分布

　　好开心呀~果然只有不看题解做出来的题目才会真正的有一种骄傲与满足吧ヾ(๑╹◡╹)ﾉ"

　　实际上这题只要顺藤摸瓜就可以了。首先按照数位dp的套路，有两维想必是省不掉：1.当前dp到到的位数；2.0/1状态表示是否受限制（这一条是因为有数字上限）。然后根据这两个维度来接着往下想。第二个维度先撇开不看，我们只考虑如何从第 \(i - 1\) 位dp到第 \(i\) 位。在这里其实卡了有点久，因为如果除数与被除数都在改变，那么两维的转移是非常凉凉的。

　　这个时候联想题目的特殊性质 ----- 当感觉无法优化转移 / 转移方式的时候，考虑状态的重新设计 & 题目的特别要求。然后很开心的发现：\(1e18\) 实际上各位数字的和最大都只有 \(162\)。那么岂不是乱搞也可以？所以我们固定除数 \(Q\) 为 \(\left ( 1, 162 \right )\) 当中的任意一个数，分别进行dp即可。此时的转移就简单了，因为除数固定，自然地追加一维表示余数。状态固定为 \(f[i][j][k][L]\)，表示dp到第 \(i\) 位，要求第 \(\left ( 1, i \right )\) 位的数字之和加起来为 \(j\)，且原数除以 \(Q\) 的余数为 \(k\)，限制为\(L\left ( 0, 1 \right )\)的总个数。

　　感觉这份代码写的还行，跑得也还行……能看。

#include <bits/stdc++.h>
using namespace std;
#define int long long
int a[], Res, mul[];
int f[][][][];

int read()
{
    int x = , k = ;
    char c;
    c = getchar();
    while(c < '' || c > '') { if(c == '-') k = -; c = getchar(); }
    while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
    return x * k;
}

#define Pre f[now][tot][rm][lim]
int DP(int now, int tot, int rm, bool lim)
{
    if(~Pre) return Pre;
    else Pre = ;
    if(tot > now * ) return ;
    if(now == )
    {
        if(tot >  || (tot > a[now] && lim)) return Pre = ;
        return Pre = ((Pre = tot % Res == rm) ?  : );
    }
    for(int i = ; i <= ; i ++)
    {
        if(i > a[now] && lim) break;
        if(tot < i) break;
        int q = (mul[now] * i) % Res; q = q % Res;
        int L = (i == a[now] && lim);
        f[now][tot][rm][lim] += DP(now - , tot - i, (rm - q + Res) % Res, L);
    }
    return f[now][tot][rm][lim];
}
#undef Pre

int Solve(int x)
{
    int k = x, ans = , num = ;
    while(k) num ++, a[num] = k % , k /= ;
    for(int i = ; i <= ; i ++)
    {
        Res = i; if(i > num * ) continue;
        memset(f, -, sizeof(f));
        ans += DP(num, i, , );
    }
    return ans;
}

signed main()
{
    int a = read(), b = read(); mul[] = ;
    for(int i = ; i <= ; i ++) mul[i] = mul[i - ] * ;
    printf("%lld\n", Solve(b) - Solve(a - ));
    return ;
}