[hihoCoder]#1033 : 交错和

描述
给定一个数 x，设它十进制展从高位到低位上的数位依次是 a0, a1, …, an - 1，定义交错和函数：

例如：

给定 l, r, k，求在 [l, r] 区间中，所有 f(x) = k 的 x 的和，即：

输入
输入数据仅一行包含三个整数，l, r, k(0 ≤ l ≤ r ≤ $10^{18}$, |k| ≤ 100)。
输出
输出一行一个整数表示结果，考虑到答案可能很大，输出结果模 $10^9$+ 7。
提示
对于样例 ，满足条件的数有 110 和 121，所以结果是 231 = 110 + 121。
更多样例：

样例输入
100 121 0
样例输出
231
数位DP，表示我一开始看错题目了，打了一个错误的代码。但还是有不少东西是一样的，所以我就直接在源程序上改了一下，因此代码很丑，一般人根本不知道在干什么（有几个数组和操作完全是为了填坑Q$\omega$Q）。反正只是为了水一下，那么我就不具体讲了。。。

#include <cstdio>
using namespace std;
#define ll long long 
const ll Mod=1000000007;
ll l,r,k,dp[20][20][420],limits[20],ans1,ans2,op[19],b[20][20][420],t[20][20][420];
ll dfs(int time,ll&ans,int len,int now,bool start,bool limit,int w)
{
    if(len==1&&w==200)
        return ans=(ans+now)%Mod,1;
    else 
        if(len==1)
            return 0; 
    if(!start)
    {
        ll sum=0,x;
        if(limit)
        {
            x=dfs(time,ans,len-1,limits[len-1],1,1,400+limits[len-1]-w);
            sum=(sum+x)%Mod;
            for (int i=1;i<limits[len-1];++i)
            {
                x=dfs(time,ans,len-1,i,1,0,400+i-w);
                sum=(sum+x)%Mod;
            }
            sum=(sum+dfs(time,ans,len-1,0,0,0,w))%Mod;
        }
        else 
        {
            for (int i=1;i<=9;++i)
                sum=(sum+dfs(time,ans,len-1,i,1,0,400+i-w))%Mod;
            sum=(sum+dfs(time,ans,len-1,0,0,0,w))%Mod;
        }
        return sum;
    }
    if(limit)
    {
        ll sum=0;
        for (int i=0;i<limits[len-1];++i)
            sum=(sum+dfs(time,ans,len-1,i,1,0,400+i-w))%Mod;
        sum=(sum+dfs(time,ans,len-1,limits[len-1],1,1,400+limits[len-1]-w))%Mod;
        ans=(ans+(now*op[len-1]%Mod*sum%Mod)%Mod)%Mod;
        return sum;
    }
    else 
    {
        if(b[len][now][w]!=time)
        {
            dp[len][now][w]=0;
            int x=ans;
            for (int i=0;i<=9;++i)
                dp[len][now][w]=(dp[len][now][w]+dfs(time,ans,len-1,i,1,0,400+i-w))%Mod;
            t[len][now][w]=((ans-x)%Mod+Mod)%Mod;
            b[len][now][w]=time;
        }
        else 
            ans=(ans+t[len][now][w])%Mod;
        ans=(ans+(now*op[len-1]%Mod*dp[len][now][w]%Mod)%Mod)%Mod;
        return dp[len][now][w];
    }
}
void solve(ll &ans,ll x,int time)
{
    ll a=x,len=0;
    while(a)
        limits[++len]=a%10,a/=10;   
    dfs(time,ans,len+1,0,0,1,k+200);
    return ;
}
int main(void)
{
    scanf("%lld%lld%d",&l,&r,&k);
    op[0]=1;
    for (int i=1;i<=18;++i)
        op[i]=op[i-1]*10%Mod;
    solve(ans1,r,1),solve(ans2,l-1,2);
    printf("%lld",(((ans1-ans2)%Mod)+Mod)%Mod);
    return 0;
}