結果
問題 | No.2351 Butterfly in Summer |
ユーザー | planes |
提出日時 | 2023-06-16 21:32:28 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 7 ms / 2,000 ms |
コード長 | 1,325 bytes |
コンパイル時間 | 1,625 ms |
コンパイル使用メモリ | 171,692 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-06-24 13:09:02 |
合計ジャッジ時間 | 2,451 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll =long long; #define all(v) v.begin(),v.end() #define rep(i,a,b) for(int i=a;i<b;i++) #define rrep(i,a,b) for(int i=a;i>=b;i--) ll INF=2e18; const ll MOD = 998244353; vector<long long> fact, fact_inv, inv; void init_nCk(int SIZE) { fact.resize(SIZE + 5); fact_inv.resize(SIZE + 5); inv.resize(SIZE + 5); fact[0] = fact[1] = 1; fact_inv[0] = fact_inv[1] = 1; inv[1] = 1; for (int i = 2; i < SIZE + 5; i++) { fact[i] = fact[i - 1] * i % MOD; inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD; fact_inv[i] = fact_inv[i - 1] * inv[i] % MOD; } } long long nCk(int n, int k) { if((n < k)) return 0; if ((n < 0 || k < 0)) return 0; return fact[n] * (fact_inv[k] * fact_inv[n - k] % MOD) % MOD; } long long nPk(int n, int k) { if((n < k)) return 0; if ((n < 0 || k < 0)) return 0; return fact[n] * fact_inv[n - k] % MOD; } ll mod_pow(ll x,ll n,ll mod) { ll res=1; while(n>0) { if(n&1) { res=res*x%mod; } x=x*x%mod; n>>=1; } return res; } int main() { ios::sync_with_stdio(false); cin.tie(0); ll N,K;cin>>N>>K; init_nCk(K+100); ll ans=nCk(K,2)*N*2%MOD; ll x=mod_pow(K,MOD-2,MOD); for(ll i=0;i<N;i++) { ans*=x; ans%=MOD; } if(ans<0) ans+=MOD; cout<<ans<<endl; }