結果
問題 | No.1887 K Consecutive Ks (Easy) |
ユーザー | niuez |
提出日時 | 2022-03-25 22:04:54 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 162 ms / 2,000 ms |
コード長 | 3,232 bytes |
コンパイル時間 | 2,115 ms |
コンパイル使用メモリ | 202,664 KB |
実行使用メモリ | 73,592 KB |
最終ジャッジ日時 | 2023-08-04 08:58:49 |
合計ジャッジ時間 | 4,648 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge15 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 161 ms
73,592 KB |
testcase_03 | AC | 1 ms
4,380 KB |
testcase_04 | AC | 1 ms
4,376 KB |
testcase_05 | AC | 1 ms
4,376 KB |
testcase_06 | AC | 1 ms
4,380 KB |
testcase_07 | AC | 1 ms
4,376 KB |
testcase_08 | AC | 2 ms
4,376 KB |
testcase_09 | AC | 2 ms
4,380 KB |
testcase_10 | AC | 2 ms
4,376 KB |
testcase_11 | AC | 23 ms
16,280 KB |
testcase_12 | AC | 12 ms
8,232 KB |
testcase_13 | AC | 153 ms
70,380 KB |
testcase_14 | AC | 162 ms
73,068 KB |
testcase_15 | AC | 160 ms
72,964 KB |
testcase_16 | AC | 129 ms
51,056 KB |
testcase_17 | AC | 91 ms
51,048 KB |
testcase_18 | AC | 161 ms
73,224 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; using i64 = long long; #define rep(i,s,e) for(i64 (i) = (s);(i) < (e);(i)++) #define all(x) x.begin(),x.end() #define STRINGIFY(n) #n #define TOSTRING(n) STRINGIFY(n) #define PREFIX "#" TOSTRING(__LINE__) "| " #define debug(x) \ { \ std::cout << PREFIX << #x << " = " << x << std::endl; \ } std::ostream& output_indent(std::ostream& os, int ind) { for(int i = 0; i < ind; i++) os << " "; return os; } template<class S, class T> std::ostream& operator<<(std::ostream& os, const std::pair<S, T>& p); template<class T> std::ostream& operator<<(std::ostream& os, const std::vector<T>& v); template<class S, class T> std::ostream& operator<<(std::ostream& os, const std::pair<S, T>& p) { return (os << "(" << p.first << ", " << p.second << ")"); } template<class T> std::ostream& operator<<(std::ostream& os, const std::vector<T>& v) { os << "["; for(int i = 0;i < v.size();i++) os << v[i] << ", "; return (os << "]"); } template<class T> static inline std::vector<T> ndvec(size_t&& n, T val) { return std::vector<T>(n, std::forward<T>(val)); } template<class... Tail> static inline auto ndvec(size_t&& n, Tail&&... tail) { return std::vector<decltype(ndvec(std::forward<Tail>(tail)...))>(n, ndvec(std::forward<Tail>(tail)...)); } i64 gcd(i64 a, i64 b) { if(b == 0) return a; return gcd(b, a % b); } #include <iostream> using i64 = long long; template<i64 M> struct modint { i64 a; constexpr modint(const i64 x = 0): a((x%M+M)%M){} constexpr i64 value() const { return a; } constexpr modint inv() const { return this->pow(M-2); } constexpr modint pow(i64 r) const { modint ans(1); modint aa = *this; while(r) { if(r & 1) ans *= aa; aa *= aa; r >>= 1; } return ans; } constexpr modint& operator=(const i64 r) { a = (r % M + M) % M; return *this; } constexpr modint& operator+=(const modint r) { a += r.a; if(a >= M) a -= M; return *this; } constexpr modint& operator-=(const modint r) { a -= r.a; if(a < 0) a += M; return *this; } constexpr modint& operator*=(const modint r) { a = a * r.a % M; return *this; } constexpr modint& operator/=(const modint r) { (*this) *= r.inv(); return *this; } constexpr modint operator+(const modint r) const { return modint(*this) += r; } constexpr modint operator-(const modint r) const { return modint(*this) -= r; } constexpr modint operator*(const modint r) const { return modint(*this) *= r; } constexpr modint operator/(const modint r) const { return modint(*this) /= r; } constexpr bool operator!=(const modint r) const { return this->value() != r.value(); } }; template<const i64 M> std::ostream& operator<<(std::ostream& os, const modint<M>& m) { os << m.value(); return os; } using fp = modint<998244353>; int main() { i64 N, M; cin >> N >> M; auto dp = ndvec(N + 1, M + 1, fp(0)); dp[0][0] = fp(1); rep(i,0,N) { fp sum; rep(j,0,M + 1) { sum += dp[i][j]; } rep(j,1,M + 1) { dp[i + 1][j] += sum - dp[i][j]; dp[i + 1][j] += dp[i][j]; // not equal assign if(i + j <= N) { dp[i + j][j] -= sum - dp[i][j]; } } } fp ans; rep(j,1,M + 1) { ans += dp[N][j]; } cout << fp(M).pow(N) - ans << endl; }