結果

問題 No.260 世界のなんとか3
ユーザー pythontanukipythontanuki
提出日時 2022-12-27 15:20:43
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 47 ms / 2,000 ms
コード長 5,458 bytes
コンパイル時間 4,295 ms
コンパイル使用メモリ 267,984 KB
実行使用メモリ 7,768 KB
最終ジャッジ日時 2024-05-01 19:02:45
合計ジャッジ時間 5,461 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
7,600 KB
testcase_01 AC 4 ms
7,400 KB
testcase_02 AC 3 ms
7,648 KB
testcase_03 AC 43 ms
7,744 KB
testcase_04 AC 42 ms
7,612 KB
testcase_05 AC 10 ms
7,668 KB
testcase_06 AC 8 ms
7,492 KB
testcase_07 AC 29 ms
7,588 KB
testcase_08 AC 22 ms
7,616 KB
testcase_09 AC 14 ms
7,572 KB
testcase_10 AC 32 ms
7,520 KB
testcase_11 AC 32 ms
7,564 KB
testcase_12 AC 20 ms
7,616 KB
testcase_13 AC 8 ms
7,568 KB
testcase_14 AC 29 ms
7,472 KB
testcase_15 AC 10 ms
7,592 KB
testcase_16 AC 25 ms
7,712 KB
testcase_17 AC 20 ms
7,572 KB
testcase_18 AC 20 ms
7,652 KB
testcase_19 AC 24 ms
7,680 KB
testcase_20 AC 18 ms
7,736 KB
testcase_21 AC 17 ms
7,756 KB
testcase_22 AC 28 ms
7,768 KB
testcase_23 AC 6 ms
7,700 KB
testcase_24 AC 23 ms
7,588 KB
testcase_25 AC 27 ms
7,688 KB
testcase_26 AC 18 ms
7,720 KB
testcase_27 AC 3 ms
7,544 KB
testcase_28 AC 34 ms
7,720 KB
testcase_29 AC 47 ms
7,464 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using namespace std;
using mint = modint1000000007;
using C = complex<double>;
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 <typename T, class... U> void read(T &t, U &...u) { cin >> t; read(u...); }
void writeln() { cout << endl; }
template <typename T, class... U, char sep = ' '> void writeln(const T &t, const U &...u) { cout << t; if (sizeof...(u)) cout << sep; writeln(u...); }
# define Pll pair<ll, ll>
# define P pair<int,int>
# define bit(x,i) (((x) >> (i)) & 1)
# define equals(a, b) (fabs((a) - (b)) < EPS) // 誤差を考慮した同値判定
template<class T>bool chmax(T &a, const T &b) { if (a < b) { a = b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b < a) { a = b; return 1; } return 0; }
template<typename T>istream& operator>>(istream&i,vector<T>&v){rep(j,v.size())i>>v[j];return i;}
template<typename T>string join(vector<T>&v){stringstream s;rep(i,v.size())s<<' '<<v[i];return s.str().substr(1);}
template<typename T>ostream& operator<<(ostream&o,vector<T>&v){if(v.size())o<<join(v);return o;}
template<typename T>ostream& operator<<(ostream&o,vector<vector<T>>&vv){if(vv.size())o<<join(vv);return o;}
template <class T, class U> ostream& operator<<(ostream& os, const pair<T, U>& p) { return os << "P(" << p.first << " " << p.second << ")";}
template<typename T> using vc = vector<T>;
template<typename T> using vv = vc<vc<T>>;
using vi = vc<int>; using vvi = vv<int>;
using vl = vc<ll>; using vvl = vv<ll>;
using vm = vc<mint>; using vvm = vv<mint>;

mint f(string s) {
    int n = sz(s);
    ll tot = 0;
    int flag = 0;
    ll mod8 = 0;
    rep(i,n) {
        int nd = s[i]-'0';
        if(nd == 3) flag |= 1;
        tot += nd;
        tot %= 3;
        mod8 = (mod8*10+nd)%8;
    }
    return (mint)((tot%3 == 0 || flag) && mod8 != 0);
}

mint solve(string s) {
    //桁数, mod8, mod3, 3があるかないか, 以下フラグ
    mint dp[11000][8][3][2][2] = {};
    dp[0][0][0][0][0] = 1;
    int n = sz(s);
    rep(i,n) {
        int nd = s[i]-'0';
        rep(j,2) {
            if(j == 0) {
                rep(mod8,8) {
                    rep(mod3,3) {
                        if(nd == 3) {
                            dp[i+1][(mod8*10+nd)%8][(mod3+nd)%3][1][0] += dp[i][mod8][mod3][0][0];
                            dp[i+1][(mod8*10+nd)%8][(mod3+nd)%3][1][0] += dp[i][mod8][mod3][1][0];
                        } else {
                            dp[i+1][(mod8*10+nd)%8][(mod3+nd)%3][0][0] += dp[i][mod8][mod3][0][0];
                            dp[i+1][(mod8*10+nd)%8][(mod3+nd)%3][1][0] += dp[i][mod8][mod3][1][0];
                        }
                    }
                }
            } else {
                rep(k,nd) {
                    rep(mod8,8) {
                        rep(mod3,3) {
                            if(k == 3) {
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][1][1] += dp[i][mod8][mod3][0][0];
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][1][1] += dp[i][mod8][mod3][1][0];
                            } else {
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][0][1] += dp[i][mod8][mod3][0][0];
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][1][1] += dp[i][mod8][mod3][1][0];
                            }
                        }
                    }
                }
                rep(k,10) {
                    rep(mod8,8) {
                        rep(mod3,3) {
                            if(k == 3) {
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][1][1] += dp[i][mod8][mod3][0][1];
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][1][1] += dp[i][mod8][mod3][1][1];
                            } else {
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][0][1] += dp[i][mod8][mod3][0][1];
                                dp[i+1][(mod8*10+k)%8][(mod3+k)%3][1][1] += dp[i][mod8][mod3][1][1];
                            }
                        }
                    }
                }
            }
        }
    }
    mint res = 0;
    rep(mod8,8) rep(mod3,3) rep(flag,2) rep(t,2) {
        if((mod8 != 0) && (mod3 == 0 || flag == 1)) res += dp[n][mod8][mod3][flag][t];
    }
    return res;
}

int main() {
    string s,t;
    read(s,t);
    mint res = solve(t)-solve(s)+f(s);
    writeln(res.val());
}
0