結果
| 問題 |
No.2383 Naphthol
|
| コンテスト | |
| ユーザー |
matsup10
|
| 提出日時 | 2023-07-15 01:10:22 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,727 bytes |
| コンパイル時間 | 5,308 ms |
| コンパイル使用メモリ | 303,996 KB |
| 最終ジャッジ日時 | 2025-02-15 14:53:45 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 13 WA * 5 |
ソースコード
#include <bits/extc++.h>
#include <atcoder/all>
using namespace std;
using ll = long long;
#define REP(i,n) for(int i=0;i<int(n);i++)
#define FOR(i,a,b) for(int i=a;i<=int(b);i++)
#define ALL(x) x.begin(),x.end()
#define INF INT_MAX
#define INFL LLONG_MAX / 4
using namespace atcoder;
using mint = modint998244353;
template<typename T> void chmin(T& a, T b) { a = min(a, b); }
template<typename T> void chmax(T& a, T b) { a = max(a, b); }
#define PR(x) cerr << #x << "=" << x << endl
using i128 = __int128_t;
struct nCkmod {
vector<long long> fanc;
vector<long long> fancinv;
long long mod;
nCkmod(long long n, long long mod_ = 998244353): fanc(n + 1), fancinv(n + 1), mod(mod_) {
fanc[0] = 1;
fancinv[0] = 1;
for(long long i = 1; i <= n; i++) {
fanc[i] = (fanc[i - 1] * i) % mod;
fancinv[i] = (fancinv[i - 1] * modinv(i)) % mod;
}
}
long long modinv(long long a) {
long long b = mod, u = 1, v = 0;
while (b) {
long long t = a / b;
a -= t * b; swap(a, b);
u -= t * v; swap(u, v);
}
u %= mod;
if (u < 0) u += mod;
return u;
}
long long nck(int n, int k) {
if(n < k) return 0;
return fanc[n] * fancinv[k] % mod * fancinv[n - k] % mod;
}
};
int main() {
ll n, k;
cin >> n >> k;
nCkmod nck(2 * n + 10);
mint ans = nck.nck(2 * n + 4, k);
if(k % 2 == 0) {
ans += nck.nck(n + 2, k / 2);
if(n % 2 == 0) {
ans += nck.nck(n + 2, k / 2);
} else {
ans += nck.nck(n+1, k / 2);
ans += nck.nck(n+1, k / 2 - 1);
}
ans += nck.nck(n + 2, k / 2);
} else {
if(n % 2 == 0) {
} else {
ans += 2 * nck.nck(n+1, k / 2);
}
}
ans /= 4;
cout << ans.val();
return 0;
}
matsup10