#include #include using namespace atcoder; using namespace std; using mint = modint1000000007; using C = complex; const int mod = 998244353; const long long LINF = 1001002003004005006; const int INF = 1001001001; const double PI = acos(-1); const double EPS = 1e-10; const int di[4] = {-1,0,1,0}; const int dj[4] = {0,-1,0,1}; const int dx[8] = {1,1,1,0,0,-1,-1,-1}; const int dy[8] = {1,0,-1,1,-1,1,0,-1}; # define sz(x) (int)(x).size() # define rsz(x,n) x.resize(n) # define yosupo(x) {cout << (x) << endl; return 0;} # define ll long long # define fi first # define se second # define pb push_back # define pf push_front # define eb emplace_back # define ef emplace_front # define pob pop_back # define pof pop_front # define GET_MACRO(_1, _2, _3, NAME, ...) NAME # define _rep(i, n) _rep2(i, 0, n) # define _rep2(i, a, b) for(int i = (int)(a); i < (int)(b); i++) # define rep(...) GET_MACRO(__VA_ARGS__, _rep2, _rep)(__VA_ARGS__) # define srep(i, a, b) for(int i = a; i <= b; ++i) # define all(obj) (obj).begin(), (obj).end() # define rall(obj) (obj).rbegin(), (obj).rend() inline void YesNo(bool f) { std::cout << (f? "Yes": "No") << std::endl; } void read() {} template void read(T &t, U &...u) { cin >> t; read(u...); } void writeln() { cout << endl; } template void writeln(const T &t, const U &...u) { cout << t; if (sizeof...(u)) cout << sep; writeln(u...); } # define Pll pair # define P pair # define bit(x,i) (((x) >> (i)) & 1) # define equals(a, b) (fabs((a) - (b)) < EPS) // 誤差を考慮した同値判定 templatebool chmax(T &a, const T &b) { if (a < b) { a = b; return 1; } return 0; } templatebool chmin(T &a, const T &b) { if (b < a) { a = b; return 1; } return 0; } templateistream& operator>>(istream&i,vector&v){rep(j,v.size())i>>v[j];return i;} templatestring join(vector&v){stringstream s;rep(i,v.size())s<<' '<ostream& operator<<(ostream&o,vector&v){if(v.size())o<ostream& operator<<(ostream&o,vector>&vv){if(vv.size())o< ostream& operator<<(ostream& os, const pair& p) { return os << "P(" << p.first << " " << p.second << ")";} template using vc = vector; template using vv = vc>; using vi = vc; using vvi = vv; using vl = vc; using vvl = vv; using vm = vc; using vvm = vv; mint f(string a, int modP) { int n = sz(a); bool is3 = false; bool ok3 = false; ll tot = 0; rep(i,n) { if(a[i] == '3') is3 = true; tot += (a[i]-'0'); } if(tot %3 == 0) ok3 = true; bool isaho = false; string tt; reverse(all(a)); rep(i,min(sz(a),5)) tt += a[i]; reverse(all(a)); reverse(all(tt)); ll kk = stoll(tt); bool is8 = false; if(kk%800 != (modP%800)) is8 = true; return ((ok3 || is3) && is8); } mint solve(string s, int modP) { //8の倍数は下3桁が8の倍数かどうか? //80の倍数は下4桁が80の倍数かどうか? //800の倍数は下5桁が800の倍数かどうか? //よって高々下5桁を管理すればいい? //桁DPは上から桁を決めていくから, 最後の3-5桁まではdp[i][flag][is3][mod3][0]みたいにしてまとめておく //最後の3-5桁だけdp[i][flag][is3][mod3][下3-5桁のmod(8,80,800)]に遷移させる。 //inline化する mint dp[2][2][3][800] = {}; dp[0][0][0][0] = 1; while(s.size() < 10) s = '0'+s; int n = sz(s); rep(i,n-5) { mint p[2][2][3][800] = {}; swap(dp,p); int nd = s[i]-'0'; rep(j,2) { if(j == 0) { rep(mod3,3) { if(nd == 3) { dp[0][1][(mod3+nd)%3][0] += p[0][0][mod3][0]; dp[0][1][(mod3+nd)%3][0] += p[0][1][mod3][0]; } else { dp[0][0][(mod3+nd)%3][0] += p[0][0][mod3][0]; dp[0][1][(mod3+nd)%3][0] += p[0][1][mod3][0]; } } } else { rep(keta,nd) { if(keta == 3) { rep(mod3,3) { dp[1][1][(mod3+keta)%3][0] += p[0][0][mod3][0]; dp[1][1][(mod3+keta)%3][0] += p[0][1][mod3][0]; } } else { rep(mod3,3) { dp[1][0][(mod3+keta)%3][0] += p[0][0][mod3][0]; dp[1][1][(mod3+keta)%3][0] += p[0][1][mod3][0]; } } } rep(keta,10) { if(keta == 3) { rep(mod3,3) { dp[1][1][(mod3+keta)%3][0] += p[1][0][mod3][0]; dp[1][1][(mod3+keta)%3][0] += p[1][1][mod3][0]; } } else { rep(mod3,3) { dp[1][1][(mod3+keta)%3][0] += p[1][1][mod3][0]; dp[1][0][(mod3+keta)%3][0] += p[1][0][mod3][0]; } } } } } } //2遷移目を考える rep(i,n-5,n) { mint p[2][2][3][800] = {}; swap(dp,p); int nd = s[i]-'0'; rep(j,2) { if(j == 0) { if(nd == 3) { rep(mod3,3) { rep(mod800,800) { dp[0][1][(mod3+nd)%3][(mod800*10+nd)%800] += p[0][0][mod3][mod800]; dp[0][1][(mod3+nd)%3][(mod800*10+nd)%800] += p[0][1][mod3][mod800]; } } } else { rep(mod3,3) { rep(mod800,800) { dp[0][0][(mod3+nd)%3][(mod800*10+nd)%800] += p[0][0][mod3][mod800]; dp[0][1][(mod3+nd)%3][(mod800*10+nd)%800] += p[0][1][mod3][mod800]; } } } } else { rep(keta,nd) { rep(mod3,3) { rep(mod800,800) { if(keta == 3) { dp[1][1][(mod3+keta)%3][(mod800*10+keta)%800] += p[0][0][mod3][mod800]; dp[1][1][(mod3+keta)%3][(mod800*10+keta)%800] += p[0][1][mod3][mod800]; } else { dp[1][0][(mod3+keta)%3][(mod800*10+keta)%800] += p[0][0][mod3][mod800]; dp[1][1][(mod3+keta)%3][(mod800*10+keta)%800] += p[0][1][mod3][mod800]; } } } } rep(keta,10) { rep(mod3,3) { rep(mod800,800) { if(keta == 3) { dp[1][1][(mod3+keta)%3][(mod800*10+keta)%800] += p[1][0][mod3][mod800]; dp[1][1][(mod3+keta)%3][(mod800*10+keta)%800] += p[1][1][mod3][mod800]; } else { dp[1][1][(mod3+keta)%3][(mod800*10+keta)%800] += p[1][1][mod3][mod800]; dp[1][0][(mod3+keta)%3][(mod800*10+keta)%800] += p[1][0][mod3][mod800]; } } } } } } } mint ans = 0; if(modP == 8) { rep(isUnder, 2) rep(is3, 2) rep(mod3, 3) rep(mod800, 800) { if ((is3 || !mod3) && ((mod800 != 0) && mod800 % 8 != 0)) ans += dp[isUnder][is3][mod3][mod800]; } return ans; } else if(modP == 80) { rep(isUnder, 2) rep(is3, 2) rep(mod3, 3) rep(mod800, 800) { if ((is3 || !mod3) && ((mod800 != 0) && mod800 % 80 != 0)) ans += dp[isUnder][is3][mod3][mod800]; } return ans; } rep(isUnder, 2) rep(is3, 2) rep(mod3, 3) rep(mod800, 800) { if ((is3 || !mod3) && ((mod800 != 0))) ans += dp[isUnder][is3][mod3][mod800]; } return ans; } int main() { string a,b; ll p; read(a,b,p); mint ans = solve(b,p) - solve(a,p) + f(a,p); cout << ans.val() << endl; }