結果
問題 | No.2141 Enumeratest |
ユーザー | ぷら |
提出日時 | 2022-12-02 21:23:37 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 231 ms / 2,000 ms |
コード長 | 1,176 bytes |
コンパイル時間 | 2,077 ms |
コンパイル使用メモリ | 201,092 KB |
実行使用メモリ | 50,432 KB |
最終ジャッジ日時 | 2024-10-09 22:27:41 |
合計ジャッジ時間 | 8,505 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 34 |
ソースコード
#include <bits/stdc++.h>using namespace std;constexpr int mod = 998244353;long long fac[2000005], finv[2000005], inv[2000005];void COMinit() {fac[0] = fac[1] = finv[0] = finv[1] = inv[1] = 1;for (int i = 2; i < 2000005; i++) {fac[i] = fac[i - 1] * i % mod;inv[i] = mod - inv[mod % i] * (mod / i) % mod;finv[i] = finv[i - 1] * inv[i] % mod;}}long long COM(int n, int k){if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % mod) % mod;}long long choose(int n,int k) {if(n < 0 || k < 0) return 0;if(n == 0) return 1;return COM(n+k-1,k-1);}long long modpow(long long a,long long b) {long long ans = 1;while(b) {if(b & 1) {(ans *= a) %= mod;}(a *= a) %= mod;b /= 2;}return ans;}int main() {COMinit();int N,M;cin >> N >> M;int ans = fac[M];for(int i = 0; i < N; i++) {if(M%N > i) {ans = ans*modpow(fac[M/N+1],mod-2)%mod;}else {ans = ans*modpow(fac[M/N],mod-2)%mod;}}cout << ans << endl;}