結果

問題 No.2381 Gift Exchange Party
コンテスト
ユーザー random contestant
提出日時 2023-07-15 15:35:06
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 109 ms / 2,000 ms
コード長 1,506 bytes
コンパイル時間 4,156 ms
コンパイル使用メモリ 251,472 KB
最終ジャッジ日時 2025-02-15 15:04:48
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

// #define _GLIBCXX_DEBUG
#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
using ll = long long;
#define rep(i,n) for (ll i = 0; i < (n); ++i)
using vl = vector<ll>;
using vvl = vector<vl>;
using P = pair<ll,ll>;
#define pb push_back
#define int long long
#define double long double
#define INF (ll) 3e18
// Ctrl + Shift + B コンパイル
// Ctrl + C 中断
// ./m 実行


// コンビネーション計算 nCr v3 570ms
constexpr int mod = 998244353;
using mint = modint998244353;
vl kaijo(200010);
vl kaijo_inv(200010);
void kaijo_calc(){
  kaijo[0] = kaijo[1] = 1;
  for(int i = 2; i <= 200000; ++i){
    kaijo[i] = kaijo[i-1] * i;
    kaijo[i] %= mod;
  }
  for(int i = 0; i <= 200000; ++i){
    kaijo_inv[i] = inv_mod(kaijo[i], mod);
  }
}
int nCr(int n, int r){
  if (kaijo[0] == 0) kaijo_calc();
  mint res = 1;
  res *= kaijo[n];
  res *= kaijo_inv[n-r];
  res *= kaijo_inv[r];
  return res.val();
}

signed main(){
  int n; cin >> n;
  int p; cin >> p; // pは素数
  using mint = modint998244353;
  mint ans = 0;
  for(int k = 0; k <= n/p; k++){
    // 長さpのサイクルがk個あり、残りは自己ループとなる。
    mint now = 1;
    now *= nCr(n, k*p); // サイクルとするk*p個の取り出し方
    now *= kaijo[k*p]; // k+pの並べ替え方
    now /= mint(p).pow(k); // サイクルの並べ替えは(P-1)!なのでPで割る
    now /= kaijo[k];
    ans += now;
  }
  ans = kaijo[n] - ans;
  cout << ans.val() << endl;
}
0