[bzoj 1799--AHOI2009]self 同类分布

给出a,b，求出[a,b]中各位数字之和能整除原数的数的个数。

这道题一看就是数位dp，但是这题有难度，因为这题牵扯到了两个东东（各位数字之和，原数），所以要固定一个量，所以要枚举各位数字之和。定义f[i][j][k]，i为枚举到的位数（从后往前），j一开始为枚举到的各位数字之和，然后减去枚举到的数值，到最后如果为0，则为满足条件，k为当前的数值%各位数字之和的值，如果最后为0，则为满足情况。然后就是套路了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
long long a[22];
long long s[22],f[22][163][163];
long long dfs(int pos,long long ss,long long sum,long long mod,bool limt)
{
    if(pos==0)
    {
        if(sum==0 && ss==0)return 1;
        else return 0;
    }
    if(limt==false && f[pos][ss][sum]!=-1)return f[pos][ss][sum];
    int up=9;long long ans=0;
    if(limt==true)up=a[pos];
    for(int i=0;i<=up;i++)
    {
        if(ss<i)break;
        long long s1=ss,s2=sum;bool bk=false;
        s1-=i;s2=(s2+i*s[pos])%mod;
        if(limt==true && i==a[pos])bk=true;
        ans+=dfs(pos-1,s1,s2,mod,bk);
    }
    if(limt==false)f[pos][ss][sum]=ans;
    return ans;
}
long long solve(long long x)
{
    int pos=0,sum=0;
    while(x)
    {
        a[++pos]=x%10;
        sum+=x%10;
        x/=10;
    }
    long long ans=0;
    for(int i=1;i<=pos*9;i++)
    {
        memset(f,-1,sizeof(f));
        ans+=dfs(pos,i,0,i,true);
    }
    return ans;
}
int main()
{
    s[1]=1;
    for(int i=2;i<=18;i++)s[i]=s[i-1]*10;
    long long n,m;
    scanf("%lld%lld",&n,&m);
    printf("%lld\n",solve(m)-solve(n-1));
    return 0;
}