結果
問題 | No.2381 Gift Exchange Party |
ユーザー |
|
提出日時 | 2023-07-14 21:33:22 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 9 ms / 2,000 ms |
コード長 | 1,106 bytes |
コンパイル時間 | 1,942 ms |
コンパイル使用メモリ | 191,804 KB |
最終ジャッジ日時 | 2025-02-15 10:44:57 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#include<bits/stdc++.h>using namespace std;using ll=long long;using pii=pair<int,int>;#define all(a) a.begin(),a.end()#define pb push_back#define sz(a) ((int)a.size())const int N=200005,mod=998244353;int add(int x, int y){x+=y; if(x>=mod) x-=mod; return x;}int sub(int x, int y){x-=y; if(x<0) x+=mod; return x;}int mul(int x, int y){return ((ll)x)*y%mod;}int Pow(int x, ll y=mod-2){int res=1; for(; y; x=mul(x,x),y>>=1) if(y&1) res=mul(res,x); return res;}int fac[N],inv[N],ifac[N];inline int C(int n, int m){if(m<0||m>n) return 0; return mul(fac[n],mul(ifac[m],ifac[n-m]));}void init_comb(){fac[0]=inv[1]=ifac[0]=1;for(int i=1; i<N; ++i) fac[i]=mul(fac[i-1],i);for(int i=2; i<N; ++i) inv[i]=mul(inv[mod%i],mod-mod/i);for(int i=1; i<N; ++i) ifac[i]=mul(ifac[i-1],inv[i]);}int n,p,dp[N];signed main(){ios_base::sync_with_stdio(0),cin.tie(0);init_comb();cin >> n >> p;dp[0]=1;for(int i=1; i<=n; ++i){dp[i]=dp[i-1];if(i>=p) dp[i]=add(dp[i],mul(dp[i-p],mul(fac[i-1],ifac[i-p])));}cout << sub(fac[n],dp[n]) << "\n";}