1. 程式人生 > >【cf】55d beautiful numbers【精妙的數位dp+離散化】

【cf】55d beautiful numbers【精妙的數位dp+離散化】

題意:

求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;
}