結果

問題 No.1191 数え上げを愛したい（数列編）
ユーザー uzzyuzzy
提出日時 2020-08-22 17:06:53
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 11 ms / 2,000 ms
コード長 1,556 bytes
コンパイル時間 2,301 ms
コンパイル使用メモリ 182,456 KB
実行使用メモリ 8,320 KB
最終ジャッジ日時 2023-08-05 13:56:38
合計ジャッジ時間 3,605 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 9 ms
8,200 KB
testcase_01 AC 10 ms
8,320 KB
testcase_02 AC 9 ms
8,060 KB
testcase_03 AC 9 ms
8,188 KB
testcase_04 AC 11 ms
8,064 KB
testcase_05 AC 9 ms
8,124 KB
testcase_06 AC 10 ms
8,200 KB
testcase_07 AC 10 ms
8,060 KB
testcase_08 AC 10 ms
8,080 KB
testcase_09 AC 10 ms
8,072 KB
testcase_10 AC 10 ms
8,100 KB
testcase_11 AC 11 ms
8,112 KB
testcase_12 AC 9 ms
8,108 KB
testcase_13 AC 10 ms
8,244 KB
testcase_14 AC 11 ms
8,060 KB
testcase_15 AC 10 ms
8,036 KB
testcase_16 AC 10 ms
8,260 KB
testcase_17 AC 9 ms
8,056 KB
testcase_18 AC 9 ms
8,136 KB
testcase_19 AC 10 ms
8,104 KB
testcase_20 AC 9 ms
8,064 KB
testcase_21 AC 9 ms
8,056 KB
testcase_22 AC 9 ms
8,240 KB
testcase_23 AC 9 ms
8,188 KB
testcase_24 AC 9 ms
8,112 KB
testcase_25 AC 9 ms
8,060 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0