bzoj1799: [Ahoi2009]self 同类分布

传送门
因为模数最多只有9*17=163，所以我们暴力悲剧模数。
然后就是dp辣。
设f[0/1][i][j][k]表示填了前i位，和是j，数字对mo取模的值是k，前边的数字是否前部和上限相同（0）的方案数。
xjb转移一下就行了。
开心的刷到第一页。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#define ll long long
using namespace std;
ll f[2][25][170][170],x,y;
int vi[2][25][170][170],a[25],cnt,mo,tot;
ll dp(int i,int j,int k,int l){
    if (!j) return !(k|l);
    if (vi[i][j][k][l]==tot) return f[i][j][k][l];
    vi[i][j][k][l]=tot;
    ll t=0;
    int le=max(0,k-(j-1)*9),r=min((i)?9:a[j],k);
    for (int p=le;p<=r;p++) t+=dp(i|(p<a[j]),j-1,k-p,(l*10+p)%mo);
    return f[i][j][k][l]=t;
}
ll get(ll x){
    ll t=0;
    for (cnt=0;x;x/=10) a[++cnt]=x%10;
    for (mo=1;mo<=cnt*9;mo++){
        tot++;
        t+=dp(0,cnt,mo,0);
    }
    return t;
}
int main(){
    scanf("%lld%lld",&x,&y);
    printf("%lld\n",get(y)-get(x-1));
}