动态规划之数位dp

数位dp

之所以叫数位dp，是因为它在记数上经行dp操作，但数位dp是很暴力的，但是比一般的纯暴力快的多，因为它有记忆化的操作。
这篇数位dp文章还是很好的：数位dp总结 之 从入门到模板

关于数位dp的题

，通常会给你一个数，然后给你一些要求，让你求从1~n中符合要求的数有多少个。
比如这道：HDU 2089 不要62
题意就是给你两个数l,r。然后让你求l到r的区间里既没有出现4,也没有出现62（连续62）的数有多少个。
思路：
我们枚举每位上的数的时候，不枚举有4和62的情况就可以了，4很好避免，但是62有两位数，那么我们可以开一个二维的dp数组，用其中一个维度记录前一位是不是6这种状态，要是的话，不枚举2出现的情况，如果不是6，那么就无所谓了。
AC:

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
int dp[20][2];
int a[20];
int dfs(int pos,int state,int limit,int pre)
{
    if(pos==-1)return 1;
    if(!limit&&dp[pos][state]!=-1)return dp[pos][state];
    int tem=0;
    int up=limit?a[pos]:9;
    for(int i=0; i<=up; i++)
    {
        if(i==4)continue;
        if(i==2&&pre==6)continue;
        tem+=dfs(pos-1,i==6,a[pos]==i&&limit,i);
    }
    if(!limit)dp[pos][state]=tem;
    return tem;
}
int solve(int n)
{
    int p=0;
    while(n)
    {
        a[p++]=n%10;
        n/=10;
    }
    return dfs(p-1,0,1,-1);
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        if(n==0&&m==0)
        {
            break;
        }
        memset(dp,-1,sizeof(dp));

        printf("%d\n",solve(m)-solve(n-1));

    }

    return 0;
}

A - Beautiful numbers

CodeForces - 55D
题目大意
给你两个数字l,r表示左右两个端点位置，然后让你求在[l,r]这个范围内有多少个数满足这个要求（可以整除它的所有数位上的数（不包括0））。
思路
题意很简单，但是却不好做，我们先要找到所有1~9数字的最小公倍数2520；然后dp数组开三维，分别表示数位，余数，当前所有数最小公倍数。但是dp[22][2522][2522]会不出意外的爆掉内存，但是总共最小公倍数就那么48个，所以hash一下，把第三个维度降到50，就可以了。
AC

#include<stdio.h>
#include<iostream>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
using namespace std;
long long dp[20][2522][50];
long long  a[20];
long long has[2550];
int gcd(int x,int y)
{
    return y==0?x:gcd(y,x%y);
}
long long dfs(int pos,int yu,int b,int limit)
{
    if(pos==-1)
    {
        if(yu%b==0)
        {
            return 1;
        }
        else
        {
            return 0;
        }
    }
    if(!limit&&dp[pos][yu][has[b]]!=-1)
        return dp[pos][yu][has[b]];
    int up=limit?a[pos]:9;
    long long  sum=0;
    for(int i=0; i<=up; i++)
    {
        int p=b;
        if(i>=2)
        {
            p=(p*i)/gcd(p,i);
        }
        sum+=dfs(pos-1,(yu*10+i)%2520,p,limit&&i==a[pos]);
    }
    if(!limit)
    {
        dp[pos][yu][has[b]]=sum;
    }
    return sum;
}
long long solve(long long x)
{
    int p=0;
    while(x)
    {
        a[p++]=x%10;
        x/=10;
    }
   return dfs(p-1,0,1,1);
}
int main()
{
    long long n,m;
    int t;
    for(int i=1,j=0; i<=2520; i++)
    {
        if(2520%i==0)
        {
            has[i]=j++;

        }
    }
    cin>>t;
    while(t--)
    {
        cin>>n>>m;
        memset(dp,-1,sizeof(dp));
        cout<<solve(m)-solve(n-1)<<endl;
    }

    return 0;
}