bzoj 3598 [ Scoi 2014 ] 方伯伯的商场之旅 ——数位DP

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=3598

数位DP...东看西看：http://www.cnblogs.com/Artanis/p/3751644.html

　　　　　　　　    https://www.cnblogs.com/MashiroSky/p/6399095.html

好巧妙的思路啊！这样统计的东西就变得很简单了；

好美的 dfs！数位DP用 dfs 好像能变得很清楚。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll l,r,f[][],ans;
int n,a[],K;
ll dfs1(int pos,int s,bool lim)
{
    if(pos==)return s;
    if(!lim&&f[pos][s]!=-)return f[pos][s];
    int end=K-; ll ret=;
    if(lim)end=a[pos];
    for(int i=;i<=end;i++)
        ret+=dfs1(pos-,s+i*(pos-),lim&&(i==end));
    if(!lim)f[pos][s]=ret;//!lim!
    return ret;
}
ll dfs(int pos,int s,int m,bool lim)
{
    if(s<)return ;//!
    if(pos==)return s;
    if(!lim&&f[pos][s]!=-)return f[pos][s];
    int end=K-; ll ret=;
    if(lim)end=a[pos];
    for(int i=;i<=end;i++)
    {
        if(pos>=m)ret+=dfs(pos-,s+i,m,lim&&(i==end));
        else ret+=dfs(pos-,s-i,m,lim&&(i==end));
    }
    if(!lim)f[pos][s]=ret;
    return ret;
}
ll calc(ll x)
{
    int n=;
    while(x)a[++n]=x%K,x/=K;
    memset(f,-,sizeof f);
    ll ret=dfs1(n,,);
    for(int i=;i<=n;i++)
    {
        memset(f,-,sizeof f);
        ret-=dfs(n,,i,);
    }
    return ret;
}
int main()
{
    scanf("%lld%lld%d",&l,&r,&K);
    printf("%lld\n",calc(r)-calc(l-));
    return ;
}