#include <bits/stdc++.h> using namespace std; typedef unsigned uint; typedef long long ll; typedef unsigned long long ull; typedef long double ldbl; typedef pair<int, int> pii; typedef pair<uint, uint> puu; typedef pair<ll, ll> pll; typedef pair<ull, ull> pull; typedef pair<double, double> pdd; typedef vector<int> vi; typedef vector<uint> vu; typedef vector<ll> vll; typedef vector<ull> vull; typedef vector<pii> vpii; typedef vector<puu> vpuu; typedef vector<pll> vpll; typedef vector<pull> vpull; typedef vector<string> vstr; typedef vector<double> vdbl; typedef vector<ldbl> vldbl; #define pb push_back #define ppb pop_back #define pfr push_front #define ppfr pop_front #define emp emplace #define empb emplace_back #define be begin #define rbe rbegin #define all(x) (x).be(), (x).end() #define rall(x) (x).rbe(), (x).rend() #define fir first #define sec second #define mkp make_pair #define brif(cond) if (cond) break #define ctif(cond) if (cond) continue #define retif(cond) if (cond) return void canhazfast() {ios_base::sync_with_stdio(false);cin.tie(nullptr);cout.tie(nullptr);} template<typename T> T gcd(T a, T b) {return b ? gcd(b, a%b) : a;} template<typename T> T extgcd(T a, T b, T &x, T &y) { T x0 = 1, y0 = 0, x1 = 0, y1 = 1; while (b) { T q = a/b; a %= b; swap(a, b); x0 -= q*x1; swap(x0, x1); y0 -= q*y1; swap(y0, y1); } x = x0; y = y0; return a; } int ctz(uint x) {return __builtin_ctz(x);} int ctzll(ull x) {return __builtin_ctzll(x);} int clz(uint x) {return __builtin_clz(x);} int clzll(ull x) {return __builtin_clzll(x);} int popcnt(uint x) {return __builtin_popcount(x);} int popcntll(ull x) {return __builtin_popcountll(x);} int bsr(uint x) {return 31^clz(x);} int bsrll(ull x) {return 63^clzll(x);} #define N 20 #define MSK (1<<9) #define MOD 2520 int dx[N]; int msk_lcm[MSK]; ll lt[N][MSK][MOD]; ll eq[N][MSK][MOD]; ll f(ll x) { int n = 0; ll res = 0; while (x) { dx[n++] = x%10; x /= 10; } reverse(dx, dx+n); lt[0][0][0] = eq[0][0][0] = 0; for (int d = 1; d < 10; ++d) { int msk = 1<<(d-1); lt[0][msk][d] = d < dx[0]; eq[0][msk][d] = d == dx[0]; } for (int i = 1; i < n; ++i) { /// init for (int msk = 0; msk < MSK; ++msk) { fill_n(lt[i][msk], MOD, 0LL); fill_n(eq[i][msk], MOD, 0LL); } for (int d = 1; d < 10; ++d) lt[i][1<<(d-1)][d] = 1; /// step for (int src_msk = 0; src_msk < MSK; ++src_msk) { for (int src_mod = 0; src_mod < MOD; ++src_mod) { ll src_lt = lt[i-1][src_msk][src_mod]; ll src_eq = eq[i-1][src_msk][src_mod]; for (int d = 1; d < 10; ++d) { int dst_msk = src_msk | (1<<(d-1)); int dst_mod = (10*src_mod+d)%MOD; lt[i][dst_msk][dst_mod] += src_lt; if (d < dx[i]) lt[i][dst_msk][dst_mod] += src_eq; if (d == dx[i]) eq[i][dst_msk][dst_mod] += src_eq; } } } } for (int msk = 0; msk < MSK; ++msk) { for (int mod = 0; mod < MOD; mod += msk_lcm[msk]) res += lt[n-1][msk][mod]; } return res; } int main() { canhazfast(); ll l, r; msk_lcm[0] = 1; for (int i = 1; i < MSK; ++i) { int bit = i&-i; int d = ctz(bit)+1; int j = i-bit; msk_lcm[i] = msk_lcm[j]*d/gcd(msk_lcm[j], d); } cin >> l >> r; cout << f(r+1)-f(l) << '\n'; return 0; }
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 | #include <bits/stdc++.h> using namespace std; typedef unsigned uint; typedef long long ll; typedef unsigned long long ull; typedef long double ldbl; typedef pair<int, int> pii; typedef pair<uint, uint> puu; typedef pair<ll, ll> pll; typedef pair<ull, ull> pull; typedef pair<double, double> pdd; typedef vector<int> vi; typedef vector<uint> vu; typedef vector<ll> vll; typedef vector<ull> vull; typedef vector<pii> vpii; typedef vector<puu> vpuu; typedef vector<pll> vpll; typedef vector<pull> vpull; typedef vector<string> vstr; typedef vector<double> vdbl; typedef vector<ldbl> vldbl; #define pb push_back #define ppb pop_back #define pfr push_front #define ppfr pop_front #define emp emplace #define empb emplace_back #define be begin #define rbe rbegin #define all(x) (x).be(), (x).end() #define rall(x) (x).rbe(), (x).rend() #define fir first #define sec second #define mkp make_pair #define brif(cond) if (cond) break #define ctif(cond) if (cond) continue #define retif(cond) if (cond) return void canhazfast() {ios_base::sync_with_stdio(false);cin.tie(nullptr);cout.tie(nullptr);} template<typename T> T gcd(T a, T b) {return b ? gcd(b, a%b) : a;} template<typename T> T extgcd(T a, T b, T &x, T &y) { T x0 = 1, y0 = 0, x1 = 0, y1 = 1; while (b) { T q = a/b; a %= b; swap(a, b); x0 -= q*x1; swap(x0, x1); y0 -= q*y1; swap(y0, y1); } x = x0; y = y0; return a; } int ctz(uint x) {return __builtin_ctz(x);} int ctzll(ull x) {return __builtin_ctzll(x);} int clz(uint x) {return __builtin_clz(x);} int clzll(ull x) {return __builtin_clzll(x);} int popcnt(uint x) {return __builtin_popcount(x);} int popcntll(ull x) {return __builtin_popcountll(x);} int bsr(uint x) {return 31^clz(x);} int bsrll(ull x) {return 63^clzll(x);} #define N 20 #define MSK (1<<9) #define MOD 2520 int dx[N]; int msk_lcm[MSK]; ll lt[N][MSK][MOD]; ll eq[N][MSK][MOD]; ll f(ll x) { int n = 0; ll res = 0; while (x) { dx[n++] = x%10; x /= 10; } reverse(dx, dx+n); lt[0][0][0] = eq[0][0][0] = 0; for (int d = 1; d < 10; ++d) { int msk = 1<<(d-1); lt[0][msk][d] = d < dx[0]; eq[0][msk][d] = d == dx[0]; } for (int i = 1; i < n; ++i) { /// init for (int msk = 0; msk < MSK; ++msk) { fill_n(lt[i][msk], MOD, 0LL); fill_n(eq[i][msk], MOD, 0LL); } for (int d = 1; d < 10; ++d) lt[i][1<<(d-1)][d] = 1; /// step for (int src_msk = 0; src_msk < MSK; ++src_msk) { for (int src_mod = 0; src_mod < MOD; ++src_mod) { ll src_lt = lt[i-1][src_msk][src_mod]; ll src_eq = eq[i-1][src_msk][src_mod]; for (int d = 1; d < 10; ++d) { int dst_msk = src_msk | (1<<(d-1)); int dst_mod = (10*src_mod+d)%MOD; lt[i][dst_msk][dst_mod] += src_lt; if (d < dx[i]) lt[i][dst_msk][dst_mod] += src_eq; if (d == dx[i]) eq[i][dst_msk][dst_mod] += src_eq; } } } } for (int msk = 0; msk < MSK; ++msk) { for (int mod = 0; mod < MOD; mod += msk_lcm[msk]) res += lt[n-1][msk][mod]; } return res; } int main() { canhazfast(); ll l, r; msk_lcm[0] = 1; for (int i = 1; i < MSK; ++i) { int bit = i&-i; int d = ctz(bit)+1; int j = i-bit; msk_lcm[i] = msk_lcm[j]*d/gcd(msk_lcm[j], d); } cin >> l >> r; cout << f(r+1)-f(l) << '\n'; return 0; } |