結果

問題 No.2381 Gift Exchange Party
ユーザー cho435
提出日時 2023-07-14 21:55:08
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 833 bytes
コンパイル時間 4,155 ms
コンパイル使用メモリ 254,028 KB
最終ジャッジ日時 2025-02-15 13:59:06
ジャッジサーバーID
（参考情報） 		judge4 / judge6
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14 RE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using ll = long long;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
using mint = atcoder::modint998244353;

// combination MOD
const int MOD=998244353;
const int MAX_N=2e5;
vector<mint> fac(MAX_N+1,1);
vector<ll> finv(MAX_N+1,1);
vector<ll> inv(MAX_N+1,1);
void comb_setup(){
	for(int i=2;i<=MAX_N;i++){
		fac.at(i)=(fac.at(i-1)*i);
		inv.at(i)=MOD-(inv.at(MOD%i)*(MOD/i))%MOD;
		finv.at(i)=(finv.at(i-1)*inv.at(i))%MOD;
	}
}
mint comb(int n,int k){
	if(n<k) return 0;
	return fac.at(n)*finv.at(n-k)*finv.at(k);
}

int main(){
	comb_setup();
	int n,p;
	cin>>n>>p;
	mint ans=fac.at(n);
	for(ll i=0;i*p<=n;i++){
		mint tmp=comb(n,i*p);
		tmp*=fac.at(i*p);
		tmp/=fac.at(p).pow(i);
		tmp/=fac.at(i);
		tmp*=fac.at(p-1).pow(i);
		ans-=tmp;
	}
	cout<<ans.val()<<endl;
}
0