#include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using lli = long long int; using Vint = std::vector; using Vlli = std::vector; using Wint = std::vector; using Wlli = std::vector; using Vbool = std::vector; using Wbool = std::vector; using pii = std::pair; using pll = std::pair; template using Vec = std::vector; constexpr int MOD = 1e9 + 7; constexpr int INFi = 2e9 + 1; constexpr lli INFl = (lli)(9e18) + 1; const std::vector> DXDY = {{1, 0}, {-1, 0}, {0, 1}, {0, -1}}; constexpr char BR = '\n'; #define DEBUG(x) std::cerr << #x << " = " << x << '\n'; #define FOR(i, a, b) for(int (i) = (a); (i) < (b); ++(i)) #define FOReq(i, a, b) for(int (i) = (a); (i) <= (b); ++(i)) #define rFOR(i, a, b) for(int (i) = (b); (i) >= (a); --(i)) #define FORstep(i, a, b, step) for(int (i) = (a); i < (b); i += (step)) #define REP(i, n) FOR(i, 0, n) #define rREP(i, n) rFOR(i, 0, (n-1)) #define vREP(ele, vec) for(auto &(ele) : (vec)) #define vREPcopy(ele, vec) for(auto (ele) : (vec)) #define SORT(A) std::sort((A).begin(), (A).end()) #define RSORT(A) std::sort((A).rbegin(), (A).rend()) #define ALL(A) (A).begin(), (A).end() // 座標圧縮 (for vector) : ソートしてから使うのが一般的 ; SORT(A) => COORDINATE_COMPRESSION(A) #define COORDINATE_COMPRESSION(A) (A).erase(unique((A).begin(),(A).end()),(A).end()) template inline int argmin(std::vector vec){return min_element(vec.begin(), vec.end()) - vec.begin();} template inline int argmax(std::vector vec){return max_element(vec.begin(), vec.end()) - vec.begin();} template inline void chmax(S &a, T b){if(a < b) a = b;} template inline void chmin(S &a, T b){if(a > b) a = b;} template inline void reverseSORT(Vec &Array){ std::sort(Array.begin(), Array.end(), std::greater()); } inline int BitI(int k){return 1 << k;} inline lli BitL(int k){return 1LL << k;} inline void putsDouble(double d){printf("%.16lf\n", d);} template inline std::string toString(T n){ if(n == 0) return "0"; std::string res; if(n < 0){n = -n;while(n != 0){res += (char)(n % 10 + '0'); n /= 10;} std::reverse(res.begin(), res.end()); return '-' + res;} while(n != 0){res += (char)(n % 10 + '0'); n /= 10;} std::reverse(res.begin(), res.end()); return res; } namespace MyFunc{ using LLi = long long int; // GCD(a, b) ; a, bの最大公約数を求める関数 inline LLi gcd(LLi a, LLi b){ while(b != 0){ a %= b; std::swap(a, b);} return a; } // LCM(a, b) ; a, bの最小公倍数を求める関数 inline LLi lcm(LLi a, LLi b){ return (a * b) / MyFunc::gcd(a, b);} // 累乗を求める関数 inline LLi power(LLi a, LLi n){ LLi res = 1LL, waiting = a; while(n != 0LL){ if((n & 1LL) != 0LL) res *= waiting; waiting *= waiting; n >>= 1;} return res; } // 累乗の余りを求める関数 inline LLi power_mod(LLi a, LLi n, LLi mod_number___ = 1e9 + 7){ LLi res = 1LL, waiting = a; while(n != 0LL){ if((n & 1LL) != 0LL){ res *= waiting; res %= mod_number___;} waiting *= waiting; waiting %= mod_number___; n >>= 1; } return res; } // Z/pZ上の逆元を求める関数 (フェルマーの小定理) inline LLi inverse_mod(LLi a, LLi mod_number___ = 1e9 + 7){ return MyFunc::power_mod(a, mod_number___-2); } inline LLi inverse_mod_euclid(LLi a, LLi mod_number___ = 1e9+7){ LLi b = mod_number___, u = 1, v = 0; while (b != 0) { LLi t = a / b; a -= t * b; std::swap(a, b); u -= t * v; std::swap(u, v);} u %= mod_number___; if (u < 0) u += mod_number___; return u; } // 素数であるかを判定する関数 template inline bool isPrime(Integer_type n){ if(n < 2) return false; if(n == 2) return true; if(n % 2 == 0) return false; for(Integer_type x = 3; x * x <= n; ++++x) if(n % x == 0) return false; return true; } // 素数であるかの真偽表を返す : n ≥ 1 inline std::vector primeTable(int n){ std::vector res(n+1, true); res[0] = false; res[1] = false; for(int x = 2; x * x <= n; ++x) if(res[x]){ for(int i = 2 * x; i <= n; i += x){ res[i] = false; } } return std::move(res); } // 素因数分解したベクトルを返す ; {素因数, 指数} template inline std::vector> prime_factorization(Integer_type n){ std::vector> res(0); if(n <= 0) return std::move(res); // 例外処理 : nが 0 以下 if(n % 2 == 0){ n /= 2; int cnt = 1; while(n % 2 == 0){ n /= 2; cnt++;} res.emplace_back(make_pair(2, cnt)); } Integer_type x = 3; while(x * x <= n){ if(n % x == 0){ n /= x; int cnt = 1; while(n % x == 0){ n /= x; cnt++; } res.emplace_back(make_pair(x, cnt)); } ++++x; } if(n > 1) res.emplace_back(make_pair(n, 1)); return std::move(res); } } // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ constexpr int MAX_SIZE = 200 + 2; constexpr int thisMOD = 1e9+9; char M[MAX_SIZE], D[MAX_SIZE]; int main(void){ scanf("%s%s", M, D); const int nM = strlen(M), nD = strlen(D); // DP_M[i][j][k] := M以下の整数のうち, 上から i 桁まで見て各位の和が j となるもの (kは未満フラグ) // ~> 前後関係から考えて今と次がわかれば十分 DP_M[j] : 和が j となるもの // DP_D[i][j][k] := D以下 同様 const int mxM = 9 * nM, mxD = 9 * nD; Vint DP_M(mxM+1, 0), DP_D(mxD+1, 0); Vint nxtM = DP_M; int ever = 0; REP(i, nM){ // Mに関する i桁目 FOReq(j, 0, mxM) nxtM[j] = 0; // 初期化 int d = M[i] - '0'; // 未満が確定したものは自由に使える FOReq(j, 0, mxM) REP(k, 10){nxtM[j+k] += DP_M[j]; nxtM[j+k] %= thisMOD;} // 未確定のもの REP(k, d) nxtM[ever+k]++; ever += d; std::swap(DP_M, nxtM); } DP_M[ever]++; Vint nxtD = DP_D; ever = 0; REP(i, nD){ FOReq(j, 0, mxD) nxtD[j] = 0; int d = D[i] - '0'; FOReq(j, 0, mxD) REP(k, 10){nxtD[j+k] += DP_D[j]; nxtD[j+k] %= thisMOD;} REP(k, d) nxtD[ever+k]++; ever += d; std::swap(DP_D, nxtD); } DP_D[ever]++; lli res = 0; int mn = (mxM >= mxD) ? mxD : mxM; FOReq(a, 1, mn){ res += (static_cast(DP_D[a])) * DP_M[a]; res %= thisMOD; } printf("%lld\n", res); return 0; }