結果

問題 No.315 世界のなんとか3.5
コンテスト
ユーザー yuppe19 😺
提出日時 2019-01-25 16:12:08
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,316 ms / 2,000 ms
コード長 2,083 bytes
コンパイル時間 1,359 ms
コンパイル使用メモリ 95,392 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-16 04:48:48
合計ジャッジ時間 15,465 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge6
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <cmath>
#include <iostream>
#include <vector>
using namespace std;

constexpr int mod = static_cast<int>(powl(10, 9)) + 7;

int f(const string &s, const int P) {
  const int length = static_cast<int>(log10(P)) + 1;
  const int n = static_cast<int>(s.size());
  // dp[未満フラグ][has3][mod 3][mod P] := パターン数
  vector<vector<vector<vector<int>>>> dp(2, vector<vector<vector<int>>>(2, vector<vector<int>>(3, vector<int>(P, 0))));
  vector<vector<vector<vector<int>>>> ndp;
  dp[0][0][0][0] = 1;
  for(int i=0; i<n; ++i) {
    const bool upper = i < n-(2+length);
    const int LIM = upper ? 1 : P;
    ndp.assign(2, vector<vector<vector<int>>>(2, vector<vector<int>>(3, vector<int>(P, 0))));
    for(int less=0; less<2; ++less) {
      for(int has3=0; has3<2; ++has3) {
        for(int rem3=0; rem3<3; ++rem3) {
          for(int remP=0; remP<LIM; ++remP) {
            if(!dp[less][has3][rem3][remP]) { continue; }
            for(int d=0; d<=(less ? 9 : s[i]-'0'); ++d) {
              int nless = less || d < s[i] - '0',
                  nhas3 = has3 || d == 3,
                  nrem3 = (rem3 + d) % 3,
                  nremP = remP;
              if(!upper) { nremP = (remP * 10 + d) % P; }
              ndp[nless][nhas3][nrem3][nremP] += dp[less][has3][rem3][remP];
              ndp[nless][nhas3][nrem3][nremP] %= mod;
            }
          }
        }
      }
    }
    swap(dp, ndp);
  }
  int res = 0;
  for(int less=0; less<2; ++less) {
    for(int remP=1; remP<P; ++remP) {
      res += dp[less][0][0][remP];
      res %= mod;
      for(int rem3=0; rem3<3; ++rem3) {
        res += dp[less][1][rem3][remP];
        res %= mod;
      }
    }
  }
  return res;
}

void minus1(string *s) {
  int n = static_cast<int>(s->size());
  for(int i=n-1; i>=0; --i) {
    if(s->at(i) == '0') {
      s->at(i) = '9';
    } else {
      --s->at(i);
      break;
    }
  }
}

int main(void) {
  string a, b; cin >> a >> b;
  minus1(&a);
  int P; scanf("%d", &P);
  int res = f(b, P) - f(a, P);
  res += mod;
  res %= mod;
  printf("%d\n", res);
  return 0;
}
0