【cf】55d beautiful numbers【精妙的數位dp+離散化】
阿新 • • 發佈:2019-01-04
題意:
求1-n中,能被所有組成他的非0數整除的數的個數
題解:
題意很簡單,但是考慮到數位dp的狀態轉移,dfs(pos,pre,status,limit)如果要記錄除以他所有數那麼勢必要用陣列來存,但這樣一來,很難用dp陣列來記錄當前狀態,這題非常巧妙的採用了lcm的方法,避免了該類操作,2-9的lcm是2520,這樣我們只要記錄前面數字的lcm和前面數字模2520剩下的數,那麼既保證了不會漏掉前面的狀態,又保證了當前狀態的可記錄性,但是這樣dp陣列要開到[20][2520][2520]還是會接近爆記憶體,我們考慮到可以把lcm離散化,因為2-9的lcm在2520區間中總共只有48個,所以dp陣列可以壓縮成dp[20][50][2530],這樣這個問題就完美的解決了。
#include<cstdio> #include<cstring> using namespace std; #define ll __int64 ll dp[30][50][2530]; int Hash[2530],bit[30]; ll b,a; int T; ll gcd(ll a,ll b) { return b==0?a:gcd(b,a%b); } ll lcm(ll a,ll b) { return a/gcd(a,b)*b; } ll dfs(int pos,int prelcm,int prenum,int limit) { if(pos==0) return prenum%prelcm==0; if(!limit&&dp[pos][Hash[prelcm]][prenum]!=-1)return dp[pos][Hash[prelcm]][prenum]; int num=limit?bit[pos]:9; ll ret=0; for(int i=0;i<=num;i++){ int prelcm1=prelcm; if(i)prelcm1=lcm(prelcm,i); int prenum1=(prenum*10+i)%2520; ret+=dfs(pos-1,prelcm1,prenum1,limit&&(i==num)); } if(!limit)dp[pos][Hash[prelcm]][prenum]=ret; return ret; } ll solve(ll n) { int len=0; while(n){ bit[++len]=n%10; n/=10; } return dfs(len,1,0,1); } int main() {//freopen("C:\\Users\\Administrator\\Desktop\\input.txt","r",stdin); scanf("%d",&T); int cnt=0; for(int i=1;i<=2520;i++)if(2520%i==0)Hash[i]=++cnt; memset(dp,-1,sizeof(dp)); while(T--){ scanf("%I64d%I64d",&a,&b); printf("%I64d\n",solve(b)-solve(a-1)); } return 0; }