#include #include using namespace std; using ll = long long; using ld = long double; using mint = atcoder::modint998244353; using Pair = pair; using Tuple = tuple; using VI1 = vector; using VI2 = vector; using VL1 = vector; using VL2 = vector; using VD1 = vector; using VD2 = vector; using VB1 = vector; using VB2 = vector; using VP1 = vector; using VP2 = vector; using VT1 = vector; using VT2 = vector; using VM1 = vector; using VM2 = vector; using Queue = queue; using DQ = deque; using PQ = priority_queue, greater>; using Table = VI2; using Graph = VI2; /** io */ template std::vector input_vec(int N); template void output_row(std::vector &row); template void output_col(std::vector &col); void outputYesNo(bool yes, const string &Yes = "Yes", const string &No = "No"); /** minmax */ template bool chmin(T &a, T b); template bool chmax(T &a, T b); /** * 組み合わせと順列を計算するクラス * * @tparam T 組み合わせ・順列を計算する体 (field) * * T の例: atcoder::modint * * 参考文献 * https://github.com/atcoder/live_library/blob/master/comb.cpp */ template struct modcombi { /** 組み合わせ・順列の引数の上限 */ int sup; /** 階乗。factorial の略 */ vector fac; /** 階乗の逆元。inverse の略 */ vector inv; /** * コンストラクタ * * @param[in] sup 組み合わせ・順列の引数の上限 * * 階乗と階乗の逆元の計算を行う。 * 計算量 O(sup) */ modcombi(int sup) : sup(sup) { assert(1 <= sup && sup <= int(1e+7)); fac.resize(sup + 1); inv.resize(sup + 1); fac.at(0) = 1; for (int i = 0; i < sup; ++i) fac.at(i + 1) = fac.at(i) * (i + 1); inv.at(sup) = 1 / fac.at(sup); for (int i = sup; i > 0; --i) { inv.at(i - 1) = inv.at(i) * i; assert(fac.at(i) * inv.at(i) == 1); } } /** * コピーコンストラクタ * * 計算量 O(sup) */ modcombi(const modcombi &org) : sup(org.sup), fac(org.fac), inv(org.inv) {} /** * 組み合わせ (combination) の計算 * * @param[in] n nCk の n * @param[in] k nCk の k * @return nCk * * 計算量 O(1) */ T combi(int n, int k) { if (n < 0 || sup < n) return 0; if (k < 0 || n < k) return 0; return fac.at(n) * inv.at(k) * inv.at(n - k); } /** * 順列 (permutation) の計算 * * @param[in] n nPk の n * @param[in] k nPk の k * @return nPk * * 計算量 O(1) */ T perm(int n, int k) { if (n < 0 || sup < n) return 0; if (k < 0 || n < k) return 0; return fac.at(n) * inv.at(n - k); } }; using MC = modcombi; auto solve() { int X, Y, Z, W; cin >> X >> Y >> Z >> W; MC mc(X + Y + Z + W); auto V = X + Y - Z - W - 1; if (Z == 0) { swap(X, Y); swap(Z, W); } mint cnt = mc.combi(V, Y - 1); cnt *= mc.perm(Y, Y); cnt *= mc.perm(X, X - Z); return cnt.val(); } int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); int t = 1; // cin >> t; while (t--) { auto result = solve(); cout << result << '\n'; // output_row(result); // output_col(result); // outputYesNo(result, "Yes", "No"); } } /** @note 使用頻度が高く毎回貼付するのが面倒なライブラリを実装しておく。*/ template std::vector input_vec(int N) { std::vector v(N); for (auto &vi : v) { cin >> vi; } return v; } void outputYesNo(bool yes, const string &Yes, const string &No) { if (yes) cout << Yes << '\n'; else cout << No << '\n'; } template void output_row(std::vector &row) { int N = row.size(); for (int i = 0; i < N; ++i) { if (i > 0) cout << ' '; cout << row.at(i); } cout << '\n'; } template void output_col(std::vector &col) { int N = col.size(); for (int i = 0; i < N; ++i) { cout << col.at(i) << '\n'; } } template bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; } template bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; }