結果
問題 |
No.2048 L(I+D)S
|
ユーザー |
|
提出日時 | 2023-01-19 12:58:32 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 6 ms / 2,000 ms |
コード長 | 1,049 bytes |
コンパイル時間 | 1,695 ms |
コンパイル使用メモリ | 193,416 KB |
最終ジャッジ日時 | 2025-02-10 04:11:45 |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 16 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define fi first #define se second #define pb push_back using vi = vector <int>; using ll = long long; using pii = pair <int, int>; const ll mod = 998244353; //~ const ll mod = 1e9 + 7; ll qpow(ll a, ll b, ll m = mod) { ll r = 1, t = a; for(; b; b /= 2) { if(b & 1) r = r * t % m; t = t * t % m; } return r; } const int N = 1e5 + 11; ll inv[N], fac[N], ifac[N]; int main() { inv[1] = 1; for(int i = 2; i < N; i ++) inv[i] = (mod - mod / i) * inv[mod % i] % mod; fac[0] = ifac[0] = 1; for(int i = 1; i < N; i ++) { fac[i] = i * fac[i - 1] % mod; ifac[i] = inv[i] * ifac[i - 1] % mod; } ios :: sync_with_stdio(0); int n; cin >> n; ll ans = 0; if(n >= 4) { for(int i = 2; i <= n - 2; i ++) { int a = i, b = n - i; ll t = fac[n]; t = t * ifac[a - 2] % mod; t = t * ifac[b - 2] % mod; t = t * inv[n - 1] % mod * inv[a] % mod * inv[b] % mod; //~ cout << a << ' ' << b << ' ' << t << '\n'; ans += t * t % mod; ans %= mod; } } cout << ans << '\n'; }