結果
| 問題 |
No.1191 数え上げを愛したい(数列編)
|
| コンテスト | |
| ユーザー |
Shibuyap
|
| 提出日時 | 2020-08-22 16:23:31 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 49 ms / 2,000 ms |
| コード長 | 1,265 bytes |
| コンパイル時間 | 2,236 ms |
| コンパイル使用メモリ | 192,256 KB |
| 最終ジャッジ日時 | 2025-01-13 10:09:51 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 26 |
ソースコード
#include <bits/stdc++.h>
#define rep(i,n) for(ll i = 0; i < (n); ++i)
#define srep(i,s,t) for(ll i = s; i < t; ++i)
#define drep(i,n) for(ll i = (n)-1; i >= 0; --i)
using namespace std;
typedef long long int ll;
typedef pair<int,int> P;
#define yn {puts("Yes");}else{puts("No");}
#define MAX_N 200005
const int MAX = 1110000;
const int MOD = 998244353;
long long fac[MAX], finv[MAX], inv[MAX];
// テーブルを作る前処理
void COMinit() {
fac[0] = fac[1] = 1;
finv[0] = finv[1] = 1;
inv[1] = 1;
for (int i = 2; i < MAX; 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 FINV(int n){
if (n < 0) return 0;
return finv[n];
}
int main() {
COMinit();
ll n, m, a, b;
cin >> n >> m >> a >> b;
ll ans = 0;
srep(i,a*(n-1),b+1){
ll amari = (ll)i - a*(n-1);
ll tmp = COM(amari+n-2,n-2) * (m-i) % MOD;
ans += tmp;
ans %= MOD;
}
ans *= fac[n];
ans %= MOD;
cout << ans << endl;
return 0;
}
Shibuyap