// 传统计数强省的 998244353 震撼
// 异或卷积呐。
// 集合幂级数呐。

/*
Bell Number:
1
1
2
5
15
52
203
877
4140
21147
115975
678570
4213597
27644437
190899322
1382958545
*/

// 每种等值联通块划分方案的容斥系数为其各部分子图的各联通选法带符号权值积之和之积。
// 先对每个图算出有其带符号权值积之和，容易发现为 1 当且仅当无边。
// 然后枚举最小点所在联通块，容斥掉不连通的情况，得到各个块的各联通选法带符号权值积之和。
// 然后把每个联通块转换成异或卷积形式，显然和联通块中最小值一致，或不对异或值做贡献。
// 可以在联通块最小值处标记该联通块，不对异或值做贡献的可以单独做子集合并 dp。
// 对联通块值为最小值的异或卷积信息如何合并？
// 因为 [0,r) 类型的异或卷积虽支持快速乘法，但不支持快速加法，我们考虑怎么办。
// 容易发现没得办。
// 因此我们写一个 Bell 数级别爆搜即可。
// 判掉 m=0，大概能拿 80pts 左右。
// 哦好像可以做到 O(3^n) 之类的，那没事了（
// 其实就是，把每个联通块最小元素也记入状态，算出方案数。
// 然后把这些最小元素选若干来暴力合并出全集即可。

#include <bits/stdc++.h>
typedef unsigned uint;
typedef unsigned long long ullt;typedef long long llt;
typedef bool bol;typedef char chr;typedef void voi;
typedef double dbl;
template<typename T>bol _min(T&a,T b){return b<a?(a=b,true):false;}
template<typename T>bol _max(T&a,T b){return a<b?(a=b,true):false;}
const ullt Mod=998244353;
struct modint
{
    ullt v;
    modint():v(0){}
    modint(ullt v):v(v%Mod){}
    ullt&operator()(){return v;}
    friend modint operator+(modint a,modint b){return a()+b();}
    friend modint operator-(modint a,modint b){return a()+Mod-b();}
    friend modint operator*(modint a,modint b){return a()*b();}
    modint operator-(){return Mod-v;}
    modint _power(ullt index){
        modint ans(1),base(*this);
        while(index){
            if(index&1)ans*=base;
            base*=base,index>>=1;
        }
        return ans;
    }
    friend modint operator^(modint a,ullt b){return a._power(b);}
    modint inv(){return _power(Mod-2);}
    modint&operator+=(modint b){return*this=*this+b;}
    modint&operator-=(modint b){return*this=*this-b;}
    modint&operator*=(modint b){return*this=*this*b;}
    voi read(){scanf("%llu",&v);}
    voi print(){printf("%llu",v);}
    voi println(){print(),putchar('\n');}
};
using modvec=std::vector<modint>;
using vec=std::vector<std::pair<bol,modint> >;
const ullt Lim=1llu<<60;
const uint Bit=60;
ullt Len[Bit+5];
vec get(ullt r){
    vec Ans(Bit+1);
    for(uint i=0;i<Bit;i++)if(r>=Len[i])Ans[i]={0,1},r-=Len[i];else Ans[i]={1,0};
    Ans[Bit]={0,r?1:0};
    return Ans;
}
vec mul(vec A,vec B){
    modvec UserA(Bit+1),UserB(Bit+1);
    for(uint i=0;i<=Bit;i++)UserA[i]=A[i].second*Len[i],UserB[i]=B[i].second*Len[i];
    for(uint i=Bit;i;i--)UserA[i-1]+=UserA[i],UserB[i-1]+=UserB[i];
    vec Ans(Bit+1);
    modint v;
    for(uint i=0;i<Bit;i++){
        Ans[i]={A[i].first==B[i].first,A[i].second*UserB[i+1]+B[i].second*UserA[i+1]+v};
        v+=A[i].second*B[i].second*Len[i];
    }
    Ans[Bit]={0,A[Bit].second*B[Bit].second+v};
    return Ans;
}
modint query(vec A,ullt p){
    for(uint i=0;i<Bit;i++)if((p>=Len[i])==A[i].first)return A[i].second;else if(p>=Len[i])p-=Len[i];
    return A[Bit].second;
}
ullt A[25],c;
uint E[25],Log2[1u<<15|1],Pop2[1u<<15|1];
modint Dp[1u<<15|1],Dp2[1u<<15|1];
inline uint lowbit(uint v){return v&-v;}
modint T[15][16][1u<<15|1],User[1u<<15|1];
modint Q[16][1u<<15|1];
modint G[1u<<15|1],W[1u<<15|1];
int main()
{
#ifdef MYEE
    freopen("QAQ.in","r",stdin);
    freopen("QAQ.out","w",stdout);
#else
    freopen("graph.in","r",stdin);
    freopen("graph.out","w",stdout);
#endif
    Len[Bit]=1;for(uint i=0;i<Bit;i++)Len[i]=Lim>>(i+1);
    uint n,m;scanf("%u%u%llu",&n,&m,&c);
    {
        static uint P[25];for(uint i=0;i<n;i++)scanf("%llu",A+i),P[i]=i,A[i]++;
        std::sort(P,P+n,[&](uint a,uint b){return A[a]<A[b];}),std::sort(A,A+n);
        std::vector<uint>User(n);
        for(uint i=0;i<n;i++)User[P[i]]=i;
        for(uint i=0;i<n;i++)P[i]=User[i];
        for(uint i=0,u,v;i<m;i++)scanf("%u%u",&u,&v),u=P[u-1],v=P[v-1],E[u]|=1u<<v,E[v]|=1u<<u;
        // for(uint i=0;i<n;i++)
        //     printf("%u%c",P[i]+1," \n"[i==n-1]);
        // for(uint i=0;i<n;i++)for(uint j=i+1;j<n;j++)
        //     if(E[i]>>j&1)printf("%u %u\n",i+1,j+1);
        // for(uint i=0;i<n;i++)
        //     printf("%llu\n",A[i]);
    }
    if(!m){
        vec a=get(1);
        for(uint i=0;i<n;i++)a=mul(a,get(A[i]));
        query(a,c).println();
        return 0;
    }
    for(uint i=0;i<=n;i++)Log2[1u<<i]=i;
    for(uint i=0;i<(1u<<n);i++){
        bol ok=1;for(uint j=0;j<n;j++)if(i>>j&1)ok&=!(i&E[j]);
        if(ok&&i)Dp[i]=1;
        Dp2[i]=Dp[i];
        for(uint j=(i-1)&i;j;j=(j-1)&i)if(lowbit(j)==lowbit(i))Dp2[i]-=Dp2[j]*Dp[i^j];
        Pop2[i]=Pop2[i>>1]+(i&1);
        if(Pop2[i]&1)T[Log2[lowbit(i)]][Pop2[i]-1][i^lowbit(i)]=Dp2[i];
        else User[i]=Dp2[i]*A[Log2[lowbit(i)]];
    }
    // for(uint i=0;i<(1u<<n);i++)Dp[i].print(),putchar(" \n"[i==(1u<<n)-1]);
    // for(uint i=0;i<(1u<<n);i++)Dp2[i].print(),putchar(" \n"[i==(1u<<n)-1]);
    for(uint i=0;i<n;i++)for(uint j=0;j<n;j++)
        for(uint k=i+1;k<n;k++)for(uint l=0;l<(1u<<n);l+=1u<<(i+1))if(l>>k&1)
            T[i][j][l]+=T[i][j][l^(1u<<k)];
    modint ans;
    for(uint i=0;i<(1u<<n);i++){
        vec a=get(1);
        for(uint j=0;j<n;j++)if(i>>j&1)a=mul(a,get(A[j]));
        G[i]=query(a,c);
    }
    W[0]=1;
    for(uint i=1;i<(1u<<n);i++)for(uint j=i;j;j=(j-1)&i)
        if(lowbit(i)==lowbit(j))W[i]+=W[i^j]*User[j];
    for(uint i=0;i<(1u<<n);i++){
        uint t=(1u<<n)-i-1;
        uint cnt=Pop2[t];
        for(uint j=t;~j;j=j-1){
            Q[0][j&=t]=1;for(uint k=1;k<=cnt;k++)Q[k][j]=0;
            for(uint p=0;p<n;p++)if((i>>p&1)&&!(j&((1u<<(p+1))-1)))for(uint k=cnt;~k;k--){
                Q[k][j]*=T[p][0][j];for(uint q=1;q<=k;q++)Q[k][j]+=T[p][q][j]*Q[k-q][j];
            }
        }
        for(uint k=0;k<=cnt;k++)for(uint a=0;a<n;a++)if(t>>a&1)
            for(uint j=t;j;j=(j-1)&t)if(j>>a&1)
                Q[k][j]-=Q[k][j^(1u<<a)];
        modint wil;
        for(uint j=t;~j;j=j-1)j&=t,wil+=W[j^t]*Q[Pop2[j]][j];
        ans+=wil*G[i];
    }
    ans.println();
    return 0;
}

// g++ graph.cpp -o graph -std=c++11 -DMYEE -Wall