【数位DP】CF55D Beautiful numbers

$dp[x][p][pp]$表示第x位，当前已有数字mod 2520（1~9数字的lcm）为p，当前各位数字的lcm为pp

观察到数组太大，考虑压缩，第三维lcm最多只有9个数字，打表发现最多只有48个状态，压掉第三维即可

打表用一个状压然后set维护（广搜也可以）即可

有一个坑点：题目里似乎没有说关于0的事情（即数字里出现0）但是有人在CF上打这个比赛的时候问了出题人，碰到0不要管即可！！！

打表代码：

 set<int>s;
 inline void Make(int x){
     int ans=;
     for(int i=;i<;i++){
         if(((<<i)&x)) ans=lcm(ans,i+);
     }s.insert(ans);
 }
 inline void States_Maker(){
     int ans=;
     for(int i=;i<(<<);i++)Make(i);
     while(!s.empty()){
         cout<<*s.begin()<<",";
         s.erase(s.begin());
     }
 }

解题代码：

 #include<bits/stdc++.h>
 #define int long long
 #define writeln(x)  write(x),puts("")
 #define writep(x)   write(x),putchar(' ')
 using namespace std;
 inline int read(){
     int ans=,f=;char chr=getchar();
     while(!isdigit(chr)){if(chr=='-') f=-;chr=getchar();}
     while(isdigit(chr)){ans=(ans<<)+(ans<<)+chr-;chr=getchar();}
     return ans*f;
 }void write(int x){
     if(x<) putchar('-'),x=-x;
     if(x>) write(x/);
     putchar(x%+'');
 }const int M = ,State = ;
 int f[M][][State],T,n,a[M],l,r,Pos[],S[State]={,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,};
 inline int gcd(int x,int y){return __gcd(x,y);}
 inline int lcm(int x,int y){if(x==)return y;if(y==)return x;return x*y/gcd(x,y);}
 int dfs(int x,int p,int pp,bool lim){
     if(!x)return (p%pp==)?():();
     if(!lim&&p!=-&&pp!=-&&f[x][p][Pos[pp]]!=-)return f[x][p][Pos[pp]];
     int ans=,up=lim?a[x]:;
     for(int i=;i<=up;i++)ans+=dfs(x-,(p*+i)%,lcm(pp,i),lim&&a[x]==i);
     if(!lim)f[x][p][Pos[pp]]=ans;
     return ans;
 }
 inline int Solve(int x){
     int len=;
     while(x)a[++len]=x%,x/=;
     return dfs(len,,,);
 }
 signed main(){
     memset(f,-,sizeof(f));
     for(int i=;i<=;i++)Pos[S[i]]=i;
     T=read();
     while(T--)l=read(),r=read(),writeln(Solve(r)-Solve(l-));
     return ;
 }