結果

問題 No.260 世界のなんとか3
コンテスト
ユーザー pekempey
提出日時 2015-07-31 23:52:48
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 84 ms / 2,000 ms
コード長 2,671 bytes
コンパイル時間 1,407 ms
コンパイル使用メモリ 162,424 KB
実行使用メモリ 11,136 KB
最終ジャッジ日時 2024-11-07 18:01:19
合計ジャッジ時間 3,338 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i, a) for (int i = 0; i < (a); i++)
#define rep2(i, a, b) for (int i = (a); i < (b); i++)
#define repr(i, a) for (int i = (a) - 1; i >= 0; i--)
#define repr2(i, a, b) for (int i = (b) - 1; i >= (a); i--)
using namespace std;
typedef long long ll;
const ll inf = 1e9;
const ll mod = 1e9 + 7;

ll dp[10002][2][2][3][8]; // ind, below, use3, mod3, mod8

ll pow8[10002];

ll solve(string N) {
    memset(dp, 0, sizeof(dp));

    dp[0][0][0][0][0] = 1;
    rep (i, N.length()) {
        rep (j, 2) rep (k, 2) rep (l, 3) rep (m, 8) {
            if (j == 0) {
                int ub = N[i] - '0';

                rep (n, ub + 1) {
                    if (n == ub) {
                        if (n == 3) {
                            (dp[i + 1][0][1][(l + n) % 3][(m + n * pow8[N.length() - i - 1]) % 8] += dp[i][j][k][l][m]) %= mod;
                        } else {
                            (dp[i + 1][0][k][(l + n) % 3][(m + n * pow8[N.length() - i - 1]) % 8] += dp[i][j][k][l][m]) %= mod;
                        }
                    } else {
                        if (n == 3) {
                            (dp[i + 1][1][1][(l + n) % 3][(m + n * pow8[N.length() - i - 1]) % 8] += dp[i][j][k][l][m]) %= mod;
                        } else {
                            (dp[i + 1][1][k][(l + n) % 3][(m + n * pow8[N.length() - i - 1]) % 8] += dp[i][j][k][l][m]) %= mod;
                        }
                    }
                }
            } else {
                rep (n, 10) {
                    if (n == 3) {
                        (dp[i + 1][1][1][(l + n) % 3][(m + n * pow8[N.length() - i - 1]) % 8] += dp[i][j][k][l][m]) %= mod;
                    } else {
                        (dp[i + 1][1][k][(l + n) % 3][(m + n * pow8[N.length() - i - 1]) % 8] += dp[i][j][k][l][m]) %= mod;
                    }
                }
            }
        }
    }

    ll ans = 0;
    rep (j, 2) rep (k, 2) rep (l, 3) rep2 (m, 1, 8) {
        if (k || l == 0) {
            (ans += dp[N.length()][j][k][l][m]) %= mod;
        }
    }
    return ans;
}

ll modulo(ll a) {
    a %= mod; a += mod; a %= mod;
    return a;
}

string decrement(string S) {
    int bollow = 1;
    repr (i, S.length()) {
        int p = S[i] - '0';
        if (p - bollow < 0) {
            S[i] = '9';
            bollow = 1;
        } else {
            S[i] = (p - bollow) + '0';
            break;
        }
    }

    return S;
}

int main() {
    pow8[0] = 1;
    rep2 (i, 1, 10002) (pow8[i] = pow8[i - 1] * 10) %= 8;

    string A, B;
    cin >> A >> B;

    ll b = solve(B);
    ll a = solve(decrement(A));

    cout << modulo(b - a) << endl;
}
0