結果
| 問題 |
No.1191 数え上げを愛したい(数列編)
|
| コンテスト | |
| ユーザー |
uzzy
|
| 提出日時 | 2020-08-22 17:06:53 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 11 ms / 2,000 ms |
| コード長 | 1,556 bytes |
| コンパイル時間 | 2,333 ms |
| コンパイル使用メモリ | 184,752 KB |
| 実行使用メモリ | 8,192 KB |
| 最終ジャッジ日時 | 2024-10-15 11:03:13 |
| 合計ジャッジ時間 | 3,377 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 26 |
ソースコード
#pragma GCC optimize ("O2")
#pragma GCC target ("avx2")
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i = 0; i < (n); i++)
#define rep1(i, n) for(int i = 1; i <= (n); i++)
#define co(x) cout << (x) << "\n"
#define cosp(x) cout << (x) << " "
#define ce(x) cerr << (x) << "\n"
#define cesp(x) cerr << (x) << " "
#define pb push_back
#define mp make_pair
#define chmin(x, y) x = min(x, y)
#define chmax(x, y) x = max(x, y)
#define Would
#define you
#define please
const int mod = 998244353;
const int IMAX = 400000;
constexpr ll modpow(ll A, ll B) {
ll kotae = 1;
while (B > 0) {
if (B & 1) kotae = kotae * A % mod;
A = A * A % mod;
B >>= 1;
}
return kotae;
}
class unko {
public:
int fac[IMAX + 1], ifac[IMAX + 1], inv[IMAX + 1];
unko() : fac(), ifac(), inv() {
fac[0] = 1;
rep1(i, IMAX) fac[i] = ll(fac[i - 1]) * i % mod;
ifac[IMAX] = modpow(fac[IMAX], mod - 2);
for (int i = IMAX; i > 0; i--) ifac[i - 1] = ll(ifac[i]) * i % mod;
inv[0] = 1;
rep1(i, IMAX) inv[i] = ll(ifac[i]) * fac[i - 1] % mod;
}
int comb(int A, int B) {
return ll(fac[A]) * ifac[B] % mod * ifac[A - B] % mod;
}
} f;
int main() {
cin.tie(0);
ios::sync_with_stdio(false);
ll N, M, A, B;
cin >> N >> M >> A >> B;
ll kotae = 0;
for (int D = 1; D <= B; D++) {
ll are = D - A * (N - 1);
if (are >= 0) {
ll tmp = f.comb(are + N - 2, N - 2);
tmp = tmp * (M - D) % mod;
kotae += tmp;
}
}
kotae %= mod;
rep1(i, N) {
kotae = kotae * i % mod;
}
co(kotae);
Would you please return 0;
}
uzzy