#include<bits/stdc++.h>
using namespace std;
const int maxn=2005; 
struct edge{
	int to,w,nex;
}e[maxn*50];
int tot=-1;
int head[maxn],now[maxn],d[maxn];
int s,t;
inline void add(int x,int y,int w){
	e[++tot].to=y,e[tot].w=w,e[tot].nex=head[x],head[x]=tot;
	e[++tot].to=x,e[tot].w=0,e[tot].nex=head[y],head[y]=tot;
}
inline bool bfs(){
	memset(now,-1,sizeof(now));
	memset(d,0,sizeof(d));
	queue<int> q;
	d[s]=1,now[s]=head[s];
	q.push(s);
	while(!q.empty()){
		int x=q.front();
		q.pop(); 
		for(int i=head[x];~i;i=e[i].nex){
			int y=e[i].to;
			if(e[i].w && !d[y]){
				d[y]=d[x]+1;
				now[y]=head[y];
				q.push(y);
				if(y==t) return 1;
			}
		}
	}
	return 0;
}
const int INF=1e9;
inline int dinic(int x,int flow){
	if(x==t) return flow;
	int rest=flow;
	for(int i=now[x];~i && rest;i=e[i].nex){
		now[x]=i;
		int y=e[i].to;
		if(e[i].w && d[y]==d[x]+1){
			int k=dinic(y,min(e[i].w,rest));
			e[i].w-=k,e[i^1].w+=k,rest-=k;
		}
	}
	return flow-rest;
}
int n,m,k;
int st,ed;
int main(){
	memset(head,-1,sizeof(head));
	scanf("%d%d%d",&n,&m,&k);
	scanf("%d%d",&st,&ed);
	for(int i=1,c;i<=n;i++){
		scanf("%d",&c);
		int x=i*2-1,y=i*2;
		for(int j=1;j<=k;j++){
			add((x-1)*k+j,(y-1)*k+j,c);
			if(j<k) add((x-1)*k+j,(y-1)*k+j+1,INF);
		}
		if(i==ed){
			for(int j=1;j<k;j++) add((y-1)*k+j,(y-1)*k+j+1,INF);
		}
	}
	for(int i=1,x,y;i<=m;i++){
		scanf("%d%d",&x,&y);
		x=x*2,y=y*2-1;
		for(int j=1;j<=k;j++) add((x-1)*k+j,(y-1)*k+j,INF);
	}
	s=(st*2-2)*k+1,t=(ed*2-1)*k+k;
	int ans=0,flow=0;
	while(bfs()) while((flow=dinic(s,INF))) ans+=flow;
	if(ans>=INF){
		puts("-1");
		return 0;
	}
	vector<int> czc;
	for(int i=1;i<=n;i++){
		int x=i*2-1,y=i*2; 
		for(int j=1;j<=k;j++){
			if(d[(x-1)*k+j] && !d[(y-1)*k+j]){
				czc.push_back(i);
				break;
			}
		}
	}
	printf("%d\n",(int)czc.size());
	for(auto x:czc) printf("%d ",x); 
	return 0;
}