void main() {
	auto ip = readAs!(ulong[]), N = ip[0], A = ip[1], B = ip[2], K = ip[3];
	ulong res;
	foreach(i; 0..ceil(K/A.to!real).to!ulong) {
		foreach(k; 0..ceil(K/B.to!real).to!ulong) {
			foreach(j; 0..ceil(K/(A+B).to!real).to!ulong) {
				if(i*A + j*(A+B) + k*B == K) {
					//writefln("R: %d, G: %d, B: %d", i, j, k);
					ulong tmp = nCr(N, i) * nCr(N-i, j) * nCr(N-i-j, k);
					res = (res + tmp) % 998244353;
				}
			}
		}
	}
	writeln(max(1, res));
}

ulong nCr(ulong n, ulong r)
{
    ulong x = 1;
    foreach (i; n - r + 1 .. n + 1)
        x *= i;
    foreach (i; 1 .. r + 1)
        x /= i;
    return x;
}


// ===================================

import std.stdio;
import std.string;
import std.conv;
import std.algorithm;
import std.range;
import std.traits;
import std.math;
import std.container;
import std.bigint;
import std.numeric;
import std.conv;
import std.typecons;
import std.uni;
import std.ascii;
import std.bitmanip;
import core.bitop;

T readAs(T)() if (isBasicType!T) {
	return readln.chomp.to!T;
}
T readAs(T)() if (isArray!T) {
	return readln.split.to!T;
}

T[][] readMatrix(T)(uint height, uint width) if (!isSomeChar!T) {
	auto res = new T[][](height, width);
	foreach(i; 0..height) {
		res[i] = readAs!(T[]);
	}
	return res;
}

T[][] readMatrix(T)(uint height, uint width) if (isSomeChar!T) {
	auto res = new T[][](height, width);
	foreach(i; 0..height) {
		auto s = rs;
		foreach(j; 0..width) res[i][j] = s[j].to!T;
	}
	return res;
}

int ri() {
	return readAs!int;
}

double rd() {
	return readAs!double;
}

string rs() {
	return readln.chomp;
}