結果

問題 No.315 世界のなんとか3.5
ユーザー yuppe19 😺
提出日時 2019-01-25 18:58:02
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 486 ms / 2,000 ms
コード長 2,136 bytes
コンパイル時間 877 ms
コンパイル使用メモリ 83,304 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-16 04:49:44
合計ジャッジ時間 6,440 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

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

vector<vector<int>> dp, ndp;

inline int calc_state(int less, int has3, int rem3) {
  return (rem3 << 2) | (has3 << 1) | less;
}

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] := パターン数
  dp.assign(12, vector<int>(P, 0));
  dp[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(12, vector<int>(P, 0));
    for(int less=0; less<2; ++less) {
      const int D = less ? 9 : s[i] - '0';
      for(int has3=0; has3<2; ++has3) {
        for(int rem3=0; rem3<3; ++rem3) {
          int state = calc_state(less, has3, rem3);
          for(int remP=0; remP<LIM; ++remP) {
            if(!dp[state][remP]) { continue; }
            for(int d=0; d<=D; ++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; }
              int nstate = calc_state(nless, nhas3, nrem3);
              ndp[nstate][nremP] += dp[state][remP];
              ndp[nstate][nremP] %= mod;
            }
          }
        }
      }
    }
    swap(dp, ndp);
  }
  int res = 0;
  for(int less=0; less<2; ++less) {
    for(int remP=1; remP<P; ++remP) {
      res += dp[calc_state(less, 0, 0)][remP];
      res %= mod;
      for(int rem3=0; rem3<3; ++rem3) {
        res += dp[calc_state(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