結果
問題 | No.1960 Guruguru Permutation |
ユーザー |
|
提出日時 | 2022-05-28 14:30:36 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 33 ms / 2,000 ms |
コード長 | 1,332 bytes |
コンパイル時間 | 1,941 ms |
コンパイル使用メモリ | 197,520 KB |
最終ジャッジ日時 | 2025-01-29 16:44:13 |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/modint>using namespace std;using namespace atcoder;using ll = long long;using mint = modint998244353;#define rep(i,n) for (int i = 0; i < n; i++)#define rrep(i,n) for (int i = n-1; i >= 0; i--)#define rep2(i,a,b) for (int i = a; i < b; i++)#define rrep2(i,a,b) for (int i = a-1; i >= b; i--)#define rep3(i,a,b,c) for (int i = a; i < b; i+=c)#define rrep3(i,a,b,c) for (int i = a-1; i >= b; i-=c)#define all(v) v.begin(),v.end()#define rall(v) v.rbegin(),v.rend()template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;}template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;}const int MAX=200010;vector<mint>fac(MAX),finv(MAX);void set_fac(){fac[0]=fac[1]=1;finv[0]=finv[1]=1;for (int i=2;i<MAX;i++){fac[i]=fac[i-1]*i;finv[i]=finv[i-1]/i;}}mint cmb_mint(int n,int r){if (r<0||n<r)return 0;return fac[n]*finv[n-r]*finv[r];}int n,m,k;int main(){ios::sync_with_stdio(false);cin.tie(nullptr);set_fac();cin>>n>>m>>k;mint ans=0;if(m+k>n){m=n-m;k=n-k;n=m+k;}mint f=fac[n]/fac[m+k];rep(i,min(m,k)+1){ans+=cmb_mint(m,i)*cmb_mint(k,i)*fac[i]*f;}cout<<ans.val()<<"\n";}