标签：Beautiful const int hs LL Codeforces hcnt numbers define

55D - Beautiful numbers

把lcm离散化一下就能过了。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ull unsigned long long
using namespace std;

const int N = 1e6 + 7;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 + 7;
const double eps = 1e-8;

int T, s[20], tot, hs[N], hcnt, who[N];
int Lcm[2521][2521];
LL L, R, f[20][2521][50];

inline int getId(int x) {
    return lower_bound(hs + 1, hs + 1 + hcnt, x) - hs;
}

LL dp(int len, int mo, int lcm, bool limit) {
    if(len == -1) return mo % hs[lcm] ? 0 : 1;
    if(!limit && ~f[len][mo][lcm]) return f[len][mo][lcm];
    LL ans = 0;
    int up = limit ? s[len] : 9;
    for(int i = 0; i <= up; i++) {
        ans += dp(len - 1, (mo*10+i)%2520, i ? who[Lcm[i][hs[lcm]]] : lcm, limit && up == i);
    }
    if(!limit) f[len][mo][lcm] = ans;
    return ans;
}

LL solve(LL x) {
    tot = 0;
    for(LL i = x; i; i /= 10) s[tot++] = i % 10;
    return dp(tot - 1, 0, 1, 1);
}

int main() {
    for(int i = 0; i < (1 << 9); i++) {
        int tmp = 1;
        for(int j = 0; j < 9; j++)
            if(i >> j & 1) tmp = tmp / __gcd(tmp, j + 1) * (j + 1);
        hs[++hcnt] = tmp;
    }
    sort(hs + 1, hs + 1 + hcnt);
    hcnt = unique(hs + 1, hs + 1 + hcnt) - hs - 1;
    for(int i = 1; i <= 2520; i++)
        who[i] = getId(i);
    memset(f, -1, sizeof(f));
    for(int i = 1; i <= 2520; i++)
        for(int j = 1; j <= 2520; j++)
            Lcm[i][j] = i / __gcd(i, j) * j;
    scanf("%d", &T);
    while(T--) {
        scanf("%lld%lld", &L, &R);
        printf("%lld\n", (solve(R) - solve(L - 1)));
    }
    return 0;
}

/*
*/

标签：Beautiful,const,int,hs,LL,Codeforces,hcnt,numbers,define
来源： https://www.cnblogs.com/CJLHY/p/10356173.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。