結果
問題 |
No.2409 Strange Werewolves
|
ユーザー |
![]() |
提出日時 | 2023-11-29 00:48:06 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 26 ms / 2,000 ms |
コード長 | 1,209 bytes |
コンパイル時間 | 689 ms |
コンパイル使用メモリ | 72,404 KB |
実行使用メモリ | 18,944 KB |
最終ジャッジ日時 | 2024-09-26 13:13:30 |
合計ジャッジ時間 | 1,775 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 16 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> using namespace std; using ll = long long; constexpr ll MOD = 998244353; struct Combination{ using ll = long long; const int n; // const int MOD = 998244353; vector<ll> fac, ifac; ll mp(ll a, ll x) const { ll res = 1; for(; x > 0; x >>= 1){ if(x&1){ res *= a; res %= MOD; } a *= a; a %= MOD; } return res; } Combination(int _n = 1000000): n(_n), fac(vector<ll>(_n+1)), ifac(vector<ll>(_n+1)){ fac[0] = 1; for(int i = 1; i <= n; i++) fac[i] = (fac[i-1]*i)%MOD; ifac[n] = mp(fac[n], MOD-2); for(int i = n; i > 0; i--) ifac[i-1] = (ifac[i]*i)%MOD; } ll operator()(int n, int k) const { if (k < 0 || k > n) return 0; return fac[n]*ifac[k]%MOD*ifac[n-k]%MOD; } }; int main(){ int x, y, z, w; cin >> x >> y >> z >> w; Combination comb{}; ll ans = (z == 0 ? x : y); if(z == 0) ans *= comb(y, w); else ans *= comb(x, z); ans %= MOD; for(int i = x-z+y-w-1; i > 1; i--){ ans *= i; ans %= MOD; } cout << ans << endl; return 0; }