結果
問題 | No.1417 100の倍数かつ正整数(2) |
ユーザー |
|
提出日時 | 2021-03-05 23:06:58 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 15 ms / 3,000 ms |
コード長 | 2,109 bytes |
コンパイル時間 | 2,130 ms |
コンパイル使用メモリ | 177,552 KB |
実行使用メモリ | 7,744 KB |
最終ジャッジ日時 | 2024-10-07 04:37:29 |
合計ジャッジ時間 | 2,724 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 |
ソースコード
#include <bits/stdc++.h> #define rep(i, n) for (int i = 0; i < (int)(n); i++) using namespace std; typedef long long ll; int main() { string N; cin >> N; ll MOD = 1000000007; int l = N.size(); vector<vector<vector<vector<ll>>>> dp(l + 1, vector<vector<vector<ll>>>(2, vector<vector<ll>>(3, vector<ll>(3, 0)))); dp[0][0][0][0] = 1; rep(i, l) dp[i + 1][1][0][0] = 1; // dp[i + 1][smaller][j][k]: i桁までで2^j5^kの数 rep(i, l) { //cout << "i=" << i << " "; int n = N[i] - '0'; rep(s, 2) { //cout << "s=" << s << " "; rep(j, 3) { //cout << "j=" << j << " "; rep(k, 3) { //cout << "k=" << k << " "; rep(x, (s ? 10 : n + 1)) { //cout << "x=" << x << " "; if (x == 0) continue; if (x % 4 == 0) { dp[i + 1][s || x < n][min(2, j + 2)][k] += dp[i][s][j][k]; dp[i + 1][s || x < n][min(2, j + 2)][k] %= MOD; } else if (x % 2 == 0) { dp[i + 1][s || x < n][min(2, j + 1)][k] += dp[i][s][j][k]; dp[i + 1][s || x < n][min(2, j + 1)][k] %= MOD; } else if (x == 5) { dp[i + 1][s || x < n][j][min(2, k + 1)] += dp[i][s][j][k]; dp[i + 1][s || x < n][j][min(2, k + 1)] %= MOD; } else { dp[i + 1][s || x < n][j][k] += dp[i][s][j][k]; dp[i + 1][s || x < n][j][k] %= MOD; } } } } } } cout << (dp[l][0][2][2] + dp[l][1][2][2]) % MOD << "\n"; /* rep(i, l) { rep(s, 2) { rep(j, 3) { rep(k, 3) { cout << "dp[" << i + 1 << "][" << s << "][" << j << "][" << k << "]=" << dp[i + 1][s][j][k] << "\n"; } } } } */ }