#include<map>
#include<set>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<vector>
#include<string>
#include<bitset>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN=205;
const int MAXS=2e3+10;
const int MAXM=2e6+10;
const int INF=0x3f3f3f3f;
const long long LINF=0x3f3f3f3f3f3f3f3f;
int n,m,k,S,T;
int s,t,cnt=0;
int idx=1;
int head[MAXS],nxt[MAXM],edg[MAXM],w[MAXM];
inline void add_edge(int x,int y,int __w){
    idx++;
    nxt[idx]=head[x];
    head[x]=idx;
    edg[idx]=y;
    w[idx]=__w;
}
#define edge(x,y,__w) add_edge(x,y,__w);add_edge(y,x,0);
int dep[MAXS],cur[MAXS];
inline bool bfs(){
    memset(dep,0,sizeof(int)*(cnt+1));
    queue <int> q;
    q.push(s);
    dep[s]=1;
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        cur[x]=head[x];
        for(int i=head[x];i;i=nxt[i])
        {
            if(!w[i]){
                continue;
            }
            int to=edg[i];
            if(!dep[to]){
                dep[to]=dep[x]+1;
                q.push(to);
            }
        }
    }
    return dep[t];
}
int dfs(int x,int flow){
    if(x==t){
        return flow;
    }
    int res=0;
    for(int i=cur[x];i;i=nxt[i])
    {
        cur[x]=i;
        int to=edg[i];
        if((!w[i])||(dep[to]^(dep[x]+1))){
            continue;
        }
        int f=dfs(to,min(flow,w[i]));
        w[i]-=f;
        w[i^1]+=f;
        res+=f;
        flow-=f;
        if(!flow){
            break;
        }
    }
    return res;
}
inline int dinic(){
    int flow=0;
    while(bfs())
    {
        flow+=dfs(s,INF);
    }
    return flow;
}
int a[MAXN];
int in[MAXN][5],out[MAXN][5];
int ans=0;
int Ans[MAXN];
signed main(){
    scanf("%d%d%d",&n,&m,&k);
    scanf("%d%d",&S,&T);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<k;j++)
        {
            cnt++;
            in[i][j]=cnt;
            cnt++;
            out[i][j]=cnt;
            edge(out[i][j],in[i][j],a[i]);
        }
        for(int j=0;j<k-1;j++)
        {
            edge(in[i][j+1],in[i][j],INF);
            edge(out[i][j+1],out[i][j],INF);
            edge(out[i][j+1],in[i][j],INF);
        }
    }
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        for(int j=0;j<k;j++)
        {
            edge(in[y][j],out[x][j],INF);
        }
    }
    s=out[T][k-1];
    t=in[S][0];
    if(dinic()>=INF){
        puts("-1");
        return 0;
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<k;j++)
        {
            if((!dep[in[i][j]])&&dep[out[i][j]]){
                ans++;
                Ans[ans]=i;
            }
        }
    }
    printf("%d\n",ans);
    for(int i=1;i<=ans;i++)
    {
        printf("%d ",Ans[i]);
    }
    putchar('\n');
    return 0;
}