#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a,i##E=b;i<=i##E;i++)
#define UF(i,a,b) for(int i=a,i##E=b;i>=i##E;i--)
#define sz(x) int(x.size())
using namespace std;typedef long long ll;
#define vec vector
#define In(x) freopen(x".in","r",stdin)
#define Out(x) freopen(x".out","w",stdout)
#define File(x) (In(x),Out(x))
#define all(x) begin(x),end(x)
#define ran(x,l,r) begin(x)+l,begin(x)+(r+1)
template<typename T>
void operator+=(vec<T> &a,const T&b) {a.push_back(b);}
#define gc() getchar()
//#define gc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
//char *p1,*p2,buf[1<<21];
int read() {
	int s=0,w=0;char ch=gc();
	while(ch<'0'||ch>'9') w|=(ch=='-'),ch=gc();
	while(ch>='0'&&ch<='9') s=(s<<3)+(s<<1)+(ch^48),ch=gc();
	return w?-s:s;
} const int p=998244353;
using calc_t=int;
calc_t sol(int A,int B,vec<pair<int,int>> E) {
    if(!sz(E)) return 1;
    if(__gcd(A,B)>1) {
        int g=__gcd(A,B);
        A/=g,B/=g;
        calc_t res=1;
        F(i,0,g-1) {
            vec<pair<int,int>> nE;
            for(auto[u,v]:E) if(u%g==i) nE+=make_pair(u/g,v/g);
            res=1ll*res*sol(A,B,nE)%p;
        }
        return res;
    }
    int n=0;
    for(auto[u,v]:E) n=max(n,max(u,v));
    ++n;
    int mx=max(n,(n+B-2)/B*B);//mx=mx*mx;
    vec<int> pA(mx),pB(mx);
    for(auto[u,v]:E) {
        assert(v-u==A||v-u==B);
        if(v-u==A) pA[v]=1;
        else pB[v]=1;
    }
    // cerr<<"SOL:"<<n<<' '<<A<<' '<<B<<endl;
    // for(auto[a,b]:E) cerr<<a<<' '<<b<<endl;
    if(B<=sqrt(2*n)) {
        // cerr<<"mode 1;"<<endl;
        vec<calc_t> f(1<<B,0);
        const int U=(1<<B)-1;
        f[U]=1;
        #define Z(x,y) (x+=y)%=p
        F(i,0,n-1) {
            vec<calc_t> g(1<<B,0);
            if(pA[i]) {
                F(j,0,U) if(!(j>>A-1&1)) Z(g[(j<<1|1|1<<A)&U],f[j]);
            } 
            if(pB[i]) {
                F(j,0,U) if(!(j>>B-1&1)) Z(g[(j<<1|1)&U],f[j]);
            }
            F(j,0,U) Z(g[(j<<1)&U],f[j]);
            swap(f,g);
        }
        calc_t ans=0;
        F(j,0,U) Z(ans,f[j]);
        return ans;
    } 
    // cerr<<"mode 2;"<<endl;
    int G=(n+B-2)/B;
    const int U=(1<<G)-1;
    int ans=0;
    F(C,0,(1<<G)-1) {
        vec<calc_t> f(1<<G,0);
        // if(C&1) continue;
        // cerr<<"ST:"<<C<<endl;
        f[C]=1;
        int x=0;
        F(_,1,B) {
            x+=A;
            if(x>=B) {
                x-=B;
                vec<calc_t> g(1<<G,0);
                F(j,0,U) Z(g[(j<<1)&U],f[j]);
                swap(f,g);
            }
            F(i,0,G-1) {
                vec<calc_t> g(1<<G,0);
                int k=x+i*B;
                if(pA[k]) {
                    F(j,0,U) if(!(j>>i&1)) Z(g[j|1<<i],f[j]);
                }
                if(pB[k]) {
                    assert(i);
                    F(j,0,U) if(!(j>>i-1&1)) Z(g[j|1<<i-1|1<<i],f[j]);
                }
                F(j,0,U) Z(g[j&~(1<<i)],f[j]);
                swap(f,g);
                // cerr<<"AFT_K:"<<k<<' '<<endl;
                // F(j,0,U) cerr<<"R:"<<j<<' '<<f[j]<<endl;
            }
        }
        assert(!x);
        Z(ans,f[C]);
    } assert(ans>0);
    return ans;
}
int main() {
    // File("match");
    // In("match_ex3");
    int n,m,A,B;
    n=read();m=read();A=read();B=read();
    vec<pair<int,int>> E;
    F(i,1,m) {
        int u=read()-1,v=read()-1;
        assert(v<n);
        E+=make_pair(u,v);
    }
    cout<<sol(A,B,E)<<endl;
	return 0;
}

/*
B<=20 好做，需要解决 B>20
按 B 分块，块数 <=10，块内按 mod A 同余分类。
考虑按照 (x+=A)%=B 的顺序处理所有 mod B 的同余类，此时类很多但是每类大小很小
复杂度 2^20.
*/