BZOJ 1799 [Ahoi2009]self 同类分布【数位dp

枚举数字总和为多少（设为mod），f[i][j][k][0..1]表示算到第i位，前i位的和为j，这个数的数值%mod=k，是否顶上界

转移就是f[i+1][j+x][(k*10+x)%mod] += f[i][j][k]

这样的（顶上界的特判一下嘛，上面那个公式就省略0..1那一维了【滑稽

手癌把i打成了x，调了一年2333

#include<bits/stdc++.h>
using namespace std; long long L,R;
long long f[20][200][200][2];
long long Lim = 162;

char s[25];
inline long long dp(long long a, int mod){
if(a<=0) return 0;
memset(f,0,sizeof f);
memset(s,0,sizeof s);
sprintf(s+1,"%lld",a);
f[0][0][0][1] = 1;
int lth = strlen(s+1);
for(register int i=0;i<lth;++i){
for(register int j=0;j<=mod;++j){
for(register int k = 0;k<mod;++k){
if(!f[i][j][k][0] && !f[i][j][k][1]) continue;
if(f[i][j][k][1]){
f[i+1][j+s[i+1]-'0'][(k*10 + s[i+1] - '0')% mod][1] = 1;
for(register int x=0;x<s[i+1] - '0' && j+x<=mod;++x)
++f[i+1][j+x][(k*10 + x)%mod][0];
}
for(register int x = 0;x<=9&&j+x<=mod;++x)
f[i+1][j+x][(k*10 +x)%mod][0] += f[i][j][k][0];
}
}
}

return f[lth][mod][0][0] + f[lth][mod][0][1];
}

long long ans = 0;

int main(){
// freopen("1.in","r",stdin);
scanf("%lld%lld",&L,&R);
if(R==1000000000000000000ll) ++ans , --R;
for(register int i = 1;i<=Lim;++i)
ans += dp(R,i) - dp(L-1,i);

printf("%lld",ans);
return 0;
}