#include #define rep(i, n) for(int(i) = 0; (i) < (n); (i)++) #define FOR(i, m, n) for(int(i) = (m); (i) < (n); (i)++) #define ALL(v) (v).begin(), (v).end() #define LLA(v) (v).rbegin(), (v).rend() #define SZ(v) (v).size() #define INT(...) \ int __VA_ARGS__; \ read(__VA_ARGS__) #define LL(...) \ ll __VA_ARGS__; \ read(__VA_ARGS__) #define DOUBLE(...) \ double __VA_ARGS__; \ read(__VA_ARGS__) #define CHAR(...) \ char __VA_ARGS__; \ read(__VA_ARGS__) #define STRING(...) \ string __VA_ARGS__; \ read(__VA_ARGS__) #define VEC(type, name, size) \ vector name(size); \ read(name) using namespace std; using ll = long long; using pii = pair; using pll = pair; using Graph = vector>; template struct edge { int to; T cost; edge(int t, T c) : to(t), cost(c) {} }; template using WGraph = vector>>; const int INF = 1 << 30; const ll LINF = 1LL << 60; const int MOD = 1e9 + 7; template inline vector make_vec(size_t a, T val) { return vector(a, val); } template inline auto make_vec(size_t a, Ts... ts) { return vector(a, make_vec(ts...)); } void read() {} template inline void read(T &a) { cin >> a; } template inline void read(pair &p) { read(p.first), read(p.second); } template inline void read(vector &v) { for(auto &a : v) read(a); } template inline void read(Head &head, Tail &...tail) { read(head), read(tail...); } template inline void chmax(T &a, T b) { (a < b ? a = b : a); } template inline void chmin(T &a, T b) { (a > b ? a = b : a); } template class ModInt { using u64 = uint_fast64_t; public: u64 val; ModInt(const u64 x = 0) : val((x + MOD) % MOD) {} constexpr u64 &value() { return val; } constexpr ModInt operator-() { return val ? MOD - val : 0; } constexpr ModInt operator+(const ModInt &rhs) const { return ModInt(*this) += rhs; } constexpr ModInt operator-(const ModInt &rhs) const { return ModInt(*this) -= rhs; } constexpr ModInt operator*(const ModInt &rhs) const { return ModInt(*this) *= rhs; } constexpr ModInt operator/(const ModInt &rhs) const { return ModInt(*this) /= rhs; } constexpr ModInt &operator+=(const ModInt &rhs) { val += rhs.val; if(val >= MOD) { val -= MOD; } return *this; } constexpr ModInt &operator-=(const ModInt &rhs) { if(val < rhs.val) { val += MOD; } val -= rhs.val; return *this; } constexpr ModInt &operator*=(const ModInt &rhs) { val = val * rhs.val % MOD; return *this; } constexpr ModInt &operator/=(const ModInt &rhs) { *this *= rhs.inv(); return *this; } constexpr bool operator==(const ModInt &rhs) { return this->val == rhs.val; } constexpr bool operator!=(const ModInt &rhs) { return this->val != rhs.val; } friend constexpr ostream &operator<<(ostream &os, const ModInt &x) { return os << x.val; } friend constexpr istream &operator>>(istream &is, ModInt &x) { return is >> x.val; } constexpr ModInt inv() const { return ModInt(*this).pow(MOD - 2); } constexpr ModInt pow(ll e) const { u64 x = 1, p = val; while(e > 0) { if(e % 2 == 0) { p = (p * p) % MOD; e /= 2; } else { x = (x * p) % MOD; e--; } } return ModInt(x); } }; using mint = ModInt; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); STRING(n); int N = n.size(); auto dp = make_vec(N + 1, 2, 2, 3, 3, mint(0)); dp[0][0][0][0][0] = 1; for(int i = 0; i < N; i++) { int now = n[i] - '0'; for(int smaller = 0; smaller < 2; smaller++) { for(int start = 0; start < 2; start++) { for(int j = 0; j < 3; j++) { for(int k = 0; k < 3; k++) { for(int x = 0; x <= (smaller ? 9 : now); x++) { if(start && x == 0) continue; if(x == 2) dp[i + 1][smaller || x < now][start || x] [min(2, j + 1)][k] += dp[i][smaller][start][j][k]; else if(x == 4) dp[i + 1][smaller || x < now][start || x] [min(2, j + 2)][k] += dp[i][smaller][start][j][k]; else if(x == 5) dp[i + 1][smaller || x < now][start || x][j] [min(2, k + 1)] += dp[i][smaller][start][j][k]; else if(x == 6) dp[i + 1][smaller || x < now][start || x] [min(2, j + 1)][k] += dp[i][smaller][start][j][k]; else if(x == 8) dp[i + 1][smaller || x < now][start || x] [min(2, j + 3)][k] += dp[i][smaller][start][j][k]; else dp[i + 1][smaller || x < now][start || x][j] [k] += dp[i][smaller][start][j][k]; } } } } } } cout << dp[N][0][1][2][2] + dp[N][1][1][2][2] << endl; }