結果
問題 | No.1417 100の倍数かつ正整数(2) |
ユーザー |
👑 |
提出日時 | 2023-11-21 15:44:26 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 4 ms / 3,000 ms |
コード長 | 5,518 bytes |
コンパイル時間 | 1,005 ms |
コンパイル使用メモリ | 111,068 KB |
最終ジャッジ日時 | 2025-02-17 22:44:41 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 |
ソースコード
#line 2 "/home/sakflat/CP/_library/cpp/template/template.cpp" //yukicoder@cpp17 #include <iostream> #include <cmath> #include <algorithm> #include <iomanip> #include <string> #include <vector> #include <set> #include <stack> #include <queue> #include <map> #include <bitset> #include <complex> #include <cctype> #include <utility> #include <climits> #include <numeric> using namespace std; using ll = long long; using P = pair<ll,ll>; const ll MOD = 998244353; const ll MODx = 1000000007; const int INF = (1<<30)-1; const ll LINF = (1LL<<62LL)-1; const double EPS = (1e-10); P ar4[4] = {{0,1},{0,-1},{1,0},{-1,0}}; P ar8[8] = {{-1,-1},{-1,0},{-1,1},{0,-1},{0,1},{1,-1},{1,0},{1,1}}; template <typename T> vector<T> make_vector(size_t a, T b) { return vector<T>(a, b); } template <typename... Ts> auto make_vector(size_t a, Ts... ts) { return vector<decltype(make_vector(ts...))>(a, make_vector(ts...)); } /* 確認ポイント cout << fixed << setprecision(n) << 小数計算//n桁の小数表記になる 計算量は変わらないが楽できるシリーズ min(max)_element(iter,iter)で一番小さい(大きい)値のポインタが帰ってくる count(iter,iter,int)でintがiterからiterの間にいくつあったかを取得できる */ /* function corner below */ /* Function corner above */ /* comment outed because can cause bugs __attribute__((constructor)) void initial() { cin.tie(0); ios::sync_with_stdio(false); } */ template <int mod> struct ModInt{ int n; ModInt():n(0){} ModInt(int n_):n(n_ >= 0 ? n_%mod : mod - ((-n_)%mod) ){} ModInt &operator+=(const ModInt &p){ if((n+=p.n) >= mod)n-=mod; return *this; } ModInt &operator-=(const ModInt &p){ n+=mod-p.n; if(n >= mod)n-=mod; return *this; } ModInt &operator*=(const ModInt &p){ n = (int) ((1LL*n*p.n)%mod); return *this; } ModInt &operator/=(const ModInt &p){ *this *= p.inverse(); return *this; } ModInt operator-() const {return ModInt(-n);} ModInt operator+(const ModInt &p) const {return ModInt(*this) += p;} ModInt operator-(const ModInt &p) const {return ModInt(*this) -= p;} ModInt operator*(const ModInt &p) const {return ModInt(*this) *= p;} ModInt operator/(const ModInt &p) const {return ModInt(*this) /= p;} bool operator==(const ModInt &p) const {return n==p.n;} bool operator<(const ModInt &p) const {return n<p.n;} bool operator>(const ModInt &p) const {return n>p.n;} bool operator>=(const ModInt &p) const {return n>=p.n;} bool operator<=(const ModInt &p) const {return n<=p.n;} bool operator!=(const ModInt &p) const {return n!=p.n;} ModInt inverse() const { int a = n,b = mod,u = 1,v = 0; while(b){ int t = a/b; a -= t*b; swap(a,b); u -= t*v; swap(u,v); } return ModInt(u); } ModInt pow(int64_t z) const { ModInt ret(1),mul(n); while(z > 0){ if(z & 1) ret *= mul; mul *= mul; z >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const ModInt &p){ return os << p.n; } friend istream &operator>>(istream &is, ModInt &a){ int64_t t; is >> t; a = ModInt<mod> (t); return (is); } }; using mint = ModInt<MODx>; mint dp[10010][3][3][2]; int main(){ string n;cin>>n; dp[0][0][0][1] = 1; for(int i = 0; n.size() > i; i++){ for(int j = 1; 10 > j; j++){ bool equal = j == n[i]-'0'; if(j == 2 || j == 6){ for(int k = 0; 3 > k; k++){ for(int l = 0; 3 > l; l++){ if(equal){ dp[i+1][min(2, k+1)][l][1] += dp[i][k][l][1]; dp[i+1][min(2, k+1)][l][0] += dp[i][k][l][0]; }else{ if(j > n[i]-'0'){ dp[i+1][min(2, k+1)][l][0] += dp[i][k][l][0]; }else{ dp[i+1][min(2, k+1)][l][0] += dp[i][k][l][0] + dp[i][k][l][1]; } } } } }else if(j == 4 || j == 8){ for(int k = 0; 3 > k; k++){ for(int l = 0; 3 > l; l++){ if(equal){ dp[i+1][min(2, k+2)][l][1] += dp[i][k][l][1]; dp[i+1][min(2, k+2)][l][0] += dp[i][k][l][0]; }else{ if(j > n[i]-'0'){ dp[i+1][min(2, k+2)][l][0] += dp[i][k][l][0]; }else{ dp[i+1][min(2, k+2)][l][0] += dp[i][k][l][0] + dp[i][k][l][1]; } } } } }else if(j == 5){ for(int k = 0; 3 > k; k++){ for(int l = 0; 3 > l; l++){ if(equal){ dp[i+1][k][min(2, l+1)][1] += dp[i][k][l][1]; dp[i+1][k][min(2, l+1)][0] += dp[i][k][l][0]; }else{ if(j > n[i]-'0'){ dp[i+1][k][min(2, l+1)][0] += dp[i][k][l][0]; }else{ dp[i+1][k][min(2, l+1)][0] += dp[i][k][l][0] + dp[i][k][l][1]; } } } } }else{ for(int k = 0; 3 > k; k++){ for(int l = 0; 3 > l; l++){ if(equal){ dp[i+1][k][l][1] += dp[i][k][l][1]; dp[i+1][k][l][0] += dp[i][k][l][0]; }else{ if(j > n[i]-'0'){ dp[i+1][k][l][0] += dp[i][k][l][0]; }else{ dp[i+1][k][l][0] += dp[i][k][l][0] + dp[i][k][l][1]; } } } } } } dp[i+1][0][0][0]+=1; } cout << dp[n.size()][2][2][1] + dp[n.size()][2][2][0] << endl; return 0; }