結果
| 問題 | 
                            No.260 世界のなんとか3
                             | 
                    
| コンテスト | |
| ユーザー | 
                             risujiroh
                         | 
                    
| 提出日時 | 2018-09-08 16:28:25 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 122 ms / 2,000 ms | 
| コード長 | 1,496 bytes | 
| コンパイル時間 | 2,312 ms | 
| コンパイル使用メモリ | 183,264 KB | 
| 実行使用メモリ | 18,892 KB | 
| 最終ジャッジ日時 | 2024-12-16 09:03:34 | 
| 合計ジャッジ時間 | 4,789 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 27 | 
ソースコード
#include <bits/stdc++.h>
using namespace std;
using lint = long long int;
template<class T = int> using V = vector<T>;
template<class T = int> using VV = V< V<T> >;
template<class T> void assign(V<T>& v, int n, const T& a = T()) { v.assign(n, a); }
template<class T, class... U> void assign(V<T>& v, int n, const U&... u) { v.resize(n); for (auto&& i : v) assign(i, u...); }
const lint mod = 1e9 + 7;
inline lint emod(lint a, lint p = mod) { return (a % p + p) % p; }
int main() {
  cin.tie(NULL); ios::sync_with_stdio(false);
  string a, b; cin >> a >> b;
  auto calc = [](string x, int in) {
    int n = x.size();
    VV<VV<V<lint>>> dp; assign(dp, n + 1, 2, 3, 2, 8);
    dp[0][1][0][0][0] = 1;
    for (int i = 0; i < n; i++) for (int r = 0; r < 3; r++) for (int s = 0; s < 8; s++) for (int d = 0; d < 10; d++) {
      (dp[i + 1][0][(r + d) % 3][d == 3][2 * s + d & 7] += dp[i][0][r][0][s]) %= mod;
      (dp[i + 1][0][(r + d) % 3][1][2 * s + d & 7] += dp[i][0][r][1][s]) %= mod;
      if (d > x[i] - '0') continue;
      (dp[i + 1][d == x[i] - '0'][(r + d) % 3][d == 3][2 * s + d & 7] += dp[i][1][r][0][s]) %= mod;
      (dp[i + 1][d == x[i] - '0'][(r + d) % 3][1][2 * s + d & 7] += dp[i][1][r][1][s]) %= mod;
    }
    lint res = 0;
    for (int f = 0; f < 2; f++) for (int r = 0; r < 3; r++) for (int g = 0; g < 2; g++) for (int s = 0; s < 8; s++) if (f <= in and (!r or g) and s) (res += dp[n][f][r][g][s]) %= mod;
    return res;
  };
  cout << emod(calc(b, 1) - calc(a, 0)) << '\n';
}
            
            
            
        
            
risujiroh