#include<bits/stdc++.h>
using namespace std;
const int inf=1e9;
namespace maxflow{
	const int NS=1e5+10,inf=1e9,MS=1e6+10;
	int s,t;
	struct as{
		int v,w;
	}e[MS];
	int k=-1;
	vector<int>f[NS];
	int now[NS],dep[NS];
	bool bfs(){
		memset(now,0,sizeof(now));
		memset(dep,0,sizeof(dep));
		queue<int>q;q.push(s);dep[s]=1;
		while(!q.empty()){
			int u=q.front();q.pop();
			for(int i=0,fs=f[u].size();i<fs;i++){
				int v=e[f[u][i]].v,w=e[f[u][i]].w;
				if(w&&!dep[v]){
					dep[v]=dep[u]+1;
					if(v==t)return 1;
					q.push(v);
				}
			}
		}
		return 0;
	}
	int dfs(int u,int sum){
		if(u==t)return sum;
		int x,fl=0;
		for(int i=now[u],fs=f[u].size();i<fs&&sum;i++){
			now[u]=i;
			int v=e[f[u][i]].v,w=e[f[u][i]].w;
			if(dep[v]==dep[u]+1&&w){
				x=dfs(v,min(sum,w));
				if(!x)dep[v]=0;
				e[f[u][i]].w-=x;e[f[u][i]^1].w+=x;
				fl+=x,sum-=x;
			}
		}
		return fl;
	}
	void add(int x,int y,int z){
		// cout<<x<<" "<<y<<" "<<z<<endl;
		e[++k]=(as){y,z};f[x].push_back(k);
		e[++k]=(as){x,0};f[y].push_back(k);
	}
	int ZDL(int ss,int tt){
		s=ss,t=tt;int ans=0;
		while(bfs())ans+=dfs(s,inf);
		return ans;
	}
	void init(){
		for(int i=0;i<=k;i++)e[i].w=e[i].v=0;
		k=-1;s=t=0;
		for(int i=0;i<NS;i++)f[i].clear();
		memset(dep,0,sizeof(dep)),memset(now,0,sizeof(now));
	}
}
const int N=1e5+10;
int n,m,k,c[N];
int s,t;
int main(){
	scanf("%d%d%d",&n,&m,&k);
	scanf("%d%d",&s,&t);
	int sum=0;
	for(int i=1;i<=n;i++){
		scanf("%d",&c[i]);sum+=c[i];
		int x=i*2,y=i*2+1;
		for(int j=1;j<=k;j++){
			maxflow::add((x-1)*k+j,(y-1)*k+j,c[i]);
			if(j<k){
				maxflow::add((x-1)*k+j,(y-1)*k+j+1,inf);
				// maxflow::add((x-1)*k+j,(x-1)*k+j+1,inf);
				if(i==n)maxflow::add((y-1)*k+j,(y-1)*k+j+1,inf);
			}
		}
	}
	for(int i=1;i<=m;i++){
		int u,v;scanf("%d%d",&u,&v);
		int x=u*2+1,y=v*2;
		for(int j=1;j<=k;j++){
			maxflow::add((x-1)*k+j,(y-1)*k+j,inf);
		}
	}
	s=(s*2-1)*k+1,t=k*(t*2)+k;
	// cout<<s<<" "<<t<<endl;
	int ans=maxflow::ZDL(s,t);
	// cout<<ans<<endl;
	if(ans>sum){
		cout<<-1;return 0;
	}
	maxflow::bfs();
	vector<int>ve;
	for(int i=1;i<=n;i++){
		int x=i*2,y=i*2+1;
		for(int j=1;j<=k;j++){
			if(maxflow::dep[(x-1)*k+j]&&!maxflow::dep[(y-1)*k+j]){
				// cout<<i<<" "<<j<<endl;
				ve.push_back(i);
				// break;
			}
		}
	}
	cout<<ve.size()<<endl;
	for(int i:ve)cout<<i<<" ";
	return 0;
}