#include <bits/stdc++.h>
#define int long long
#define pii pair<int,int>
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(c) (c).begin(),(c).end()
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
#define MINF(a) memset(a,0x3f,sizeof(a))
#define POW(n) (1LL<<(n))
#define IN(i,a,b) (a <= i && i <= b)
using namespace std;
template <typename T> inline bool CHMIN(T& a,T b) { if(a>b) { a=b; return 1; } return 0; }
template <typename T> inline bool CHMAX(T& a,T b) { if(a<b) { a=b; return 1; } return 0; }
template <typename T> inline void SORT(T& a) { sort(ALL(a)); }
template <typename T> inline void REV(T& a) { reverse(ALL(a)); }
template <typename T> inline void UNI(T& a) { sort(ALL(a)); a.erase(unique(ALL(a)),a.end()); }
template <typename T> inline void LB(vector<T>& v, T& a) { return *lower_bound(ALL(v),a); }
template <typename T> inline void LBP(vector<T>& v, T& a) { return lower_bound(ALL(v),a) - v.begin(); }
template <typename T> inline void UB(vector<T>& v, T& a) { return *upper_bound(ALL(v),a); }
template <typename T> inline void UBP(vector<T>& v, T& a) { return upper_bound(ALL(v),a) - v.begin(); }
template <typename T1, typename T2> ostream& operator<< (ostream& os, const pair<T1,T2>& p) { os << p.first << " " << p.second; return os; }
template <typename T1, typename T2> istream& operator>> (istream& is, pair<T1,T2>& p) { is >> p.first >> p.second; return is; }
template <typename T> ostream& operator<< (ostream& os, const vector<T>& v) { REP(i,v.size()) { if (i) os << " "; os << v[i]; } return os; }
template <typename T> istream& operator>> (istream& is, vector<T>& v) { for(T& in : v) is >> in; return is; }
template <typename T = int> vector<T> make_v(size_t a) { return vector<T>(a); }
template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts) { return vector<decltype(make_v<T>(ts...))>(a,make_v<T>(ts...)); }
template< typename T, typename V > typename enable_if< is_class< T >::value == 0 >::type fill_v(T &t, const V &v) { t = v; }
template< typename T, typename V > typename enable_if< is_class< T >::value != 0 >::type fill_v(T &t, const V &v) { for(auto &e : t) fill_v(e, v); }
const int MOD = 998244353;
const int INF = 0x3f3f3f3f3f3f3f3f;
const double EPS = 1e-10;

int inv_mod(int a, int m = MOD) {
  int b = m, x = 1, y = 0, t;
  while (b > 0) {
    t = a / b;
    swap(a -= t * b, b);
    swap(x -= t * y, y);
  }
  x %= m;
  if (x < 0) x += m;
  return x;
}

const int MAX_N = 3000010;
int fac[MAX_N];
int ifac[MAX_N];
int inv[MAX_N];

void init(int n = MAX_N, int m = MOD) {
  fac[0] = 1;
  for (int i = 1; i < n; i++) fac[i] = fac[i - 1] * i % m;
  ifac[n - 1] = inv_mod(fac[n - 1]);
  for (int i = n - 2; i >= 0; i--) ifac[i] = ifac[i+1] * (i + 1) % m;
  for (int i = 1; i < n; i++) inv[i] = ifac[i] * fac[i - 1] % m;
}

int nPr(int n, int r, int m = MOD) {
  if (r < 0 || n < r) return 0;
  return fac[n] * ifac[n - r] % m;
}

int nCr(int n, int r, int m = MOD) {
  if (r < 0 || n < r) return 0;
  return fac[n] * ifac[r] % m * ifac[n - r] % m;
}

int nHr(int n, int r, int m = MOD) {
  if (n < 0 || r < 0) return 0;
  return r == 0 ? 1 : nCr(n + r - 1, r);
}

int N,A,B,K;

signed main() {
  cin.tie(nullptr);
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(10);

  init();
  cin >> N >> A >> B >> K;
  if (K == 0) {
    cout << 1 << endl;
    return 0;
  }
  int ans = 0;
  REP(i,N+1) {
    int suma = A*i;
    if (suma > K) break;
    int rem = K - suma;
    if (rem % B) continue;
    int num = rem / B;
    int tmp = nCr(N,i);
    if (num > 0) (tmp *= nCr(N,num)) %= MOD;
    (ans += tmp) %= MOD;
  }
  cout << ans << endl;

  return 0;
}

Measured memory_bytes: 70528