結果

問題 No.315 世界のなんとか3.5
ユーザー OIer1048576OIer1048576
提出日時 2024-09-22 22:24:48
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 2,869 bytes
コンパイル時間 1,695 ms
コンパイル使用メモリ 176,504 KB
実行使用メモリ 153,520 KB
最終ジャッジ日時 2024-09-22 22:25:05
合計ジャッジ時間 15,996 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

#define int long long
const int M = 1e9 + 7;

using state = tuple<int, int, int, int, int>;

enum { ind, ful, ex3, md3, md8 };

state add(const state &a, const int &b, string &s) {
    return {get<ind>(a) + 1, get<ful>(a) && (b + '0' == s[get<ind>(a)]),
            get<ex3>(a) || (b == 3), (get<md3>(a) * 10 + b) % 3,
            (get<md8>(a) * 10 + b) % 8};
}

int mp[] = {3,  13, 23, 30, 31, 32, 33, 34, 35, 36,
            37, 38, 39, 43, 53, 63, 73, 83, 93};

int fun(int n, int k) {
    int ret = 0;
    for (int i = 1; i <= n; i++)
        if (i % 3 == k % 3)
            ret++;
        else if (binary_search(mp, mp + 19, i))
            ret++;
    return ret;
}

int mem[200000][2][2][3][8];

int dp(string s, int lvl) {
    memset(mem, 0, sizeof(mem));
    int n = s.size() - lvl, ldg = 0;
    for (int i = n; i < n + lvl; i++) ldg = (ldg * 10 + (s[i] - '0'));
    auto get_st = [&](state st) -> int & {
        return mem[get<ind>(st)][get<ful>(st)][get<ex3>(st)][get<md3>(st)]
                  [get<md8>(st)];
    };

    get_st({0, 1, 0, 0, 0}) = 1;
    for (int i = 0; i < n; i++)
        for (int f = 0; f < 2; f++)
            for (int e3 = 0; e3 < 2; e3++)
                for (int m3 = 0; m3 < 3; m3++)
                    for (int m8 = 0; m8 < 8; m8++)
                        for (int nxt = 0; nxt < 10; nxt++) {
                            if (f && nxt > s[i] - '0') break;
                            get_st(add({i, f, e3, m3, m8}, nxt, s)) +=
                                get_st({i, f, e3, m3, m8});
                            get_st(add({i, f, e3, m3, m8}, nxt, s)) %= M;
                        }
    int ret = 0;
    for (int f = 0; f < 2; f++)
        for (int e3 = 0; e3 < 2; e3++)
            for (int m3 = 0; m3 < 3; m3++)
                if ((e3 || (m3 == 0)))
                    for (int m8 = 1; m8 < 8; m8++)
                        ret += get_st(state{n, f, e3, m3, m8}), ret %= M;
    for (int m3 = 0; m3 < 3; m3++)
        for (int m8 = 0; m8 < 8; m8++)
            ret +=
                get_st(state{n, 0, 0, m3, m8}) *
                    fun((int)pow(10, lvl) - 1, 3 - m3) % M +
                get_st(state{n, 1, 0, m3, m8}) * fun(ldg, 3 - m3) % M +
                get_st(state{n, 0, 1, m3, m8}) * ((int)pow(10, lvl) - 1) % M +
                get_st(state{n, 1, 1, m3, m8}) * ldg % M,
                ret %= M;
    return ret;
}

string minus1(string s) {
    for (int i = s.size() - 1; i >= 0; i--) {
        if (s[i] == '0')
            s[i] = '9';
        else {
            s[i]--;
            break;
        }
    }
    return s;
}

signed main() {
    freopen("vr.in", "r", stdin);
    freopen("vr.out", "w", stdout);
    string A, B;
    int p;
    cin >> A >> B >> p;
    cout << (dp(B, floor(log10(p))) - dp(minus1(A), floor(log10(p))) + M) % M
         << endl;
}
0