/////////////////////////////////////////////////////////////////////////////// #include #include #include #include using namespace std; /////////////////////////////////////////////////////////////////////////////// #define DEBUG 0 #define pb push_back #define V vector #define S static #define rep(i,n) for(ll i=0LL; i=0LL; --i) #define ALL(a) (a).begin(),(a).end() #define RALL(a) (a).rbegin(),(a).rend() #define CIN(x) do { \ assert(!cin.eof()); \ cin >> x; \ assert(!cin.fail()); \ } while(0); #if DEBUG #define debug_print(...) _debug_print(__VA_ARGS__) #define debug_printf(...) printf(__VA_ARGS__) #define debug_print_time _debug_print_time #else // DEBUG #define debug_print(...) #define debug_printf(...) #define debug_print_time #endif // DEBUG typedef long long ll; typedef unsigned long long ull; typedef __int128_t ll128; typedef tuple t2; typedef tuple t3; typedef tuple t4; typedef tuple t5; template using priority_queue_incr = priority_queue, greater >; /////////////////////////////////////////////////////////////////////////////// template void incr_m(map &m, TT k) { if (m.find(k) == m.end()) m[k] = 0; m[k]++; } void llin(ll &a) { CIN(a); } void llinl1(V &v, ll count) { for (ll i = 0LL; i < count; ++i) { ll a; CIN(a); v.push_back(a); } } void llinl2(V &v, ll count) { for (ll i = 0LL; i < count; ++i) { ll a, b; CIN(a >> b); v.push_back(t2(a, b)); } } void llinl3(V &v, ll count) { for (ll i = 0LL; i < count; ++i) { ll a, b, c; CIN(a >> b >> c); v.push_back(t3(a, b, c)); } } void llinl4(V &v, ll count) { for (ll i = 0LL; i < count; ++i) { ll a, b, c, d; CIN(a >> b >> c >> d); v.push_back(t4(a, b, c, d)); } } void llina(V &v, ll count) { llinl1(v, count); } template void sort(V &v) { sort(v.begin(), v.end()); } template void sort_reverse(V &v) { sort(v.begin(), v.end(), greater()); } t2 _ext_gcd(ll a, ll b, ll g) { if (!b) return t2(1, 0); ll q = a / b; ll r = a % b; t2 sol = _ext_gcd(b, r, g); ll sx = get<0>(sol); ll sy = get<1>(sol); ll x = sy; ll y = sx - q * sy; return t2(x, y); } t2 ext_gcd(ll a, ll b) { return _ext_gcd(a, b, gcd(a, b)); } // x and mod must be coprime ll mod_inv(ll x, ll mod) { t2 result = ext_gcd(x, mod); ll ret = get<0>(result); while (ret < 0) ret += mod; ret %= mod; return ret; } ll msec() { struct timeval tv; gettimeofday(&tv, NULL); ll ret = 0; ret += tv.tv_sec * 1000LL; ret += tv.tv_usec / 1000LL; return ret; } template void _debug_print(T x) { cout << x << " "; } void _debug_print(const t2 &x) { ll x1 = get<0>(x); ll x2 = get<1>(x); cout << "(" << x1 << ", " << x2 << ") "; } void _debug_print(const t3 &x) { ll x1 = get<0>(x); ll x2 = get<1>(x); ll x3 = get<2>(x); cout << "(" << x1 << ", " << x2 << ", " << x3 << ") "; } void _debug_print(const t4 &x) { ll x1 = get<0>(x); ll x2 = get<1>(x); ll x3 = get<2>(x); ll x4 = get<3>(x); cout << "(" << x1 << ", " << x2 << ", " << x3 << ", " << x4 << ") "; } template void _debug_print(T xarray[], ll n) { rep (i, n) _debug_print(xarray[i]); cout << endl; } template void _debug_print(const V &xlist) { for (auto x : xlist) _debug_print(x); cout << endl; } template void _debug_print(const set &xset) { for (auto x : xset) _debug_print(x); cout << endl; } template void _debug_print(const map &xlist) { for (auto x : xlist) { TT k = x.first; T v = x.second; cout << "K="; _debug_print(k); cout << " V="; _debug_print(v); cout << endl; } } #define MOD (1000000007LL) // #define MOD (998244353LL) // O(log(exp)) ll mod_pow(ll base, ll exp, ll mod) { ll ret = 1LL; for ( ; exp; ) { if (exp & 1LL) { ret *= base; ret %= mod; } base = (base * base) % mod; exp >>= 1; } return ret; } ll mod_mlt(ll x, ll y, ll mod) { ll ret = 0LL; x %= mod; while (y) { if (y & 1LL) { ret += x; ret %= mod; } y >>= 1; x <<= 1; x %= mod; } return ret; } // returns t2(solution, mod) t2 chinese_remainder(ll a1, ll m1, ll a2, ll m2) { assert(a1 >= 0); assert(m1 > 0); assert(a2 >= 0); assert(m2 > 0); ll mgcd = gcd(m1, m2); if (a1 % mgcd != a2 % mgcd) return t2(0, 0); ll mlcm = m1 * m2 / mgcd; t2 z = ext_gcd(m1, m2); ll z1 = get<0>(z); ll z2 = get<1>(z); // ll x = a1 + ((a2 - a1) / mgcd) * m1 * z1; ll x = z1; while (x < 0) x += mlcm; x = mod_mlt(x, m1, mlcm); ll coef = (a2 - a1) / mgcd; while (coef < 0) coef += mlcm; x = mod_mlt(x, coef, mlcm); x += a1; x %= mlcm; return t2(x, mlcm); } void get_divisors(V &retlist, ll x) { for (ll i = 1LL; i < sqrt(x) + 3LL; ++i) { if (x % i == 0LL) { retlist.push_back(i); retlist.push_back(x / i); } } } // returns factors and 1 void get_factors(V &retlist, ll x) { retlist.pb(1LL); for (ll i = 2LL; i < (ll)(sqrt(x)) + 3LL; ++i) { while (x % i == 0LL) { retlist.pb(i); x /= i; } } retlist.pb(x); } bool is_prime(ll x) { V factors, factors2; get_factors(factors, x); for (auto factor : factors) { if (factor > 1) factors2.pb(factor); } return factors2.size() == 1 && x == factors2[0]; } void eratosthenes(set &primes, ll n) { bool *is_not_prime = new bool[n+3LL]; memset(is_not_prime, 0, sizeof(bool) * (n+3LL)); srep (v, 2LL, (ll)sqrt(n)+10LL) { if (is_not_prime[v]) continue; for (ll vv = v * 2LL; vv <= n; vv += v) { is_not_prime[vv] = true; } } srep (v, 2LL, n+1LL) if (!is_not_prime[v]) primes.insert(v); delete [] is_not_prime; } // p must be prime ll mod_root(ll p) { if (p == 2) return 1; if (p == 3) return 2; V flist; get_factors(flist, p-1LL); set fs; for (auto f : flist) if (f > 1) fs.insert(f); srep (a, 2, p) { bool ok = true; for (auto f : fs) { if (mod_pow(a, (p-1LL) / f, p) == 1) { ok = false; break; } } if (ok) return a; } assert(false); } ull combination(ll x, ll y) { if (y > x / 2LL) y = x - y; ull ret = 1LL; for (ll i = 0LL; i < y; ++i) { ret *= x--; ret /= (i + 1LL); } return ret; } void make_factorial(ll factorials[], ll inv_factorials[], ll size, ll mod) { factorials[0] = 1LL; srep (i, 1, size) { factorials[i] = factorials[i-1] * i; factorials[i] %= mod; } inv_factorials[size-1] = mod_inv(factorials[size-1], mod); rrep (i, size-1) { inv_factorials[i] = inv_factorials[i+1LL] * (i+1LL); inv_factorials[i] %= mod; } } #if 0 S ll factorials[200010]; S ll inv_factorials[200010]; make_factorial(factorials, inv_factorials, 200010, MOD); #endif /////////////////////////////////////////////////////////////////////////////// void _main(); int main() { cout << setprecision(12); #if !DEBUG ios::sync_with_stdio(false); cin.tie(0); #endif _main(); return 0; } void _main() { string sn; CIN(sn); ll ans = 0; S ll dp1[10010][3][3]; memset(dp1, 0, sizeof(dp1)); dp1[0][0][0] = 1; srep (len, 1, sn.size()) { dp1[len][0][0] += dp1[len-1][0][0] * 4LL; dp1[len][0][0] %= MOD; dp1[len][0][1] += dp1[len-1][0][0]; dp1[len][0][1] += dp1[len-1][0][1] * 4LL; dp1[len][0][1] %= MOD; dp1[len][0][2] += dp1[len-1][0][1]; dp1[len][0][2] += dp1[len-1][0][2] * 5LL; dp1[len][0][2] %= MOD; dp1[len][1][0] += dp1[len-1][0][0] * 2LL; dp1[len][1][0] += dp1[len-1][1][0] * 4LL; dp1[len][1][0] %= MOD; dp1[len][1][1] += dp1[len-1][1][0]; dp1[len][1][1] += dp1[len-1][0][1] * 2LL; dp1[len][1][1] += dp1[len-1][1][1] * 4LL; dp1[len][1][1] %= MOD; dp1[len][1][2] += dp1[len-1][0][2] * 2LL; dp1[len][1][2] += dp1[len-1][1][1]; dp1[len][1][2] += dp1[len-1][1][2] * 5LL; dp1[len][1][2] %= MOD; dp1[len][2][0] += dp1[len-1][0][0] * 2LL; dp1[len][2][0] += dp1[len-1][1][0] * 4LL; dp1[len][2][0] += dp1[len-1][2][0] * 8LL; dp1[len][2][0] %= MOD; dp1[len][2][1] += dp1[len-1][0][1] * 2LL; dp1[len][2][1] += dp1[len-1][1][1] * 4LL; dp1[len][2][1] += dp1[len-1][2][1] * 8LL; dp1[len][2][1] += dp1[len-1][2][0]; dp1[len][2][1] %= MOD; dp1[len][2][2] += dp1[len-1][0][2] * 2LL; dp1[len][2][2] += dp1[len-1][1][2] * 4LL; dp1[len][2][2] += dp1[len-1][2][1]; dp1[len][2][2] += dp1[len-1][2][2] * 9LL; dp1[len][2][2] %= MOD; } rep (len, 4) { ll cnt = 0; rep (i, 3) rep (j, 3) cnt += dp1[len][i][j]; debug_printf("---- len=%lld cnt=%lld\n", len, cnt); } rep (len, sn.size()) { ans += dp1[len][2][2]; ans %= MOD; } debug_printf("---- ans 1 = %lld\n", ans); S ll dp2[3][3]; // (just, 2cnt, 5cnt) memset(dp2, 0, sizeof(dp2)); ll just2 = 0; ll just5 = 0; bool just_valid = true; for (auto c : sn) { ll dp2_prev[3][3]; memcpy(dp2_prev, dp2, sizeof(dp2_prev)); memset(dp2, 0, sizeof(dp2)); if (just_valid) { srep (d, '1', c) { if (d == '2' || d == '6') { dp2[min(2LL, just2+1LL)][just5]++; dp2[min(2LL, just2+1LL)][just5] %= MOD; } else if (d == '4' || d == '8') { dp2[2][just5]++; dp2[2][just5] %= MOD; } else if (d == '5') { dp2[just2][min(2LL, just5+1LL)]++; dp2[just2][min(2LL, just5+1LL)] %= MOD; } else { dp2[just2][just5]++; dp2[just2][just5] %= MOD; } } } srep (d, '1', '9'+1) rep (cnt2, 3) rep (cnt5, 3) { if (d == '2' || d == '6') { dp2[min(2LL, cnt2+1LL)][cnt5] += dp2_prev[cnt2][cnt5]; dp2[min(2LL, cnt2+1LL)][cnt5] %= MOD; } else if (d == '4' || d == '8') { dp2[2][cnt5] += dp2_prev[cnt2][cnt5]; dp2[2][cnt5] %= MOD; } else if (d == '5') { dp2[cnt2][min(2LL, cnt5+1LL)] += dp2_prev[cnt2][cnt5]; dp2[cnt2][min(2LL, cnt5+1LL)] %= MOD; } else { dp2[cnt2][cnt5] += dp2_prev[cnt2][cnt5]; dp2[cnt2][cnt5] %= MOD; } } if (c == '2' || c == '6') just2 = min(2LL, just2+1LL); if (c == '4' || c == '8') just2 = 2; if (c == '5') just5 = min(2LL, just5+1LL); if (c == '0') just_valid = false; debug_printf("------------------\n"); rep (x2, 3) rep (x5, 3) debug_printf("dp2[%lld][%lld] = %lld\n", x2, x5, dp2[x2][x5]); } ans += dp2[2][2]; if (just_valid && just2 == 2 && just5 == 2) ans++; ans %= MOD; debug_printf("---- ans\n"); cout << ans << endl; } ///////////////////////////////////////////////////////////////////////////////