bzoj1833: [ZJOI2010]count 数字计数&&USACO37 Cow Queueing 数数的梦（数位DP）

难受啊，怎么又遇到我不会的题了（捂脸）

如题，这是一道数位DP，随便找了个博客居然就是我们大YZ的……果然nb，然后就是改改模版++注释就好的了，直接看注释吧，就是用1~B - 1~A-1而已，枚举全部位然后判一下是不是上限边缘和前导零就OK

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
LL bit[];

int len,a[];
LL as[],ret[];//ret表示假如是上限边缘，后面有多少种填的方式
LL f[][][];
//f[k][fg][z]表示当前这个数x，搞到第k位的时候的数目
//fg表示该位是不是在上限边缘 z表示目前是否还是前导零
LL dfs(int x,int k,int fg,int z)
{
    if(k==len+)return ;
    if(f[k][fg][z]!=-)return f[k][fg][z];

    LL ans=;
    if(fg==)//如果前面的部分处于上限边缘
    {
        for(int i=;i<a[k];i++)//不填上限边缘，枚举填什么
        {
            if(z==&&i==)ans+=dfs(x,k+,,);
            else
            {
                ans+=dfs(x,k+,,);
                if(i==x)ans+=bit[len-k];
            }
            //如果i是要统计次数的那个数字，而且不是在前导零的时候，而且还不是在上限的边缘
            //那么后面几位的数字有多少种填法x就出现了几次
        }
        //该位填上限边缘
        if(z==&&a[k]==)ans+=dfs(x,k+,,);
        else
        {
            ans+=dfs(x,k+,,);
            if(a[k]==x)ans+=ret[k];
        }
        //如果填的是要统计次数的那个数字 而且不是在前导零的时候
        //那么后面的数字最多有多少种填法x就出现了几次
    }
    else//没有限制
    {
        for(int i=;i<=;i++)//同上
        {
            if(z==&&i==)ans+=dfs(x,k+,,);
            else
            {
                ans+=dfs(x,k+,,);
                if(i==x)ans+=bit[len-k];
            }
        }
    }
    f[k][fg][z]=ans;
    return ans;
}
char s[];
void cl()//把A减1
{
    int t=len;
    while(t>&&s[t]==''){s[t]='';t--;}
    s[t]--;
    if(s[t]==)
    {
        len--;
        for(int i=;i<=len;i++)s[i]=s[i+]-'';
    }
}
int main()
{
    freopen("dream.in","r",stdin);
    freopen("dream.out","w",stdout);
    bit[]=;for(int i=;i<=;i++)bit[i]=bit[i-]*;
    memset(as,,sizeof(as));

    scanf("%s",s+);len=strlen(s+);cl();
    for(int i=;i<=len;i++)a[i]=s[i]-'';
    //0的情况，所以填的方案都要+1
    ret[len]=;for(int i=len-;i>=;i--)ret[i]=a[i+]*bit[len-i-]+ret[i+];
    for(int i=;i<=;i++)
    {
        memset(f,-,sizeof(f));
        as[i]-=dfs(i,,,);
    }

    scanf("%s",s+);len=strlen(s+);
    for(int i=;i<=len;i++)a[i]=s[i]-'';

    ret[len]=;for(int i=len-;i>=;i--)ret[i]=a[i+]*bit[len-i-]+ret[i+];
    for(int i=;i<=;i++)
    {
        memset(f,-,sizeof(f));
        as[i]+=dfs(i,,,);
    }

    for(int i=;i<;i++)printf("%d ",as[i]);
    printf("%d\n",as[]);
    return ;
}