結果
問題 | No.986 Present |
ユーザー |
|
提出日時 | 2020-02-11 15:00:18 |
言語 | C++17(clang) (17.0.6 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 59 ms / 2,000 ms |
コード長 | 3,147 bytes |
コンパイル時間 | 640 ms |
コンパイル使用メモリ | 120,192 KB |
実行使用メモリ | 7,168 KB |
最終ジャッジ日時 | 2024-11-30 16:22:42 |
合計ジャッジ時間 | 2,538 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <iostream>#include <vector>template <int MOD>struct ModInt {using lint = long long;int val;// constructorModInt(lint v = 0) : val(v % MOD) {if (val < 0) val += MOD;};// unary operatorModInt operator+() const { return ModInt(val); }ModInt operator-() const { return ModInt(MOD - val); }ModInt inv() const { return this->pow(MOD - 2); }// arithmeticModInt operator+(const ModInt& x) const { return ModInt(*this) += x; }ModInt operator-(const ModInt& x) const { return ModInt(*this) -= x; }ModInt operator*(const ModInt& x) const { return ModInt(*this) *= x; }ModInt operator/(const ModInt& x) const { return ModInt(*this) /= x; }ModInt pow(lint n) const {auto x = ModInt(1);auto b = *this;while (n > 0) {if (n & 1) x *= b;n >>= 1;b *= b;}return x;}// compound assignmentModInt& operator+=(const ModInt& x) {if ((val += x.val) >= MOD) val -= MOD;return *this;}ModInt& operator-=(const ModInt& x) {if ((val -= x.val) < 0) val += MOD;return *this;}ModInt& operator*=(const ModInt& x) {val = lint(val) * x.val % MOD;return *this;}ModInt& operator/=(const ModInt& x) { return *this *= x.inv(); }// comparebool operator==(const ModInt& b) const { return val == b.val; }bool operator!=(const ModInt& b) const { return val != b.val; }// I/Ofriend std::istream& operator>>(std::istream& is, ModInt& x) noexcept { return is >> x.val; }friend std::ostream& operator<<(std::ostream& os, const ModInt& x) noexcept { return os << x.val; }};constexpr int MOD = 998244353;using mint = ModInt<MOD>;template <class T>struct Combination {int max_n;std::vector<T> f, invf;explicit Combination(int n): max_n(n), f(n + 1), invf(n + 1) {f[0] = 1;for (int i = 1; i <= n; ++i) {f[i] = f[i - 1] * i;}invf[max_n] = f[max_n].inv();for (int i = max_n - 1; i >= 0; --i) {invf[i] = invf[i + 1] * (i + 1);}}T fact(int n) const { return f[n]; }T invfact(int n) const { return invf[n]; }T perm(int a, int b) const {return a < b || b < 0 ? T(0) : f[a] * invf[a - b];}T comb(int a, int b) const {return a < b || b < 0 ? T(0) : f[a] * invf[a - b] * invf[b];}};const Combination<mint> C(500000);void solve() {int n, m;std::cin >> n >> m;// task 1std::cout << mint(2).pow(n) << " ";// task 2mint ans = 1;for (int k = 0; k < n; ++k) {ans *= mint(2).pow(m) - mint(2).pow(k);}ans /= C.fact(n);std::cout << ans << " ";// task 3mint pat = 1;for (int k = 0; k < n; ++k) {pat *= mint(2).pow(n) - mint(2).pow(k);}pat /= C.fact(n);std::cout << ans / pat << std::endl;}int main() {std::cin.tie(nullptr);std::cout.tie(nullptr);std::ios::sync_with_stdio(false);solve();return 0;}