Codeforces - 55 D Beautiful numbers (数位 dp +状态压缩)

链接 ：Beautiful numbers

题意：
求 a到 b 中有多少个数满足 - 这个数能被它除0以外的每一个数位整除。
思路：

1. 感觉这题挺不错的，状态有点难搞，首先 9 个数 的 lcm 为 2520，一开始想的状态是到第 i 位 该数对 2520 的余数mod和当前所选数的 lcm 。如果最后 mod%lcm=0，就说明这个数是满足的。但这样要开的空间是 20 * 2520 * 2520.然后就 。。。。。
   
2. 然后右想到可以用 二进制位 表示选了哪些数，把每种状态的 lcm 预处理出来就好了。这样就变成了 20 * 2520 *1024,emmmm,这不是差不多，还是一样的结局 （哭了）。
3. 然后又发现 0 和 1 好像对状态没什么影响（选它们和没选是一样的），那么就可以只要 8 个二进制位 就可以把所有的状态存下来，就是 20 * 2520 * 256 ,这下算是可以了。
4. 知道了状态就好写了。

代码：

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e4+7;
ll dp[20][2522][260],T,Lcm[2000],x[20];
ll a,b;
int gcd(int a,int b){
    if(b==0) return a;
    return gcd(b,a%b);
}
int lcm(int a,int b){
     return a*b/gcd(a,b);
}
void init(){
     memset(dp,-1,sizeof(dp));
     for(int i=0;i<(1<<8);i++){        //预处理 每个状态的 lcm
         int s=1;
         for(int j=0;j<8;j++){
             if((i>>j)&1) s=lcm(s,j+2);
         }
         Lcm[i]=s;
     }
}
ll dfs(int pos,int mod,int bit,int limit){
    if(pos<0){
       return mod%Lcm[bit]==0;
    }
    if(!limit&&dp[pos][mod][bit]!=-1) return dp[pos][mod][bit];
    int up=limit ? x[pos] : 9;
    ll tmp=0;
    for(int i=0;i<=up;i++){
        int bitt=bit,modd;
        if(i>=2) bitt=bit|(1<<(i-2));         //把第 i-2位变成 1.
        modd=(mod*10+i)%2520;
        tmp+=dfs(pos-1,modd,bitt,limit&&x[pos]==i);
    }
    if(!limit) dp[pos][mod][bit]=tmp;
    return tmp;
}
ll solve(ll s){
    int pos=0;
    while(s>0){
        x[pos++]=s%10;
        s/=10;
    }
    return dfs(pos-1,0,0,1);

}
int main (){
    cin>>T;
    init();
    while(T--){
          cin>>a>>b;
          cout<<solve(b)-solve(a-1)<<endl;
    }
}