CodeForces 55D

题意略。

思路:

本题可以说是醉翁之意不在酒了。要解开本题有几个关键点:

1.意识到数X = An An-1 An-2 An-3 .... A2 A1能被{An,An-1,An-2,....,A1}这n个数整除的充要条件是lcm(An,An-1,An-2,....,A1) | X。

2.要知道1~9这9个数进行组合最大的lcm是2520,而且组合出的lcm个数有限,最多48个。

3.我们为了判断这个数能否被它的所有数位整除,我们还需要这个数的值,显然要记录值是不可能的,其实我们只需记录它对2520的模即可。

因为所有的lcm都是2520的因子,所以我们用mod 2520的值来表示当前已经摆好组成的值,到了最后pos == -1的时候,我们再把当前的值mod 当前若干数位组合好的lcm,看看mod lcm 是不是等于0。

如果等于0,则是合格的解。

4.不用考虑前导0的影响。因为在 !limit 的前提下 1111[ ][ ] 与 [ ][ ]是一样的,对当前组合好的lcm不会产生影响。

代码附上:

 

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int maxn0 = 20;
const int maxn1 = 3000;
const int maxn2 = 50;
const int mod = 2520;

int a[maxn0],tail,T;
LL dp[maxn0][maxn2][maxn1],l,r;
int mp[maxn1],store[maxn1],t;

int gcd(int a,int b){
    return b == 0 ? a : gcd(b,a % b);
} 
int lcm(int a,int b){
    int d = gcd(a,b);
    return a / d * b;
}
void process(LL x){
    tail = 0;
    while(x){
        a[tail++] = int(x % 10);
        x /= 10;
    }
}
LL dfs(int pos,int Lcm,int remainder,bool limit){
    if(pos == -1 && remainder % store[Lcm] == 0) return 1;
    if(pos == -1 && remainder % store[Lcm] != 0) return 0;
    if(!limit && dp[pos][Lcm][remainder] != -1) return dp[pos][Lcm][remainder];
    
    int up = limit ? a[pos] : 9;
    LL ret = 0;
    for(int i = 0;i <= up;++i){
        int truelcm = store[Lcm];
        int ntruelcm = i == 0 ? truelcm : lcm(truelcm,i);
        int nlcm = mp[ntruelcm];
        int nremainder = (remainder * 10 + i) % mod;
        ret += dfs(pos - 1,nlcm,nremainder,limit && i == up);
    }
    if(!limit) dp[pos][Lcm][remainder] = ret;
    return ret;
}
void prepare(){
    int total = (1<<9);
    for(int i = 1;i < total;++i){
        int temp = 1;
        for(int j = 0;j < 9;++j){
            if(((i>>j) & 1) == 0) continue;
            temp = lcm(temp,(j + 1));
        }
        store[t++] = temp;
    }
    sort(store,store + t);
    t = unique(store,store + t) - store;
    for(int i = 0;i < t;++i){
        mp[store[i]] = i;
    }
}

int main(){
    prepare();
    memset(dp,-1,sizeof(dp));
    scanf("%d",&T);
    while(T--){
        scanf("%I64d%I64d",&l,&r);
        
        process(l - 1);
        LL ans0 = dfs(tail - 1,0,0,true);
        
        process(r);
        LL ans1 = dfs(tail - 1,0,0,true);
        
        printf("%I64d\n",ans1 - ans0);
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/tiberius/p/11437488.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值